[petsc-users] Optimised vs Debug

2011-12-14 Thread Jed Brown
An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20111214/989e2d17/attachment.htm>

[petsc-users] run successfully on 20 processors but failed on 24 processors

2011-12-14 Thread Xiangdong Liang
I use MatNorm and VecNorm on A and b before I call kspsolve, and both of them are finite. Then I use fp_trap to catch where the nan comes from. However, it traces down to unlikely place pthread_join. The gdb where information is given below. Can you give me some help? Thanks. 0x7faeb6d8cbe5 i

[petsc-users] multigrid as preconditioner--message 8

2011-12-14 Thread M.Younas
: Warning not preallocating matrix storage > > What this means? Did I correctly preallocate the matrix? > > 2- I need to test MatSeqBAIJSetPreallocation() too but I little confused > with descriptions of the command in man-page: > > what are "nz" and "nnz[]" i

[petsc-users] run successfully on 20 processors but failed on 24 processors

2011-12-14 Thread Matthew Knepley
> >> > >> Xiangdong > > > > -- > > What most experimenters take for granted before they begin their > experiments > > is infinitely more interesting than any results to which their > experiments > > lead. > > -- Norbert Wiener > -- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20111214/4567210f/attachment.htm>

[petsc-users] save vector to bin-file at each time step

2011-12-14 Thread robert
On 12/14/2011 03:32 PM, Barry Smith wrote: > On Dec 14, 2011, at 8:00 AM, robert wrote: > >> Hello, >> >> I am quite new to petsc and therefore my question might be quite trivial. >> >> I have a 1D finite difference code and want to store several vectors in a >> binary: >> >> ierr = PetscView

[petsc-users] AMG as preconditioner

2011-12-14 Thread M.Younas
n the liner system first (once only) and then GMRES(m) is being used for stopping criterion? Kind regards, Younas RUG, Groningen, The Netherlands. -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/

[petsc-users] save vector to bin-file at each time step

2011-12-14 Thread robert
- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20111214/5a31167d/attachment.htm>

[petsc-users] Multigrid as preconditioner

2011-12-14 Thread M.Younas
n the liner system first (once only) and then GMRES(m) is being used for stopping criterion? Kind regards, Younas RUG, Groningen, The Netherlands. -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/

[petsc-users] AMG as preconditioner

2011-12-14 Thread Matthew Knepley
-- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20111214/ae73c9d3/attachment.htm>

[petsc-users] Two Questions about Matrix Preallocation

2011-12-14 Thread behzad baghapour
at are "nz" and "nnz[]" in this case? what about number of blocks per row? Thanks bothering too much :-) BehZad -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20111214/abec449b/attachment.htm>

[petsc-users] multigrid as preconditioner--message 8

2011-12-14 Thread Barry Smith
gt; I am facing with two questions about Block-Sparse Preallocation: >> >> 1- I set array "nnz[]" and did: >> >> MatSetType ( JAC, MATSEQBAIJ ); >> MatSetSizes( JAC, PETSC_DECIDE, PETSC_DECIDE, nt, nt ); >> MatSeqAIJSetPreallocation( JAC, 0

[petsc-users] save vector to bin-file at each time step

2011-12-14 Thread Barry Smith
On Dec 14, 2011, at 8:00 AM, robert wrote: > Hello, > > I am quite new to petsc and therefore my question might be quite trivial. > > I have a 1D finite difference code and want to store several vectors in a > binary: > > ierr = PetscViewerCreate(PETSC_COMM_WORLD,&viewer);CHKERRQ(ierr);

[petsc-users] Multigrid as preconditioner

2011-12-14 Thread Barry Smith
On Dec 14, 2011, at 7:07 AM, M.Younas wrote: > Hi, PETSc, Team, > > I have solved a non-symmetric Poisson problem using GMRES(m) as solver and > BoomerAMG (from Hypre) as preconditioner. But I don't understand how the > BoomerAMG works as preconditioner with GMRES(m). Can anyone explain the

[petsc-users] Two Questions about Matrix Preallocation

2011-12-14 Thread Barry Smith
On Dec 14, 2011, at 2:02 AM, behzad baghapour wrote: > Dear Developers, > > I am facing with two questions about Block-Sparse Preallocation: > > 1- I set array "nnz[]" and did: > > MatSetType ( JAC, MATSEQBAIJ ); > MatSetSizes( JAC, PETSC_DECIDE, PETSC_DECIDE, nt, nt ); > MatSe