4070   PetscInt n_constraints, n_R, old_size;  in PCBDDCSetUpLocalWorkVectors()  local
4073   n_constraints = pcbddc->local_primal_size - pcbddc->benign_n - pcbddc->n_vertices;  in PCBDDCSetUpLocalWorkVectors()
4100   if (n_constraints && n_constraints != old_size) {  in PCBDDCSetUpLocalWorkVectors()
4103     PetscCall(VecSetSizes(pcbddc->vec1_C, PETSC_DECIDE, n_constraints));  in PCBDDCSetUpLocalWorkVectors()
4297   PetscInt        lda_rhs, n_vertices, n_constraints, *p0_lidx_I;  in PCBDDCSetUpCorrection()  local
4316   n_constraints = pcbddc->local_primal_size - pcbddc->benign_n - n_vertices;  in PCBDDCSetUpCorrection()
4347     PetscCall(PetscMalloc1(n_constraints, &C_to_eff_C));  in PCBDDCSetUpCorrection()
4354     for (PetscInt i = 0; i < n_constraints; i++) {  in PCBDDCSetUpCorrection()
4362     PetscCall(PetscMalloc1(n_vertices + n_constraints, &nnz));  in PCBDDCSetUpCorrection()
4368     for (PetscInt i = 0; i < n_constraints; i++) {  in PCBDDCSetUpCorrection()
4395     for (PetscInt i = 0; i < n_constraints; i++) {  in PCBDDCSetUpCorrection()
4436     PetscCall(PetscMalloc1(n_constraints, &idxp));  in PCBDDCSetUpCorrection()
4447     PetscCall(ISCreateGeneral(PETSC_COMM_SELF, n_constraints, idxp, PETSC_OWN_POINTER, &is_C_perm));  in PCBDDCSetUpCorrection()
4453     n_eff_constraints = n_constraints;  in PCBDDCSetUpCorrection()
4502   PetscCall(ISCreateStride(PETSC_COMM_SELF, n_constraints, n_vertices, 1, &is_C));  in PCBDDCSetUpCorrection()
4507   if (n_constraints) {  in PCBDDCSetUpCorrection()
4521       for (i = 0; i < n_constraints; i++) {  in PCBDDCSetUpCorrection()
4543 … PetscCall(MatCreateSeqAIJWithArrays(PETSC_COMM_SELF, n_constraints, lda_rhs, ii, jj, aa, &tC_CR));  in PCBDDCSetUpCorrection()
4609         PetscCall(MatCreateSeqAIJFromDenseExpand(local_auxmat2_R, n_constraints, R_eff_C_J, &T));  in PCBDDCSetUpCorrection()
4612 …    PetscCall(MatCreateSeqAIJFromDenseExpand(pcbddc->local_auxmat2, n_constraints, B_eff_C_J, &T));  in PCBDDCSetUpCorrection()
4621         PetscCall(MatCreateSeqAIJFromDenseExpand(local_auxmat2_R, n_constraints, R_eff_C_J, &T));  in PCBDDCSetUpCorrection()
4659 …PetscCall(MatSetValuesSubMat(*coarse_submat, S_CC, n_constraints, idx_C, n_constraints, idx_C, INS…  in PCBDDCSetUpCorrection()
4907       if (n_constraints) {  in PCBDDCSetUpCorrection()
4937 …PetscCall(MatSetValuesSubMat(*coarse_submat, S_CV, n_constraints, idx_C, n_vertices, idx_V, INSERT…  in PCBDDCSetUpCorrection()
5037   if (n_constraints) {  in PCBDDCSetUpCorrection()
5060 …PetscCall(MatSetValuesSubMat(*coarse_submat, S_VC, n_vertices, idx_V, n_constraints, idx_C, INSERT…  in PCBDDCSetUpCorrection()
5067 …ix(pcbddc->coarse_phi_B, PETSC_DECIDE, PETSC_DECIDE, n_vertices, n_vertices + n_constraints, &B2));  in PCBDDCSetUpCorrection()
5071 …ix(pcbddc->coarse_phi_D, PETSC_DECIDE, PETSC_DECIDE, n_vertices, n_vertices + n_constraints, &B2));  in PCBDDCSetUpCorrection()
5074 …for (i = 0; i < n_constraints; i++) PetscCall(MatSetValues(B2, pcbddc->benign_n, p0_lidx_I, 1, &i,…  in PCBDDCSetUpCorrection()
5126     if (n_constraints) {  in PCBDDCSetUpCorrection()
5183       for (i = n_vertices; i < n_constraints + n_vertices; i++) {  in PCBDDCSetUpCorrection()
5248   if (n_constraints) PetscCall(MatDestroy(&C_CR));  in PCBDDCSetUpCorrection()