and from your email I got off-list (please try to use the mailinglist):
> [0]PETSC ERROR: #1 PetscCommDuplicate() line 137 in
> /home/anna/petsc-3.6.4/src/sys/objects/tagm.c
> An error occurred in line <724> of file
> </home/anna/dealii-8.5.0/source/lac/petsc_solver.cc> in function
> void dealii::PETScWrappers::SparseDirectMUMPS::solve(const
> dealii::PETScWrappers::MatrixBase&, dealii::PETScWrappers::VectorBase&,
> const dealii::PETScWrappers::VectorBase&)
I assume that you either did not set the MPI communicator in the mumps
object or you set it to something invalid.
On Mon, Sep 25, 2017 at 8:29 AM, Timo Heister <heis...@clemson.edu> wrote:
> Anna,
>
>> The main reason I do this is that I do not understand how to reuse this
>> decomposition in deal.ii.
>> I am relatively new to deal.ii and C++, and I have never used MUMPS before.
>
> Well, this has nothing to do with MUMPS or deal.II. It sounds like you
> are struggling because you are not familiar with c++ templates and
> template specialization and how we use it in step-22. I also think you
> are hung up on a detail ("how can I switch between solvers based on
> the dimension?"), while you could just always use MUMPS.
>
>> since PETScWrappers::SparseDirectMUMPS does not have member function vmult.
>
> You could create a specialized InverseMatrix class that does something
> different depending on the type PreconditionerType. I would just
> unconditionally replace the vmult() inside InverseMatrix to call
> ::solve() of SparseDirectMUMPS.
>
>> I would appreciate any advice on this. Maybe there is some example in
>> deal.ii that use MUMPS to construct a preconditioner.
>
> for example:
> https://urldefense.proofpoint.com/v2/url?u=https-3A__github.com_dealii_dealii_blob_master_tests_mpi_step-2D40-5Fdirect-5Fsolver.cc&d=DwIBaQ&c=Ngd-ta5yRYsqeUsEDgxhcqsYYY1Xs5ogLxWPA_2Wlc4&r=4k7iKXbjGC8LfYxVJJXiaYVu6FRWmEjX38S7JmlS9Vw&m=1A8jBef7oOJePKduQoE3_poqbYEfwhW0Egy3MO3dVSk&s=g_sD1-HBOWyhL4OywtSIRMXOm4UH7tzq3rvYXsKc4rg&e=
>
> https://urldefense.proofpoint.com/v2/url?u=https-3A__github.com_dealii_dealii_blob_master_tests_petsc_sparse-5Fdirect-5Fmumps.cc&d=DwIBaQ&c=Ngd-ta5yRYsqeUsEDgxhcqsYYY1Xs5ogLxWPA_2Wlc4&r=4k7iKXbjGC8LfYxVJJXiaYVu6FRWmEjX38S7JmlS9Vw&m=1A8jBef7oOJePKduQoE3_poqbYEfwhW0Egy3MO3dVSk&s=t5b9JVDBAWGDxrQI-5oI-gBFq7pfzKrWpXjhXN-HUhA&e=
>
>
>
>
>
> --
> Timo Heister
> http://www.math.clemson.edu/~heister/
--
Timo Heister
http://www.math.clemson.edu/~heister/
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to dealii+unsubscr...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.
