[petsc-users] snessetjacobian question

Milan Mitrovic milanm at student.ethz.ch
Thu Sep 22 15:10:21 CDT 2011


ok... one more error... concerning this line
call MatSetValuesLocal(Jac,1,(/i-1/),Particles%nvlist(i)+1,ind(:),val(:),INSERT_VALUES,info)

when I try to deallocate val and ind I get a glibc error saying memory
is corrupt and that it is probably a double dealloc...

testing val and ind with isalloc gives true right before deallocation...

commenting out the above line of code, and the call to snes solve,
removes the error...

any ideas?


More information about the petsc-users mailing list