ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/rt/o_instance.c
Revision: 2.7
Committed: Mon Jul 21 22:30:19 2003 UTC (20 years, 9 months ago) by schorsch
Content type: text/plain
Branch: MAIN
Changes since 2.6: +2 -2 lines
Log Message:
Eliminated copystruct() macro, which is unnecessary in ANSI.
Reduced ambiguity warnings for nested if/if/else clauses.

File Contents

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