16 …actSubMatrix_Private(Mat pmat, Mat *submat, MatReuse reuse, PetscInt component, PetscInt blocksize)  in PCHMGExtractSubMatrix_Private()  argument
24 …ent < blocksize, comm, PETSC_ERR_ARG_INCOMP, "Component %" PetscInt_FMT " should be less than bloc…  in PCHMGExtractSubMatrix_Private()
26 …) % blocksize == 0, comm, PETSC_ERR_ARG_INCOMP, "Block size %" PetscInt_FMT " is inconsistent for …  in PCHMGExtractSubMatrix_Private()
27 …PetscCall(ISCreateStride(comm, (rend - rstart) / blocksize, rstart + component, blocksize, &isrow)…  in PCHMGExtractSubMatrix_Private()
33 …tic PetscErrorCode PCHMGExpandInterpolation_Private(Mat subinterp, Mat *interp, PetscInt blocksize)  in PCHMGExpandInterpolation_Private()  argument
45   rowsize    = subrowsize * blocksize;  in PCHMGExpandInterpolation_Private()
49   colsize    = subcolsize * blocksize;  in PCHMGExpandInterpolation_Private()
60     for (i = 0; i < blocksize; i++) {  in PCHMGExpandInterpolation_Private()
61       d_nnz[(subrow - subrstart) * blocksize + i] = dnz;  in PCHMGExpandInterpolation_Private()
62       o_nnz[(subrow - subrstart) * blocksize + i] = onz;  in PCHMGExpandInterpolation_Private()
77     for (i = 0; i < blocksize; i++) {  in PCHMGExpandInterpolation_Private()
78       row = subrow * blocksize + i;  in PCHMGExpandInterpolation_Private()
79       for (j = 0; j < nz; j++) indices[j] = idx[j] * blocksize + i;  in PCHMGExpandInterpolation_Private()
99   PetscInt         blocksize;  in PCSetUp_HMG()  local
137   PetscCall(MatGetBlockSize(PA, &blocksize));  in PCSetUp_HMG()
138   if (blocksize <= 1) hmg->subcoarsening = PETSC_FALSE;  in PCSetUp_HMG()
141 …tscCall(PCHMGExtractSubMatrix_Private(PA, &submat, MAT_INITIAL_MATRIX, hmg->component, blocksize));  in PCSetUp_HMG()
168         PetscCall(MatCreateMAIJ(interpolations[level - 1], blocksize, &P));  in PCSetUp_HMG()
172         PetscCall(PCHMGExpandInterpolation_Private(interpolations[level - 1], &P, blocksize));  in PCSetUp_HMG()