Lines Matching refs:am
2056 PetscInt reallocs = 0, i, *ai = a->i, *aj = a->j, am = A->rmap->n, *ui, *udiag; in MatICCFactorSymbolic_SeqAIJ() local
2075 PetscCall(PetscShmgetAllocateArray(am + 1, sizeof(PetscInt), (void **)&ui)); in MatICCFactorSymbolic_SeqAIJ()
2076 PetscCall(PetscMalloc1(am + 1, &udiag)); in MatICCFactorSymbolic_SeqAIJ()
2081 for (i = 0; i < am; i++) { in MatICCFactorSymbolic_SeqAIJ()
2086 PetscCall(PetscMalloc1(ui[am] + 1, &uj)); in MatICCFactorSymbolic_SeqAIJ()
2088 for (i = 0; i < am; i++) { in MatICCFactorSymbolic_SeqAIJ()
2099 PetscCall(PetscMalloc1(am + 1, &ajtmp)); in MatICCFactorSymbolic_SeqAIJ()
2103 PetscCall(PetscMalloc4(am, &uj_ptr, am, &uj_lvl_ptr, am, &jl, am, &il)); in MatICCFactorSymbolic_SeqAIJ()
2104 for (i = 0; i < am; i++) { in MatICCFactorSymbolic_SeqAIJ()
2105 jl[i] = am; in MatICCFactorSymbolic_SeqAIJ()
2110 nlnk = am + 1; in MatICCFactorSymbolic_SeqAIJ()
2111 PetscCall(PetscIncompleteLLCreate(am, am, nlnk, lnk, lnk_lvl, lnkbt)); in MatICCFactorSymbolic_SeqAIJ()
2114 PetscCall(PetscFreeSpaceGet(PetscRealIntMultTruncate(fill, (ai[am] + am) / 2), &free_space)); in MatICCFactorSymbolic_SeqAIJ()
2116 … PetscCall(PetscFreeSpaceGet(PetscRealIntMultTruncate(fill, (ai[am] + am) / 2), &free_space_lvl)); in MatICCFactorSymbolic_SeqAIJ()
2119 for (k = 0; k < am; k++) { /* for each active row k */ in MatICCFactorSymbolic_SeqAIJ()
2132 PetscCall(PetscIncompleteLLInit(ncols_upper, ajtmp, am, riip, &nlnk, lnk, lnk_lvl, lnkbt)); in MatICCFactorSymbolic_SeqAIJ()
2149 PetscCall(PetscICCLLAddSorted(ncols, cols, levels, uj, am, &nlnk, lnk, lnk_lvl, lnkbt, j)); in MatICCFactorSymbolic_SeqAIJ()
2164 i = am - k + 1; /* num of unfactored rows */ in MatICCFactorSymbolic_SeqAIJ()
2173 …PetscCall(PetscIncompleteLLClean(am, am, nzk, lnk, lnk_lvl, current_space->array, current_space_lv… in MatICCFactorSymbolic_SeqAIJ()
2202 PetscCall(PetscShmgetAllocateArray(ui[am] + 1, sizeof(PetscInt), (void **)&uj)); in MatICCFactorSymbolic_SeqAIJ()
2203 …PetscCall(PetscFreeSpaceContiguous_Cholesky(&free_space, uj, am, ui, udiag)); /* store matrix fact… in MatICCFactorSymbolic_SeqAIJ()
2212 PetscCall(PetscShmgetAllocateArray(ui[am], sizeof(PetscScalar), (void **)&b->a)); in MatICCFactorSymbolic_SeqAIJ()
2226 PetscCall(PetscMalloc1(am, &b->solve_work)); in MatICCFactorSymbolic_SeqAIJ()
2228 b->maxnz = b->nz = ui[am]; in MatICCFactorSymbolic_SeqAIJ()
2232 if (ai[am] != 0) { in MatICCFactorSymbolic_SeqAIJ()
2234 fact->info.fill_ratio_needed = ((PetscReal)2 * ui[am]) / (ai[am] + am); in MatICCFactorSymbolic_SeqAIJ()
2239 if (ai[am] != 0) { in MatICCFactorSymbolic_SeqAIJ()
2259 PetscInt i, am = A->rmap->n, *ai = a->i, *aj = a->j, reallocs = 0, prow; in MatCholeskyFactorSymbolic_SeqAIJ() local
2279 PetscCall(PetscShmgetAllocateArray(am + 1, sizeof(PetscInt), (void **)&ui)); in MatCholeskyFactorSymbolic_SeqAIJ()
2280 PetscCall(PetscMalloc1(am + 1, &udiag)); in MatCholeskyFactorSymbolic_SeqAIJ()
2285 PetscCall(PetscMalloc4(am, &ui_ptr, am, &jl, am, &il, am, &cols)); in MatCholeskyFactorSymbolic_SeqAIJ()
2286 for (i = 0; i < am; i++) { in MatCholeskyFactorSymbolic_SeqAIJ()
2287 jl[i] = am; in MatCholeskyFactorSymbolic_SeqAIJ()
2292 nlnk = am + 1; in MatCholeskyFactorSymbolic_SeqAIJ()
2293 PetscCall(PetscLLCreate(am, am, nlnk, lnk, lnkbt)); in MatCholeskyFactorSymbolic_SeqAIJ()
2296 PetscCall(PetscFreeSpaceGet(PetscRealIntMultTruncate(fill, (ai[am] + am) / 2), &free_space)); in MatCholeskyFactorSymbolic_SeqAIJ()
2299 for (k = 0; k < am; k++) { /* for each active row k */ in MatCholeskyFactorSymbolic_SeqAIJ()
2312 PetscCall(PetscLLAdd(ncols_upper, cols, am, &nlnk, lnk, lnkbt)); in MatCholeskyFactorSymbolic_SeqAIJ()
2325 PetscCall(PetscLLAddSorted(ncols, uj_ptr, am, &nlnk, lnk, lnkbt)); in MatCholeskyFactorSymbolic_SeqAIJ()
2340 i = am - k + 1; /* num of unfactored rows */ in MatCholeskyFactorSymbolic_SeqAIJ()
2347 PetscCall(PetscLLClean(am, am, nzk, lnk, current_space->array, lnkbt)); in MatCholeskyFactorSymbolic_SeqAIJ()
2370 PetscCall(PetscShmgetAllocateArray(ui[am], sizeof(PetscInt), (void **)&uj)); in MatCholeskyFactorSymbolic_SeqAIJ()
2371 …PetscCall(PetscFreeSpaceContiguous_Cholesky(&free_space, uj, am, ui, udiag)); /* store matrix fact… in MatCholeskyFactorSymbolic_SeqAIJ()
2377 PetscCall(PetscShmgetAllocateArray(ui[am], sizeof(PetscScalar), (void **)&b->a)); in MatCholeskyFactorSymbolic_SeqAIJ()
2393 PetscCall(PetscMalloc1(am, &b->solve_work)); in MatCholeskyFactorSymbolic_SeqAIJ()
2395 b->maxnz = b->nz = ui[am]; in MatCholeskyFactorSymbolic_SeqAIJ()
2399 if (ai[am] != 0) { in MatCholeskyFactorSymbolic_SeqAIJ()
2401 fact->info.fill_ratio_needed = ((PetscReal)2 * ui[am]) / (ai[am] + am); in MatCholeskyFactorSymbolic_SeqAIJ()
2406 if (ai[am] != 0) { in MatCholeskyFactorSymbolic_SeqAIJ()