7050     PetscInt dual_dofs, primal_dofs;  in PCBDDCConstraintsSetUp()  local
7185         primal_dofs = constraints_n[total_counts];  in PCBDDCConstraintsSetUp()
7186         dual_dofs   = size_of_constraint - primal_dofs;  in PCBDDCConstraintsSetUp()
7188 …MT " need a change of basis (size %" PetscInt_FMT ")\n", total_counts, primal_dofs, size_of_constr…  in PCBDDCConstraintsSetUp()
7193 …bg_work, &constraints_data[constraints_data_ptr[total_counts]], size_of_constraint * primal_dofs));  in PCBDDCConstraintsSetUp()
7195 …r_basis, &constraints_data[constraints_data_ptr[total_counts]], size_of_constraint * primal_dofs));  in PCBDDCConstraintsSetUp()
7199           PetscCall(PetscBLASIntCast(primal_dofs, &Blas_N));  in PCBDDCConstraintsSetUp()
7207           PetscCall(PetscArrayzero(trs_rhs, primal_dofs * primal_dofs));  in PCBDDCConstraintsSetUp()
7208           for (j = 0; j < primal_dofs; j++) trs_rhs[j * (primal_dofs + 1)] = 1.0;  in PCBDDCConstraintsSetUp()
7209           PetscCall(PetscBLASIntCast(primal_dofs, &Blas_N));  in PCBDDCConstraintsSetUp()
7210           PetscCall(PetscBLASIntCast(primal_dofs, &Blas_NRHS));  in PCBDDCConstraintsSetUp()
7212           PetscCall(PetscBLASIntCast(primal_dofs, &Blas_LDB));  in PCBDDCConstraintsSetUp()
7221           PetscCall(PetscBLASIntCast(primal_dofs, &Blas_K));  in PCBDDCConstraintsSetUp()
7232           PetscCall(PetscBLASIntCast(primal_dofs, &Blas_N));  in PCBDDCConstraintsSetUp()
7233           PetscCall(PetscBLASIntCast(primal_dofs, &Blas_K));  in PCBDDCConstraintsSetUp()
7235           PetscCall(PetscBLASIntCast(primal_dofs, &Blas_LDB));  in PCBDDCConstraintsSetUp()
7240 …r_basis, &constraints_data[constraints_data_ptr[total_counts]], size_of_constraint * primal_dofs));  in PCBDDCConstraintsSetUp()
7245           for (j = 0; j < primal_dofs; j++) {  in PCBDDCConstraintsSetUp()
7253               start_vals = &qr_basis[(primal_dofs + j) * size_of_constraint];  in PCBDDCConstraintsSetUp()
7264             PetscCall(PetscBLASIntCast(primal_dofs, &Blas_M));  in PCBDDCConstraintsSetUp()
7269             PetscCall(PetscBLASIntCast(primal_dofs, &Blas_LDC));  in PCBDDCConstraintsSetUp()
7271 …k, &Blas_LDA, qr_basis, &Blas_LDB, &zero, &dbg_work[size_of_constraint * primal_dofs], &Blas_LDC));  in PCBDDCConstraintsSetUp()
7274               for (ii = 0; ii < primal_dofs; ii++) {  in PCBDDCConstraintsSetUp()
7275 …if (ii != jj && PetscAbsScalar(dbg_work[size_of_constraint * primal_dofs + jj * primal_dofs + ii])…  in PCBDDCConstraintsSetUp()
7276 …if (ii == jj && PetscAbsScalar(dbg_work[size_of_constraint * primal_dofs + jj * primal_dofs + ii] …  in PCBDDCConstraintsSetUp()
7282                 for (ii = 0; ii < primal_dofs; ii++) {  in PCBDDCConstraintsSetUp()
7283 …if (ii != jj && PetscAbsScalar(dbg_work[size_of_constraint * primal_dofs + jj * primal_dofs + ii])…  in PCBDDCConstraintsSetUp()
7284 …jj, ii, (double)PetscAbsScalar(dbg_work[size_of_constraint * primal_dofs + jj * primal_dofs + ii])…  in PCBDDCConstraintsSetUp()
7286 …if (ii == jj && PetscAbsScalar(dbg_work[size_of_constraint * primal_dofs + jj * primal_dofs + ii] …  in PCBDDCConstraintsSetUp()
7287 …jj, ii, (double)PetscAbsScalar(dbg_work[size_of_constraint * primal_dofs + jj * primal_dofs + ii])…  in PCBDDCConstraintsSetUp()