Re: [deal.II] MUMPS with PETScWrappers::MPI::BlockSparseMatrix

2019-09-09 Thread Richard Schussnig
Hi Konrad! You can use the parallel direct solver in the schur complement, for orientation, take a look at step-57 (should be Navier-Stokes with direct solver for the A-block, if im not mistaken). However, my inferior C++ knowledge did not allow me to do the factorization in the constructor of

Re: [deal.II] MUMPS with PETScWrappers::MPI::BlockSparseMatrix

2019-09-08 Thread Konrad
Thank you, David, so I will see if I can maybe still use it in the Schur complement somehow. Best, Konrad On Friday, September 6, 2019 at 7:46:23 PM UTC+2, David Wells wrote: > > Hi Konrad, > > I don't think that it is possible to use MUMPS with a block matrix for > exactly this reason. I

Re: [deal.II] MUMPS with PETScWrappers::MPI::BlockSparseMatrix

2019-09-06 Thread David Wells
Hi Konrad, I don't think that it is possible to use MUMPS with a block matrix for exactly this reason. I think that if you want to use MUMPS you will need to copy the block matrix into a non-block matrix. Thanks, David On Fri, Sep 6, 2019 at 12:25 PM Konrad wrote: > Dear deal.ii community, >