8 |
|
#include "copyright.h" |
9 |
|
|
10 |
|
#include "ray.h" |
11 |
– |
|
11 |
|
#include "otypes.h" |
12 |
+ |
#include "rtotypes.h" |
13 |
|
|
14 |
|
#ifdef DISPERSE |
15 |
|
#include "source.h" |
16 |
< |
static disperse(); |
17 |
< |
static int lambda(); |
16 |
> |
static int disperse(OBJREC *m,RAY *r,FVECT vt,double tr,COLOR cet,COLOR abt); |
17 |
> |
static int lambda(OBJREC *m, FVECT v2, FVECT dv, FVECT lr); |
18 |
|
#endif |
19 |
|
|
20 |
+ |
static double mylog(double x); |
21 |
+ |
|
22 |
+ |
|
23 |
|
/* |
24 |
|
* Explicit calculations for Fresnel's equation are performed, |
25 |
|
* but only one square root computation is necessary. |
52 |
|
|
53 |
|
|
54 |
|
static double |
55 |
< |
mylog(x) /* special log for extinction coefficients */ |
56 |
< |
double x; |
55 |
> |
mylog( /* special log for extinction coefficients */ |
56 |
> |
double x |
57 |
> |
) |
58 |
|
{ |
59 |
|
if (x < 1e-40) |
60 |
|
return(-100.); |
64 |
|
} |
65 |
|
|
66 |
|
|
67 |
< |
m_dielectric(m, r) /* color a ray which hit a dielectric interface */ |
68 |
< |
OBJREC *m; |
69 |
< |
register RAY *r; |
67 |
> |
extern int |
68 |
> |
m_dielectric( /* color a ray which hit a dielectric interface */ |
69 |
> |
OBJREC *m, |
70 |
> |
register RAY *r |
71 |
> |
) |
72 |
|
{ |
73 |
|
double cos1, cos2, nratio; |
74 |
|
COLOR ctrans; |
75 |
|
COLOR talb; |
76 |
|
int hastexture; |
77 |
+ |
double transdist, transtest=0; |
78 |
+ |
double mirdist, mirtest=0; |
79 |
+ |
int flatsurface; |
80 |
|
double refl, trans; |
81 |
|
FVECT dnorm; |
82 |
|
double d1, d2; |
88 |
|
|
89 |
|
raytexture(r, m->omod); /* get modifiers */ |
90 |
|
|
91 |
< |
if (hastexture = DOT(r->pert,r->pert) > FTINY*FTINY) |
91 |
> |
if ( (hastexture = DOT(r->pert,r->pert) > FTINY*FTINY) ) |
92 |
|
cos1 = raynormal(dnorm, r); /* perturb normal */ |
93 |
|
else { |
94 |
|
VCOPY(dnorm, r->ron); |
95 |
|
cos1 = r->rod; |
96 |
|
} |
97 |
+ |
flatsurface = !hastexture && r->ro != NULL && isflat(r->ro->otype); |
98 |
+ |
|
99 |
|
/* index of refraction */ |
100 |
|
if (m->otype == MAT_DIELECTRIC) |
101 |
|
nratio = m->oargs.farg[3] + m->oargs.farg[4]/MLAMBDA; |
194 |
|
rayvalue(&p); |
195 |
|
scalecolor(p.rcol, trans); |
196 |
|
addcolor(r->rcol, p.rcol); |
197 |
< |
if (nratio >= 1.0-FTINY && nratio <= 1.0+FTINY) |
198 |
< |
r->rt = r->rot + p.rt; |
197 |
> |
/* virtual distance */ |
198 |
> |
if (flatsurface || |
199 |
> |
(1.-FTINY <= nratio && |
200 |
> |
nratio <= 1.+FTINY)) { |
201 |
> |
transtest = 2*bright(p.rcol); |
202 |
> |
transdist = r->rot + p.rt; |
203 |
> |
} |
204 |
|
} |
205 |
|
} |
206 |
|
} |
220 |
|
|
221 |
|
scalecolor(p.rcol, refl); /* color contribution */ |
222 |
|
addcolor(r->rcol, p.rcol); |
223 |
+ |
/* virtual distance */ |
224 |
+ |
if (flatsurface) { |
225 |
+ |
mirtest = 2*bright(p.rcol); |
226 |
+ |
mirdist = r->rot + p.rt; |
227 |
+ |
} |
228 |
|
} |
229 |
+ |
/* check distance to return */ |
230 |
+ |
d1 = bright(r->rcol); |
231 |
+ |
if (transtest > d1) |
232 |
+ |
r->rt = transdist; |
233 |
+ |
else if (mirtest > d1) |
234 |
+ |
r->rt = mirdist; |
235 |
|
/* rayvalue() computes absorption */ |
236 |
|
return(1); |
237 |
|
} |
239 |
|
|
240 |
|
#ifdef DISPERSE |
241 |
|
|
242 |
< |
static |
243 |
< |
disperse(m, r, vt, tr, cet, abt) /* check light sources for dispersion */ |
244 |
< |
OBJREC *m; |
245 |
< |
RAY *r; |
246 |
< |
FVECT vt; |
247 |
< |
double tr; |
248 |
< |
COLOR cet, abt; |
242 |
> |
static int |
243 |
> |
disperse( /* check light sources for dispersion */ |
244 |
> |
OBJREC *m, |
245 |
> |
RAY *r, |
246 |
> |
FVECT vt, |
247 |
> |
double tr, |
248 |
> |
COLOR cet, |
249 |
> |
COLOR abt |
250 |
> |
) |
251 |
|
{ |
252 |
|
RAY sray, *entray; |
253 |
|
FVECT v1, v2, n1, n2; |
378 |
|
|
379 |
|
|
380 |
|
static int |
381 |
< |
lambda(m, v2, dv, lr) /* compute lambda for material */ |
382 |
< |
register OBJREC *m; |
383 |
< |
FVECT v2, dv, lr; |
381 |
> |
lambda( /* compute lambda for material */ |
382 |
> |
register OBJREC *m, |
383 |
> |
FVECT v2, |
384 |
> |
FVECT dv, |
385 |
> |
FVECT lr |
386 |
> |
) |
387 |
|
{ |
388 |
|
FVECT lrXdv, v2Xlr; |
389 |
|
double dtmp, denom; |