[ViewVC] Diff of: radiance/ray/src/common/fvect.c

Comparing ray/src/common/fvect.c (file contents):
Revision 2.18 by greg, Wed Apr 3 00:22:12 2013 UTC vs.
Revision 2.21 by greg, Mon Dec 8 23:51:12 2014 UTC

#	Line 7 \| Line 7 \| static const char RCSid[] = "$Id$";
7
8		#include "copyright.h"
9
10	+	#define _USE_MATH_DEFINES
11		#include <math.h>
12		#include "fvect.h"
13	+	#include "random.h"
14
15	+	double
16	+	Acos(double x) /* insurance for touchy math library */
17	+	{
18	+	if (x <= -1.+FTINY*FTINY)
19	+	return(M_PI);
20	+	if (x >= 1.-FTINY*FTINY)
21	+	return(.0);
22	+	return(acos(x));
23	+	}
24
25		double
26	+	Asin(double x) /* insurance for touchy math library */
27	+	{
28	+	if (x <= -1.+FTINY*FTINY)
29	+	return(-M_PI/2.);
30	+	if (x >= 1.-FTINY*FTINY)
31	+	return(M_PI/2);
32	+	return(asin(x));
33	+	}
34	+
35	+	double
36		fdot( /* return the dot product of two vectors */
37		const FVECT v1,
38		const FVECT v2
#	Line 85 \| Line 106 \| const FVECT v1,
106		const FVECT v2
107		)
108		{
109	+	if ((vres == v1) \| (vres == v2)) {
110	+	FVECT vtmp;
111	+	VCROSS(vtmp, v1, v2);
112	+	VCOPY(vres, vtmp);
113	+	return;
114	+	}
115		VCROSS(vres, v1, v2);
116		}
117
#	Line 125 \| Line 152 \| FVECT v
152		v[2] *= d;
153
154		return(len);
155	+	}
156	+
157	+
158	+	int
159	+	getperpendicular( /* choose random perpedicular direction */
160	+	FVECT vp, /* returns normalized */
161	+	const FVECT v /* input vector must be normalized */
162	+	)
163	+	{
164	+	FVECT v1;
165	+	int i;
166	+	/* randomize other coordinates */
167	+	v1[0] = 0.5 - frandom();
168	+	v1[1] = 0.5 - frandom();
169	+	v1[2] = 0.5 - frandom();
170	+	for (i = 3; i--; )
171	+	if ((-0.6 < v[i]) & (v[i] < 0.6))
172	+	break;
173	+	if (i < 0)
174	+	return(0);
175	+	v1[i] = 1.0;
176	+	fcross(vp, v1, v);
177	+	return(normalize(vp) > 0.0);
178		}
179
180

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