--- ray/src/rt/rayinit.cal 1996/02/29 11:13:34 2.10 +++ ray/src/rt/rayinit.cal 1996/10/31 11:05:35 2.11 @@ -50,8 +50,8 @@ hermite(p0,p1,r0,r1,t) - 1-dimensional hermite polynomial - noise3(x,y,z), noise3a(x,y,z), - noise3b(x,y,z), noise3c(x,y,z) - noise function with gradient (-1 to 1) + noise3(x,y,z), noise3x(x,y,z), + noise3y(x,y,z), noise3z(x,y,z) - noise function with gradient (-1 to 1) fnoise3(x,y,z) - fractal noise function (-1 to 1) } @@ -61,11 +61,15 @@ AC = arg(0); A1 = arg(1); A2 = arg(2); A3 = arg(3); A4 = arg(4); A5 = arg(5); A6 = arg(6); A7 = arg(7); A8 = arg(8); A9 = arg(9); A10 = arg(10); +noise3a(x,y,z) : noise3x(x,y,z); +noise3b(x,y,z) : noise3y(x,y,z); +noise3c(x,y,z) : noise3z(x,y,z); + { Forward compatibility (?) } D(i) = select(i, Dx, Dy, Dz); N(i) = select(i, Nx, Ny, Nz); P(i) = select(i, Px, Py, Pz); -noise3d(i,x,y,z) = select(i, noise3a(x,y,z), noise3b(x,y,z), noise3c(x,y,z)); +noise3d(i,x,y,z) : select(i, noise3x(x,y,z), noise3y(x,y,z), noise3z(x,y,z)); { More robust versions of library functions } bound(a,x,b) : if(a-x, a, if(x-b, b, x)); @@ -128,13 +132,13 @@ bspline(pp, p0, p1, pn, t) = pp * (1/6+t*(-.5+t*(.5-1/ turbulence(x,y,z,s) = if( s-1.01, 0, abs(noise3(x/s,y/s,z/s)*s) + turbulence(x,y,z,2*s) ); turbulencea(x,y,z,s) = if( s-1.01, 0, - sgn(noise3(x/s,y/s,z/s))*noise3a(x/s,y/s,z/s) + + sgn(noise3(x/s,y/s,z/s))*noise3x(x/s,y/s,z/s) + turbulencea(x,y,z,2*s) ); turbulenceb(x,y,z,s) = if( s-1.01, 0, - sgn(noise3(x/s,y/s,z/s))*noise3b(x/s,y/s,z/s) + + sgn(noise3(x/s,y/s,z/s))*noise3y(x/s,y/s,z/s) + turbulenceb(x,y,z,2*s) ); turbulencec(x,y,z,s) = if( s-1.01, 0, - sgn(noise3(x/s,y/s,z/s))*noise3c(x/s,y/s,z/s) + + sgn(noise3(x/s,y/s,z/s))*noise3z(x/s,y/s,z/s) + turbulencec(x,y,z,2*s) ); { Normal distribution from uniform range (0,1) }