Lines Matching refs:elemMat
2940 PetscScalar *elemMat; in DMPlexComputeInterpolatorNested() local
3003 PetscCall(PetscMalloc1(rTotDim * cTotDim, &elemMat)); in DMPlexComputeInterpolatorNested()
3004 PetscCall(PetscArrayzero(elemMat, rTotDim * cTotDim)); in DMPlexComputeInterpolatorNested()
3066 …for (c = 0; c < Nc; ++c) elemMat[(offsetI + i) * cTotDim + offsetJ + j] += T->T[0][k * cpdim * NcJ… in DMPlexComputeInterpolatorNested()
3087 …for (c = 0; c < Nc; ++c) elemMat[(offsetI + i) * cTotDim + offsetJ + j] += 1.0 * qweights[p * qNc … in DMPlexComputeInterpolatorNested()
3098 if (mesh->printFEM > 1) PetscCall(DMPrintCellMatrix(0, name, rTotDim, cTotDim, elemMat)); in DMPlexComputeInterpolatorNested()
3132 …ined(dmf, fsection, fglobalSection, dmc, csection, cglobalSection, In, c, elemMat, INSERT_VALUES)); in DMPlexComputeInterpolatorNested()
3134 …lSection, PETSC_FALSE, dmc, csection, cglobalSection, PETSC_FALSE, In, c, elemMat, INSERT_VALUES)); in DMPlexComputeInterpolatorNested()
3139 PetscCall(PetscFree(elemMat)); in DMPlexComputeInterpolatorNested()
3182 PetscScalar *elemMat; in DMPlexComputeInterpolatorGeneral() local
3199 PetscCall(PetscMalloc1(totDim, &elemMat)); in DMPlexComputeInterpolatorGeneral()
3286 PetscCall(PetscArrayzero(elemMat, cpdim)); in DMPlexComputeInterpolatorGeneral()
3289 … for (c = 0; c < Nc; ++c) elemMat[j] += T->T[0][j * Nc + c] * qweights[ccell * qNc + c]; in DMPlexComputeInterpolatorGeneral()
3293 for (c = 0; c < Nc; ++c) elemMat[j] += 1.0 * qweights[ccell * qNc + c]; in DMPlexComputeInterpolatorGeneral()
3296 … if (mesh->printFEM > 1) PetscCall(DMPrintCellMatrix(cell, name, 1, numCIndices, elemMat)); in DMPlexComputeInterpolatorGeneral()
3300 PetscCall(MatSetValues(interp, 1, &row, cpdim, &cindices[off], elemMat, INSERT_VALUES)); in DMPlexComputeInterpolatorGeneral()
3321 PetscCall(PetscFree(elemMat)); in DMPlexComputeInterpolatorGeneral()
3355 PetscScalar *elemMat; in DMPlexComputeMassMatrixGeneral() local
3371 PetscCall(PetscMalloc1(totDim, &elemMat)); in DMPlexComputeMassMatrixGeneral()
3523 PetscCall(PetscArrayzero(elemMat, cpdim)); in DMPlexComputeMassMatrixGeneral()
3525 …for (c = 0; c < Nc; ++c) elemMat[j] += T->T[0][j * Nc + c] * Tfine->T[0][(ccell * numFIndices + i)… in DMPlexComputeMassMatrixGeneral()
3528 … if (mesh->printFEM > 1) PetscCall(DMPrintCellMatrix(cell, name, 1, numCIndices, elemMat)); in DMPlexComputeMassMatrixGeneral()
3530 … PetscCall(MatSetValues(mass, 1, &findices[i], numCIndices, cindices, elemMat, ADD_VALUES)); in DMPlexComputeMassMatrixGeneral()
3535 PetscCall(PetscArrayzero(elemMat, cpdim)); in DMPlexComputeMassMatrixGeneral()
3537 for (c = 0; c < Nc; ++c) elemMat[j] += 1.0 * 1.0 * qweights[ccell * Nc + c] * detJ; in DMPlexComputeMassMatrixGeneral()
3540 … if (mesh->printFEM > 1) PetscCall(DMPrintCellMatrix(cell, name, 1, numCIndices, elemMat)); in DMPlexComputeMassMatrixGeneral()
3543 … PetscCall(MatSetValues(mass, 1, &findices[i], numCIndices, cindices, elemMat, ADD_VALUES)); in DMPlexComputeMassMatrixGeneral()
3557 PetscCall(PetscFree(elemMat)); in DMPlexComputeMassMatrixGeneral()
4630 …PetscScalar *work, *u = NULL, *u_t = NULL, *a = NULL, *elemMat = NULL, *elemMatP = NULL, *elemM… in DMPlexComputeJacobian_Patch_Internal() local
4698 elemMat = hasJac ? (sz = chunkSize * totDim * totDim, off += sz, work + off - sz) : NULL; in DMPlexComputeJacobian_Patch_Internal()
4751 …ian(prob, prob, PETSCFE_JACOBIAN, key, Ncell, cgeomFEM, u, u_t, probAux, a, t, X_tShift, elemMat)); in DMPlexComputeJacobian_Patch_Internal()
4766 if (hasJac) elemMat[eOffset + i * totDim + i] = 1.0; in DMPlexComputeJacobian_Patch_Internal()
4774 for (c = 0; c < chunkSize * totDim * totDim; ++c) elemMat[c] += X_tShift * elemMatD[c]; in DMPlexComputeJacobian_Patch_Internal()
4780 …if (hasJac) PetscCall(DMPrintCellMatrix(cell, name, totDim, totDim, &elemMat[(c - cStart) * totDim… in DMPlexComputeJacobian_Patch_Internal()
4783 …ure_Internal(dm, section, globalSection, mesh->useMatClPerm, Jac, cell, &elemMat[(c - cStart) * to… in DMPlexComputeJacobian_Patch_Internal()
4784 …sure_Internal(dm, section, globalSection, mesh->useMatClPerm, JP, cell, &elemMat[(c - cStart) * to… in DMPlexComputeJacobian_Patch_Internal()
5915 PetscScalar *u = NULL, *u_t = NULL, *a = NULL, *elemMat = NULL, *elemMatP = NULL; in DMPlexComputeBdJacobianSingleByLabel() local
5969 …* totDim : 0), &u_t, (hasJac ? (size_t)numFaces * totDim * totDim : 0), &elemMat, (hasPrec ? (size… in DMPlexComputeBdJacobianSingleByLabel()
6003 if (elemMat) PetscCall(PetscArrayzero(elemMat, numFaces * totDim * totDim)); in DMPlexComputeBdJacobianSingleByLabel()
6027 …eBdJacobian(ds, wf, PETSCFE_JACOBIAN, key, Ne, chunkGeom, u, u_t, dsAux, a, t, X_tShift, elemMat)); in DMPlexComputeBdJacobianSingleByLabel()
6034 … dsAux, PetscSafePointerPlusOffset(a, offset * totDimAux), t, X_tShift, &elemMat[offset * totDim *… in DMPlexComputeBdJacobianSingleByLabel()
6045 …nsformPointTensor_Internal(dm, tdm, tv, support[0], PETSC_TRUE, totDim, &elemMat[face * totDim * t… in DMPlexComputeBdJacobianSingleByLabel()
6049 …M > 1) PetscCall(DMPrintCellMatrix(point, "BdJacobian", totDim, totDim, &elemMat[face * totDim * t… in DMPlexComputeBdJacobianSingleByLabel()
6050 …rnal(plex, section, globalSection, mesh->useMatClPerm, Jac, support[0], &elemMat[face * totDim * t… in DMPlexComputeBdJacobianSingleByLabel()
6056 …M > 1) PetscCall(DMPrintCellMatrix(point, "BdJacobian", totDim, totDim, &elemMat[face * totDim * t… in DMPlexComputeBdJacobianSingleByLabel()
6057 …rnal(plex, section, globalSection, mesh->useMatClPerm, Jac, support[0], &elemMat[face * totDim * t… in DMPlexComputeBdJacobianSingleByLabel()
6065 PetscCall(PetscFree5(u, u_t, elemMat, elemMatP, a)); in DMPlexComputeBdJacobianSingleByLabel()
6182 PetscScalar *elemMat, *elemMatP, *elemMatD, *u, *u_t, *a = NULL; in DMPlexComputeJacobianByKey() local
6219 …* totDim : 0), &u_t, (hasJac ? (size_t)numCells * totDim * totDim : 0), &elemMat, (hasPrec ? (size… in DMPlexComputeJacobianByKey()
6243 if (hasJac) PetscCall(PetscArrayzero(elemMat, numCells * totDim * totDim)); in DMPlexComputeJacobianByKey()
6286 …obian(prob, prob, PETSCFE_JACOBIAN, key, Ne, chunkGeom, u, u_t, probAux, a, t, X_tShift, elemMat)); in DMPlexComputeJacobianByKey()
6287 …robAux, PetscSafePointerPlusOffset(a, offset * totDimAux), t, X_tShift, &elemMat[offset * totDim *… in DMPlexComputeJacobianByKey()
6305 for (c = 0; c < numCells * totDim * totDim; ++c) elemMat[c] += X_tShift * elemMatD[c]; in DMPlexComputeJacobianByKey()
6329 if (hasJac) elemMat[eOffset + i * totDim + i] = vol; in DMPlexComputeJacobianByKey()
6332 elemMat[eOffset + i * totDim + i] = vol; in DMPlexComputeJacobianByKey()
6347 …sisTransformPointTensor_Internal(dm, tdm, tv, cell, PETSC_TRUE, totDim, &elemMat[cind * totDim * t… in DMPlexComputeJacobianByKey()
6350 …if (mesh->printFEM > 1) PetscCall(DMPrintCellMatrix(cell, name, totDim, totDim, &elemMat[cind * to… in DMPlexComputeJacobianByKey()
6351 …ure_Internal(dm, section, globalSection, mesh->useMatClPerm, Jac, cell, &elemMat[cind * totDim * t… in DMPlexComputeJacobianByKey()
6357 …if (mesh->printFEM > 1) PetscCall(DMPrintCellMatrix(cell, name, totDim, totDim, &elemMat[cind * to… in DMPlexComputeJacobianByKey()
6358 …re_Internal(dm, section, globalSection, mesh->useMatClPerm, JacP, cell, &elemMat[cind * totDim * t… in DMPlexComputeJacobianByKey()
6364 PetscCall(PetscFree5(u, u_t, elemMat, elemMatP, elemMatD)); in DMPlexComputeJacobianByKey()
6397 PetscScalar *elemMat, *elemMatP, *elemMatD, *u, *u_t, *a = NULL; in DMPlexComputeJacobianByKeyGeneral() local
6442 … totDim : 0), &u_t, (hasJac ? (size_t)numCells * totDim * ctotDim : 0), &elemMat, (hasPrec ? (size… in DMPlexComputeJacobianByKeyGeneral()
6466 if (hasJac) PetscCall(PetscArrayzero(elemMat, numCells * totDim * ctotDim)); in DMPlexComputeJacobianByKeyGeneral()
6509 …eJacobian(rds, cds, PETSCFE_JACOBIAN, key, Ne, chunkGeom, u, u_t, dsAux, a, t, X_tShift, elemMat)); in DMPlexComputeJacobianByKeyGeneral()
6510 … dsAux, PetscSafePointerPlusOffset(a, offset * totDimAux), t, X_tShift, &elemMat[offset * totDim *… in DMPlexComputeJacobianByKeyGeneral()
6528 for (PetscInt c = 0; c < numCells * totDim * ctotDim; ++c) elemMat[c] += X_tShift * elemMatD[c]; in DMPlexComputeJacobianByKeyGeneral()
6552 if (hasJac) elemMat[eOffset + i * ctotDim + i] = vol; in DMPlexComputeJacobianByKeyGeneral()
6555 elemMat[eOffset + i * ctotDim + i] = vol; in DMPlexComputeJacobianByKeyGeneral()
6570 …isTransformPointTensor_Internal(dmr, tdm, tv, cell, PETSC_TRUE, totDim, &elemMat[cind * totDim * c… in DMPlexComputeJacobianByKeyGeneral()
6573 …if (printFEM > 1) PetscCall(DMPrintCellMatrix(cell, name, totDim, ctotDim, &elemMat[cind * totDim … in DMPlexComputeJacobianByKeyGeneral()
6574 …lobalSection, clPerm, dmc, csection, cglobalSection, clPerm, Jac, cell, &elemMat[cind * totDim * c… in DMPlexComputeJacobianByKeyGeneral()
6580 …if (printFEM > 1) PetscCall(DMPrintCellMatrix(cell, name, totDim, ctotDim, &elemMat[cind * totDim … in DMPlexComputeJacobianByKeyGeneral()
6581 …obalSection, clPerm, dmc, csection, cglobalSection, clPerm, JacP, cell, &elemMat[cind * totDim * c… in DMPlexComputeJacobianByKeyGeneral()
6587 PetscCall(PetscFree5(u, u_t, elemMat, elemMatP, elemMatD)); in DMPlexComputeJacobianByKeyGeneral()
7013 PetscScalar *elemMat, *elemMatD, *u, *u_t, *a = NULL, *y, *z; in DMPlexComputeJacobianActionByKey() local
7041 …X_t ? (size_t)numCells * totDim : 0), &u_t, numCells * totDim * totDim, &elemMat, (hasDyn ? (size_… in DMPlexComputeJacobianActionByKey()
7069 PetscCall(PetscArrayzero(elemMat, numCells * totDim * totDim)); in DMPlexComputeJacobianActionByKey()
7105 …obian(prob, prob, PETSCFE_JACOBIAN, key, Ne, chunkGeom, u, u_t, probAux, a, t, X_tShift, elemMat)); in DMPlexComputeJacobianActionByKey()
7106 …robAux, PetscSafePointerPlusOffset(a, offset * totDimAux), t, X_tShift, &elemMat[offset * totDim *… in DMPlexComputeJacobianActionByKey()
7118 for (c = 0; c < numCells * totDim * totDim; ++c) elemMat[c] += X_tShift * elemMatD[c]; in DMPlexComputeJacobianActionByKey()
7128 …PetscCallBLAS("BLASgemv", BLASgemv_("N", &M, &M, &a, &elemMat[cind * totDim * totDim], &M, &y[cind… in DMPlexComputeJacobianActionByKey()
7130 PetscCall(DMPrintCellMatrix(c, name, totDim, totDim, &elemMat[cind * totDim * totDim])); in DMPlexComputeJacobianActionByKey()
7136 PetscCall(PetscFree6(u, u_t, elemMat, elemMatD, y, z)); in DMPlexComputeJacobianActionByKey()