ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/common/cone.c
(Generate patch)

Comparing ray/src/common/cone.c (file contents):
Revision 1.5 by greg, Wed Jul 17 14:10:16 1991 UTC vs.
Revision 2.11 by greg, Fri Sep 16 15:09:21 2016 UTC

# Line 1 | Line 1
1 – /* Copyright (c) 1991 Regents of the University of California */
2 –
1   #ifndef lint
2 < static char SCCSid[] = "$SunId$ LBL";
2 > static const char       RCSid[] = "$Id$";
3   #endif
6 –
4   /*
5   *  cone.c - routines for making cones
9 – *
10 – *     2/12/86
6   */
7  
8 + #include "copyright.h"
9 +
10   #include  "standard.h"
11  
12   #include  "object.h"
# Line 44 | Line 41 | static char SCCSid[] = "$SunId$ LBL";
41  
42  
43   CONE *
44 < getcone(o, getxf)                       /* get cone structure */
45 < register OBJREC  *o;
46 < int  getxf;
44 > getcone(                                /* get cone structure */
45 >        OBJREC  *o,
46 >        int  getxf
47 > )
48   {
51 –        double  fabs(), sqrt();
49          int  sgn0, sgn1;
50 <        register CONE  *co;
50 >        CONE  *co;
51  
52          if ((co = (CONE *)o->os) == NULL) {
53  
# Line 60 | Line 57 | int  getxf;
57  
58                  co->ca = o->oargs.farg;
59                                                  /* get radii */
60 <                if (o->otype == OBJ_CYLINDER || o->otype == OBJ_TUBE) {
60 >                if ((o->otype == OBJ_CYLINDER) | (o->otype == OBJ_TUBE)) {
61                          if (o->oargs.nfargs != 7)
62                                  goto argcerr;
63                          if (co->ca[6] < -FTINY) {
# Line 70 | Line 67 | int  getxf;
67                                  co->ca[6] = -co->ca[6];
68                          } else if (co->ca[6] <= FTINY)
69                                  goto raderr;
70 +                        co->p0 = 0; co->p1 = 3;
71                          co->r0 = co->r1 = 6;
72                  } else {
73                          if (o->oargs.nfargs != 8)
# Line 82 | Line 80 | int  getxf;
80                          else sgn1 = 0;
81                          if (sgn0+sgn1 == 0)
82                                  goto raderr;
83 <                        if (sgn0 < 0 || sgn1 < 0) {
84 <                                objerror(o, o->otype==OBJ_RING?USER:WARNING,
85 <                                        "negative radii");
83 >                        if ((sgn0 < 0) | (sgn1 < 0)) {
84 >                                if (o->otype == OBJ_RING)
85 >                                        goto raderr;
86 >                                objerror(o, WARNING, "negative radii");
87                                  o->otype = o->otype == OBJ_CONE ?
88                                                  OBJ_CUP : OBJ_CONE;
89                          }
90                          co->ca[6] = co->ca[6]*sgn0;
91                          co->ca[7] = co->ca[7]*sgn1;
92 <                        co->r0 = 6;
94 <                        co->r1 = 7;
95 <                        if (fabs(co->ca[7] - co->ca[6]) <= FTINY) {
92 >                        if (co->ca[7] - co->ca[6] > FTINY) {
93                                  if (o->otype == OBJ_RING)
94 +                                        co->p0 = co->p1 = 0;
95 +                                else {
96 +                                        co->p0 = 0; co->p1 = 3;
97 +                                }
98 +                                co->r0 = 6; co->r1 = 7;
99 +                        } else if (co->ca[6] - co->ca[7] > FTINY) {
100 +                                if (o->otype == OBJ_RING)
101 +                                        co->p0 = co->p1 = 0;
102 +                                else {
103 +                                        co->p0 = 3; co->p1 = 0;
104 +                                }
105 +                                co->r0 = 7; co->r1 = 6;
106 +                        } else {
107 +                                if (o->otype == OBJ_RING)
108                                          goto raderr;
109                                  o->otype = o->otype == OBJ_CONE ?
110                                                  OBJ_CYLINDER : OBJ_TUBE;
111                                  o->oargs.nfargs = 7;
112 <                                co->r1 = 6;
112 >                                co->p0 = 0; co->p1 = 3;
113 >                                co->r0 = co->r1 = 6;
114                          }
115                  }
116                                                  /* get axis orientation */
117 <                co->p0 = 0;
106 <                if (o->otype == OBJ_RING) {
107 <                        if (co->ca[6] > co->ca[7]) {    /* make r0 smaller */
108 <                                co->r0 = 7;
109 <                                co->r1 = 6;
110 <                        }
111 <                        co->p1 = 0;
117 >                if (o->otype == OBJ_RING)
118                          VCOPY(co->ad, o->oargs.farg+3);
119 <                } else {
120 <                        co->p1 = 3;
121 <                        co->ad[0] = co->ca[3] - co->ca[0];
122 <                        co->ad[1] = co->ca[4] - co->ca[1];
117 <                        co->ad[2] = co->ca[5] - co->ca[2];
119 >                else {
120 >                        co->ad[0] = CO_P1(co)[0] - CO_P0(co)[0];
121 >                        co->ad[1] = CO_P1(co)[1] - CO_P0(co)[1];
122 >                        co->ad[2] = CO_P1(co)[2] - CO_P0(co)[2];
123                  }
124                  co->al = normalize(co->ad);
125 <                if (co->al == 0.0)
126 <                        objerror(o, USER, "zero orientation");
125 >                if (co->al == 0.0) {
126 >                        objerror(o, WARNING, "unknown orientation");
127 >                        free(co);
128 >                        return(NULL);
129 >                }
130                                          /* compute axis and side lengths */
131                  if (o->otype == OBJ_RING) {
132                          co->al = 0.0;
133 <                        co->sl = co->ca[co->r1] - co->ca[co->r0];
134 <                } else if (o->otype == OBJ_CONE || o->otype == OBJ_CUP) {
133 >                        co->sl = CO_R1(co) - CO_R0(co);
134 >                } else if ((o->otype == OBJ_CONE) | (o->otype == OBJ_CUP)) {
135                          co->sl = co->ca[7] - co->ca[6];
136                          co->sl = sqrt(co->sl*co->sl + co->al*co->al);
137 <                } else { /* OBJ_CYLINDER || OBJ_TUBE */
137 >                } else { /* OBJ_CYLINDER or OBJ_TUBE */
138                          co->sl = co->al;
139                  }
140                  co->tm = NULL;
# Line 139 | Line 147 | int  getxf;
147   argcerr:
148          objerror(o, USER, "bad # arguments");
149   raderr:
150 <        objerror(o, USER, "illegal radii");
150 >        objerror(o, WARNING, "illegal radii");
151 >        free(co);
152 >        return(NULL);
153   }
154  
155  
156 < freecone(o)                     /* free memory associated with cone */
157 < OBJREC  *o;
156 > void
157 > freecone(OBJREC *o)             /* free memory associated with cone */
158   {
159 <        register CONE  *co = (CONE *)o->os;
159 >        CONE  *co = (CONE *)o->os;
160  
161 <        if (o->os == NULL)
161 >        if (co == NULL)
162                  return;
163          if (co->tm != NULL)
164 <                free((char *)co->tm);
165 <        free(o->os);
164 >                free((void *)co->tm);
165 >        free((void *)co);
166          o->os = NULL;
167   }
168  
169  
170 < conexform(co)                   /* get cone transformation matrix */
171 < register CONE  *co;
170 > void
171 > conexform(CONE *co)                     /* get cone transformation matrix */
172   {
173 <        double  sqrt(), fabs();
174 <        double  m4[4][4];
175 <        register double  d;
166 <        register int  i;
173 >        MAT4  m4;
174 >        double  d;
175 >        int  i;
176  
177 <        co->tm = (double (*)[4])malloc(sizeof(m4));
177 >        co->tm = (RREAL (*)[4])malloc(sizeof(MAT4));
178          if (co->tm == NULL)
179                  error(SYSTEM, "out of memory in conexform");
180  
# Line 174 | Line 183 | register CONE  *co;
183          if (co->r0 == co->r1)
184                  d = 0.0;
185          else
186 <                d = co->ca[co->r0] / (co->ca[co->r1] - co->ca[co->r0]);
186 >                d = CO_R0(co) / (CO_R1(co) - CO_R0(co));
187          for (i = 0; i < 3; i++)
188 <                co->tm[3][i] = d*(co->ca[co->p1+i] - co->ca[co->p0+i])
189 <                                - co->ca[co->p0+i];
188 >                co->tm[3][i] = d*(CO_P1(co)[i] - CO_P0(co)[i])
189 >                                - CO_P0(co)[i];
190          
191                                  /* rotate to positive z-axis */
192          setident4(m4);
# Line 201 | Line 210 | register CONE  *co;
210          multmat4(co->tm, co->tm, m4);
211  
212                                  /* scale z-axis */
213 <        setident4(m4);
214 <        if (co->p0 != co->p1 && co->r0 != co->r1) {
215 <                d = fabs(co->ca[co->r1] - co->ca[co->r0]);
216 <                m4[2][2] = d/co->al;
213 >        if ((co->p0 != co->p1) & (co->r0 != co->r1)) {
214 >                setident4(m4);
215 >                m4[2][2] = (CO_R1(co) - CO_R0(co)) / co->al;
216 >                multmat4(co->tm, co->tm, m4);
217          }
209 –        multmat4(co->tm, co->tm, m4);
218   }

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines