Re: [petsc-users] PETSC ERROR: Can only handle MPIU_REAL or MPIU_COMPLEX with ML preconditioner

2014-10-01 Thread Matthew Knepley
On Wed, Oct 1, 2014 at 3:51 AM, Chris Richardson ch...@bpi.cam.ac.uk wrote: I have just run into this error after pulling the master branch of petsc, when using the ML preconditioner in FEniCS. It looks like a type error, but I guess FEniCS uses PETScInt, so should be safe? PETSC ERROR:

Re: [petsc-users] PETSC ERROR: Can only handle MPIU_REAL or MPIU_COMPLEX with ML preconditioner

2014-10-01 Thread Chris Richardson
On 01/10/2014 11:33, Matthew Knepley wrote: On Wed, Oct 1, 2014 at 3:51 AM, Chris Richardson ch...@bpi.cam.ac.uk wrote: I have just run into this error after pulling the master branch of petsc, when using the ML preconditioner in FEniCS. It looks like a type error, but I guess FEniCS uses

Re: [petsc-users] Scaling/Preconditioners for Poisson equation

2014-10-01 Thread Jed Brown
Filippo Leonardi filippo.leona...@sam.math.ethz.ch writes: I am actually having hard time figuring out where I am spending my time. Reading the report I am spending time on KSPSolve and PCApply (e+02). Since the number of those operations is well under control. I guess is some

Re: [petsc-users] Scaling/Preconditioners for Poisson equation

2014-10-01 Thread Filippo Leonardi
On Wednesday 01 October 2014 07:23:03 Jed Brown wrote: Filippo Leonardi filippo.leona...@sam.math.ethz.ch writes: I am actually having hard time figuring out where I am spending my time. Reading the report I am spending time on KSPSolve and PCApply (e+02). Since the number of those

Re: [petsc-users] PETSC ERROR: Can only handle MPIU_REAL or MPIU_COMPLEX with ML preconditioner

2014-10-01 Thread Barry Smith
Recommend running under valgrind. http://www.mcs.anl.gov/petsc/documentation/faq.html#valgrind On Oct 1, 2014, at 7:03 AM, Chris Richardson ch...@bpi.cam.ac.uk wrote: On 01/10/2014 11:33, Matthew Knepley wrote: On Wed, Oct 1, 2014 at 3:51 AM, Chris Richardson ch...@bpi.cam.ac.uk wrote:

Re: [petsc-users] PETSC ERROR: Can only handle MPIU_REAL or MPIU_COMPLEX with ML preconditioner

2014-10-01 Thread Matthew Knepley
On Wed, Oct 1, 2014 at 9:32 AM, Barry Smith bsm...@mcs.anl.gov wrote: Recommend running under valgrind. http://www.mcs.anl.gov/petsc/documentation/faq.html#valgrind Barry, it looks like this commit introduced it

[petsc-users] examples of automatic differentiation with PETSc

2014-10-01 Thread Xiangdong
Hello everyone, Could you please point me to examples of using automatic differentiation for Jacobian calculation in PETSc? Chapter 9 of the manual briefly talked about adifor, but seems not quite complete. Thank you. Xiangdong

Re: [petsc-users] examples of automatic differentiation with PETSc

2014-10-01 Thread Barry Smith
Sorry that functionality is no longer in PETSc. On Oct 1, 2014, at 2:51 PM, Xiangdong epsco...@gmail.com wrote: Hello everyone, Could you please point me to examples of using automatic differentiation for Jacobian calculation in PETSc? Chapter 9 of the manual briefly talked about

Re: [petsc-users] examples of automatic differentiation with PETSc

2014-10-01 Thread Xiangdong
Is there any reason that AD is not supported in PETSc? Is it because the ADIFOR/ADIC library itself? Thank you. Xiangdong On Wed, Oct 1, 2014 at 4:09 PM, Barry Smith bsm...@mcs.anl.gov wrote: Sorry that functionality is no longer in PETSc. On Oct 1, 2014, at 2:51 PM, Xiangdong

Re: [petsc-users] Default solver type for ksp?

2014-10-01 Thread Barry Smith
Never guess or rely on possible wrong documentation or other peoples responses -ksp_view or KSPView(ksp,0); Barry On Oct 1, 2014, at 3:06 PM, Saurabh Chawdhary chaw0...@umn.edu wrote: I understand that we can chose on of many solvers from KSPType form here. But if do not set any

Re: [petsc-users] Default solver type for ksp?

2014-10-01 Thread Jed Brown
Saurabh Chawdhary chaw0...@umn.edu writes: I understand that we can chose on of many solvers from KSPType form here http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/KSP/KSPType.html. But if do not set any option to the type which one does it chose by default? GMRES I could

Re: [petsc-users] examples of automatic differentiation with PETSc

2014-10-01 Thread Barry Smith
On Oct 1, 2014, at 3:16 PM, Jed Brown j...@jedbrown.org wrote: Xiangdong epsco...@gmail.com writes: Is there any reason that AD is not supported in PETSc? Is it because the ADIFOR/ADIC library itself? Those libraries are not maintained and have license problems. You're welcome to use

Re: [petsc-users] Default solver type for ksp?

2014-10-01 Thread Saurabh Chawdhary
Thanks a lot guys. I also did -ksp_view and saw it was GMRES. ___ KSP Object: 32 MPI processes type: gmres GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement GMRES: happy breakdown tolerance 1e-30

Re: [petsc-users] examples of automatic differentiation with PETSc

2014-10-01 Thread Gaetan Kenway
HI Xiangdong We have use the Tapenade AD package to compute sparse matrix jacobians for CFD problems. We use the result from the AD to populate a PETSc sparse matrix and then use PETSc to solve the resulting linear system. Hope that helps. Gaetan On Wed, Oct 1, 2014 at 4:23 PM, Barry Smith

Re: [petsc-users] examples of automatic differentiation with PETSc

2014-10-01 Thread Barry Smith
Thanks. Note that us removing the “simple minded” AD stuff from PETSc doesn’t mean we think AD is a bad idea. It can be an effective way to generate Jacobians. Barry On Oct 1, 2014, at 4:09 PM, Gaetan Kenway gaet...@gmail.com wrote: HI Xiangdong We have use the Tapenade AD

Re: [petsc-users] include/finclude/petscsysdef.h and daimag

2014-10-01 Thread Barry Smith
Randy, Thanks for the report. I am not sure why we use the d.. version of any of those functions. I have eliminated them in the branch barry/remove-dreal and merged that to next. If all the testing looks good on it we will put it into maint (the next patch) and master. Barry On Sep