Re: [Libmesh-users] sparse direct solver

2017-03-14 Thread David Knezevic
On Tue, Mar 14, 2017 at 7:48 PM, Kathrin Smetana wrote: > Hi John, > > Thanks very much for your answer and sorry for forgetting to let you know > the command line options. > > I am using the following command line options: > > ./example-opt -ksp_type preonly -pctype lu -pc_factor_mat_solver_pack

Re: [Libmesh-users] sparse direct solver

2017-03-14 Thread Kathrin Smetana
Hi David, Thanks very much! That solved the problem. Thanks, Kathrin Am 14.03.2017 um 18:52 schrieb David Knezevic mailto:[email protected]>>: On Tue, Mar 14, 2017 at 7:48 PM, Kathrin Smetana mailto:[email protected]>> wrote: Hi John, Thanks very much for your answer and sorry for f

Re: [Libmesh-users] sparse direct solver

2017-03-14 Thread Kathrin Smetana
Hi John, Thanks very much for your answer and sorry for forgetting to let you know the command line options. I am using the following command line options: ./example-opt -ksp_type preonly -pctype lu -pc_factor_mat_solver_package mumps So, if I understood everything correctly this should use a

Re: [Libmesh-users] sparse direct solver

2017-03-14 Thread John Peterson
On Tue, Mar 14, 2017 at 4:59 PM, Kathrin Smetana wrote: I believe this error: > [0]PETSC ERROR: - Error Message > -- > [0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/ > linearsolvertable.html for pos

[Libmesh-users] sparse direct solver

2017-03-14 Thread Kathrin Smetana
Dear libmesh users/developers, I am trying to solve several linear systems of equations for the same sparse matrix and different right-hand sides. To that end I am doing the following in my code: UniquePtr< LinearSolver > linear_solver = LinearSolver::build(init.comm(),libMesh::default_solver_