Searched refs:sqrtM1 (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()
1516 sqrtM1[i * dim + k] = 0.0; in DMPlexMetricIntersection_Private()
1520 sqrtM1[i * dim + k] += evecs[j * dim + i] * sqrtj * evecs[j * dim + k]; in DMPlexMetricIntersection_Private()
1571 …0; l < dim; ++l) M2[i * dim + m] += sqrtM1[j * dim + i] * evecs[j * dim + k] * evals[k] * evecs[l … in DMPlexMetricIntersection_Private()
1579 PetscCall(PetscFree4(evecs, sqrtM1, isqrtM1, evals)); in DMPlexMetricIntersection_Private()
|