VecDestroy and memory leak

2009-11-28 Thread nicolas aunai
seems to refer to sequential vectors... well I think I'm freeing everything I've created unless I misunderstood something ? Nico 2009/11/28 Jed Brown jed at 59a2.org: On Sat, 28 Nov 2009 10:18:04 +0100, nicolas aunai nicolas.aunai at gmail.com wrote: I've looked at the function VecDestroy

bug ? number of processor, coarse and fine DA for DMMG

2009-09-07 Thread nicolas aunai
way to extend the code for these other configurations. ? Barry On Sep 4, 2009, at 7:39 AM, nicolas aunai wrote: Hello, I have a code that uses the DMMG objet for a linear problem. The DA I create for the coarse grid is DA_XPERIODIC and 2D. I let petsc decide for the number

bug ? number of processor, coarse and fine DA for DMMG

2009-09-04 Thread nicolas aunai
Hello, I have a code that uses the DMMG objet for a linear problem. The DA I create for the coarse grid is DA_XPERIODIC and 2D. I let petsc decide for the number of processors in each direction. I have a problem which depends on the number of processors I use. If I want a 1024x1025 grid with 4

VecScatterCreateToZero

2009-08-28 Thread nicolas aunai
size. 2009/8/27 Matthew Knepley knepley at gmail.com: This should be good enough http://www.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-current/docs/manualpages/DA/DAGlobalToNaturalAllCreate.html Matt On Thu, Aug 27, 2009 at 12:51 PM, nicolas aunai nicolas.aunai at gmail.com wrote

VecScatterCreateToZero

2009-08-27 Thread nicolas aunai
Hello, I have a problem using 'VecScatterCreateToZero'. I attach to this email a test code that illustrates the problem. I create a 2D DA then I ask for the coordinates (that have been set just before), which are returned in a global vector. I want to write the coordinates in a file, so I use

linear elliptic vector equation

2009-08-20 Thread nicolas aunai
within the domain when I'm on periodic borders, i.e. is col[1].i = i-1 is ok when i=0 (left border) ? Ok that's it for the moment, Thx Nico 2009/8/19 Barry Smith bsmith at mcs.anl.gov: On Aug 19, 2009, at 2:04 AM, nicolas aunai wrote: Hi, OK then, lets go for multigrid. When I'm looking

linear elliptic vector equation

2009-08-19 Thread nicolas aunai
the residual ? I guess this is it for the moment. Thx Nico 2009/8/19 Barry Smith bsmith at mcs.anl.gov: On Aug 18, 2009, at 11:42 AM, nicolas aunai wrote: Hello, Yes it is a bit overkill indeed :-) I was looking at ex29 actually... which is still linear and with DA's but scalar. But I'd

linear elliptic vector equation

2009-08-18 Thread nicolas aunai
Dear all, I have a simulation code (particle in cell) in which I have to solve a linear vector equation 4 times per time step. A typical run consists of 50 000 time steps. The grid is rectangular and uniform of size Lx, Ly, with nx+1 and ny+1 points in the x and y direction respectively., (nx,

linear elliptic vector equation

2009-08-18 Thread nicolas aunai
. A On Tue, Aug 18, 2009 at 3:30 AM, nicolas aunai nicolas.aunai at gmail.com wrote: Dear all, I have a simulation code (particle in cell) in which I have to solve a linear vector equation 4 times per time step. A typical run consists of 50 000 time steps. The grid is rectangular and uniform