96   Nest_Dense        *contents;  in MatProductNumeric_Nest_Dense()  local
113   contents = (Nest_Dense *)C->product->data;  in MatProductNumeric_Nest_Dense()
114   PetscCheck(contents, PetscObjectComm((PetscObject)C), PETSC_ERR_PLIB, "Product data empty");  in MatProductNumeric_Nest_Dense()
125 …ectComm((PetscObject)A), contents->dm[i + 1] - contents->dm[i], PETSC_DECIDE, M, N, PetscSafePoint…  in MatProductNumeric_Nest_Dense()
130 …((PetscObject)A), contents->dn[j + 1] - contents->dn[j], PETSC_DECIDE, M, N, PetscSafePointerPlusO…  in MatProductNumeric_Nest_Dense()
134       workC             = contents->workC[i * nc + j];  in MatProductNumeric_Nest_Dense()
142       PetscCall(MatAXPY(viewC, 1.0, contents->workC[i * nc + j], SAME_NONZERO_PATTERN));  in MatProductNumeric_Nest_Dense()
157   Nest_Dense *contents = *(Nest_Dense **)ctx;  in MatNest_DenseDestroy()  local
161   PetscCall(PetscFree(contents->tarray));  in MatNest_DenseDestroy()
162   for (i = 0; i < contents->k; i++) PetscCall(MatDestroy(contents->workC + i));  in MatNest_DenseDestroy()
163   PetscCall(PetscFree3(contents->dm, contents->dn, contents->workC));  in MatNest_DenseDestroy()
164   PetscCall(PetscFree(contents));  in MatNest_DenseDestroy()
174   Nest_Dense        *contents = NULL;  in MatProductSymbolic_Nest_Dense()  local
205   PetscCall(PetscNew(&contents));  in MatProductSymbolic_Nest_Dense()
206   C->product->data    = contents;  in MatProductSymbolic_Nest_Dense()
208   PetscCall(PetscCalloc3(nr + 1, &contents->dm, nc + 1, &contents->dn, nr * nc, &contents->workC));  in MatProductSymbolic_Nest_Dense()
209   contents->k = nr * nc;  in MatProductSymbolic_Nest_Dense()
211     PetscCall(ISGetLocalSize(bA->isglobal.row[i], contents->dm + i + 1));  in MatProductSymbolic_Nest_Dense()
212     maxm = PetscMax(maxm, contents->dm[i + 1]);  in MatProductSymbolic_Nest_Dense()
213     contents->dm[i + 1] += contents->dm[i];  in MatProductSymbolic_Nest_Dense()
216     PetscCall(ISGetLocalSize(bA->isglobal.col[i], contents->dn + i + 1));  in MatProductSymbolic_Nest_Dense()
217     contents->dn[i + 1] += contents->dn[i];  in MatProductSymbolic_Nest_Dense()
219   PetscCall(PetscMalloc1(maxm * N, &contents->tarray));  in MatProductSymbolic_Nest_Dense()
226 …((PetscObject)A), contents->dn[j + 1] - contents->dn[j], PETSC_DECIDE, M, N, PetscSafePointerPlusO…  in MatProductSymbolic_Nest_Dense()
232       PetscCall(MatProductCreate(bA->m[i][j], viewB, NULL, &contents->workC[i * nc + j]));  in MatProductSymbolic_Nest_Dense()
233       workC = contents->workC[i * nc + j];  in MatProductSymbolic_Nest_Dense()
241       PetscCall(MatSeqDenseSetPreallocation(workC, contents->tarray));  in MatProductSymbolic_Nest_Dense()
242       PetscCall(MatMPIDenseSetPreallocation(workC, contents->tarray));  in MatProductSymbolic_Nest_Dense()