6553   PetscInt  size_of_constraint, max_size_of_constraint = 0, max_constraints, temp_constraints;  in PCBDDCConstraintsSetUp()  local
6767       temp_constraints = 0; /* zero the number of constraints I have on this conn comp */  in PCBDDCConstraintsSetUp()
6784         temp_constraints++;  in PCBDDCConstraintsSetUp()
6792 …ptr_to_data = &constraints_data[constraints_data_ptr[total_counts_cc] + temp_constraints * size_of…  in PCBDDCConstraintsSetUp()
6799           temp_constraints++;  in PCBDDCConstraintsSetUp()
6808       valid_constraints = temp_constraints;  in PCBDDCConstraintsSetUp()
6809       if (!pcbddc->use_nnsp_true && temp_constraints) {  in PCBDDCConstraintsSetUp()
6810         if (temp_constraints == 1) { /* just normalize the constraint */  in PCBDDCConstraintsSetUp()
6829             PetscCall(PetscArrayzero(correlation_mat, temp_constraints * temp_constraints));  in PCBDDCConstraintsSetUp()
6833             for (j = 0; j < temp_constraints; j++) {  in PCBDDCConstraintsSetUp()
6834 …k = 0; k < j + 1; k++) PetscCallBLAS("BLASdot", correlation_mat[j * temp_constraints + k] = BLASdo…  in PCBDDCConstraintsSetUp()
6837             PetscCall(PetscBLASIntCast(temp_constraints, &Blas_N));  in PCBDDCConstraintsSetUp()
6838             PetscCall(PetscBLASIntCast(temp_constraints, &Blas_LDA));  in PCBDDCConstraintsSetUp()
6848 …  while (j < temp_constraints && singular_vals[j] / singular_vals[temp_constraints - 1] < tol) j++;  in PCBDDCConstraintsSetUp()
6850             valid_constraints = temp_constraints - j;  in PCBDDCConstraintsSetUp()
6853             PetscCall(PetscBLASIntCast(temp_constraints, &Blas_N));  in PCBDDCConstraintsSetUp()
6854             PetscCall(PetscBLASIntCast(temp_constraints, &Blas_K));  in PCBDDCConstraintsSetUp()
6856             PetscCall(PetscBLASIntCast(temp_constraints, &Blas_LDB));  in PCBDDCConstraintsSetUp()
6858             if (j < temp_constraints) {  in PCBDDCConstraintsSetUp()
6860 …   for (k = j; k < temp_constraints; k++) singular_vals[k] = 1.0 / PetscSqrtReal(singular_vals[k]);  in PCBDDCConstraintsSetUp()
6864               for (k = 0; k < temp_constraints - j; k++) {  in PCBDDCConstraintsSetUp()
6865 …ta[k * size_of_constraint + ii] = singular_vals[temp_constraints - 1 - k] * temp_basis[(temp_const…  in PCBDDCConstraintsSetUp()
6871             PetscCall(PetscBLASIntCast(temp_constraints, &Blas_N));  in PCBDDCConstraintsSetUp()
6882             k = temp_constraints;  in PCBDDCConstraintsSetUp()
6887             total_counts      = total_counts - temp_constraints + valid_constraints;  in PCBDDCConstraintsSetUp()