1352   Mat                 intMat;  in PetscDualSpaceLagrangeCreateSimplexNodeMat()  local
1414   PetscCall(MatCreateSeqAIJ(PETSC_COMM_SELF, nNodes * Nk, nNodes * Nk, Nk, NULL, &intMat));  in PetscDualSpaceLagrangeCreateSimplexNodeMat()
1415   PetscCall(PetscObjectSetOptionsPrefix((PetscObject)intMat, "lag_"));  in PetscDualSpaceLagrangeCreateSimplexNodeMat()
1416   PetscCall(MatSetOption(intMat, MAT_IGNORE_ZERO_ENTRIES, PETSC_FALSE));  in PetscDualSpaceLagrangeCreateSimplexNodeMat()
1422 …for (a = aprev; a < anext; a++) PetscCall(MatSetValue(intMat, j, a, (a == j) ? 1. : 0., INSERT_VAL…  in PetscDualSpaceLagrangeCreateSimplexNodeMat()
1424   PetscCall(MatAssemblyBegin(intMat, MAT_FINAL_ASSEMBLY));  in PetscDualSpaceLagrangeCreateSimplexNodeMat()
1425   PetscCall(MatAssemblyEnd(intMat, MAT_FINAL_ASSEMBLY));  in PetscDualSpaceLagrangeCreateSimplexNodeMat()
1427   *iMat        = intMat;  in PetscDualSpaceLagrangeCreateSimplexNodeMat()
1516     Mat             intMat;  in PetscDualSpaceCreateAllDataFromInteriorData()  local
1524     PetscCall(PetscDualSpaceGetInteriorData(psp, &intNodes, &intMat));  in PetscDualSpaceCreateAllDataFromInteriorData()
1531     if (intMat) {  in PetscDualSpaceCreateAllDataFromInteriorData()
1535       PetscCall(MatSeqAIJGetMaxRowNonzeros(intMat, &maxNzsp));  in PetscDualSpaceCreateAllDataFromInteriorData()
1554     Mat             intMat;  in PetscDualSpaceCreateAllDataFromInteriorData()  local
1561     PetscCall(PetscDualSpaceGetInteriorData(psp, &intNodes, &intMat));  in PetscDualSpaceCreateAllDataFromInteriorData()
1562     if (intNodes == NULL && intMat == NULL) continue;  in PetscDualSpaceCreateAllDataFromInteriorData()
1596     if (intMat) {  in PetscDualSpaceCreateAllDataFromInteriorData()
1609         PetscCall(MatGetRow(intMat, j, &ncols, &cols, &vals));  in PetscDualSpaceCreateAllDataFromInteriorData()
1626         PetscCall(MatRestoreRow(intMat, j, &ncols, &cols, &vals));  in PetscDualSpaceCreateAllDataFromInteriorData()
2188     sp->allMat = sp->intMat = allMat;  in PetscDualSpaceSetUp_Lagrange()
2213 …implexNodeMat(nodeFamily, dim, sum, Nk, numNodeSkip, &sp->intNodes, &sp->intMat, &lag->intNodeIndi…  in PetscDualSpaceSetUp_Lagrange()
2214         PetscCall(MatGetSize(sp->intMat, &nDofs, NULL));  in PetscDualSpaceSetUp_Lagrange()
2227 …implexNodeMat(nodeFamily, dim, sum, Nk, numNodeSkip, &sp->intNodes, &sp->intMat, &lag->intNodeIndi…  in PetscDualSpaceSetUp_Lagrange()
2228         PetscCall(MatGetSize(sp->intMat, &nDofs, NULL));  in PetscDualSpaceSetUp_Lagrange()
2233         PetscCall(PetscObjectReference((PetscObject)sp->intMat));  in PetscDualSpaceSetUp_Lagrange()
2234         sp->allMat = sp->intMat;  in PetscDualSpaceSetUp_Lagrange()
2246           Mat                 intMat;  in PetscDualSpaceSetUp_Lagrange()  local
2255           PetscCall(PetscDualSpaceGetAllData(trimmedsp, &intNodes, &intMat));  in PetscDualSpaceSetUp_Lagrange()
2260           PetscCall(PetscObjectReference((PetscObject)intMat));  in PetscDualSpaceSetUp_Lagrange()
2267             PetscCall(MatDuplicate(intMat, MAT_COPY_VALUES, &intMatT));  in PetscDualSpaceSetUp_Lagrange()
2268             PetscCall(MatGetSize(intMat, &nRows, &nCols));  in PetscDualSpaceSetUp_Lagrange()
2276               PetscCall(MatGetRow(intMat, row, &nrCols, &rCols, &rVals));  in PetscDualSpaceSetUp_Lagrange()
2287               PetscCall(MatRestoreRow(intMat, row, &nrCols, &rCols, &rVals));  in PetscDualSpaceSetUp_Lagrange()
2291             PetscCall(MatDestroy(&intMat));  in PetscDualSpaceSetUp_Lagrange()
2292             intMat = intMatT;  in PetscDualSpaceSetUp_Lagrange()
2312           sp->intMat = intMat;  in PetscDualSpaceSetUp_Lagrange()
2313           PetscCall(MatGetSize(sp->intMat, &nDofs, NULL));  in PetscDualSpaceSetUp_Lagrange()
2327     Mat                 intMat         = NULL;  in PetscDualSpaceSetUp_Lagrange()  local
2349         PetscCall(MatTensorAltV(intMatTrace, intMatFiber, dim - 1, formDegree, 1, 0, &intMat));  in PetscDualSpaceSetUp_Lagrange()
2378         if (!intMat) {  in PetscDualSpaceSetUp_Lagrange()
2379           intMat         = intMat2;  in PetscDualSpaceSetUp_Lagrange()
2391           PetscCall(MatGetSize(intMat, &nDof, NULL));  in PetscDualSpaceSetUp_Lagrange()
2395 …   PetscCall(MatricesMerge(intMat, intMat2, dim, formDegree, nM, toMerged, toMerged2, &matMerged));  in PetscDualSpaceSetUp_Lagrange()
2399           PetscCall(MatDestroy(&intMat));  in PetscDualSpaceSetUp_Lagrange()
2406           intMat         = matMerged;  in PetscDualSpaceSetUp_Lagrange()
2415             PetscCall(MatPermuteByNodeIdx(intMat, intNodeIndices, &intMatPerm));  in PetscDualSpaceSetUp_Lagrange()
2416             PetscCall(MatDestroy(&intMat));  in PetscDualSpaceSetUp_Lagrange()
2417             intMat = intMatPerm;  in PetscDualSpaceSetUp_Lagrange()
2427     sp->intMat          = intMat;  in PetscDualSpaceSetUp_Lagrange()
2432       if (intMat) PetscCall(MatGetSize(intMat, &nDofs, NULL));  in PetscDualSpaceSetUp_Lagrange()
2442       PetscCall(PetscObjectReference((PetscObject)intMat));  in PetscDualSpaceSetUp_Lagrange()
2443       sp->allMat = intMat;  in PetscDualSpaceSetUp_Lagrange()