Re: [petsc-users] solving saddle point problem in distributed memory

2020-07-29 Thread Bin Liu
Yes, I did not wait until it aborts. I will try a small problem. From: Matthew Knepley [mailto:knep...@gmail.com] Sent: Wednesday, 29 July 2020 7:54 PM To: Bin Liu Cc: Barry Smith ; petsc-users@mcs.anl.gov Subject: Re: [petsc-users] solving saddle point problem in distributed memory On Wed, Jul

Re: [petsc-users] solving saddle point problem in distributed memory

2020-07-29 Thread Matthew Knepley
petsc-users@mcs.anl.gov > *Subject:* Re: [petsc-users] solving saddle point problem in distributed > memory > > > > On Tue, Jul 28, 2020 at 11:33 PM Bin Liu wrote: > > Dear Barry, > > > > Thanks for your suggestions. I just tested it. I still get the same error.

Re: [petsc-users] solving saddle point problem in distributed memory

2020-07-29 Thread Bin Liu
Knepley [mailto:knep...@gmail.com] Sent: Wednesday, 29 July 2020 5:21 PM To: Bin Liu Cc: Barry Smith ; petsc-users@mcs.anl.gov Subject: Re: [petsc-users] solving saddle point problem in distributed memory On Tue, Jul 28, 2020 at 11:33 PM Bin Liu mailto:lbllm2...@hotmail.com>> wrote: Dear

Re: [petsc-users] solving saddle point problem in distributed memory

2020-07-29 Thread Matthew Knepley
; -- > *From:* Barry Smith > *Sent:* Wednesday, July 29, 2020 10:49 AM > *To:* Bin Liu > *Cc:* petsc-users@mcs.anl.gov > *Subject:* Re: [petsc-users] solving saddle point problem in distributed > memory > > >No, PETSc direct solv

Re: [petsc-users] solving saddle point problem in distributed memory

2020-07-28 Thread Bin Liu
: Wednesday, July 29, 2020 10:49 AM To: Bin Liu Cc: petsc-users@mcs.anl.gov Subject: Re: [petsc-users] solving saddle point problem in distributed memory No, PETSc direct solvers do not work in parallel. Only SuperLU_DIST, MUMPS, and PasTix. -pc_type lu (or cholesky depending on the package

Re: [petsc-users] solving saddle point problem in distributed memory

2020-07-28 Thread Barry Smith
  <> > From: Barry Smith [mailto:bsm...@petsc.dev] > Sent: Wednesday, 29 July 2020 9:48 AM > To: Bin Liu > Cc: petsc-users@mcs.anl.gov > Subject: Re: [petsc-users] solving saddle point problem in distributed memory > > >SuperLU_DIST won't "magically"

Re: [petsc-users] solving saddle point problem in distributed memory

2020-07-28 Thread Bin Liu
2020 9:48 AM To: Bin Liu Cc: petsc-users@mcs.anl.gov Subject: Re: [petsc-users] solving saddle point problem in distributed memory SuperLU_DIST won't "magically" run on saddle point problems. It only does limited pivoting, realistically a parallel LU cannot always do complete

Re: [petsc-users] solving saddle point problem in distributed memory

2020-07-28 Thread Bin Liu
: Re: [petsc-users] solving saddle point problem in distributed memory On Tue, Jul 28, 2020 at 9:00 PM Bin Liu mailto:lbllm2...@hotmail.com>> wrote: Thanks for your tutorials. Yes. I tried PCFIELDSPLIT. However, it only works for sequential runs. When I run the code in distributed memo

Re: [petsc-users] solving saddle point problem in distributed memory

2020-07-28 Thread Barry Smith
: Stefano Zampini [mailto:stefano.zamp...@gmail.com > <mailto:stefano.zamp...@gmail.com>] > Sent: Tuesday, 28 July 2020 6:55 PM > To: Bin Liu mailto:lbllm2...@hotmail.com>> > Cc: petsc-users@mcs.anl.gov <mailto:petsc-users@mcs.anl.gov> > Subject: Re: [petsc-users] so

Re: [petsc-users] solving saddle point problem in distributed memory

2020-07-28 Thread Matthew Knepley
o:stefano.zamp...@gmail.com] > *Sent:* Tuesday, 28 July 2020 6:55 PM > *To:* Bin Liu > *Cc:* petsc-users@mcs.anl.gov > *Subject:* Re: [petsc-users] solving saddle point problem in distributed > memory > > > > If you want advice you should post the error trace PETSc reports. &

Re: [petsc-users] solving saddle point problem in distributed memory

2020-07-28 Thread Bin Liu
) does the direct solvers in petsc can run in distributed memory for solving saddle point problem? From: Stefano Zampini [mailto:stefano.zamp...@gmail.com] Sent: Tuesday, 28 July 2020 6:55 PM To: Bin Liu Cc: petsc-users@mcs.anl.gov Subject: Re: [petsc-users] solving saddle point problem in

Re: [petsc-users] solving saddle point problem in distributed memory

2020-07-28 Thread Stefano Zampini
If you want advice you should post the error trace PETSc reports. Anyway, solving Stokes is not so trivial (without direct solvers, you may need mesh dependent information), but we have examples for it https://gitlab.com/petsc/petsc/-/blob/master/src/ksp/ksp/tutorials/ex42.c

[petsc-users] solving saddle point problem in distributed memory

2020-07-28 Thread Bin Liu
I would like to solve a saddle point problem arising from the stokes equation. I got successful to use the direct solvers in sequential runs. However, I would like to extend it for distributed memory computation. I tried to use superlu_dist, but the program returns errors. Is it possible to solv