[petsc-users] false-positive leak report in log_view?

2016-08-02 Thread Mohammad Mirzadeh
I often use the memory usage information in log_view as a way to check on memory leaks and so far it has worked perfect. However, I had long noticed a false-positive report in memory leak for Viewers, i.e. destruction count is always one less than creation. Today, I noticed what seems to be a

Re: [petsc-users] block matrix without MatCreateNest

2016-08-02 Thread Matthew Knepley
On Tue, Aug 2, 2016 at 2:25 AM, Klaij, Christiaan wrote: > Thanks for your help! Going from individual blocks to a whole > matrix makes perfect sense if the blocks are readily available or > needed as fully functional matrices. Don't change that! Maybe add > the opposite? > >

[petsc-users] TS and petscFE

2016-08-02 Thread Maximilian Hartig
Hello all, I would like to run a transient problem with PetscFE. Example ex11.c seems relevant since it uses the PestcFV context to create boundary conditions and RHS Functions for the TS. Is there an easy way to do transient analysis with TS and petscFE or do I have to code my own

Re: [petsc-users] vec norm for local portion of a vector

2016-08-02 Thread Patrick Sanan
There is also VecGetLocalVector() and VecRestoreLocalVector() . You should be able to use those to (copy-free, hopefully) obtain the local entries of a global vector as a local vector which you could call the usual VecNorm() on. On Fri, Jul 29, 2016 at 4:41 PM, Barry Smith

Re: [petsc-users] block matrix without MatCreateNest

2016-08-02 Thread Klaij, Christiaan
Thanks for your help! Going from individual blocks to a whole matrix makes perfect sense if the blocks are readily available or needed as fully functional matrices. Don't change that! Maybe add the opposite? I'm surprised it's broken though: on this mailing list several petsc developers have