9 |
|
*/ |
10 |
|
|
11 |
|
#include <stdio.h> |
12 |
+ |
#include <math.h> |
13 |
|
#include <string.h> |
14 |
|
#include "fvect.h" |
15 |
|
#include "object.h" |
18 |
|
|
19 |
|
#define PI 3.14159265358979323846 |
20 |
|
|
21 |
+ |
#define C_1SIDEDTHICK 0.005 |
22 |
+ |
|
23 |
|
int o_face(), o_cone(), o_sphere(), o_ring(), o_cylinder(); |
24 |
|
int o_instance(), o_source(), o_illum(); |
25 |
< |
int o_plastic(), o_metal(), o_glass(), o_mirror(), o_trans(), o_light(); |
25 |
> |
int o_plastic(), o_metal(), o_glass(), o_dielectric(), |
26 |
> |
o_mirror(), o_trans(), o_light(); |
27 |
|
|
28 |
|
extern void free(); |
29 |
|
extern char *malloc(); |
288 |
|
add2dispatch("metal", o_metal); |
289 |
|
add2dispatch("metal2", o_metal); |
290 |
|
add2dispatch("glass", o_glass); |
291 |
+ |
add2dispatch("dielectric", o_dielectric); |
292 |
|
add2dispatch("trans", o_trans); |
293 |
|
add2dispatch("trans2", o_trans); |
294 |
|
add2dispatch("mirror", o_mirror); |
629 |
|
newmat(id, NULL); |
630 |
|
if (fa->nfargs == 4) |
631 |
|
nrfr = fa->farg[3]; |
632 |
+ |
printf("\tir %f 0\n", nrfr); |
633 |
|
F = (1. - nrfr)/(1. + nrfr); /* use normal incidence */ |
634 |
|
F *= F; |
635 |
|
for (i = 0; i < 3; i++) { |
644 |
|
if (d > FTINY) |
645 |
|
printf("\t\tcxy %.4f %.4f\n", cxyz[0]/d, cxyz[1]/d); |
646 |
|
printf("\trs %.4f 0\n", cxyz[1]); |
647 |
+ |
rgb_cie(cxyz, trgb); /* put transmitted component */ |
648 |
+ |
puts("\tc"); |
649 |
+ |
d = cxyz[0] + cxyz[1] + cxyz[2]; |
650 |
+ |
if (d > FTINY) |
651 |
+ |
printf("\t\tcxy %.4f %.4f\n", cxyz[0]/d, cxyz[1]/d); |
652 |
+ |
printf("\tts %.4f 0\n", cxyz[1]); |
653 |
+ |
return(0); |
654 |
+ |
} |
655 |
+ |
|
656 |
+ |
|
657 |
+ |
int |
658 |
+ |
o_dielectric(mod, typ, id, fa) /* convert a dielectric material */ |
659 |
+ |
char *mod, *typ, *id; |
660 |
+ |
register FUNARGS *fa; |
661 |
+ |
{ |
662 |
+ |
COLOR cxyz, trgb; |
663 |
+ |
double F, d; |
664 |
+ |
register int i; |
665 |
+ |
|
666 |
+ |
if (fa->nfargs != 5) |
667 |
+ |
return(-1); |
668 |
+ |
newmat(id, NULL); |
669 |
+ |
F = (1. - fa->farg[3])/(1. + fa->farg[3]); /* normal incidence */ |
670 |
+ |
F *= F; |
671 |
+ |
for (i = 0; i < 3; i++) |
672 |
+ |
trgb[i] = (1. - F)*pow(fa->farg[i], C_1SIDEDTHICK/unit_mult); |
673 |
+ |
printf("\tir %f 0\n", fa->farg[3]); /* put index of refraction */ |
674 |
+ |
printf("\tsides 1\n"); |
675 |
+ |
puts("\tc"); /* put reflected component */ |
676 |
+ |
printf("\trs %.4f 0\n", F); |
677 |
|
rgb_cie(cxyz, trgb); /* put transmitted component */ |
678 |
|
puts("\tc"); |
679 |
|
d = cxyz[0] + cxyz[1] + cxyz[2]; |