[petsc-users] Number of nonzeros > 2^31-1

Thanks for your answer, Matt. The problem is that I am using ML as a preconditioner and it does not support 64 bit indices (information reported in the PETSc configure phase). I am afraid I am stuck with this wrong nonzeros number. Nicolas Nicolas TARDIEU Ing. Chercheur EDF - R&D Dpt

[petsc-users] Number of nonzeros > 2^31-1

) 8 Is there a workaround? Should I worry about the result of my simulation? Thanks in advance, Nicolas Nicolas TARDIEU Ing. Chercheur EDF - R&D Dpt AMA nicolas.tardieu at edf.fr T?l. : 01 47 65 39 05 U

The MatMatSolve mystery

; I will add an error check telling you the B matrix is not dense so no one > will fall for this > again. > > Barry > > > > On Sep 10, 2008, at 3:39 PM, Nicolas Tardieu wrote: > > Dear PETSc users, >> >> The attached file is a naive test of MatMatSolve :

The MatMatSolve mystery

Dear PETSc users, The attached file is a naive test of MatMatSolve : I create a matrix,say A, I duplicate it, say BB, and I would like to compute A^{-1}*BB. Since A is the identity matrix, I would like to get the identity matrix. But I don't! Here is what I obtain by doing "./niko1f -n 3" : =

Installation with Intel or PGI compilers

Hi, I have some troubles using PETSc compiled with Intel compilers (version 10.1) in Fortran language in parallel on a 64 bits machine. The PetscInitialize always fails. In order to make it work, I have to make the following changes in petscconf.h. 320,321c320,321 < #ifndef PETSC_HAVE_IPXFARGC_

Block-smoothers with variable block size

Barry, Satish, Thank you very much for these precious informations! Nicolas 2008/3/21, Barry Smith : > > > On Mar 21, 2008, at 8:03 AM, Nicolas Tardieu wrote: > > > Dear Barry, > > > > > I have a few questions about inodes : > > - can the user specify the

Block-smoothers with variable block size

> http://www-unix.mcs.anl.gov/petsc/petsc-as/developers/index.html > > it is in src/mat/impls/aij/seq/inode.c > > > > Barry > > > > On Mar 20, 2008, at 1:52 PM, Nicolas Tardieu wrote: > > > Dear Barry, > > > > Thank you very much for your q

Block-smoothers with variable block size

t; If you do not have inodes then you can mimic the MatRelax_Inode() > code to make a new > routine MatPBRelax_SeqAIJ(). If you do implement this we'd love to add > it to PETSc. > > > Please let us know if you have additional questions. > > > Barry > >

Block-smoothers with variable block size

Dear PETSc experts, I am currently designing a multigrid PC that needs to operate on variable block size matrices. I know PETSc does not support this feature. So, I will do it "by hand". Neverthless I would like to be able to use some smoothers like Jacobi or SOR in their block version. What is

Parallel use of Prometheus

... > call KSPSetOperators (... > call KSPGetPC (... > call PCSetCoordinates(... ! I'm not fully certain > but your problem may have to be defined over R^3 for this to work well > call KSPSetFromOptions(... > > call KSPSolve (... > > > Nicol

Parallel use of Prometheus

Hi Sanjay, Have you used it on simple examples of the PETSc distribution? Nicolas 2008/1/16, Sanjay Govindjee : > > We have used it extensively with good success in the past. > -sanjay > > Nicolas Tardieu wrote: > > Hi everyone, > > > > I have been using

Parallel use of Prometheus

Hi everyone, I have been using Prometheus for scalar and vector problems with great success on sequential computers. But when trying to use it in parallel, Prometheus crashes. Then I tried it on some PETSc tests (src/ksp/ksp/examples/tests/ex3.c for instance) on parallel and they also crash... Doe