1 |
greg |
1.1 |
#ifndef lint |
2 |
greg |
2.11 |
static const char RCSid[] = "$Id: t_func.c,v 2.10 2014/07/08 18:25:00 greg Exp $"; |
3 |
greg |
1.1 |
#endif |
4 |
|
|
/* |
5 |
|
|
* t_func.c - routine for procedural textures. |
6 |
greg |
2.4 |
*/ |
7 |
|
|
|
8 |
greg |
2.5 |
#include "copyright.h" |
9 |
greg |
1.1 |
|
10 |
|
|
#include "ray.h" |
11 |
greg |
2.2 |
#include "func.h" |
12 |
schorsch |
2.7 |
#include "rtotypes.h" |
13 |
greg |
2.2 |
|
14 |
greg |
1.1 |
/* |
15 |
|
|
* A procedural texture perturbs the surface normal |
16 |
|
|
* at the point of intersection with an object. It has |
17 |
|
|
* the form: |
18 |
|
|
* |
19 |
|
|
* modifier texfunc name |
20 |
|
|
* 4+ xvarname yvarname zvarname filename xf |
21 |
|
|
* 0 |
22 |
|
|
* n A1 A2 .. |
23 |
|
|
* |
24 |
|
|
* Filename is the name of the file where the variable definitions |
25 |
|
|
* can be found. The list of real arguments can be accessed by |
26 |
|
|
* definitions in the file. The xf is a transformation to get |
27 |
|
|
* from the original coordinates to the current coordinates. |
28 |
|
|
*/ |
29 |
|
|
|
30 |
|
|
|
31 |
greg |
2.10 |
int |
32 |
schorsch |
2.7 |
t_func( /* compute texture for ray */ |
33 |
greg |
2.10 |
OBJREC *m, |
34 |
|
|
RAY *r |
35 |
schorsch |
2.7 |
) |
36 |
greg |
1.1 |
{ |
37 |
|
|
FVECT disp; |
38 |
greg |
1.6 |
double d; |
39 |
greg |
2.10 |
MFUNC *mf; |
40 |
|
|
int i; |
41 |
greg |
1.1 |
|
42 |
|
|
if (m->oargs.nsargs < 4) |
43 |
|
|
objerror(m, USER, "bad # arguments"); |
44 |
greg |
2.2 |
mf = getfunc(m, 3, 0x7, 1); |
45 |
|
|
setfunc(m, r); |
46 |
greg |
1.1 |
errno = 0; |
47 |
greg |
2.2 |
for (i = 0; i < 3; i++) { |
48 |
|
|
disp[i] = evalue(mf->ep[i]); |
49 |
greg |
2.11 |
if ((errno == EDOM) | (errno == ERANGE)) { |
50 |
greg |
2.2 |
objerror(m, WARNING, "compute error"); |
51 |
greg |
2.3 |
return(0); |
52 |
greg |
2.2 |
} |
53 |
greg |
1.1 |
} |
54 |
greg |
2.9 |
if (mf->fxp != &unitxf) |
55 |
|
|
multv3(disp, disp, mf->fxp->xfm); |
56 |
greg |
1.8 |
if (r->rox != NULL) { |
57 |
|
|
multv3(disp, disp, r->rox->f.xfm); |
58 |
greg |
2.9 |
d = 1.0 / (mf->fxp->sca * r->rox->f.sca); |
59 |
greg |
1.8 |
} else |
60 |
greg |
2.9 |
d = 1.0 / mf->fxp->sca; |
61 |
greg |
2.8 |
VSUM(r->pert, r->pert, disp, d); |
62 |
greg |
2.3 |
return(0); |
63 |
greg |
1.1 |
} |