348     uint32_t i, j, k, l, m;  in mat33GetEigenbasis()  local
369         l = ind[m];  in mat33GetEigenbasis()
370         float p = S->elem[k][l];  in mat33GetEigenbasis()
376         float y = (_eigenvals[l] - _eigenvals[k]) * 0.5f;  in mat33GetEigenbasis()
389         S->elem[k][l] = 0.0f;  in mat33GetEigenbasis()
392         _eigenvals[l] += t;  in mat33GetEigenbasis()
395             mat33Rotate(S, c, s, i, k, i, l);  in mat33GetEigenbasis()
398         for (i = k + 1; i < l; ++i) {  in mat33GetEigenbasis()
399             mat33Rotate(S, c, s, k, i, i, l);  in mat33GetEigenbasis()
402         for (i = l + 1; i < N; ++i) {  in mat33GetEigenbasis()
403             mat33Rotate(S, c, s, k, i, l, i);  in mat33GetEigenbasis()
407             float tmp = c * eigenvecs->elem[k][i] - s * eigenvecs->elem[l][i];  in mat33GetEigenbasis()
408             eigenvecs->elem[l][i] = s * eigenvecs->elem[k][i] + c * eigenvecs->elem[l][i];  in mat33GetEigenbasis()
413         ind[l] = mat33Maxind(S, l);  in mat33GetEigenbasis()
429         for (l = k + 1; l < N; ++l) {  in mat33GetEigenbasis()
430             if (_eigenvals[l] > _eigenvals[m]) {  in mat33GetEigenbasis()
431                 m = l;  in mat33GetEigenbasis()
466 void mat33Rotate(struct Mat33 *A, float c, float s, uint32_t k, uint32_t l, uint32_t i, uint32_t j)  in mat33Rotate()  argument
468     float tmp = c * A->elem[k][l] - s * A->elem[i][j];  in mat33Rotate()
469     A->elem[i][j] = s * A->elem[k][l] + c * A->elem[i][j];  in mat33Rotate()
470     A->elem[k][l] = tmp;  in mat33Rotate()