Searched refs:invVscalar (Results 1 – 1 of 1) sorted by relevance
| /petsc/src/dm/dt/fe/impls/composite/ |
| H A D | fecomposite.c | 25 PetscScalar *work, *invVscalar; in PetscFESetUp_Composite() local
66 PetscCall(PetscMalloc1(cmp->numSubelements * spdim * spdim, &invVscalar)); in PetscFESetUp_Composite()
68 invVscalar = fem->invV; in PetscFESetUp_Composite()
83 invVscalar[(s * spdim + j) * spdim + k] = 0.0; in PetscFESetUp_Composite()
84 …for (q = 0; q < Nq; ++q) invVscalar[(s * spdim + j) * spdim + k] += Bf[q * spdim + k] * weights[q]; in PetscFESetUp_Composite()
89 …PetscCallBLAS("LAPACKgetrf", LAPACKgetrf_(&n, &n, &invVscalar[s * spdim * spdim], &n, pivots, &inf… in PetscFESetUp_Composite()
90 …PetscCallBLAS("LAPACKgetri", LAPACKgetri_(&n, &invVscalar[s * spdim * spdim], &n, pivots, work, &n… in PetscFESetUp_Composite()
93 …for (s = 0; s < cmp->numSubelements * spdim * spdim; s++) fem->invV[s] = PetscRealPart(invVscalar[… in PetscFESetUp_Composite()
94 PetscCall(PetscFree(invVscalar)); in PetscFESetUp_Composite()
|