Re: [petsc-users] gear tools being molded

2013-10-30 Thread Mila
Dear Sir or Madam, My name is Bill from WISDOM, which is specialized in product designing , mold making, plastic silicone parts and assembly. Although you are super in this line, if you need any help please feel free to tell me , WISDOM entire team will do the best to quote our best price and

Re: [petsc-users] ILUTP

2013-10-30 Thread Torquil Macdonald Sørensen
Thanks Matthew! The problem I'm working on is the Dirac equation, in various number of dimensions. I'm going to take a stab at multigrid preconditioning. Best regards Torquil Sørensen On 27 October 2013 14:48, Matthew Knepley knep...@gmail.com wrote: On Sat, Oct 26, 2013 at 3:05 PM, Torquil

Re: [petsc-users] ILUTP

2013-10-30 Thread Matthew Knepley
On Wed, Oct 30, 2013 at 8:58 AM, Torquil Macdonald Sørensen torq...@gmail.com wrote: Thanks Matthew! The problem I'm working on is the Dirac equation, in various number of dimensions. I'm going to take a stab at multigrid preconditioning. My limited knowledge of the Dirac Equation tells me

Re: [petsc-users] DMPlex example with manual grid layout

2013-10-30 Thread Bernhard Reinhardt
On 28.10.2013 16:43, Matthew Knepley wrote:: On Mon, Oct 28, 2013 at 6:17 AM, Bernhard Reinhardt b.reinha...@physik.uni-muenchen.de mailto:b.reinha...@physik.uni-muenchen.de wrote: Hi there, I am trying to familiarize myself with unstructured grids. For this I started with the

Re: [petsc-users] DMPlex example with manual grid layout

2013-10-30 Thread Matthew Knepley
On Wed, Oct 30, 2013 at 11:13 AM, Bernhard Reinhardt b.reinha...@physik.uni-muenchen.de wrote: On 28.10.2013 16:43, Matthew Knepley wrote:: On Mon, Oct 28, 2013 at 6:17 AM, Bernhard Reinhardt b.reinha...@physik.uni-**muenchen.deb.reinha...@physik.uni-muenchen.de

Re: [petsc-users] ILUTP

2013-10-30 Thread Xiaoye S. Li
The sequential SuperLU has ILUTP implementation, not in parallel versions. PETSc already supports the option of using SuperLU. In SuperLU distribution: EXAMPLE/zitersol.c : an example to use GMRES with ILUTP preconditioner (returned from driver SRC/zgsisx.c) SRC/zgsitrf.c : the actual ILUTP

[petsc-users] Compiling PETSc 3.4.3 with ICC 14.0 and OpenMPI 1.6.5

2013-10-30 Thread Maxime Boissonneault
Hi, I am currently trying to compile PETSc 3.4.3 with ICC 14.0 and OpenMPI 1.6.5 and it fails SuperLU_Dist. I am getting the following error : 68303 /software6/mpi/openmpi/1.6.5_intel/bin/mpicc -I/software6/mpi/openmpi/1.6.5_intel/include -I/software6/mpi/openmpi/1.6.5_intel/include -O2

[petsc-users] Seeking an example for TS

2013-10-30 Thread Gautam Bisht
Hi, I'm wondering if there is a TS example for a similar system: d/dt [f(x)] = g(x) For example, in PFLOTRAN, we have the following governing equation for subsurface flow: d/dt [ sat(P)*den(P)*por] = - g(P) where, sat = saturation, den = density, por = porosity, P = Pressure PFLOTRAN uses

Re: [petsc-users] Compiling PETSc 3.4.3 with ICC 14.0 and OpenMPI 1.6.5

2013-10-30 Thread Matthew Knepley
On Wed, Oct 30, 2013 at 3:39 PM, Maxime Boissonneault maxime.boissonnea...@calculquebec.ca wrote: Hi, I am currently trying to compile PETSc 3.4.3 with ICC 14.0 and OpenMPI 1.6.5 and it fails SuperLU_Dist. I am getting the following error : 68303

Re: [petsc-users] Seeking an example for TS

2013-10-30 Thread Jed Brown
Gautam Bisht gbi...@lbl.gov writes: Hi, I'm wondering if there is a TS example for a similar system: d/dt [f(x)] = g(x) You can apply chain rule to produce f'(x) dx/dt = g(x) and then implement TSSetIFunction and (optionally) TSSetIJacobian. pgpvD80y26QKQ.pgp Description: PGP

Re: [petsc-users] Seeking an example for TS

2013-10-30 Thread Matthew Knepley
On Wed, Oct 30, 2013 at 3:47 PM, Gautam Bisht gbi...@lbl.gov wrote: Hi, I'm wondering if there is a TS example for a similar system: d/dt [f(x)] = g(x) For example, in PFLOTRAN, we have the following governing equation for subsurface flow: d/dt [ sat(P)*den(P)*por] = - g(P) where,

Re: [petsc-users] Compiling PETSc 3.4.3 with ICC 14.0 and OpenMPI 1.6.5

2013-10-30 Thread Jed Brown
Maxime Boissonneault maxime.boissonnea...@calculquebec.ca writes: Hi, I am currently trying to compile PETSc 3.4.3 with ICC 14.0 and OpenMPI 1.6.5 and it fails SuperLU_Dist. I am getting the following error : 68303 /software6/mpi/openmpi/1.6.5_intel/bin/mpicc

Re: [petsc-users] Compiling PETSc 3.4.3 with ICC 14.0 and OpenMPI 1.6.5

2013-10-30 Thread Matthew Knepley
On Wed, Oct 30, 2013 at 4:41 PM, Maxime Boissonneault maxime.boissonnea...@calculquebec.ca wrote: Hi, Here is the attached configure.log file. I will try the flag --with-debugging tomorrow. Get rid of -O2: CFLAGS=-O2 -xHost -mkl -fPIC -m64 FFLAGS=-O2 -xHost -mkl -fPIC -m64

Re: [petsc-users] Compiling PETSc 3.4.3 with ICC 14.0 and OpenMPI 1.6.5

2013-10-30 Thread Maxime Boissonneault
Le 2013-10-30 17:45, Matthew Knepley a écrit : On Wed, Oct 30, 2013 at 4:41 PM, Maxime Boissonneault maxime.boissonnea...@calculquebec.ca mailto:maxime.boissonnea...@calculquebec.ca wrote: Hi, Here is the attached configure.log file. I will try the flag --with-debugging

Re: [petsc-users] Compiling PETSc 3.4.3 with ICC 14.0 and OpenMPI 1.6.5

2013-10-30 Thread Matthew Knepley
On Wed, Oct 30, 2013 at 4:48 PM, Maxime Boissonneault maxime.boissonnea...@calculquebec.ca wrote: Le 2013-10-30 17:45, Matthew Knepley a écrit : On Wed, Oct 30, 2013 at 4:41 PM, Maxime Boissonneault maxime.boissonnea...@calculquebec.ca wrote: Hi, Here is the attached configure.log

Re: [petsc-users] Compiling PETSc 3.4.3 with ICC 14.0 and OpenMPI 1.6.5

2013-10-30 Thread Maxime Boissonneault
Le 2013-10-30 17:50, Matthew Knepley a écrit : On Wed, Oct 30, 2013 at 4:48 PM, Maxime Boissonneault maxime.boissonnea...@calculquebec.ca mailto:maxime.boissonnea...@calculquebec.ca wrote: Le 2013-10-30 17:45, Matthew Knepley a écrit : On Wed, Oct 30, 2013 at 4:41 PM, Maxime

Re: [petsc-users] Compiling PETSc 3.4.3 with ICC 14.0 and OpenMPI 1.6.5

2013-10-30 Thread Matthew Knepley
On Wed, Oct 30, 2013 at 4:55 PM, Maxime Boissonneault maxime.boissonnea...@calculquebec.ca wrote: Le 2013-10-30 17:50, Matthew Knepley a écrit : On Wed, Oct 30, 2013 at 4:48 PM, Maxime Boissonneault maxime.boissonnea...@calculquebec.ca wrote: Le 2013-10-30 17:45, Matthew Knepley a écrit

Re: [petsc-users] Could not find make when building PETSc in Windows 7

2013-10-30 Thread Satish Balay
On Mon, 28 Oct 2013, Satish Balay wrote: /usr/bin/ar: petscvecmod.o: No such file or directory Mostlikely this is related to using 'winzip' or '7zip' to untar petsc sources. Try using 'tar' from cygwin. Or you can ignore this issue if not using f90 modules from your code. Could you

[petsc-users] How to force the code run with a given number of iterations?

2013-10-30 Thread Fande Kong
Hi all, How to force the code run with a given number of iterations? I use the following script to run code, but the code stop after doing *two*iterations: * mpirun -n 2 ./linearElasticity -ksp_type gmres -pc_type asm -pc_asm_overlap 1 -sub_pc_type ilu -sub_pc_factor_levels 2 -sub_ksp_type

Re: [petsc-users] How to force the code run with a given number of iterations?

2013-10-30 Thread Jed Brown
Fande Kong fd.k...@siat.ac.cn writes: Hi all, How to force the code run with a given number of iterations? I use the following script to run code, but the code stop after doing *two*iterations: * mpirun -n 2 ./linearElasticity -ksp_type gmres -pc_type asm -pc_asm_overlap 1 -sub_pc_type

Re: [petsc-users] How to force the code run with a given number of iterations?

2013-10-30 Thread Fande Kong
Thanks, I got this message *Linear solve converged due to CONVERGED_RTOL iterations 2*. What is the definition of RTOL in petsc? I doubt the RTOL is small.* * On Wed, Oct 30, 2013 at 6:38 PM, Jed Brown jedbr...@mcs.anl.gov wrote: Fande Kong fd.k...@siat.ac.cn writes: Hi all, How

Re: [petsc-users] How to force the code run with a given number of iterations?

2013-10-30 Thread Jed Brown
Fande Kong fd.k...@siat.ac.cn writes: Thanks, I got this message *Linear solve converged due to CONVERGED_RTOL iterations 2*. What is the definition of RTOL in petsc? I doubt the RTOL is small.* * On Wed, Oct 30, 2013 at 6:38 PM, Jed Brown jedbr...@mcs.anl.gov wrote: Fande Kong

[petsc-users] MatMatSolve and mumps

2013-10-30 Thread Jose David Bermeol
Hi, will it work if I call MatMatSolve and the solver is mumps?? Thanks

Re: [petsc-users] [petsc-dev] MatMatSolve and mumps

2013-10-30 Thread Jed Brown
Jose David Bermeol jberm...@purdue.edu writes: Hi, will it work if I call MatMatSolve and the solver is mumps?? Yes pgp2OjtId8Rok.pgp Description: PGP signature

Re: [petsc-users] [petsc-dev] MatMatSolve and mumps

2013-10-30 Thread Jose David Bermeol
So for which solvers MatMatSolve does not work?? Thanks - Original Message - From: Jed Brown jedbr...@mcs.anl.gov To: Jose David Bermeol jberm...@purdue.edu, petsc-dev petsc-...@mcs.anl.gov, petsc-users petsc-users@mcs.anl.gov Sent: Wednesday, October 30, 2013 9:01:48 PM Subject: Re:

Re: [petsc-users] Could not find make when building PETSc in Windows 7

2013-10-30 Thread Qin Lu
Satish, I just got chance to try your two suggestions and the errors are gone! Many thanks, Qin On Wednesday, October 30, 2013 6:59 PM, Satish Balay ba...@mcs.anl.gov wrote: On Mon, 28 Oct 2013, Satish Balay wrote:   /usr/bin/ar: petscvecmod.o: No such file or directory Mostlikely

Re: [petsc-users] Could not find make when building PETSc in Windows 7

2013-10-30 Thread Satish Balay
Thanks for the confirmation. We'll try to add in a check for this in configure. Satish On Wed, 30 Oct 2013, Qin Lu wrote: Satish, I just got chance to try your two suggestions and the errors are gone! Many thanks, Qin On Wednesday, October 30, 2013 6:59 PM, Satish Balay

[petsc-users] Fwd: MatZeroRowsColumns for mpibaij

2013-10-30 Thread Rhodri Davies
Dear all, Would you mind giving us an update on your thoughts re: this question? We're very keen to be able to use such functionally and it would be good to know if you're looking into this / have any suggestions? Best wishes, Rhodri