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.

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

PETScWrappers::MPI::SparseMatrix: use SEQAIJ with one process #14683

Merged
merged 2 commits into from Jan 18, 2023

Conversation

stefanozampini
Copy link
Contributor

With one process, it is more efficient to use MATSEQAIJ

So, use MATAIJ to create the object (which is a proxy for SEQAIJ in sequential and MPIAIJ in parallel) and add missing preallocators

@stefanozampini stefanozampini force-pushed the stefanozampini/petsc-usematseqaij branch from 640cd0d to 8439b7c Compare January 17, 2023 07:07
@stefanozampini
Copy link
Contributor Author

@bangerth This is ready to go. I have implemented the comments you requested

Comment on lines +439 to +442
// Here we both call the specific API since this is how
// PETSc polymorphism works. If the matrix is of type MPIAIJ,
// the second call is dummy. If the matrix is of type SEQAIJ,
// the first call is dummy.
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Ah, fascinating. It's not a design I would have chosen, but ok :-)

@bangerth
Copy link
Member

/rebuild

@tamiko
Copy link
Member

tamiko commented Jan 18, 2023

We have an unrelated intermittend failure of mpi/petsc_step_27.mpirun=2. I observe this on a number of configurations at the moment.

@tamiko tamiko merged commit 8bd7332 into dealii:master Jan 18, 2023
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

Successfully merging this pull request may close these issues.

None yet

3 participants