[petsc-users] No speedup of tutorial in parallel

2010-04-07 Thread Ivo Roghair
Hi, I am trying to build PETSc into my CFD code. I just built the 3.1-p0 against a version of open-mpi 1.4.1. I thought src/ksp/ksp/ex/tut/ex2.c seemed like a comparable case (though much smaller and fewer matrix bands) of what I need. I tried running the default ex2.c example on a single-core

[petsc-users] PaStiX does not work

2010-04-07 Thread Gong Ding
Dear Petsc developer, I found a problem that PASTIX solver can not be loaded from PETSC, even for 3.1 version. When I use following code, the ksp_view told me that PC is ilu. And the snes solver can not convergence at all. ierr = KSPSetType (ksp, (char*) KSPPREONLY); assert(!ierr);//it

[petsc-users] configure failed on MPI and Intel Compiler environment

2010-04-07 Thread Takuya Sekikawa
Dear Petsc developpers, I tried to configure PETSc (ver 3.0) on Linux (64bit), MPI-enabled, Intel Compiler (icc/icpc) environment, but it failed. ./config/configure.py --with-mpi=1 --with-mpi-dir=/opt/home/sekikawa/bin-personal/mpich2 \ --with-x=0 --with-fc=0 --with-debugging=0

[petsc-users] sqrt for a PetscScalar

2010-04-07 Thread Umut Tabak
Dear all, I am getting a compilation error from the sqrt(cmath header) function with PETSc 3.1, I configured PETSc with complex scalar type so I am suspecting that this is the source of the problem. Is there a simple workaround for this? And Is there an sqrt function for a 'PetscScalar', I

[petsc-users] sqrt for a PetscScalar

2010-04-07 Thread Matthew Knepley
was scrubbed... URL: http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20100407/a00f8e73/attachment.htm

[petsc-users] sqrt for a PetscScalar

2010-04-07 Thread Umut Tabak
Matthew Knepley wrote: My guess is that you want PetscReal xTMx; Dear Matthew, Thanks for the quick reply, Not sure, here is the code for the related part, I used a VecTDot before that so I need PetscScalar I guess. PetscErrorCode ierr; Vec x, Bx; int rSz, cSz; PetscScalar

[petsc-users] sqrt for a PetscScalar

2010-04-07 Thread Matthew Knepley
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/20100407/675a643a/attachment-0001.htm

[petsc-users] PaStiX does not work

2010-04-07 Thread Gong Ding
Dear Barry, Very sorry. It is my mistake. I set PC to ILU again in another function when PasTIX is used (but MUMPS is right). Now it works, thanks. Yours Gong Ding If you use MUMPS then you should be able to use PaStiX just by changing that one line. It should not suddenly say that

[petsc-users] sqrt for a PetscScalar

2010-04-07 Thread Barry Smith
no match for ?operator/? in ?1 / std::sqrt [with It is complaining about the operator / not the sqrt. Try 1.0 / sqrt(xTMx) instead. Sometimes the complex class is real braindead about casting from ints. Barry On Apr 7, 2010, at 5:47

[petsc-users] configure failed on MPI and Intel Compiler environment

2010-04-07 Thread Satish Balay
petsc-3.0 [with latest patches] configure should be able to handle MLK-10. Don't remember MLK-11. without MPI (--with-mpi=0), this configure succeeded with no problem. Claims MLK-10 was detected properly without mpi. This is weird - but again configure.log for both builds will clarify the

[petsc-users] minres: PETScMATLAB

2010-04-07 Thread Barry Smith
Umut, Send one of the matrices with this behavior to petsc-maint at mcs.anl.gov and we'll see what is going on. Also in that email tell use exactly how you run it in Matlab. Barry On Apr 7, 2010, at 5:58 PM, Umut Tabak wrote: Dear all, I have tried on some