Skip to content

Hypre iterative ILU features #4818

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
merged 4 commits into from
Dec 3, 2025
Merged

Hypre iterative ILU features #4818

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
4 commits
Select commit Hold shift + click to select a range
3dada17
Iterative ILU as a BoomerAMG smoother
itopcuoglu Aug 29, 2025
0ad3084
ILU as a stand-alone preconditioner, including iterative ILU
itopcuoglu Aug 29, 2025
50666bf
Formatting
itopcuoglu Nov 25, 2025
e8fdb29
Adding HYPRE_RELEASE_NUMBER check
itopcuoglu Nov 25, 2025
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
1 change: 1 addition & 0 deletions Src/Extern/HYPRE/AMReX_HypreIJIface.H
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -60,6 +60,7 @@ private:
// Preconditioners
void boomeramg_precond_configure(const std::string& prefix);
void euclid_precond_configure(const std::string& prefix);
void ilu_precond_configure(const std::string& prefix);

// Solvers
void boomeramg_solver_configure(const std::string& prefix);
Expand Down
46 changes: 46 additions & 0 deletions Src/Extern/HYPRE/AMReX_HypreIJIface.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -203,6 +203,8 @@ void HypreIJIface::init_preconditioner (
boomeramg_precond_configure(prefix);
} else if (name == "euclid") {
euclid_precond_configure(prefix);
} else if (name == "ILU") {
ilu_precond_configure(prefix);
} else {
amrex::Abort("Invalid HYPRE preconditioner specified: " + name);
}
Expand Down Expand Up @@ -324,6 +326,10 @@ void HypreIJIface::boomeramg_precond_configure (const std::string& prefix)
hpp("bamg_ilu_level", HYPRE_BoomerAMGSetILULevel);
hpp("bamg_ilu_max_iter", HYPRE_BoomerAMGSetILUMaxIter);
#if defined(HYPRE_RELEASE_NUMBER) && (HYPRE_RELEASE_NUMBER >= 22900)
hpp("bamg_ilu_iterative_algorithm_type", HYPRE_BoomerAMGSetILUIterSetupType);
hpp("bamg_ilu_iterative_setup_type", HYPRE_BoomerAMGSetILUIterSetupOption);
hpp("bamg_ilu_iterative_max_iter", HYPRE_BoomerAMGSetILUIterSetupMaxIter);
hpp("bamg_ilu_iterative_tolerance", HYPRE_BoomerAMGSetILUIterSetupTolerance);
hpp("bamg_ilu_reordering_type", HYPRE_BoomerAMGSetILULocalReordering);
hpp("bamg_ilu_tri_solve", HYPRE_BoomerAMGSetILUTriSolve);
hpp("bamg_ilu_lower_jacobi_iters", HYPRE_BoomerAMGSetILULowerJacobiIters);
Expand Down Expand Up @@ -373,6 +379,46 @@ void HypreIJIface::euclid_precond_configure (const std::string& prefix)
hpp("euclid_mem", HYPRE_EuclidSetMem, 0);
}

void HypreIJIface::ilu_precond_configure (const std::string& prefix)
{
if (m_verbose > 2) {
amrex::Print() << "Creating ILU preconditioner" << '\n';
}
HYPRE_ILUCreate(&m_precond);

// Setup the pointers
m_precondDestroyPtr = &HYPRE_ILUDestroy;
m_precondSetupPtr = &HYPRE_ILUSetup;
m_precondSolvePtr = &HYPRE_ILUSolve;


HypreOptParse hpp(prefix, m_precond);
#if defined(HYPRE_RELEASE_NUMBER) && (HYPRE_RELEASE_NUMBER >= 22100)
// Process ILU smoother parameters
// ParILUK
hpp("ilu_type", HYPRE_ILUSetType);
hpp("ilu_max_iter", HYPRE_ILUSetMaxIter);
hpp("ilu_tolerance", HYPRE_ILUSetTol);
hpp("ilu_reordering_type", HYPRE_ILUSetLocalReordering);
hpp("ilu_print_level", HYPRE_ILUSetPrintLevel);

// ILUK
hpp("ilu_fill", HYPRE_ILUSetLevelOfFill);
// ILUT
hpp("ilu_max_nnz_per_row", HYPRE_ILUSetMaxNnzPerRow);
hpp("ilu_drop_threshold", HYPRE_ILUSetDropThreshold);
#if defined(HYPRE_RELEASE_NUMBER) && (HYPRE_RELEASE_NUMBER >= 22900)
hpp("ilu_iterative_algorithm_type", HYPRE_ILUSetIterativeSetupType);
hpp("ilu_iterative_setup_type", HYPRE_ILUSetIterativeSetupOption);
hpp("ilu_iterative_max_iter", HYPRE_ILUSetIterativeSetupMaxIter);
hpp("ilu_iterative_tolerance", HYPRE_ILUSetIterativeSetupTolerance);
hpp("ilu_tri_solve", HYPRE_ILUSetTriSolve);
hpp("ilu_lower_jacobi_iters", HYPRE_ILUSetLowerJacobiIters);
hpp("ilu_upper_jacobi_iters", HYPRE_ILUSetUpperJacobiIters);
#endif
#endif
}

void HypreIJIface::boomeramg_solver_configure (const std::string& prefix)
{
if (m_has_preconditioner) {
Expand Down
Loading