842     const PetscScalar *rS_data;  in PCBDDCSubSchursSetUp()  local
1450     if (!multi_element) PetscCall(MatDenseGetArrayRead(S_all, &rS_data));  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 …for (PetscInt j = 0; j < subset_size; j++) work[k * subset_size + j] = rS_data[cum2 + k * size_sch…  in PCBDDCSubSchursSetUp()
1568     if (!multi_element) PetscCall(MatDenseRestoreArrayRead(S_all, &rS_data));  in PCBDDCSubSchursSetUp()
1790         rS_data    = NULL;  in PCBDDCSubSchursSetUp()
1791         if (S_all_inv && !multi_element) PetscCall(MatDenseGetArrayRead(S_all_inv, &rS_data));  in PCBDDCSubSchursSetUp()
1813           } else if (rS_data) {  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()
1841 …PetscCall(MatCreateSeqDense(PETSC_COMM_SELF, subset_size, subset_size, (rS_data || multi_element) …  in PCBDDCSubSchursSetUp()
1856 …if (rS_data || multi_element) PetscCall(PetscArraycpy(SEjinv_arr, work, subset_size * subset_size)…  in PCBDDCSubSchursSetUp()
1863           if (!multi_element) PetscCall(MatDenseRestoreArrayRead(S_all_inv, &rS_data));  in PCBDDCSubSchursSetUp()