Re: [petsc-dev] implementation header incluing ../src

2015-02-19 Thread Matthew Knepley
On Thu, Feb 19, 2015 at 1:42 AM, Lisandro Dalcin dalc...@gmail.com wrote: Please take look to the last line of the following output: $ grep include $PETSC_DIR/include/petsc-private/dmpleximpl.h #include petscmat.h /*I petscmat.h I*/ #include petscdmplex.h /*I

Re: [petsc-dev] [petsc-users] SNESSetFunctionDomainError

2015-02-19 Thread Dmitry Karpeyev
On Thu Feb 19 2015 at 12:59:12 PM Barry Smith bsm...@mcs.anl.gov wrote: On Feb 19, 2015, at 1:56 PM, Dmitry Karpeyev karp...@mcs.anl.gov wrote: On Thu Feb 19 2015 at 12:41:59 PM Barry Smith bsm...@mcs.anl.gov wrote: Yeah, that sounds like a good fix, except for this: we have to

Re: [petsc-dev] [petsc-users] SNESSetFunctionDomainError

2015-02-19 Thread Barry Smith
All the Krylov methods are supposed to use the routines KSP_MatMult(), KSP_PCApply() etc and not the MatMult(), PCApply() etc directly (I'm fixing some that don't as we speak). I propose that these KSP_xxx() routines have code like MatGetFailure(mat,failure); if (failure)

Re: [petsc-dev] CUDA guys, stuff still broken in master/next

2015-02-19 Thread Dominic Meiser
Hi Karl, Have you had a chance to look into this? Cheers, Dominic On 02/19/2015 12:46 PM, Barry Smith wrote: http://ftp.mcs.anl.gov/pub/petsc/nightlylogs/archive/2015/02/18/examples_master_arch-cuda-double_bb-proxy.log = 1,2c1,150 Linear stokes_

Re: [petsc-dev] [petsc-users] SNESSetFunctionDomainError

2015-02-19 Thread Barry Smith
On Feb 19, 2015, at 1:56 PM, Dmitry Karpeyev karp...@mcs.anl.gov wrote: On Thu Feb 19 2015 at 12:41:59 PM Barry Smith bsm...@mcs.anl.gov wrote: Yeah, that sounds like a good fix, except for this: we have to make sure all ranks diverge with this failure so that the user can retry the

[petsc-dev] CUDA guys, stuff still broken in master/next

2015-02-19 Thread Barry Smith
http://ftp.mcs.anl.gov/pub/petsc/nightlylogs/archive/2015/02/18/examples_master_arch-cuda-double_bb-proxy.log = 1,2c1,150 Linear stokes_ solve converged due to CONVERGED_RTOL iterations 16 9.0909e-02 1.9273e-03 2.9842e-04 4.8732e-03 --- Linear

Re: [petsc-dev] implementation header incluing ../src

2015-02-19 Thread Jed Brown
Matthew Knepley knep...@gmail.com writes: On Thu, Feb 19, 2015 at 1:42 AM, Lisandro Dalcin dalc...@gmail.com wrote: Please take look to the last line of the following output: $ grep include $PETSC_DIR/include/petsc-private/dmpleximpl.h #include petscmat.h /*I petscmat.h

Re: [petsc-dev] [petsc-users] SNESSetFunctionDomainError

2015-02-19 Thread Barry Smith
On Feb 19, 2015, at 2:06 PM, Dmitry Karpeyev karp...@mcs.anl.gov wrote: On Thu Feb 19 2015 at 12:59:12 PM Barry Smith bsm...@mcs.anl.gov wrote: On Feb 19, 2015, at 1:56 PM, Dmitry Karpeyev karp...@mcs.anl.gov wrote: On Thu Feb 19 2015 at 12:41:59 PM Barry Smith

[petsc-dev] Naming convention of examples

2015-02-19 Thread Dave May
Is there a naming convention one should follow when introducing new examples? Specifically, in src/vec/vec/examples/tutorials I see the following ex1.c ex10.c ex11.c ex12.c ex15.c ex16.c ex18.c ex19.c ex2.c ex21.c ex3.c ex42.c ex42a.c ex5.c ex6.c ex7.c ex8.c ex9.c If I

Re: [petsc-dev] Naming convention of examples

2015-02-19 Thread Barry Smith
On Feb 19, 2015, at 4:26 AM, Dave May dave.mayhe...@gmail.com wrote: Is there a naming convention one should follow when introducing new examples? Specifically, in src/vec/vec/examples/tutorials I see the following ex1.c ex10.c ex11.c ex12.c ex15.c ex16.c ex18.c ex19.c

Re: [petsc-dev] [petsc-users] SNESSetFunctionDomainError

2015-02-19 Thread Dmitry Karpeyev
I wanted to revive this thread and move it to petsc-dev. This problem seems to be harder than I realized. Suppose MatMult inside KSPSolve() inside SNESSolve() cannot compute a valid output vector. For example, it's a MatMFFD and as part of its function evaluation it has to evaluate an

Re: [petsc-dev] fix for MatGetRowIJ_SeqSBAIJ

2015-02-19 Thread Hong
Stefano : why didn't you merge it to master and main yet? is the wrong patch? Just merged. Sorry, I thought it was merged. Thanks once again, Hong 2015-02-12 20:26 GMT+03:00 h...@aspiritech.org h...@aspiritech.org: Barry: On Feb 12, 2015, at 9:27 AM, h...@aspiritech.org wrote:

Re: [petsc-dev] implementation header incluing ../src

2015-02-19 Thread Matthew Knepley
On Thu, Feb 19, 2015 at 2:39 PM, Jed Brown j...@jedbrown.org wrote: Matthew Knepley knep...@gmail.com writes: On Thu, Feb 19, 2015 at 1:42 AM, Lisandro Dalcin dalc...@gmail.com wrote: Please take look to the last line of the following output: $ grep include

Re: [petsc-dev] implementation header incluing ../src

2015-02-19 Thread Barry Smith
Is it intended for endusers - petschash.h (requires decent documentation and error checking) or for use by PETSc - petsc-private/hashimpl.h On Feb 19, 2015, at 3:47 PM, Matthew Knepley knep...@gmail.com wrote: On Thu, Feb 19, 2015 at 2:39 PM, Jed Brown j...@jedbrown.org wrote: Matthew

Re: [petsc-dev] implementation header incluing ../src

2015-02-19 Thread Matthew Knepley
On Thu, Feb 19, 2015 at 3:54 PM, Barry Smith bsm...@mcs.anl.gov wrote: Is it intended for endusers - petschash.h (requires decent documentation and error checking) or for use by PETSc - petsc-private/hashimpl.h I am fine with the impl version. I will move all the includes to source files

Re: [petsc-dev] CUDA guys, stuff still broken in master/next

2015-02-19 Thread Dominic Meiser
On 02/19/2015 02:31 PM, Karl Rupp wrote: Hi, Have you had a chance to look into this? sorry, no, I was busy with a project proposal. I'll look into it tomorrow. I've been really swamped, too. I'm afraid these failures might be related to the changes we made for asm and bjacobi with GPU

Re: [petsc-dev] implementation header incluing ../src

2015-02-19 Thread Jed Brown
Matthew Knepley knep...@gmail.com writes: On Thu, Feb 19, 2015 at 2:39 PM, Jed Brown j...@jedbrown.org wrote: Matthew Knepley knep...@gmail.com writes: It does not look like anything else uses hash.h. Can we just move it to include/petschash.h? Why does it need to be included by a

Re: [petsc-dev] [petsc-users] SNESSetFunctionDomainError

2015-02-19 Thread Derek Gaston
No problem on our end with MatSetFailure() being collective... we can guarantee that easily. Thanks for working on this! This has been a high priority feature request for a while. Derek On Thu, Feb 19, 2015 at 3:41 PM Barry Smith bsm...@mcs.anl.gov wrote: On Feb 19, 2015, at 2:06 PM, Dmitry