Re: [petsc-users] TSSolve problems

2015-07-23 Thread Abhyankar, Shrirang G.
Xinya, Take a look at the example $PETSC_DIR/src/snes/examples/tutorials/network/pflow/pf.c This example models power grid steady-state equations (power flow) using DMNetwork. Shri -Original Message- From: Li, Xinya xinya...@pnnl.gov Date: Thursday, July 23, 2015 at 3:18 PM To: Shri

Re: [petsc-users] TSSolve problems

2015-07-23 Thread Li, Xinya
Shri, Thank you for your suggestion. In the code, DMDACreat1D was used to create distributed array to manage parallel grid and vectors: DMDACreate1d(PETSC_COMM_WORLD, DM_BOUNDARY_NONE, 4*ngen, 1, 1, NULL, da); How to use DMPlex/DMnetwork instead of DMDA? Which function should I call? Thanks

Re: [petsc-users] Double free in KSPDestroy

2015-07-23 Thread Florian Lindner
Am Mittwoch, 22. Juli 2015, 13:05:57 schrieben Sie: On Jul 22, 2015, at 11:33 AM, Florian Lindner mailingli...@xgm.de wrote: Am Dienstag, 21. Juli 2015, 18:32:02 schrieben Sie: Try putting a breakpoint in KSPSetUp_GMRES and check the values of all the pointers immediately after

Re: [petsc-users] Double free in KSPDestroy

2015-07-23 Thread Matthew Knepley
On Thu, Jul 23, 2015 at 10:07 AM, Florian Lindner mailingli...@xgm.de wrote: Am Mittwoch, 22. Juli 2015, 13:05:57 schrieben Sie: On Jul 22, 2015, at 11:33 AM, Florian Lindner mailingli...@xgm.de wrote: Am Dienstag, 21. Juli 2015, 18:32:02 schrieben Sie: Try putting a

Re: [petsc-users] Null space/Near null space

2015-07-23 Thread Barry Smith
On Jul 23, 2015, at 3:42 AM, Lawrence Mitchell lawrence.mitch...@imperial.ac.uk wrote: -BEGIN PGP SIGNED MESSAGE- Hash: SHA1 On 22/07/15 23:57, Barry Smith wrote: If your matrix has a null space you should always be setting MatSetNullSpace. If you know the near null

Re: [petsc-users] Null space/Near null space

2015-07-23 Thread Lawrence Mitchell
-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 On 22/07/15 23:57, Barry Smith wrote: If your matrix has a null space you should always be setting MatSetNullSpace. If you know the near null space then you should also always set it (it cannot do any harm) because some preconditioners can