home *** CD-ROM | disk | FTP | other *** search
/ Stars of Shareware: Raytrace & Morphing / SOS-RAYTRACE.ISO / programm / raytrace / radiance / lib / perezsky.cal < prev    next >
Encoding:
Text File  |  1992-07-02  |  617 b   |  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.