Re: [Libmesh-users] Bug in PetscVector::create_subvector

2016-07-15 Thread Jonas Ballani
users Betreff: Re: [Libmesh-users] Bug in PetscVector::create_subvector On Thu, Jul 14, 2016 at 9:08 AM, Jonas Ballani mailto:[email protected]>> wrote: Hi John, Thanks for your quick reply. If you don't initialize y, say by just calling PetscVector< Number > y( mesh.

Re: [Libmesh-users] Bug in PetscVector::create_subvector

2016-07-14 Thread John Peterson
On Thu, Jul 14, 2016 at 9:08 AM, Jonas Ballani wrote: > Hi John, > > > Thanks for your quick reply. If you don't initialize y, say by just calling > > > PetscVector< Number > y( mesh.comm() ); > > > you get a memory leak. > > [1]PETSC ERROR: - Error Message > -

[Libmesh-users] Bug in PetscVector::create_subvector

2016-07-14 Thread Jonas Ballani
Hi all, I would like to extract a subvector of a distributed PetscVector using the create_subvector routine. My example programs looks as follows: #include #include #include #include "libmesh/libmesh.h" #include "libmesh/mesh.h" #include "libmesh/mesh_generation.h" #include "libmesh/linear_i