PaStiX crash

Hi Barry, here again with line numbers: http://pastebin.com/m630324e I noticed, the Pastix-version with only '-g' gives no errors. Hope this output now helps for debugging. Cheers, ando Barry Smith schrieb: I think if you compile all the code (including Scotch) with the -g option as

GMRES performance

Hello, I followed the suggestion in PETSc previous threads by adding options -pc_type asm -sub_pc_type lu Now the solver is really fast! Thanks you Jarunan Quoting Barry Smith bsmith at mcs.anl.gov: On Oct 20, 2009, at 9:52 AM, jarunan at ascomp.ch wrote: Hello, I would like

KSP different results with the default and direct LU

Dear all, I am trying get myself more acquainted with PETSc. I am trying to solve a linear equation system which is conditioned badly. If I use the approach as given in ex1.c of the KSP section(by default, I read from the manual that it uses restarted GMRES , preconditioned with JACOBI in

KSP different results with the default and direct LU

On Oct 22, 2009, at 6:08 AM, Umut Tabak wrote: Dear all, I am trying get myself more acquainted with PETSc. I am trying to solve a linear equation system which is conditioned badly. If I use the approach as given in ex1.c of the KSP section(by default, I read from the manual that it

KSP different results with the default and direct LU

Barry Smith wrote: When testing always run with -ksp_converged_reason or call KSPGetConvergedReason() after KSP solve to determine if PETSc thinks it has actually solved the system. Also since iterative solvers only solve to some tolerance the answer may not be wrong it may just be

Manually setting DA min and max indexes

Hello: Is it possible to set the min/max m,n,p indexes when creating a 3D distributed array? I'm writing a 3D code on a structured mesh, and the mesh is partitioned among the processes elsewhere. Assuming that each process already knows what chunk of the mesh it owns, is it possible to give Petsc

Manually setting DA min and max indexes

... Name: signature.asc Type: application/pgp-signature Size: 261 bytes Desc: OpenPGP digital signature URL: http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20091022/31960152/attachment.pgp

One simple question about complex number

: http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20091022/7910e4e9/attachment.htm

One simple question about complex number

we use 'std:complex' with --with-clanguage=cxx [and with c - we use c99 complex support]. I'm not sure what stl::complex is. Satish On Thu, 22 Oct 2009, Yujie wrote: Hi, PETSc Developers, What is the difference between PETSc complex number and stl::complex. If I define a variable var using

One simple question about complex number

. Regards, Yujie -- next part -- An HTML attachment was scrubbed... URL: http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20091022/5a637134/attachment.htm

One simple question about complex number

-- A non-text attachment was scrubbed... Name: signature.asc Type: application/pgp-signature Size: 261 bytes Desc: OpenPGP digital signature URL: http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20091022/91c0d10d/attachment.pgp

MatZeroRows ignore negative indices?

/attachments/20091022/1d88a75f/attachment.htm

MatZeroRows ignore negative indices?

No, there is no support for this. You should pass in an array that only contains the rows you want zeroed. You cannot pass in negative entries for row numbers. Barry On Oct 22, 2009, at 5:38 PM, Pedro Juan Torres Lopez wrote: Hello, In MatZeroRows() if I put negative value in the