Searched refs:evecs (Results 1 – 4 of 4) sorted by relevance
| /petsc/src/mat/tests/ |
| H A D | ex118.c | 20 Vec *evecs; in main()
80 PetscCall(PetscMalloc1(nevs + 1, &evecs)); in main()
82 PetscCall(VecCreate(PETSC_COMM_SELF, &evecs[i])); in main()
83 PetscCall(VecSetSizes(evecs[i], PETSC_DECIDE, n)); in main()
84 PetscCall(VecSetFromOptions(evecs[i])); in main()
85 PetscCall(VecPlaceArray(evecs[i], evecs_array + i * n)); in main()
90 PetscCall(CkEigenSolutions(cklvl, T, il - 1, iu - 1, evals, evecs, tols)); in main()
92 for (i = 0; i < nevs; i++) PetscCall(VecResetArray(evecs[i])); in main()
98 for (i = 0; i < nevs; i++) PetscCall(VecDestroy(&evecs[i])); in main()
99 PetscCall(PetscFree(evecs)); in main()
|
| H A D | ex116.c | 18 Vec *evecs; in main() local
104 PetscCall(PetscMalloc1(nevs + 1, &evecs)); in main()
106 PetscCall(VecCreate(PETSC_COMM_SELF, &evecs[i])); in main()
107 PetscCall(VecSetSizes(evecs[i], PETSC_DECIDE, n)); in main()
108 PetscCall(VecSetFromOptions(evecs[i])); in main()
109 PetscCall(VecPlaceArray(evecs[i], evecs_array + i * n)); in main()
113 PetscCall(CkEigenSolutions(cklvl, A, il - 1, iu - 1, evals, evecs, tols)); in main()
116 for (i = 0; i < nevs; i++) PetscCall(VecDestroy(&evecs[i])); in main()
117 PetscCall(PetscFree(evecs)); in main()
|
| H A D | ex120.c | 12 Vec *evecs; in main() local
208 PetscCall(PetscMalloc1(nevs + 1, &evecs)); in main()
210 PetscCall(VecCreate(PETSC_COMM_SELF, &evecs[i])); in main()
211 PetscCall(VecSetSizes(evecs[i], PETSC_DECIDE, n)); in main()
212 PetscCall(VecSetFromOptions(evecs[i])); in main()
213 PetscCall(VecPlaceArray(evecs[i], evecs_array + i * n)); in main()
218 PetscCall(CkEigenSolutions(cklvl, A, il - 1, iu - 1, evals, evecs, tols)); in main()
219 for (i = 0; i < nevs; i++) PetscCall(VecDestroy(&evecs[i])); in main()
220 PetscCall(PetscFree(evecs)); in main()
|
| /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()
1482 for (j = 0; j < dim; ++j) evecs[i * dim + j] = M1[i * dim + j]; in DMPlexMetricIntersection_Private()
1501 …PetscCallBLAS("LAPACKsyev", LAPACKsyev_("V", "U", &nb, evecs, &nb, evals, work, &lwork, rwork, &li… in DMPlexMetricIntersection_Private()
1505 … PetscCallBLAS("LAPACKsyev", LAPACKsyev_("V", "U", &nb, evecs, &nb, evals, work, &lwork, &lierr)); in DMPlexMetricIntersection_Private()
1520 sqrtM1[i * dim + k] += evecs[j * dim + i] * sqrtj * evecs[j * dim + k]; in DMPlexMetricIntersection_Private()
1521 isqrtM1[i * dim + k] += evecs[j * dim + i] * (1.0 / sqrtj) * evecs[j * dim + k]; in DMPlexMetricIntersection_Private()
1529 evecs[i * dim + l] = 0.0; 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()
1542 …PetscCallBLAS("LAPACKsyev", LAPACKsyev_("V", "U", &nb, evecs, &nb, evals, work, &lwork, rwork, &li… in DMPlexMetricIntersection_Private()
[all …]
|