4070   PetscInt n_constraints, n_R, old_size;  in PCBDDCSetUpLocalWorkVectors()  local
4074   n_R           = pcis->n - pcbddc->n_vertices;  in PCBDDCSetUpLocalWorkVectors()
4080   if (n_R != old_size) {  in PCBDDCSetUpLocalWorkVectors()
4084     PetscCall(VecSetSizes(pcbddc->vec1_R, PETSC_DECIDE, n_R));  in PCBDDCSetUpLocalWorkVectors()
4299   PetscInt        i, n_R, n_D, n_B;  in PCBDDCSetUpCorrection()  local
4319   n_R           = pcis->n - n_vertices;  in PCBDDCSetUpCorrection()
4403     PetscCall(PetscMalloc1(n_R * n_eff_vertices, &R_eff_V_J));  in PCBDDCSetUpCorrection()
4404     PetscCall(PetscMalloc1(n_R * n_eff_constraints, &R_eff_C_J));  in PCBDDCSetUpCorrection()
4407     for (PetscInt i = 0; i < n_R * n_eff_vertices; i++) R_eff_V_J[i] = -1;  in PCBDDCSetUpCorrection()
4408     for (PetscInt i = 0; i < n_R * n_eff_constraints; i++) R_eff_C_J[i] = -1;  in PCBDDCSetUpCorrection()
4413     for (PetscInt i = 0; i < n_R; i++) {  in PCBDDCSetUpCorrection()
4463   lda_rhs                = n_R;  in PCBDDCSetUpCorrection()
4491   if (need_benign_correction && lda_rhs != n_R && F) {  in PCBDDCSetUpCorrection()
4500   PetscCall(ISCreateStride(PETSC_COMM_SELF, n_R, 0, 1, &is_R));  in PCBDDCSetUpCorrection()
4535       if (lda_rhs != n_R) {  in PCBDDCSetUpCorrection()
4568         if (lda_rhs != n_R) {  in PCBDDCSetUpCorrection()
4625       if (lda_rhs != n_R) {  in PCBDDCSetUpCorrection()
4734     if (n_R) {  in PCBDDCSetUpCorrection()
4755         if (lda_rhs == n_R && !multi_element) {  in PCBDDCSetUpCorrection()
4823         if (lda_rhs != n_R) {  in PCBDDCSetUpCorrection()
4849           if (lda_rhs != n_R) {  in PCBDDCSetUpCorrection()
4872           if (lda_rhs != n_R) {  in PCBDDCSetUpCorrection()
4956       if (lda_rhs != n_R) {  in PCBDDCSetUpCorrection()
4981             for (PetscInt k = 0; k < nz; k++) sums[j] += marr[idxs_zero[k] + j * n_R];  in PCBDDCSetUpCorrection()
5051         if (lda_rhs != n_R) {  in PCBDDCSetUpCorrection()
5145       PetscCall(MatCreateSeqDense(PETSC_COMM_SELF, n_R, n_vertices, NULL, &B_V));  in PCBDDCSetUpCorrection()
5147     if (n_vertices && n_R) {  in PCBDDCSetUpCorrection()
5160         for (j = xadj[i]; j < xadj[i + 1]; j++) marray[i * n_R + adjncy[j]] -= av[j];  in PCBDDCSetUpCorrection()
5168     PetscCall(PetscMalloc1(n_R * pcbddc->local_primal_size, &work));  in PCBDDCSetUpCorrection()
5172         PetscCall(VecPlaceArray(pcbddc->vec1_R, marray + i * n_R));  in PCBDDCSetUpCorrection()
5173         PetscCall(VecPlaceArray(pcbddc->vec2_R, work + i * n_R));  in PCBDDCSetUpCorrection()
5184         PetscCall(VecPlaceArray(pcbddc->vec1_R, marray + (i - n_vertices) * n_R));  in PCBDDCSetUpCorrection()
5185         PetscCall(VecPlaceArray(pcbddc->vec2_R, work + i * n_R));  in PCBDDCSetUpCorrection()
5199       PetscCall(VecPlaceArray(pcbddc->vec1_R, work + i * n_R));  in PCBDDCSetUpCorrection()
5714   PetscInt        n_vertices, i, j, n_R, n_D, n_B;  in PCBDDCSetUpLocalScatters()  local
5745     for (i = 0, n_R = 0; i < pcis->n; i++) {  in PCBDDCSetUpLocalScatters()
5746       if (!PetscBTLookup(bitmask, i)) idx_R_local[n_R++] = i;  in PCBDDCSetUpLocalScatters()
5752     PetscCall(ISGetLocalSize(reuse_solver->is_R, &n_R));  in PCBDDCSetUpLocalScatters()
5776       for (i = 0; i < n_R / bs; i++) {  in PCBDDCSetUpLocalScatters()
5788       for (i = 0; i < n_R / vbs; i++) idx_R_local[i] = idx_R_local[vbs * i] / vbs;  in PCBDDCSetUpLocalScatters()
5791 …PetscCall(ISCreateBlock(PETSC_COMM_SELF, vbs, n_R / vbs, idx_R_local, PETSC_COPY_VALUES, &pcbddc->…  in PCBDDCSetUpLocalScatters()
5810 …_FMT ", constraints = %" PetscInt_FMT ", local_primal_size = %" PetscInt_FMT "\n", n_R, n_vertices,  in PCBDDCSetUpLocalScatters()
5826     for (i = 0, j = 0; i < n_R; i++) {  in PCBDDCSetUpLocalScatters()
5842       for (i = 0, j = 0; i < n_R; i++) {  in PCBDDCSetUpLocalScatters()
5943   PetscInt     n_D, n_R;  in PCBDDCSetUpLocalSolvers()  local
6085     PetscCall(ISGetSize(pcbddc->is_R_local, &n_R));  in PCBDDCSetUpLocalSolvers()
6091       if (nn_R != n_R) { /* old ksp is not reusable, so reset it */  in PCBDDCSetUpLocalSolvers()
6200     if (!n_R) {  in PCBDDCSetUpLocalSolvers()