Skip to content
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

Ifpack2: Add Schur complements to the BlockTriDiag solver. #12344

Merged
merged 5 commits into from
Oct 24, 2023
Merged

Ifpack2: Add Schur complements to the BlockTriDiag solver. #12344

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
5 commits
Select commit Hold shift + click to select a range
4a17393
Add Schur complements to the BlockTriDiag solver.
kliegeois Sep 26, 2023
67ca800
Build warning free
kliegeois Oct 3, 2023
699b7da
Update the PR with Luc's comments
kliegeois Oct 17, 2023
89501b0
Reorder the loops to avoid branching
kliegeois Oct 18, 2023
418504c
Merge origin/develop
kliegeois Oct 18, 2023
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
6 changes: 4 additions & 2 deletions packages/ifpack2/example/BlockTriDi.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -18,7 +18,7 @@ namespace { // (anonymous)

// Values of command-line arguments.
struct CmdLineArgs {
CmdLineArgs ():blockSize(-1),numIters(10),numRepeats(1),tol(1e-12),nx(172),ny(-1),nz(-1),mx(1),my(1),mz(1),sublinesPerLine(1),useStackedTimer(false),overlapCommAndComp(false){}
CmdLineArgs ():blockSize(-1),numIters(10),numRepeats(1),tol(1e-12),nx(172),ny(-1),nz(-1),mx(1),my(1),mz(1),sublinesPerLine(1),sublinesPerLineSchur(1),useStackedTimer(false),overlapCommAndComp(false){}

std::string mapFilename;
std::string matrixFilename;
Expand All @@ -35,6 +35,7 @@ struct CmdLineArgs {
int my;
int mz;
int sublinesPerLine;
int sublinesPerLineSchur;
bool useStackedTimer;
bool overlapCommAndComp;
std::string problemName;
Expand Down Expand Up @@ -71,6 +72,7 @@ getCmdLineArgs (CmdLineArgs& args, int argc, char* argv[])
"Whether to run with overlapCommAndComp)");
cmdp.setOption("problemName", &args.problemName, "Human-readable problem name for Watchr plot");
cmdp.setOption("matrixType", &args.matrixType, "matrixType");
cmdp.setOption("sublinesPerLineSchur", &args.sublinesPerLineSchur, "sublinesPerLineSchur");
auto result = cmdp.parse (argc, argv);
return result == Teuchos::CommandLineProcessor::PARSE_SUCCESSFUL;
}
Expand Down Expand Up @@ -573,7 +575,7 @@ main (int argc, char* argv[])

{
Teuchos::TimeMonitor precSetupTimeMon (*precSetupTime);
precond = rcp(new BTDC(Ablock,parts,args.overlapCommAndComp));
precond = rcp(new BTDC(Ablock,parts,args.sublinesPerLineSchur,args.overlapCommAndComp));

if(rank0) std::cout<<"Initializing preconditioner..."<<std::endl;
precond->initialize ();
Expand Down
859 changes: 859 additions & 0 deletions packages/ifpack2/src/Ifpack2_BlockComputeResidualVector.hpp
View file
Open in desktop

Large diffs are not rendered by default.

Loading