rev |
line source |
nuclear@0
|
1 /*
|
nuclear@0
|
2 libvmath javascript port (original C version: http://gfxtools.sourceforge.net)
|
nuclear@0
|
3 Copyright (C) 2011 John Tsiombikas <nuclear@member.fsf.org>
|
nuclear@0
|
4
|
nuclear@0
|
5 This program is free software: you can redistribute it and/or modify
|
nuclear@0
|
6 it under the terms of the GNU General Public License as published by
|
nuclear@0
|
7 the Free Software Foundation, either version 3 of the License, or
|
nuclear@0
|
8 (at your option) any later version.
|
nuclear@0
|
9
|
nuclear@0
|
10 This program is distributed in the hope that it will be useful,
|
nuclear@0
|
11 but WITHOUT ANY WARRANTY; without even the implied warranty of
|
nuclear@0
|
12 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
|
nuclear@0
|
13 GNU General Public License for more details.
|
nuclear@0
|
14
|
nuclear@0
|
15 You should have received a copy of the GNU General Public License
|
nuclear@0
|
16 along with this program. If not, see <http://www.gnu.org/licenses/>.
|
nuclear@0
|
17 */
|
nuclear@0
|
18
|
nuclear@0
|
19 function v3_cons(res, x, y, z)
|
nuclear@0
|
20 {
|
nuclear@0
|
21 res[0] = x;
|
nuclear@0
|
22 res[1] = y;
|
nuclear@0
|
23 res[2] = z;
|
nuclear@0
|
24 }
|
nuclear@0
|
25
|
nuclear@0
|
26 function v3_add(res, a, b)
|
nuclear@0
|
27 {
|
nuclear@0
|
28 res[0] = a[0] + b[0];
|
nuclear@0
|
29 res[1] = a[1] + b[1];
|
nuclear@0
|
30 res[2] = a[2] + b[2];
|
nuclear@0
|
31 }
|
nuclear@0
|
32
|
nuclear@0
|
33 function v3_sub(res, a, b)
|
nuclear@0
|
34 {
|
nuclear@0
|
35 res[0] = a[0] - b[0];
|
nuclear@0
|
36 res[1] = a[1] - b[1];
|
nuclear@0
|
37 res[2] = a[2] - b[2];
|
nuclear@0
|
38 }
|
nuclear@0
|
39
|
nuclear@0
|
40 function v3_mul(res, a, b)
|
nuclear@0
|
41 {
|
nuclear@0
|
42 res[0] = a[0] * b[0];
|
nuclear@0
|
43 res[1] = a[1] * b[1];
|
nuclear@0
|
44 res[2] = a[2] * b[2];
|
nuclear@0
|
45 }
|
nuclear@0
|
46
|
nuclear@0
|
47 function v3_div(res, a, b)
|
nuclear@0
|
48 {
|
nuclear@0
|
49 res[0] = a[0] / b[0];
|
nuclear@0
|
50 res[1] = a[1] / b[1];
|
nuclear@0
|
51 res[2] = a[2] / b[2];
|
nuclear@0
|
52 }
|
nuclear@0
|
53
|
nuclear@0
|
54 function v3_dot(a, b)
|
nuclear@0
|
55 {
|
nuclear@0
|
56 return a[0] * b[0] + a[1] * b[1] + a[2] * b[2];
|
nuclear@0
|
57 }
|
nuclear@0
|
58
|
nuclear@0
|
59 function v3_cross(res, a, b)
|
nuclear@0
|
60 {
|
nuclear@0
|
61 res[0] = a[1] * b[2] - a[2] * b[1];
|
nuclear@0
|
62 res[1] = a[2] * b[0] - a[0] * b[2];
|
nuclear@0
|
63 res[2] = a[0] * b[1] - a[1] * b[0];
|
nuclear@0
|
64 }
|
nuclear@0
|
65
|
nuclear@0
|
66 function v3_length(v)
|
nuclear@0
|
67 {
|
nuclear@0
|
68 return Math.sqrt(v[0] * v[0] + v[1] * v[1] + v[2] * v[2]);
|
nuclear@0
|
69 }
|
nuclear@0
|
70
|
nuclear@0
|
71 function v3_length_sq(v)
|
nuclear@0
|
72 {
|
nuclear@0
|
73 return v[0] * v[0] + v[1] * v[1] + v[2] * v[2];
|
nuclear@0
|
74 }
|
nuclear@0
|
75
|
nuclear@0
|
76 function v3_normalize(v)
|
nuclear@0
|
77 {
|
nuclear@0
|
78 var len = v3_length(v);
|
nuclear@0
|
79 v[0] /= len;
|
nuclear@0
|
80 v[1] /= len;
|
nuclear@0
|
81 v[2] /= len;
|
nuclear@0
|
82 }
|
nuclear@0
|
83
|
nuclear@0
|
84 function v3_transform(res, v, m)
|
nuclear@0
|
85 {
|
nuclear@0
|
86 res[0] = v[0] * m[0] + v[1] * m[1] + v[2] * m[2] + m[3];
|
nuclear@0
|
87 res[1] = v[0] * m[4] + v[1] * m[5] + v[2] * m[6] + m[7];
|
nuclear@0
|
88 res[3] = v[0] * m[8] + v[1] * m[9] + v[2] * m[10] + m[11];
|
nuclear@0
|
89 }
|
nuclear@0
|
90
|
nuclear@0
|
91 /* --- matrix 4x4 --- */
|
nuclear@0
|
92
|
nuclear@0
|
93 function m4_cons(res, m0, m1, m2, m3, m4, m5, m6, m7, m8, m9, ma, mb, mc, md, me, mf)
|
nuclear@0
|
94 {
|
nuclear@0
|
95 res[0] = m0;
|
nuclear@0
|
96 res[1] = m1;
|
nuclear@0
|
97 res[2] = m2;
|
nuclear@0
|
98 res[3] = m3;
|
nuclear@0
|
99 res[4] = m4;
|
nuclear@0
|
100 res[5] = m5;
|
nuclear@0
|
101 res[6] = m6;
|
nuclear@0
|
102 res[7] = m7;
|
nuclear@0
|
103 res[8] = m8;
|
nuclear@0
|
104 res[9] = m9;
|
nuclear@0
|
105 res[10] = ma;
|
nuclear@0
|
106 res[11] = mb;
|
nuclear@0
|
107 res[12] = mc;
|
nuclear@0
|
108 res[13] = md;
|
nuclear@0
|
109 res[14] = me;
|
nuclear@0
|
110 res[15] = mf;
|
nuclear@0
|
111 }
|
nuclear@0
|
112
|
nuclear@0
|
113 function m4_identity(res)
|
nuclear@0
|
114 {
|
nuclear@0
|
115 res[0] = res[5] = res[10] = res[15] = 1.0;
|
nuclear@0
|
116 res[1] = res[2] = res[3] = res[4] = res[6] = res[7] = 0.0;
|
nuclear@0
|
117 res[8] = res[9] = res[11] = res[12] = res[13] = res[14] = 0.0;
|
nuclear@0
|
118 }
|
nuclear@0
|
119
|
nuclear@0
|
120 function m4_copy(res, m)
|
nuclear@0
|
121 {
|
nuclear@0
|
122 res[0] = m[0]; res[1] = m[1]; res[2] = m[2]; res[3] = m[3];
|
nuclear@0
|
123 res[4] = m[4]; res[5] = m[5]; res[6] = m[6]; res[7] = m[7];
|
nuclear@0
|
124 res[8] = m[8]; res[9] = m[9]; res[10] = m[10]; res[11] = m[11];
|
nuclear@0
|
125 res[12] = m[12]; res[13] = m[13]; res[14] = m[14]; res[15] = m[15];
|
nuclear@0
|
126 }
|
nuclear@0
|
127
|
nuclear@0
|
128 function m4_mul(res, a, b)
|
nuclear@0
|
129 {
|
nuclear@0
|
130 var a00 = a[0], a01 = a[1], a02 = a[2], a03 = a[3];
|
nuclear@0
|
131 var a10 = a[4], a11 = a[5], a12 = a[6], a13 = a[7];
|
nuclear@0
|
132 var a20 = a[8], a21 = a[9], a22 = a[10], a23 = a[11];
|
nuclear@0
|
133 var a30 = a[12], a31 = a[13], a32 = a[14], a33 = a[15];
|
nuclear@0
|
134 var b00 = b[0], b01 = b[1], b02 = b[2], b03 = b[3];
|
nuclear@0
|
135 var b10 = b[4], b11 = b[5], b12 = b[6], b13 = b[7];
|
nuclear@0
|
136 var b20 = b[8], b21 = b[9], b22 = b[10], b23 = b[11];
|
nuclear@0
|
137 var b30 = b[12], b31 = b[13], b32 = b[14], b33 = b[15];
|
nuclear@0
|
138
|
nuclear@0
|
139 res[0] = b00 * a00 + b01 * a10 + b02 * a20 + b03 * a30;
|
nuclear@0
|
140 res[1] = b00 * a01 + b01 * a11 + b02 * a21 + b03 * a31;
|
nuclear@0
|
141 res[2] = b00 * a02 + b01 * a12 + b02 * a22 + b03 * a32;
|
nuclear@0
|
142 res[3] = b00 * a03 + b01 * a13 + b02 * a23 + b03 * a33;
|
nuclear@0
|
143 res[4] = b10 * a00 + b11 * a10 + b12 * a20 + b13 * a30;
|
nuclear@0
|
144 res[5] = b10 * a01 + b11 * a11 + b12 * a21 + b13 * a31;
|
nuclear@0
|
145 res[6] = b10 * a02 + b11 * a12 + b12 * a22 + b13 * a32;
|
nuclear@0
|
146 res[7] = b10 * a03 + b11 * a13 + b12 * a23 + b13 * a33;
|
nuclear@0
|
147 res[8] = b20 * a00 + b21 * a10 + b22 * a20 + b23 * a30;
|
nuclear@0
|
148 res[9] = b20 * a01 + b21 * a11 + b22 * a21 + b23 * a31;
|
nuclear@0
|
149 res[10] = b20 * a02 + b21 * a12 + b22 * a22 + b23 * a32;
|
nuclear@0
|
150 res[11] = b20 * a03 + b21 * a13 + b22 * a23 + b23 * a33;
|
nuclear@0
|
151 res[12] = b30 * a00 + b31 * a10 + b32 * a20 + b33 * a30;
|
nuclear@0
|
152 res[13] = b30 * a01 + b31 * a11 + b32 * a21 + b33 * a31;
|
nuclear@0
|
153 res[14] = b30 * a02 + b31 * a12 + b32 * a22 + b33 * a32;
|
nuclear@0
|
154 res[15] = b30 * a03 + b31 * a13 + b32 * a23 + b33 * a33;
|
nuclear@0
|
155 }
|
nuclear@0
|
156
|
nuclear@0
|
157 function m4_translation(res, x, y, z)
|
nuclear@0
|
158 {
|
nuclear@0
|
159 res[0] = res[5] = res[10] = res[15] = 1.0;
|
nuclear@0
|
160 res[1] = res[2] = res[3] = res[4] = res[6] = res[7] = res[8] = res[9] = res[11] = 0.0;
|
nuclear@0
|
161 res[12] = x;
|
nuclear@0
|
162 res[13] = y;
|
nuclear@0
|
163 res[14] = z;
|
nuclear@0
|
164 }
|
nuclear@0
|
165
|
nuclear@0
|
166 function m4_rotation(res, angle, x, y, z)
|
nuclear@0
|
167 {
|
nuclear@0
|
168 var sina = Math.sin(angle);
|
nuclear@0
|
169 var cosa = Math.cos(angle);
|
nuclear@0
|
170 var one_minus_cosa = 1.0 - cosa;
|
nuclear@0
|
171 var nxsq = x * x;
|
nuclear@0
|
172 var nysq = y * y;
|
nuclear@0
|
173 var nzsq = z * z;
|
nuclear@0
|
174
|
nuclear@0
|
175 res[0] = nxsq + (1.0 - nxsq) * cosa;
|
nuclear@0
|
176 res[4] = x * y * one_minus_cosa - z * sina;
|
nuclear@0
|
177 res[8] = x * z * one_minus_cosa + y * sina;
|
nuclear@0
|
178 res[1] = x * y * one_minus_cosa + z * sina;
|
nuclear@0
|
179 res[5] = nysq + (1.0 - nysq) * cosa;
|
nuclear@0
|
180 res[9] = y * z * one_minus_cosa - x * sina;
|
nuclear@0
|
181 res[2] = x * z * one_minus_cosa - y * sina;
|
nuclear@0
|
182 res[6] = y * z * one_minus_cosa + x * sina;
|
nuclear@0
|
183 res[10] = nzsq + (1.0 - nzsq) * cosa;
|
nuclear@0
|
184
|
nuclear@0
|
185 res[3] = res[7] = res[11] = res[12] = res[13] = res[14] = 0.0;
|
nuclear@0
|
186 res[15] = 1.0;
|
nuclear@0
|
187 }
|
nuclear@0
|
188
|
nuclear@0
|
189 function m4_scale(res, x, y, z)
|
nuclear@0
|
190 {
|
nuclear@0
|
191 res[0] = x;
|
nuclear@0
|
192 res[5] = y;
|
nuclear@0
|
193 res[10] = z;
|
nuclear@0
|
194 res[15] = 1.0;
|
nuclear@0
|
195 res[1] = res[2] = res[3] = res[4] = res[6] = res[7] = 0.0;
|
nuclear@0
|
196 res[8] = res[9] = res[11] = res[12] = res[13] = res[14] = 0.0;
|
nuclear@0
|
197 }
|
nuclear@0
|
198
|
nuclear@0
|
199 function m4_string(mat)
|
nuclear@0
|
200 {
|
nuclear@0
|
201 return '[' + mat[0] + ', ' + mat[1] + ', ' + mat[2] + ', ' + mat[3] + ']\n' +
|
nuclear@0
|
202 '['+ mat[4] + ', ' + mat[5] + ', ' + mat[6] + ', ' + mat[7] + ']\n' +
|
nuclear@0
|
203 '['+ mat[8] + ', ' + mat[9] + ', ' + mat[10] + ', ' + mat[11] + ']\n' +
|
nuclear@0
|
204 '['+ mat[12] + ', ' + mat[13] + ', ' + mat[14] + ', ' + mat[15] + ']\n';
|
nuclear@0
|
205
|
nuclear@0
|
206 }
|