ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/rt/o_instance.c
Revision: 2.8
Committed: Tue Mar 30 16:13:01 2004 UTC (20 years ago) by schorsch
Content type: text/plain
Branch: MAIN
CVS Tags: rad3R7P2, rad3R7P1, rad4R1, rad4R0, rad3R6, rad3R6P1, rad3R8, rad3R9
Changes since 2.7: +7 -5 lines
Log Message:
Continued ANSIfication. There are only bits and pieces left now.

File Contents

# Content
1 #ifndef lint
2 static const char RCSid[] = "$Id: o_instance.c,v 2.7 2003/07/21 22:30:19 schorsch Exp $";
3 #endif
4 /*
5 * o_instance.c - routines for computing ray intersections with octrees.
6 */
7
8 #include "copyright.h"
9
10 #include "ray.h"
11 #include "instance.h"
12 #include "rtotypes.h"
13
14
15 extern int
16 o_instance( /* compute ray intersection with octree */
17 OBJREC *o,
18 register RAY *r
19 )
20 {
21 RAY rcont;
22 double d;
23 register INSTANCE *ins;
24 register int i;
25 /* get the octree */
26 ins = getinstance(o, IO_ALL);
27 /* copy and transform ray */
28 rcont = *r;
29 multp3(rcont.rorg, r->rorg, ins->x.b.xfm);
30 multv3(rcont.rdir, r->rdir, ins->x.b.xfm);
31 for (i = 0; i < 3; i++)
32 rcont.rdir[i] /= ins->x.b.sca;
33 rcont.rmax *= ins->x.b.sca;
34 /* clear and trace it */
35 rayclear(&rcont);
36 if (!localhit(&rcont, &ins->obj->scube))
37 return(0); /* missed */
38 if (rcont.rot * ins->x.f.sca >= r->rot)
39 return(0); /* not close enough */
40
41 if (o->omod != OVOID) { /* if we have modifier, use it */
42 r->ro = o;
43 r->rox = NULL;
44 } else { /* else use theirs */
45 r->ro = rcont.ro;
46 if (rcont.rox != NULL) {
47 newrayxf(r); /* allocate transformation */
48 /* NOTE: r->rox may equal rcont.rox! */
49 multmat4(r->rox->f.xfm, rcont.rox->f.xfm, ins->x.f.xfm);
50 r->rox->f.sca = rcont.rox->f.sca * ins->x.f.sca;
51 multmat4(r->rox->b.xfm, ins->x.b.xfm, rcont.rox->b.xfm);
52 r->rox->b.sca = ins->x.b.sca * rcont.rox->b.sca;
53 } else
54 r->rox = &ins->x;
55 }
56 /* transform it back */
57 r->rot = rcont.rot * ins->x.f.sca;
58 multp3(r->rop, rcont.rop, ins->x.f.xfm);
59 multv3(r->ron, rcont.ron, ins->x.f.xfm);
60 multv3(r->pert, rcont.pert, ins->x.f.xfm);
61 d = 1./ins->x.f.sca;
62 for (i = 0; i < 3; i++) {
63 r->ron[i] *= d;
64 r->pert[i] *= d;
65 }
66 r->rod = rcont.rod;
67 r->uv[0] = rcont.uv[0];
68 r->uv[1] = rcont.uv[1];
69 /* return hit */
70 return(1);
71 }