Skip to content
GitLab
Menu
Projects
Groups
Snippets
Help
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Contribute to GitLab
Sign in / Register
Toggle navigation
Menu
Open sidebar
dumux-repositories
dumux
Commits
fbc362e1
Commit
fbc362e1
authored
May 08, 2021
by
Timo Koch
Committed by
hanchuan
May 08, 2021
Browse files
[amg] Add warning for minimal setup without direct solver
parent
2d0145a1
Changes
1
Show whitespace changes
Inline
Side-by-side
dumux/linear/amgbackend.hh
View file @
fbc362e1
...
@@ -60,6 +60,8 @@ public:
...
@@ -60,6 +60,8 @@ public:
{
{
if
(
isParallel_
)
if
(
isParallel_
)
DUNE_THROW
(
Dune
::
InvalidStateException
,
"Using sequential constructor for parallel run. Use signature with gridView and dofMapper!"
);
DUNE_THROW
(
Dune
::
InvalidStateException
,
"Using sequential constructor for parallel run. Use signature with gridView and dofMapper!"
);
checkAvailabilityOfDirectSolver_
();
}
}
/*!
/*!
...
@@ -81,6 +83,7 @@ public:
...
@@ -81,6 +83,7 @@ public:
if
(
isParallel_
)
if
(
isParallel_
)
phelper_
=
std
::
make_unique
<
ParallelISTLHelper
<
LinearSolverTraits
>>
(
gridView
,
dofMapper
);
phelper_
=
std
::
make_unique
<
ParallelISTLHelper
<
LinearSolverTraits
>>
(
gridView
,
dofMapper
);
#endif
#endif
checkAvailabilityOfDirectSolver_
();
}
}
/*!
/*!
...
@@ -118,6 +121,16 @@ public:
...
@@ -118,6 +121,16 @@ public:
}
}
private:
private:
//! see https://gitlab.dune-project.org/core/dune-istl/-/issues/62
void
checkAvailabilityOfDirectSolver_
()
{
#if !HAVE_SUPERLU && !HAVE_UMFPACK
std
::
cout
<<
"
\n
AMGBiCGSTABBackend: No direct solver backend found. Using iterative solver as coarse grid solver.
\n
"
<<
"Note that dune-istl currently hard-codes a tolerance of 1e-2 for the iterative coarse grid solver.
\n
"
<<
"This may result in reduced accuracy or performance depending on your setup.
\n
Consider installing "
<<
"UMFPack (SuiteSparse) or SuperLU or apply the istl patch, see dumux/patches/README.md."
<<
std
::
endl
;
#endif
}
#if HAVE_MPI
#if HAVE_MPI
template
<
class
Matrix
,
class
Vector
>
template
<
class
Matrix
,
class
Vector
>
...
...
Write
Preview
Supports
Markdown
0%
Try again
or
attach a new file
.
Attach a file
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment