Lines Matching refs:cum
718 PetscInt cum, i, j, k; in PCBDDCSubSchursSetUp() local
724 cum = 0; in PCBDDCSubSchursSetUp()
728 const PetscInt row = cum + j; in PCBDDCSubSchursSetUp()
729 PetscInt col = cum; in PCBDDCSubSchursSetUp()
737 cum += subset_size; in PCBDDCSubSchursSetUp()
843 PetscInt n, n_I, size_schur, size_active_schur, cum, cum2; in PCBDDCSubSchursSetUp() local
857 PetscCall(ISGetLocalSize(sub_schurs->is_I, &cum)); in PCBDDCSubSchursSetUp()
858 if (cum != n_I) economic = PETSC_TRUE; in PCBDDCSubSchursSetUp()
883 cum = n_I + size_active_schur; in PCBDDCSubSchursSetUp()
890 PetscCall(PetscArraycpy(all_local_idx_N + cum, idxs, n_dir)); in PCBDDCSubSchursSetUp()
892 for (PetscInt j = 0; j < n_dir; j++) all_local_subid_N[j + cum] = local_subs[idxs[j]]; in PCBDDCSubSchursSetUp()
894 cum += n_dir; in PCBDDCSubSchursSetUp()
906 PetscCall(PetscArraycpy(all_local_idx_N + cum, idxs, n_v)); in PCBDDCSubSchursSetUp()
908 for (PetscInt j = 0; j < n_v; j++) all_local_subid_N[j + cum] = local_subs[idxs[j]]; in PCBDDCSubSchursSetUp()
910 cum += n_v; in PCBDDCSubSchursSetUp()
915 size_schur = cum - n_I; in PCBDDCSubSchursSetUp()
916 PetscCall(ISCreateGeneral(PETSC_COMM_SELF, cum, all_local_idx_N, PETSC_OWN_POINTER, &is_A_all)); in PCBDDCSubSchursSetUp()
921 if (cum == n) { in PCBDDCSubSchursSetUp()
1083 for (i = 0, cum = 0; i < sub_schurs->n_subs; i++) { in PCBDDCSubSchursSetUp()
1089 PetscCall(ISCreateStride(PETSC_COMM_SELF, subset_size, cum, 1, &is)); in PCBDDCSubSchursSetUp()
1104 cum += subset_size; in PCBDDCSubSchursSetUp()
1449 cum = cum2 = 0; in PCBDDCSubSchursSetUp()
1467 for (PetscInt j = 0; j < subset_size; j++) idx_work[j] = cum + j; in PCBDDCSubSchursSetUp()
1553 PetscCall(VecCreateSeqWithArray(PETSC_COMM_SELF, 1, subset_size, array + cum, &D)); in PCBDDCSubSchursSetUp()
1562 cum += subset_size; in PCBDDCSubSchursSetUp()
1633 cum = 0; in PCBDDCSubSchursSetUp()
1636 … PetscCall(PetscArraycpy(schur_factor + cum, tdata + i * (news + 1), size_active_schur - i)); in PCBDDCSubSchursSetUp()
1637 cum += size_active_schur - i; in PCBDDCSubSchursSetUp()
1639 …for (i = 0; i < nd; i++) schur_factor[cum + i] = PetscSqrtReal(PetscRealPart(data[(i + size_active… in PCBDDCSubSchursSetUp()
1748 for (i = 0, cum = 0; i < sub_schurs->n_subs; i++) { /* naive implementation */ in PCBDDCSubSchursSetUp()
1752 PetscCall(ISCreateStride(PETSC_COMM_SELF, subset_size, cum, 1, &EE)); in PCBDDCSubSchursSetUp()
1782 cum += subset_size; in PCBDDCSubSchursSetUp()
1789 cum = cum2 = 0; in PCBDDCSubSchursSetUp()
1802 for (PetscInt j = 0; j < subset_size; j++) idx_work[j] = cum + j; in PCBDDCSubSchursSetUp()
1857 cum += subset_size; in PCBDDCSubSchursSetUp()
1888 cum = 0; in PCBDDCSubSchursSetUp()
1890 … PetscCall(PetscArraycpy(data + i * (size_schur + 1), schur_factor + cum, size_active_schur - i)); in PCBDDCSubSchursSetUp()
1891 cum += size_active_schur - i; in PCBDDCSubSchursSetUp()
1898 … for (i = 0; i < nd; i++) data[(i + size_active_schur) * (size_schur + 1)] = schur_factor[cum + i]; in PCBDDCSubSchursSetUp()
1994 PetscInt cum; in PCBDDCSubSchursSetUp() local
1997 cum = 0; in PCBDDCSubSchursSetUp()
2003 PetscCallBLAS("LAPACKpotrf", LAPACKpotrf_("L", &B_N, array + cum, &B_N, &B_ierr)); in PCBDDCSubSchursSetUp()
2005 PetscCallBLAS("LAPACKpotri", LAPACKpotri_("L", &B_N, array + cum, &B_N, &B_ierr)); in PCBDDCSubSchursSetUp()
2008 …PetscCallBLAS("LAPACKsytrf", LAPACKsytrf_("L", &B_N, array + cum, &B_N, pivots, Bwork, &B_lwork, &… in PCBDDCSubSchursSetUp()
2010 … PetscCallBLAS("LAPACKsytri", LAPACKsytri_("L", &B_N, array + cum, &B_N, pivots, Bwork, &B_ierr)); in PCBDDCSubSchursSetUp()
2013 … PetscCallBLAS("LAPACKgetrf", LAPACKgetrf_(&B_N, &B_N, array + cum, &B_N, pivots, &B_ierr)); in PCBDDCSubSchursSetUp()
2015 …PetscCallBLAS("LAPACKgetri", LAPACKgetri_(&B_N, array + cum, &B_N, pivots, Bwork, &B_lwork, &B_ier… in PCBDDCSubSchursSetUp()
2020 cum += subset_size * subset_size; in PCBDDCSubSchursSetUp()