Searched refs:isqrtM1 (Results 1 – 1 of 1) sorted by relevance
| /petsc/src/dm/impls/plex/ |
| H A D | plexmetric.c | 1470 PetscScalar *evecs, *sqrtM1, *isqrtM1; in DMPlexMetricIntersection_Private() local
1480 PetscCall(PetscMalloc4(dim * dim, &evecs, dim * dim, &sqrtM1, dim * dim, &isqrtM1, dim, &evals)); in DMPlexMetricIntersection_Private()
1517 isqrtM1[i * dim + k] = 0.0; in DMPlexMetricIntersection_Private()
1521 isqrtM1[i * dim + k] += evecs[j * dim + i] * (1.0 / sqrtj) * evecs[j * dim + k]; in DMPlexMetricIntersection_Private()
1531 …for (k = 0; k < dim; ++k) evecs[i * dim + l] += isqrtM1[j * dim + i] * M2[j * dim + k] * isqrtM1[k… in DMPlexMetricIntersection_Private()
1553 …for (k = 0; k < dim; ++k) evecs[i * dim + l] += isqrtM1[j * dim + i] * M2[j * dim + k] * isqrtM1[k… in DMPlexMetricIntersection_Private()
1579 PetscCall(PetscFree4(evecs, sqrtM1, isqrtM1, evals)); in DMPlexMetricIntersection_Private()
|