Re: [petsc-users] Loading Laplacian as Module

2016-09-23 Thread Barry Smith
I don't know much about modules so can't help, but PETSc variables are just like any other variables and so should behave in the same way. Barry > On Sep 23, 2016, at 2:31 PM, Manuel Valera wrote: > > Ok, i got this: > > RHS loaded, size: 213120 /

Re: [petsc-users] Loading Laplacian as Module

2016-09-23 Thread Manuel Valera
Ok, i got this: RHS loaded, size: 213120 / 213120 CSRMAt loaded, sizes: 213120 x 213120 8.43036175 s solveP pass: 1 RHS loaded, size: 213120 / 213120 [0]PETSC ERROR: MatGetSize() line 6295 in

Re: [petsc-users] Loading Laplacian as Module

2016-09-23 Thread Barry Smith
Ok, so the problem is not memory corruption. > 0]PETSC ERROR: - Error Message > -- > [0]PETSC ERROR: Invalid argument > [0]PETSC ERROR: Wrong type of object: Parameter # 1 > [0]PETSC ERROR: #1 MatGetSize() line

Re: [petsc-users] Loading Laplacian as Module

2016-09-23 Thread Manuel Valera
Barry, that was awesome, all the valgrind error dissappeared after using the mpiexec from petsc folder, the more you know... Anyway this is my output from valgrind running with those options: Last Update: 9/23/201612: 5:12 ucmsSeamount Entering MAIN loop. RHS loaded, size: 213120 /

Re: [petsc-users] Loading Laplacian as Module

2016-09-23 Thread Barry Smith
> On Sep 23, 2016, at 1:09 PM, Manuel Valera wrote: > > Thanks Barry, for the quick reply, > > I tried doing that once recently, not for this problem though, but it looks > like the model i'm working on isn't optimized at all for memory leaks, and > valgrind stopped

Re: [petsc-users] Loading Laplacian as Module

2016-09-23 Thread Manuel Valera
Thanks Barry, for the quick reply, I tried doing that once recently, not for this problem though, but it looks like the model i'm working on isn't optimized at all for memory leaks, and valgrind stopped with thousands of errors before reaching this part of the execution. Is there maybe an

Re: [petsc-users] Loading Laplacian as Module

2016-09-23 Thread Barry Smith
Run with valgrind to find the exact location of the first memory corruption. http://www.mcs.anl.gov/petsc/documentation/faq.html#valgrind > On Sep 23, 2016, at 12:47 PM, Manuel Valera wrote: > > Hello all, > > I'm trying to load my laplacian matrix into a fortran

[petsc-users] Loading Laplacian as Module

2016-09-23 Thread Manuel Valera
Hello all, I'm trying to load my laplacian matrix into a fortran module, and i have implemented it and it works for the first iteration of laplacian solver, but when starts the second step the laplacian matrix object becomes corrupts and looks like it loses one of it's dimensions. Can you help