20     PetscInt n_I, size_schur;  in PCBDDCReuseSolversBenignAdapt()  local
23     PetscCall(MatGetSize(ctx->benign_csAIB, &size_schur, NULL));  in PCBDDCReuseSolversBenignAdapt()
25     n_I = n_I - size_schur;  in PCBDDCReuseSolversBenignAdapt()
93     PetscInt n_I, size_schur;  in PCBDDCReuseSolversBenignAdapt()  local
96     PetscCall(MatGetSize(ctx->benign_csAIB, &size_schur, NULL));  in PCBDDCReuseSolversBenignAdapt()
98     n_I = n_I - size_schur;  in PCBDDCReuseSolversBenignAdapt()
843     PetscInt           n, n_I, size_schur, size_active_schur, cum, cum2;  in PCBDDCSubSchursSetUp()  local
915     size_schur = cum - n_I;  in PCBDDCSubSchursSetUp()
933     PetscCall(ISCreateStride(PETSC_COMM_SELF, size_schur, n_I, 1, &is_schur));  in PCBDDCSubSchursSetUp()
939       PetscCall(PetscMalloc1(n + size_schur, &idx_sub));  in PCBDDCSubSchursSetUp()
946         for (PetscInt j = n_I; j < n_I + size_schur; j++)  in PCBDDCSubSchursSetUp()
975       PetscCall(MatCreateSeqDense(PETSC_COMM_SELF, size_schur, benign_n, NULL, &cs_AIB_mat));  in PCBDDCSubSchursSetUp()
994         for (j = 0; j < size_schur; j++) {  in PCBDDCSubSchursSetUp()
996 …cs_AIB[i * size_schur + j] = (PetscRealPart(array[j + n_I]) / nz + PETSC_i * (PetscImaginaryPart(a…  in PCBDDCSubSchursSetUp()
998           cs_AIB[i * size_schur + j] = array[j + n_I] / nz;  in PCBDDCSubSchursSetUp()
1077         PetscCall(ISCreateStride(PETSC_COMM_SELF, size_schur, n_I, 1, &is));  in PCBDDCSubSchursSetUp()
1178             PetscCall(VecCreateSeqWithArray(PETSC_COMM_SELF, 1, size_schur, array + n_I, &vv));  in PCBDDCSubSchursSetUp()
1186           PetscCall(PetscBLASIntCast(size_schur, &B_n));  in PCBDDCSubSchursSetUp()
1192 …allBLAS("BLASsyrk", BLASsyrk_("L", "N", &B_n, &B_k, &sum, cs_AIB + i * size_schur, &B_n, &one, S_d…  in PCBDDCSubSchursSetUp()
1193 …allBLAS("BLASsyrk", BLASsyrk_("L", "N", &B_n, &B_k, &sum, cs_AIB + i * size_schur, &B_n, &one, S3_…  in PCBDDCSubSchursSetUp()
1196 …k", BLASsyr2k_("L", "N", &B_n, &B_k, &sum, cs_AIB + k * size_schur, &B_n, cs_AIB + i * size_schur,…  in PCBDDCSubSchursSetUp()
1197 …k", BLASsyr2k_("L", "N", &B_n, &B_k, &sum, cs_AIB + k * size_schur, &B_n, cs_AIB + i * size_schur,…  in PCBDDCSubSchursSetUp()
1201 …BLASsyr2k_("L", "N", &B_n, &B_k, &sum, array + n_I, &B_n, cs_AIB + i * size_schur, &B_n, &one, S_d…  in PCBDDCSubSchursSetUp()
1202 …BLASsyr2k_("L", "N", &B_n, &B_k, &sum, array + n_I, &B_n, cs_AIB + i * size_schur, &B_n, &one, S2_…  in PCBDDCSubSchursSetUp()
1218           for (k = 0; k < size_schur; k++) {  in PCBDDCSubSchursSetUp()
1219 …for (j = k; j < size_schur; j++) S_data[j * size_schur + k] = PetscConj(S_data[k * size_schur + j]…  in PCBDDCSubSchursSetUp()
1262       PetscCall(PetscCalloc1(size_schur, &nnz));  in PCBDDCSubSchursSetUp()
1269       PetscCall(MatCreateSeqAIJ(PETSC_COMM_SELF, size_schur, size_schur, 0, nnz, &S_all));  in PCBDDCSubSchursSetUp()
1403           dual = size_schur - n_v;  in PCBDDCSubSchursSetUp()
1433 …PetscCall(VecCreateSeqWithArray(PETSC_COMM_SELF, 1, size_schur, array, &msolv_ctx->benign_dummy_sc…  in PCBDDCSubSchursSetUp()
1481               work[k * subset_size + j] = rS_data[cum2 + k * size_schur + j];  in PCBDDCSubSchursSetUp()
1482               work[j * subset_size + k] = PetscConj(rS_data[cum2 + k * size_schur + j]);  in PCBDDCSubSchursSetUp()
1487 …tscInt j = 0; j < subset_size; j++) work[k * subset_size + j] = rS_data[cum2 + k * size_schur + j];  in PCBDDCSubSchursSetUp()
1563       cum2 += subset_size * (size_schur + 1);  in PCBDDCSubSchursSetUp()
1590 …if (sub_schurs->n_subs == 1 && size_schur == size_active_schur && deluxe) { /* we already computed…  in PCBDDCSubSchursSetUp()
1592         PetscCall(PetscArraycpy(SEjinv_arr, work, size_schur * size_schur));  in PCBDDCSubSchursSetUp()
1621 …PetscCall(PetscArraycpy(tdata + i * (news + 1), data + i * (size_schur + 1), size_active_schur - i…  in PCBDDCSubSchursSetUp()
1622 …tscArraycpy(tdata + i * (news + 1) + size_active_schur - i, data + i * size_schur + size_active_sc…  in PCBDDCSubSchursSetUp()
1626 …      PetscCall(PetscArraycpy(tdata + k * (news + 1), data + (k + nd) * (size_schur + 1), nv - i));  in PCBDDCSubSchursSetUp()
1639 …r_factor[cum + i] = PetscSqrtReal(PetscRealPart(data[(i + size_active_schur) * (size_schur + 1)]));  in PCBDDCSubSchursSetUp()
1647 …for (i = 0; i < size_active_schur; i++) PetscCall(PetscArraycpy(data + i * size_schur, tdata + i *…  in PCBDDCSubSchursSetUp()
1655 …          for (i = 0; i < nd; i++) aux[i] = 1.0 / data[(i + size_active_schur) * (size_schur + 1)];  in PCBDDCSubSchursSetUp()
1657               PetscCall(MatDenseSetLDA(M, size_schur));  in PCBDDCSubSchursSetUp()
1662 …PetscCall(MatCreateSeqDense(PETSC_COMM_SELF, size_schur, nd, data + size_active_schur * size_schur…  in PCBDDCSubSchursSetUp()
1665 …       PetscCall(MatCreateSeqDense(PETSC_COMM_SELF, nd, size_schur, data + size_active_schur, &M));  in PCBDDCSubSchursSetUp()
1666               PetscCall(MatDenseSetLDA(M, size_schur));  in PCBDDCSubSchursSetUp()
1669               for (i = 0; i < nd; i++) data[(i + size_active_schur) * (size_schur + 1)] = aux[i];  in PCBDDCSubSchursSetUp()
1716             PetscCall(PetscBLASIntCast(size_schur, &B_N));  in PCBDDCSubSchursSetUp()
1734             PetscCall(PetscLogFlops(1.0 * size_schur * size_schur * size_schur));  in PCBDDCSubSchursSetUp()
1819                     work[k * subset_size + j] = rS_data[cum2 + k * size_schur + j];  in PCBDDCSubSchursSetUp()
1825 …  for (j = k; j < subset_size; j++) work[k * subset_size + j] = rS_data[cum2 + k * size_schur + j];  in PCBDDCSubSchursSetUp()
1831 …  for (j = 0; j < subset_size; j++) work[k * subset_size + j] = rS_data[cum2 + k * size_schur + j];  in PCBDDCSubSchursSetUp()
1858           cum2 += subset_size * (size_schur + 1);  in PCBDDCSubSchursSetUp()
1885         PetscCall(PetscArrayzero(data, size_schur * size_schur));  in PCBDDCSubSchursSetUp()
1890 …  PetscCall(PetscArraycpy(data + i * (size_schur + 1), schur_factor + cum, size_active_schur - i));  in PCBDDCSubSchursSetUp()
1894           PetscCall(PetscArraycpy(data, schur_factor, size_schur * size_schur));  in PCBDDCSubSchursSetUp()
1898 … for (i = 0; i < nd; i++) data[(i + size_active_schur) * (size_schur + 1)] = schur_factor[cum + i];  in PCBDDCSubSchursSetUp()
1902         for (i = size_active_schur + nd; i < size_schur; i++) data[i * (size_schur + 1)] = infty;  in PCBDDCSubSchursSetUp()
1915 …active_schur; i++) PetscCall(PetscArrayzero(data + i * size_schur + size_active_schur, size_schur …  in PCBDDCSubSchursSetUp()
1916       for (i = size_active_schur + nd; i < size_schur; i++) {  in PCBDDCSubSchursSetUp()
1917 …PetscCall(PetscArrayzero(data + i * size_schur + size_active_schur, size_schur - size_active_schur…  in PCBDDCSubSchursSetUp()
1918         data[i * (size_schur + 1)] = infty;  in PCBDDCSubSchursSetUp()