src/px/mx3.c

#ifndef lint
static const char       RCSid[] = "$Id$";
#endif

double
mx3d_adjoint(a, b)
register double a[3][3], b[3][3];
{
    b[0][0] =  ((a[1][1])*( a[2][2]) - ( a[1][2])*( a[2][1]));
    b[1][0] =  ((a[1][2])*( a[2][0]) - ( a[1][0])*( a[2][2]));
    b[2][0] =  ((a[1][0])*( a[2][1]) - ( a[1][1])*( a[2][0]));
    b[0][1] =  ((a[2][1])*( a[0][2]) - ( a[2][2])*( a[0][1]));
    b[1][1] =  ((a[2][2])*( a[0][0]) - ( a[2][0])*( a[0][2]));
    b[2][1] =  ((a[2][0])*( a[0][1]) - ( a[2][1])*( a[0][0]));
    b[0][2] =  ((a[0][1])*( a[1][2]) - ( a[0][2])*( a[1][1]));
    b[1][2] =  ((a[0][2])*( a[1][0]) - ( a[0][0])*( a[1][2]));
    b[2][2] =  ((a[0][0])*( a[1][1]) - ( a[0][1])*( a[1][0]));
    return a[0][0]*b[0][0] + a[1][0]*b[0][1] + a[2][0]*b[0][2];
}


mx3d_mul(a, b, c)
register double a[3][3], b[3][3], c[3][3];
{
    c[0][0] = a[0][0]*b[0][0] + a[0][1]*b[1][0] + a[0][2]*b[2][0];
    c[0][1] = a[0][0]*b[0][1] + a[0][1]*b[1][1] + a[0][2]*b[2][1];
    c[0][2] = a[0][0]*b[0][2] + a[0][1]*b[1][2] + a[0][2]*b[2][2];
    c[1][0] = a[1][0]*b[0][0] + a[1][1]*b[1][0] + a[1][2]*b[2][0];
    c[1][1] = a[1][0]*b[0][1] + a[1][1]*b[1][1] + a[1][2]*b[2][1];
    c[1][2] = a[1][0]*b[0][2] + a[1][1]*b[1][2] + a[1][2]*b[2][2];
    c[2][0] = a[2][0]*b[0][0] + a[2][1]*b[1][0] + a[2][2]*b[2][0];
    c[2][1] = a[2][0]*b[0][1] + a[2][1]*b[1][1] + a[2][2]*b[2][1];
    c[2][2] = a[2][0]*b[0][2] + a[2][1]*b[1][2] + a[2][2]*b[2][2];
}


mx3d_transform(p, a, q)
register double p[3], q[3];
register double a[3][3];
{
    q[0] = p[0]*a[0][0] + p[1]*a[1][0] + p[2]*a[2][0];
    q[1] = p[0]*a[0][1] + p[1]*a[1][1] + p[2]*a[2][1];
    q[2] = p[0]*a[0][2] + p[1]*a[1][2] + p[2]*a[2][2];
}


double
mx3d_transform_div(p, a, q)
register double p[3], q[2];
register double a[3][3];
{
    double q2;
    q2 = p[0]*a[0][2] + p[1]*a[1][2] + p[2]*a[2][2];
    if (q2!=0.) {
        q[0] = (p[0]*a[0][0] + p[1]*a[1][0] + p[2]*a[2][0]) / q2;
        q[1] = (p[0]*a[0][1] + p[1]*a[1][1] + p[2]*a[2][1]) / q2;
    }
    return q2;
}


mx3d_translate_mat(tx, ty, m)
double tx, ty;
register double m[3][3];
{
    m[0][0] = 1.; m[0][1] = 0.; m[0][2] = 0.;
    m[1][0] = 0.; m[1][1] = 1.; m[1][2] = 0.;
    m[2][0] = tx; m[2][1] = ty; m[2][2] = 1.;
}


mx3d_translate(m, tx, ty)
register double m[3][3];
double tx, ty;
{
    m[3][0] += tx*m[0][0]+ty*m[1][0];
    m[3][1] += tx*m[0][1]+ty*m[1][1];
    m[3][2] += tx*m[0][2]+ty*m[1][2];
}


mx3d_scale_mat(sx, sy, m)
double sx, sy;
register double m[3][3];
{
    m[0][0] = sx; m[0][1] = 0.; m[0][2] = 0.;
    m[1][0] = 0.; m[1][1] = sy; m[1][2] = 0.;
    m[2][0] = 0.; m[2][1] = 0.; m[2][2] = 1.;
}


mx3d_scale(m, sx, sy)
register double m[3][3];
double sx, sy;
{
    m[0][0] *= sx; m[0][1] *= sx; m[0][2] *= sx;
    m[1][0] *= sy; m[1][1] *= sy; m[1][2] *= sy;
}
Revision:	2.2
Committed:	Sat Feb 22 02:07:27 2003 UTC (21 years, 2 months ago) by greg
Content type:	text/plain
Branch:	MAIN
CVS Tags:	rad3R5
Changes since 2.1:	+1 -4 lines
Log Message:	Changes and check-in for 3.5 release Includes new source files and modifications not recorded for many years See ray/doc/notes/ReleaseNotes for notes between 3.1 and 3.5 release
#	Content
1	#ifndef lint
2	static const char RCSid[] = "$Id$";
3	#endif
4
5	double
6	mx3d_adjoint(a, b)
7	register double a[3][3], b[3][3];
8	{
9	b[0][0] = ((a[1][1])( a[2][2]) - ( a[1][2])( a[2][1]));
10	b[1][0] = ((a[1][2])( a[2][0]) - ( a[1][0])( a[2][2]));
11	b[2][0] = ((a[1][0])( a[2][1]) - ( a[1][1])( a[2][0]));
12	b[0][1] = ((a[2][1])( a[0][2]) - ( a[2][2])( a[0][1]));
13	b[1][1] = ((a[2][2])( a[0][0]) - ( a[2][0])( a[0][2]));
14	b[2][1] = ((a[2][0])( a[0][1]) - ( a[2][1])( a[0][0]));
15	b[0][2] = ((a[0][1])( a[1][2]) - ( a[0][2])( a[1][1]));
16	b[1][2] = ((a[0][2])( a[1][0]) - ( a[0][0])( a[1][2]));
17	b[2][2] = ((a[0][0])( a[1][1]) - ( a[0][1])( a[1][0]));
18	return a[0][0]b[0][0] + a[1][0]b[0][1] + a[2][0]*b[0][2];
19	}
20
21
22	mx3d_mul(a, b, c)
23	register double a[3][3], b[3][3], c[3][3];
24	{
25	c[0][0] = a[0][0]b[0][0] + a[0][1]b[1][0] + a[0][2]*b[2][0];
26	c[0][1] = a[0][0]b[0][1] + a[0][1]b[1][1] + a[0][2]*b[2][1];
27	c[0][2] = a[0][0]b[0][2] + a[0][1]b[1][2] + a[0][2]*b[2][2];
28	c[1][0] = a[1][0]b[0][0] + a[1][1]b[1][0] + a[1][2]*b[2][0];
29	c[1][1] = a[1][0]b[0][1] + a[1][1]b[1][1] + a[1][2]*b[2][1];
30	c[1][2] = a[1][0]b[0][2] + a[1][1]b[1][2] + a[1][2]*b[2][2];
31	c[2][0] = a[2][0]b[0][0] + a[2][1]b[1][0] + a[2][2]*b[2][0];
32	c[2][1] = a[2][0]b[0][1] + a[2][1]b[1][1] + a[2][2]*b[2][1];
33	c[2][2] = a[2][0]b[0][2] + a[2][1]b[1][2] + a[2][2]*b[2][2];
34	}
35
36
37	mx3d_transform(p, a, q)
38	register double p[3], q[3];
39	register double a[3][3];
40	{
41	q[0] = p[0]a[0][0] + p[1]a[1][0] + p[2]*a[2][0];
42	q[1] = p[0]a[0][1] + p[1]a[1][1] + p[2]*a[2][1];
43	q[2] = p[0]a[0][2] + p[1]a[1][2] + p[2]*a[2][2];
44	}
45
46
47	double
48	mx3d_transform_div(p, a, q)
49	register double p[3], q[2];
50	register double a[3][3];
51	{
52	double q2;
53	q2 = p[0]a[0][2] + p[1]a[1][2] + p[2]*a[2][2];
54	if (q2!=0.) {
55	q[0] = (p[0]a[0][0] + p[1]a[1][0] + p[2]*a[2][0]) / q2;
56	q[1] = (p[0]a[0][1] + p[1]a[1][1] + p[2]*a[2][1]) / q2;
57	}
58	return q2;
59	}
60
61
62	mx3d_translate_mat(tx, ty, m)
63	double tx, ty;
64	register double m[3][3];
65	{
66	m[0][0] = 1.; m[0][1] = 0.; m[0][2] = 0.;
67	m[1][0] = 0.; m[1][1] = 1.; m[1][2] = 0.;
68	m[2][0] = tx; m[2][1] = ty; m[2][2] = 1.;
69	}
70
71
72	mx3d_translate(m, tx, ty)
73	register double m[3][3];
74	double tx, ty;
75	{
76	m[3][0] += txm[0][0]+tym[1][0];
77	m[3][1] += txm[0][1]+tym[1][1];
78	m[3][2] += txm[0][2]+tym[1][2];
79	}
80
81
82	mx3d_scale_mat(sx, sy, m)
83	double sx, sy;
84	register double m[3][3];
85	{
86	m[0][0] = sx; m[0][1] = 0.; m[0][2] = 0.;
87	m[1][0] = 0.; m[1][1] = sy; m[1][2] = 0.;
88	m[2][0] = 0.; m[2][1] = 0.; m[2][2] = 1.;
89	}
90
91
92	mx3d_scale(m, sx, sy)
93	register double m[3][3];
94	double sx, sy;
95	{
96	m[0][0] = sx; m[0][1] = sx; m[0][2] *= sx;
97	m[1][0] = sy; m[1][1] = sy; m[1][2] *= sy;
98	}