Lines Matching refs:am
2474 PetscInt bs = A->rmap->bs, am = a->mbs, *ai = a->i, *aj = a->j; in MatICCFactorSymbolic_SeqSBAIJ() local
2498 PetscCall(PetscMalloc1(am + 1, &ui)); in MatICCFactorSymbolic_SeqSBAIJ()
2499 PetscCall(PetscMalloc1(am + 1, &udiag)); in MatICCFactorSymbolic_SeqSBAIJ()
2505 for (i = 0; i < am; i++) { in MatICCFactorSymbolic_SeqSBAIJ()
2510 PetscCall(PetscMalloc1(ui[am], &uj)); in MatICCFactorSymbolic_SeqSBAIJ()
2512 for (i = 0; i < am; i++) { in MatICCFactorSymbolic_SeqSBAIJ()
2524 PetscCall(PetscMalloc4(am, &uj_ptr, am, &uj_lvl_ptr, am, &il, am, &jl)); in MatICCFactorSymbolic_SeqSBAIJ()
2525 for (i = 0; i < am; i++) { in MatICCFactorSymbolic_SeqSBAIJ()
2526 jl[i] = am; in MatICCFactorSymbolic_SeqSBAIJ()
2531 nlnk = am + 1; in MatICCFactorSymbolic_SeqSBAIJ()
2532 PetscCall(PetscIncompleteLLCreate(am, am, nlnk, lnk, lnk_lvl, lnkbt)); in MatICCFactorSymbolic_SeqSBAIJ()
2535 PetscCall(PetscFreeSpaceGet(PetscRealIntMultTruncate(fill, ai[am] + 1), &free_space)); in MatICCFactorSymbolic_SeqSBAIJ()
2539 PetscCall(PetscFreeSpaceGet(PetscRealIntMultTruncate(fill, ai[am] + 1), &free_space_lvl)); in MatICCFactorSymbolic_SeqSBAIJ()
2543 for (k = 0; k < am; k++) { /* for each active row k */ in MatICCFactorSymbolic_SeqSBAIJ()
2549 PetscCall(PetscIncompleteLLInit(ncols, cols, am, rip, &nlnk, lnk, lnk_lvl, lnkbt)); in MatICCFactorSymbolic_SeqSBAIJ()
2567 PetscCall(PetscICCLLAddSorted(ncols, cols, levels, uj, am, &nlnk, lnk, lnk_lvl, lnkbt, j)); in MatICCFactorSymbolic_SeqSBAIJ()
2583 i = am - k + 1; /* num of unfactored rows */ in MatICCFactorSymbolic_SeqSBAIJ()
2592 …PetscCall(PetscIncompleteLLClean(am, am, nzk, lnk, lnk_lvl, current_space->array, current_space_lv… in MatICCFactorSymbolic_SeqSBAIJ()
2618 PetscCall(PetscMalloc1(ui[am], &uj)); in MatICCFactorSymbolic_SeqSBAIJ()
2619 …PetscCall(PetscFreeSpaceContiguous_Cholesky(&free_space, uj, am, ui, udiag)); /* store matrix fact… in MatICCFactorSymbolic_SeqSBAIJ()
2630 PetscCall(PetscShmgetAllocateArray(ui[am], sizeof(PetscScalar), (void **)&b->a)); in MatICCFactorSymbolic_SeqSBAIJ()
2646 PetscCall(PetscMalloc1(am + 1, &b->solve_work)); in MatICCFactorSymbolic_SeqSBAIJ()
2648 b->maxnz = b->nz = ui[am]; in MatICCFactorSymbolic_SeqSBAIJ()
2652 if (ai[am] != 0) { in MatICCFactorSymbolic_SeqSBAIJ()
2653 fact->info.fill_ratio_needed = ((PetscReal)ui[am]) / ai[am]; in MatICCFactorSymbolic_SeqSBAIJ()
2658 if (ai[am] != 0) { in MatICCFactorSymbolic_SeqSBAIJ()
2676 PetscInt bs = A->rmap->bs, am = a->mbs; in MatICCFactorSymbolic_SeqSBAIJ_inplace() local
2717 PetscCall(PetscMalloc1(am + 1, &ui)); in MatICCFactorSymbolic_SeqSBAIJ_inplace()
2722 PetscCall(PetscMalloc4(am, &uj_ptr, am, &uj_lvl_ptr, am, &il, am, &jl)); in MatICCFactorSymbolic_SeqSBAIJ_inplace()
2723 for (i = 0; i < am; i++) { in MatICCFactorSymbolic_SeqSBAIJ_inplace()
2724 jl[i] = am; in MatICCFactorSymbolic_SeqSBAIJ_inplace()
2729 nlnk = am + 1; in MatICCFactorSymbolic_SeqSBAIJ_inplace()
2730 PetscCall(PetscIncompleteLLCreate(am, am, nlnk, lnk, lnk_lvl, lnkbt)); in MatICCFactorSymbolic_SeqSBAIJ_inplace()
2733 PetscCall(PetscFreeSpaceGet(PetscRealIntMultTruncate(fill, ai[am] + 1), &free_space)); in MatICCFactorSymbolic_SeqSBAIJ_inplace()
2737 PetscCall(PetscFreeSpaceGet(PetscRealIntMultTruncate(fill, ai[am] + 1), &free_space_lvl)); in MatICCFactorSymbolic_SeqSBAIJ_inplace()
2741 for (k = 0; k < am; k++) { /* for each active row k */ in MatICCFactorSymbolic_SeqSBAIJ_inplace()
2746 PetscCall(PetscIncompleteLLInit(ncols, cols, am, rip, &nlnk, lnk, lnk_lvl, lnkbt)); in MatICCFactorSymbolic_SeqSBAIJ_inplace()
2763 … PetscCall(PetscICCLLAddSorted(ncols, cols, levels, cols_lvl, am, &nlnk, lnk, lnk_lvl, lnkbt, j)); in MatICCFactorSymbolic_SeqSBAIJ_inplace()
2778 …i = am - k + 1; /* num of unfactored r… in MatICCFactorSymbolic_SeqSBAIJ_inplace()
2786 …PetscCall(PetscIncompleteLLClean(am, am, nzk, lnk, lnk_lvl, current_space->array, current_space_lv… in MatICCFactorSymbolic_SeqSBAIJ_inplace()
2812 PetscCall(PetscMalloc1(ui[am] + 1, &uj)); in MatICCFactorSymbolic_SeqSBAIJ_inplace()
2816 if (ai[am] != 0) { in MatICCFactorSymbolic_SeqSBAIJ_inplace()
2817 ratio_needed = ((PetscReal)ui[am]) / ((PetscReal)ai[am]); in MatICCFactorSymbolic_SeqSBAIJ_inplace()
2828 PetscCall(PetscShmgetAllocateArray(ui[am], sizeof(PetscScalar), (void **)&b->a)); in MatICCFactorSymbolic_SeqSBAIJ_inplace()
2842 PetscCall(PetscMalloc1(am + 1, &b->solve_work)); in MatICCFactorSymbolic_SeqSBAIJ_inplace()
2844 b->maxnz = b->nz = ui[am]; in MatICCFactorSymbolic_SeqSBAIJ_inplace()
2850 if (ai[am] != 0) { in MatICCFactorSymbolic_SeqSBAIJ_inplace()