Searched refs:vec_D (Results 1 – 1 of 1) sorted by relevance
1448 CeedScalar *mat_C, *mat_G, *vec_D; in CeedSimultaneousDiagonalization() local1452 CeedCall(CeedCalloc(n, &vec_D)); in CeedSimultaneousDiagonalization()1456 CeedCall(CeedSymmetricSchurDecomposition(ceed, mat_G, vec_D, n)); in CeedSimultaneousDiagonalization()1461 if (fabs(vec_D[j]) > fabs(vec_D[j + 1])) { in CeedSimultaneousDiagonalization()1462 CeedScalarSwap(vec_D[j], vec_D[j + 1]); in CeedSimultaneousDiagonalization()1471 for (CeedInt i = 0; i < n; i++) vec_D[i] = 1. / sqrt(vec_D[i]); in CeedSimultaneousDiagonalization()1476 mat_G[i * n + j] *= vec_D[j]; in CeedSimultaneousDiagonalization()1505 CeedCall(CeedFree(&vec_D)); in CeedSimultaneousDiagonalization()