home *** CD-ROM | disk | FTP | other *** search
/ DP Tool Club 8 / CDASC08.ISO / NEWS / RADIANCE / LIB / PEREZSKY.CAL < prev    next >
Text File  |  1993-10-07  |  617b  |  24 lines
  1. {
  2.     Sky luminance model developed by Richard Perez, Robert Seals,
  3.     Joseph Michalsky and Pierre Ineichen
  4.  
  5.     7/1/92
  6. }
  7.  
  8. fact(n) : if(n, n*fact(n-1), 1);    { factorial -- n must be integer }
  9.  
  10. A=.3; B=6; C=.1310; D=3; E=.3200;    { Test coefficients }
  11.  
  12. Coef1 : (1-A)/fact(B)/C^(B+1);
  13. Coef2 : A/fact(D)/E^(D+1);
  14.  
  15.             { Probability density function for intensity q }
  16. P(q) = Coef1*q^B*exp(-q/C) + Coef2*q^D*exp(-q/E);
  17.  
  18.             { Integral used to generate random intensity }
  19. QsumStep : .02;
  20. Qfind(X,q,s) = if( s-X, -QsumStep/2+q,
  21.         Qfind(X, q+QsumStep, s+QsumStep*P(QsumStep/2+q)) );
  22.  
  23. Qrand(x) = Qfind(rand(x),0,0);
  24.