| 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(); |
| 43 |
|
* Stuff for tracking and reusing vertices: |
| 44 |
|
*/ |
| 45 |
|
|
| 46 |
< |
char VKFMT[] = "%+1.9e %+1.9e %+1.9e"; |
| 46 |
> |
char VKFMT[] = "%+16.9e %+16.9e %+16.9e"; |
| 47 |
|
#define VKLEN 64 |
| 48 |
|
|
| 49 |
|
#define mkvkey(k,v) sprintf(k, VKFMT, (v)[0], (v)[1], (v)[2]) |
| 50 |
|
|
| 51 |
|
#define NVERTS 256 |
| 52 |
|
|
| 53 |
< |
long clock; /* incremented at each vertex request */ |
| 53 |
> |
long vclock; /* incremented at each vertex request */ |
| 54 |
|
|
| 55 |
|
struct vert { |
| 56 |
|
long lused; /* when last used (0 if unassigned) */ |
| 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); |
| 348 |
|
char *vname; |
| 349 |
|
FVECT vp; |
| 350 |
|
{ |
| 351 |
< |
char vkey[VKLEN]; |
| 351 |
> |
static char vkey[VKLEN]; |
| 352 |
|
register LUENT *lp; |
| 353 |
|
register int i, vndx; |
| 354 |
|
|
| 355 |
< |
clock++; /* increment counter */ |
| 355 |
> |
vclock++; /* increment counter */ |
| 356 |
|
mkvkey(vkey, vp); |
| 357 |
|
if ((lp = lu_find(&vertab, vkey)) == NULL) |
| 358 |
|
goto memerr; |
| 378 |
|
lp->data = (char *)&vert[vndx]; /* set it */ |
| 379 |
|
} else |
| 380 |
|
vndx = (struct vert *)lp->data - vert; |
| 381 |
< |
vert[vndx].lused = clock; /* record this use */ |
| 381 |
> |
vert[vndx].lused = vclock; /* record this use */ |
| 382 |
|
sprintf(vname, "v%d", vndx); |
| 383 |
|
return(vname); |
| 384 |
|
memerr: |
| 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]; |
