webgl-tools

annotate vmath.js @ 8:056da157f21f

find changesets by author, revision, files, or words in the commit message
merged something ... can't remember what
author John Tsiombikas <nuclear@member.fsf.org>
date Thu, 18 Aug 2011 00:09:14 +0300
parents
children
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 }