Re: [petsc-users] Release canditate?

2015-12-11 Thread Jed Brown
Éric Chamberland writes: > Hi, > > we just waited since petsc 3.5.4 until 3.6.3 to retrieve a working petsc > for our usages (got bugs with 3.6.1, 3.6.2) > > Now, just discovered 3.6.3 have a bug with mkl_pardiso... :/ v3.6.0 is now 6 months old, but the mkl_pardiso stuff was just fixed recentl

Re: [petsc-users] Release canditate?

2015-12-11 Thread Satish Balay
We do autmoatically generate tarballs everynight - its avaliable at: http://ftp.mcs.anl.gov/pub/petsc/petsc-master.tar.gz [but we don't maintain snapshots - i.e old tarballs for this] However adding in 'release' strings is a manual process - so we do this at the release time. [so we would have to

Re: [petsc-users] Release canditate?

2015-12-11 Thread Barry Smith
> On Dec 11, 2015, at 10:35 PM, Éric Chamberland > wrote: > > > > Le 2015-12-11 23:22, Barry Smith a écrit : >> Eric, >> >>Would it be possible for the release candidates to be a git repository >> branch or does it have to be a tarball? Generating and regenerating the > for sure, if y

Re: [petsc-users] Release canditate?

2015-12-11 Thread Éric Chamberland
Le 2015-12-11 23:22, Barry Smith a écrit : Eric, Would it be possible for the release candidates to be a git repository branch or does it have to be a tarball? Generating and regenerating the for sure, if you want a maximum of RC testers, I think it should be a tarball... In my book,

Re: [petsc-users] Release canditate?

2015-12-11 Thread Barry Smith
Eric, Would it be possible for the release candidates to be a git repository branch or does it have to be a tarball? Generating and regenerating the tarball is a time consuming process which is why we don't use release candidates, but if you are able to test off the git repository we coul

[petsc-users] Release canditate?

2015-12-11 Thread Éric Chamberland
Hi, we just waited since petsc 3.5.4 until 3.6.3 to retrieve a working petsc for our usages (got bugs with 3.6.1, 3.6.2) Now, just discovered 3.6.3 have a bug with mkl_pardiso... :/ We will have to wait until 3.6.4 (or patch the source with 3f7bb31

Re: [petsc-users] How to specify how many digits we could print using PetscPrintf

2015-12-11 Thread Barry Smith
> On Dec 11, 2015, at 9:40 PM, Fande Kong wrote: > > Hi all, > > By default, PetscPrintf prints 6 digits. Not really, it uses the standard printf format commands. > How to specify how many digits we could print using PetscPrintf? For example, > I want to print 16 digits for a double value

[petsc-users] How to specify how many digits we could print using PetscPrintf

2015-12-11 Thread Fande Kong
Hi all, By default, PetscPrintf prints 6 digits. How to specify how many digits we could print using PetscPrintf? For example, I want to print 16 digits for a double value (PetscReal). Thanks, Fande,

Re: [petsc-users] SPRNG package

2015-12-11 Thread Barry Smith
I believe that we compile sprng with MPI but PETSc uses sprng only as a sequential random number generator hence I believe the comment is correct and should be left. Barry > On Dec 11, 2015, at 1:51 PM, Gideon Simpson wrote: > > It seems confused. You can compile SPRNG with/without MP

Re: [petsc-users] Solving/creating SPD systems

2015-12-11 Thread Jed Brown
Justin Chang writes: > Pardon me for my apparent lack of understanding over what may be simple > concepts, but why is div[u]*div[v] singular in the context of LSFEM? The corresponding strong form is grad(div(u)). If u has n*d entries (in d dimensions -- it is a vector), div(u) has only n entrie

Re: [petsc-users] Solving/creating SPD systems

2015-12-11 Thread Justin Chang
Pardon me for my apparent lack of understanding over what may be simple concepts, but why is div[u]*div[v] singular in the context of LSFEM? On Fri, Dec 11, 2015 at 12:15 PM, Jed Brown wrote: > Justin Chang writes: > > > Jed, > > > > 1) What exactly are the PETSc options for CGNE? > > -ksp_type

Re: [petsc-users] SPRNG package

2015-12-11 Thread Gideon Simpson
It seems confused. You can compile SPRNG with/without MPI. However, even if you compile it with MPI support, it can still be run as a serial RNG. -gideon > On Dec 11, 2015, at 2:40 PM, Hong wrote: > > Barry : > > there is a comment: > > > >This is NOT currently using a parallel random nu

Re: [petsc-users] SPRNG package

2015-12-11 Thread Hong
Barry : > > > there is a comment: > > > >This is NOT currently using a parallel random number generator. Sprng > does have > >an MPI version we should investigate. > Shall we remove this comment? Hong > > >> On Dec 11, 2015, at 11:30 AM, Hong wrote: > >> > >> Gideon: > >> I was looking a

Re: [petsc-users] SPRNG package

2015-12-11 Thread Barry Smith
To answer your question If you look in config/BuildSystem/config/packages/sprng.py you'll see cflags += ' ' + '-DSPRNG_MPI' # either using MPI or MPIUNI so it looks like --download-sprng does install the MPI version of sprng that you should be able to use the parallel sprng functional

Re: [petsc-users] Solving/creating SPD systems

2015-12-11 Thread Jed Brown
Justin Chang writes: > Jed, > > 1) What exactly are the PETSc options for CGNE? -ksp_type cgne (Conjugate Gradients on the Normal Equations) > Also, would LSQR be worth trying? I am doing all of this through > Firedrake, so I hope these things can be done directly through simply > providing c

Re: [petsc-users] Solving/creating SPD systems

2015-12-11 Thread Justin Chang
Jed, 1) What exactly are the PETSc options for CGNE? Also, would LSQR be worth trying? I am doing all of this through Firedrake, so I hope these things can be done directly through simply providing command line PETSc options :) 2) So i spoke with Matt the other day, and the primary issue I am hav

Re: [petsc-users] Solving/creating SPD systems

2015-12-11 Thread Jed Brown
Justin Chang writes: > So I am wanting to compare the performance of various FEM discretization > with their respective "best" possible solver/pre conditioner. There > are saddle-point systems which HDiv formulations like RT0 work, but then > there are others like LSFEM that are naturally SPD and

Re: [petsc-users] SPRNG package

2015-12-11 Thread Gideon Simpson
Ok, but if I look at http://www.mcs.anl.gov/petsc/petsc-dev/src/sys/classes/random/impls/sprng/sprng.c there is a comment: This is NOT currently using a parallel random number generator. Sprng does have an MPI version we should investigate. -gideon > On Dec 11, 2015, at 11:30 AM, Hong

Re: [petsc-users] SPRNG package

2015-12-11 Thread Hong
Gideon: > I was looking at the source files and noticed in the comments that when > petsc is built with sprng, the petsc interface isn’t to the parallel sprng > RNG. But is the package built with MPI, so that I could manually use the > parallel RNG? > Petsc- sprng-1.0 interface was written many y

Re: [petsc-users] user-defined snes stopping criteria

2015-12-11 Thread anton
SNESSetConvergenceTest Best, Anton On 12/11/2015 04:29 PM, Xiangdong wrote: Hello everyone, Instead of using the standard atol, rtol, stol etc, can user define a new stopping criterion for snes? For example, if I want to use a specialized norm as my snes stopping criteria, which functions sh

[petsc-users] user-defined snes stopping criteria

2015-12-11 Thread Xiangdong
Hello everyone, Instead of using the standard atol, rtol, stol etc, can user define a new stopping criterion for snes? For example, if I want to use a specialized norm as my snes stopping criteria, which functions should I call to pass this information? Thank you. Best, Xiangdong

[petsc-users] SPRNG package

2015-12-11 Thread Gideon Simpson
I was looking at the source files and noticed in the comments that when petsc is built with sprng, the petsc interface isn’t to the parallel sprng RNG. But is the package built with MPI, so that I could manually use the parallel RNG? -gideon