Lines Matching refs:elem
30 memset(A->elem, 0.0f, sizeof(A->elem));
40 A->elem[i][i] = x;
46 A->elem[0][0] = v1->x;
47 A->elem[0][1] = v2->x;
48 A->elem[0][2] = v3->x;
50 A->elem[1][0] = v1->y;
51 A->elem[1][1] = v2->y;
52 A->elem[1][2] = v3->y;
54 A->elem[2][0] = v1->z;
55 A->elem[2][1] = v2->z;
56 A->elem[2][2] = v3->z;
62 out->x = A->elem[0][0] * v->x + A->elem[0][1] * v->y + A->elem[0][2] * v->z;
63 out->y = A->elem[1][0] * v->x + A->elem[1][1] * v->y + A->elem[1][2] * v->z;
64 out->z = A->elem[2][0] * v->x + A->elem[2][1] * v->y + A->elem[2][2] * v->z;
80 sum += A->elem[i][k] * B->elem[k][j];
83 out->elem[i][j] = sum;
95 A->elem[i][j] *= c;
107 out->elem[i][j] += A->elem[i][j];
119 out->elem[i][j] -= A->elem[i][j];
129 if (A->elem[i][i] < 0.0f) {
137 if (fabsf(A->elem[i][j] - A->elem[j][i]) > tolerance) {
159 sum += A->elem[k][i] * B->elem[k][j];
162 out->elem[i][j] = sum;
180 sum += A->elem[i][k] * B->elem[j][k];
183 out->elem[i][j] = sum;
201 if (fabsf(tmp.elem[j][i]) > fabsf(tmp.elem[i][i])) {
208 t = tmp.elem[i][k];
209 tmp.elem[i][k] = tmp.elem[swap][k];
210 tmp.elem[swap][k] = t;
212 t = out->elem[i][k];
213 out->elem[i][k] = out->elem[swap][k];
214 out->elem[swap][k] = t;
218 t = 1.0f / tmp.elem[i][i];
220 tmp.elem[i][k] *= t;
221 out->elem[i][k] *= t;
226 t = tmp.elem[j][i];
228 tmp.elem[j][k] -= tmp.elem[i][k] * t;
229 out->elem[j][k] -= out->elem[i][k] * t;
243 out->elem[i][j] = A->elem[j][i];
256 pivot->elem[k] = k;
257 float max = fabsf(LU->elem[k][k]);
259 if (max < fabsf(LU->elem[j][k])) {
260 max = fabsf(LU->elem[j][k]);
261 pivot->elem[k] = j;
265 if (pivot->elem[k] != k) {
266 mat33SwapRows(LU, k, pivot->elem[k]);
269 if (fabsf(LU->elem[k][k]) < kEps) {
274 LU->elem[k][j] /= LU->elem[k][k];
279 LU->elem[i][j] -= LU->elem[i][k] * LU->elem[k][j];
297 float tmp = A->elem[i][k];
298 A->elem[i][k] = A->elem[j][k];
299 A->elem[j][k] = tmp;
317 if (pivot->elem[k] != k) {
319 bCopy[k] = bCopy[pivot->elem[k]];
320 bCopy[pivot->elem[k]] = tmp;
325 _x[k] -= _x[i] * A->elem[k][i];
327 _x[k] /= A->elem[k][k];
332 _x[k] -= _x[i] * A->elem[k][i];
355 _eigenvals[k] = S->elem[k][k];
363 if (fabsf(S->elem[k][ind[k]]) > fabsf(S->elem[m][ind[m]])) {
370 float p = S->elem[k][l];
389 S->elem[k][l] = 0.0f;
407 float tmp = c * eigenvecs->elem[k][i] - s * eigenvecs->elem[l][i];
408 eigenvecs->elem[l][i] = s * eigenvecs->elem[k][i] + c * eigenvecs->elem[l][i];
409 eigenvecs->elem[k][i] = tmp;
418 sum += fabsf(S->elem[i][j]);
458 if (fabsf(A->elem[k][i]) > fabsf(A->elem[k][m])) {
468 float tmp = c * A->elem[k][l] - s * A->elem[i][j];
469 A->elem[i][j] = s * A->elem[k][l] + c * A->elem[i][j];
470 A->elem[k][l] = tmp;
478 A->elem[0][0] * v->x
479 + A->elem[0][1] * v->y
480 + A->elem[0][2] * v->z
481 + A->elem[0][3] * v->w;
484 A->elem[1][0] * v->x
485 + A->elem[1][1] * v->y
486 + A->elem[1][2] * v->z
487 + A->elem[1][3] * v->w;
490 A->elem[2][0] * v->x
491 + A->elem[2][1] * v->y
492 + A->elem[2][2] * v->z
493 + A->elem[2][3] * v->w;
496 A->elem[3][0] * v->x
497 + A->elem[3][1] * v->y
498 + A->elem[3][2] * v->z
499 + A->elem[3][3] * v->w;
509 pivot->elem[k] = k;
510 float max = fabsf(LU->elem[k][k]);
512 if (max < fabsf(LU->elem[j][k])) {
513 max = fabsf(LU->elem[j][k]);
514 pivot->elem[k] = j;
518 if (pivot->elem[k] != k) {
519 mat44SwapRows(LU, k, pivot->elem[k]);
522 if (fabsf(LU->elem[k][k]) < kEps) {
527 LU->elem[k][j] /= LU->elem[k][k];
532 LU->elem[i][j] -= LU->elem[i][k] * LU->elem[k][j];
550 float tmp = A->elem[i][k];
551 A->elem[i][k] = A->elem[j][k];
552 A->elem[j][k] = tmp;
571 if (pivot->elem[k] != k) {
573 bCopy[k] = bCopy[pivot->elem[k]];
574 bCopy[pivot->elem[k]] = tmp;
579 _x[k] -= _x[i] * A->elem[k][i];
581 _x[k] /= A->elem[k][k];
586 _x[k] -= _x[i] * A->elem[k][i];
