6537   PetscInt i, j, k, total_counts, total_counts_cc, cum;  in PCBDDCConstraintsSetUp()  local
6645     total_counts = n_ISForFaces + n_ISForEdges;  in PCBDDCConstraintsSetUp()
6646     total_counts *= max_constraints;  in PCBDDCConstraintsSetUp()
6647     total_counts += n_vertices;  in PCBDDCConstraintsSetUp()
6648     PetscCall(PetscBTCreate(total_counts, &change_basis));  in PCBDDCConstraintsSetUp()
6650     total_counts           = 0;  in PCBDDCConstraintsSetUp()
6660       total_counts += j;  in PCBDDCConstraintsSetUp()
6663 …l(PetscMalloc3(total_counts * max_constraints + n_vertices, &constraints_data, total_counts + n_ve…  in PCBDDCConstraintsSetUp()
6737     total_counts            = 0;  in PCBDDCConstraintsSetUp()
6745         constraints_n[total_counts]            = 1;  in PCBDDCConstraintsSetUp()
6746         constraints_data[total_counts]         = 1.0;  in PCBDDCConstraintsSetUp()
6747         constraints_idxs_ptr[total_counts + 1] = constraints_idxs_ptr[total_counts] + 1;  in PCBDDCConstraintsSetUp()
6748         constraints_data_ptr[total_counts + 1] = constraints_data_ptr[total_counts] + 1;  in PCBDDCConstraintsSetUp()
6749         total_counts++;  in PCBDDCConstraintsSetUp()
6755     total_counts_cc = total_counts;  in PCBDDCConstraintsSetUp()
6785         total_counts++;  in PCBDDCConstraintsSetUp()
6800           total_counts++;  in PCBDDCConstraintsSetUp()
6849             total_counts      = total_counts - j;  in PCBDDCConstraintsSetUp()
6887             total_counts      = total_counts - temp_constraints + valid_constraints;  in PCBDDCConstraintsSetUp()
6921     total_counts = 0;  in PCBDDCConstraintsSetUp()
6927       total_counts += pcbddc->adaptive_constraints_n[i];  in PCBDDCConstraintsSetUp()
6953   pcbddc->local_primal_size = total_counts;  in PCBDDCConstraintsSetUp()
7003   total_counts = total_primal_vertices;  in PCBDDCConstraintsSetUp()
7012 … pcbddc->primal_indices_local_idxs[total_counts++] = constraints_idxs[constraints_idxs_ptr[i] + k];  in PCBDDCConstraintsSetUp()
7025   total_counts = total_primal_vertices;  in PCBDDCConstraintsSetUp()
7033         PetscInt     row = total_counts + k;  in PCBDDCConstraintsSetUp()
7039       total_counts += constraints_n[i];  in PCBDDCConstraintsSetUp()
7181     for (total_counts = n_vertices; total_counts < total_counts_cc; total_counts++) {  in PCBDDCConstraintsSetUp()
7182 …  size_of_constraint = constraints_idxs_ptr[total_counts + 1] - constraints_idxs_ptr[total_counts];  in PCBDDCConstraintsSetUp()
7183       if (PetscBTLookup(change_basis, total_counts)) {  in PCBDDCConstraintsSetUp()
7185         primal_dofs = constraints_n[total_counts];  in PCBDDCConstraintsSetUp()
7188 …%" PetscInt_FMT " need a change of basis (size %" PetscInt_FMT ")\n", total_counts, primal_dofs, s…  in PCBDDCConstraintsSetUp()
7190         if (PetscBTLookup(qr_needed_idx, total_counts)) { /* QR */  in PCBDDCConstraintsSetUp()
7193 …scCall(PetscArraycpy(dbg_work, &constraints_data[constraints_data_ptr[total_counts]], size_of_cons…  in PCBDDCConstraintsSetUp()
7195 …PetscCall(PetscArraycpy(qr_basis, &constraints_data[constraints_data_ptr[total_counts]], size_of_c…  in PCBDDCConstraintsSetUp()
7238 …LDA, trs_rhs, &Blas_LDB, &zero, constraints_data + constraints_data_ptr[total_counts], &Blas_LDC));  in PCBDDCConstraintsSetUp()
7240 …PetscCall(PetscArraycpy(qr_basis, &constraints_data[constraints_data_ptr[total_counts]], size_of_c…  in PCBDDCConstraintsSetUp()
7243           start_rows = &constraints_idxs[constraints_idxs_ptr[total_counts]];  in PCBDDCConstraintsSetUp()
7247             start_cols = &constraints_idxs[constraints_idxs_ptr[total_counts] + j];  in PCBDDCConstraintsSetUp()
7252 …       if (!PetscBTLookup(is_primal, constraints_idxs_B[constraints_idxs_ptr[total_counts] + k])) {  in PCBDDCConstraintsSetUp()
7254               start_cols = &constraints_idxs[constraints_idxs_ptr[total_counts] + k];  in PCBDDCConstraintsSetUp()
7300 …straints_data + constraints_data_ptr[total_counts], &Blas_one, constraints_data + constraints_data…  in PCBDDCConstraintsSetUp()
7302             PetscInt row_B = constraints_idxs_B[constraints_idxs_ptr[total_counts] + j];  in PCBDDCConstraintsSetUp()
7303             row            = constraints_idxs[constraints_idxs_ptr[total_counts] + j];  in PCBDDCConstraintsSetUp()
7305               col = constraints_idxs[constraints_idxs_ptr[total_counts]];  in PCBDDCConstraintsSetUp()
7307 …lChangeOfBasisMatrix, row, col, constraints_data[constraints_data_ptr[total_counts] + j] / norm, I…  in PCBDDCConstraintsSetUp()
7310                 col = constraints_idxs[constraints_idxs_ptr[total_counts] + k];  in PCBDDCConstraintsSetUp()
7312 …= -constraints_data[constraints_data_ptr[total_counts] + k] / constraints_data[constraints_data_pt…  in PCBDDCConstraintsSetUp()
7314                   val = constraints_data[constraints_data_ptr[total_counts]] / norm;  in PCBDDCConstraintsSetUp()
7323 …nt_FMT " does not need a change of basis (size %" PetscInt_FMT ")\n", total_counts, size_of_constr…  in PCBDDCConstraintsSetUp()