Re: [petsc-users] jmp_buf error while linking petsc

2017-03-21 Thread Satish Balay
On Tue, 21 Mar 2017, Bikash Kanungo wrote: > Hi Satish, > > Thanks for the suggestions. I will change make makefiles and try to make it > work. However, I would like to state that the above problem appears only > while linking to petsc-3.7.5. My code compiles well with petsc-3.6.3. The >

Re: [petsc-users] jmp_buf error while linking petsc

2017-03-21 Thread Bikash Kanungo
Hi Satish, Thanks for the suggestions. I will change make makefiles and try to make it work. However, I would like to state that the above problem appears only while linking to petsc-3.7.5. My code compiles well with petsc-3.6.3. The configuration options were same for both petsc-3.7.5 and

Re: [petsc-users] Segmentation fault due to TSDestroy

2017-03-21 Thread Barry Smith
Thanks for letting us know, I added some additional error checking (soon to be) in master that will flag the use of incorrect NULL types for function pointers and thus prevent such crashes. Barry > On Mar 20, 2017, at 7:46 AM, Praveen C wrote: > > It turns out the

Re: [petsc-users] [petsc-maint] Question about DMDA BOUNDARY_CONDITION set

2017-03-21 Thread zhaowenbo.n...@gmail.com
Barry,Thanks. It helps me a lot.Five point (stencil width is one) is enough for begin.BEST,Wenbo 原始邮件 主题:Re: [petsc-maint] Question about DMDA BOUNDARY_CONDITION set发件人:Barry Smith 收件人:Wenbo Zhao 抄送:petsc-users@mcs.anl.gov,Matthew

Re: [petsc-users] jmp_buf error while linking petsc

2017-03-21 Thread Satish Balay
Should have said: make getccompiler getincludedirs getlinklibs Satish On Tue, 21 Mar 2017, Satish Balay wrote: > Ok - Then the issue might be with your makefile. Perhaps its not a > petsc formatted makefile. > > In this case - run 'make getincludedirs getlibs' in PETSC_DIR - and > use the

Re: [petsc-users] jmp_buf error while linking petsc

2017-03-21 Thread Satish Balay
Ok - Then the issue might be with your makefile. Perhaps its not a petsc formatted makefile. In this case - run 'make getincludedirs getlibs' in PETSC_DIR - and use the compilers, compile/link options listed there. Satish On Tue, 21 Mar 2017, Bikash Kanungo wrote: > Hi Satish, > > Yes, I ran

Re: [petsc-users] jmp_buf error while linking petsc

2017-03-21 Thread Bikash Kanungo
Hi Satish, Yes, I ran 'make test' and it didn't throw any error. Here's my test.log: Running test examples to verify correct installation Using PETSC_DIR=/home/bikashk/softwares/femdft_softwares/petsc/opt_3.7.5_mvapich2_64bit_complex and PETSC_ARCH= C/C++ example src/snes/examples/tutorials/ex19

Re: [petsc-users] jmp_buf error while linking petsc

2017-03-21 Thread Satish Balay
Did you run 'make test' after insatalling PETSc? Were there errors with this test? [you can send us corresponding test.log] Satish On Tue, 21 Mar 2017, Bikash Kanungo wrote: > Hi, > > I have recently installed petsc-3.7.5. However, while linking petsc with my > code, I get the following

[petsc-users] jmp_buf error while linking petsc

2017-03-21 Thread Bikash Kanungo
Hi, I have recently installed petsc-3.7.5. However, while linking petsc with my code, I get the following compilation error: /home/bikashk/softwares/femdft_softwares/petsc/opt_3.7.5_mvapich2_64bit_complex/include/petscdraw.h(335): error: identifier "jmp_buf" is undefined PETSC_EXTERN jmp_buf

Re: [petsc-users] [petsc-maint] Question about DMDA BOUNDARY_CONDITION set

2017-03-21 Thread Barry Smith
I can cook up a simple example that communicates the needed ghost values for your boundary conditions. It may take a day or two. Are you using a five point stencil or a nine point stencil (the nine point stencil seems "weird" for that very corner point)? Barry > On Mar 21, 2017,

Re: [petsc-users] Question about DMDA BOUNDARY_CONDITION set

2017-03-21 Thread Matthew Knepley
On Tue, Mar 21, 2017 at 1:46 PM, Wenbo Zhao wrote: > Matt, > > Thanks. > > I want to solve neutron diffusion equations using finite difference method > and PETSc. > This rotation boundary condition is very common in my cases. > Though the mesh consists of ~ 10 Miliion

Re: [petsc-users] Configure nested PCFIELDSPLIT with general index sets

2017-03-21 Thread Natacha BEREUX
Thanks for your quick answers. To be honest, I am not familiar at all with DMShells and DMPlexes. But since it is what I need, I am going to try it. Thanks again for your advices, Natacha On Tue, Mar 21, 2017 at 2:27 PM, Lawrence Mitchell < lawrence.mitch...@imperial.ac.uk> wrote: > > > On 21

Re: [petsc-users] Question about DMDA BOUNDARY_CONDITION set

2017-03-21 Thread Wenbo Zhao
Matt, Thanks. I want to solve neutron diffusion equations using finite difference method and PETSc. This rotation boundary condition is very common in my cases. Though the mesh consists of ~ 10 Miliion structured hexahedron cells, the mesh is simple and could be discribed by three vectors about

Re: [petsc-users] Configure nested PCFIELDSPLIT with general index sets

2017-03-21 Thread Lawrence Mitchell
> On 21 Mar 2017, at 13:24, Matthew Knepley wrote: > > I think the remedy is as easy as specifying a DMShell that has a PetscSection > (DMSetDefaultSection) with your ordering, and > I think this is how Firedrake (http://www.firedrakeproject.org/) does it. We actually don't

Re: [petsc-users] Configure nested PCFIELDSPLIT with general index sets

2017-03-21 Thread Matthew Knepley
On Tue, Mar 21, 2017 at 12:58 PM, Natacha BEREUX wrote: > Dear PETSc user's, > I am trying to solve a poroelasticity problem with an additional > temperature equation. The problem is a 3 fields problem involving a > displacement field (u), a pressure field (p) and a

[petsc-users] Configure nested PCFIELDSPLIT with general index sets

2017-03-21 Thread Natacha BEREUX
Dear PETSc user's, I am trying to solve a poroelasticity problem with an additional temperature equation. The problem is a 3 fields problem involving a displacement field (u), a pressure field (p) and a temperature field (t). I have seen similar examples in

Re: [petsc-users] TSAdaptSetType in fortran

2017-03-21 Thread Barry Smith
> On Mar 21, 2017, at 3:39 AM, Praveen C wrote: > > Thank you. > > Will this fix be added to future versions or I have to apply this manually ? I will put it in maint and master as soon as I know it isn't buggy. So within a day or two. Barry > > Best > praveen >

Re: [petsc-users] TSAdaptSetType in fortran

2017-03-21 Thread Praveen C
Thank you. Will this fix be added to future versions or I have to apply this manually ? Best praveen On Tue, Mar 21, 2017 at 1:41 PM, Barry Smith wrote: > > Sorry about that; since it involves character strings so it is something > we currently need to do manually. > >

Re: [petsc-users] TSAdaptSetType in fortran

2017-03-21 Thread Barry Smith
Sorry about that; since it involves character strings so it is something we currently need to do manually. If you are 1) using petsc 3.7 from the tar ball I've attached a patch file, apply it in the petsc directory with patch -p1 < add-fortran-tsadaptsettype 2) if you are using petsc

[petsc-users] TSAdaptSetType in fortran

2017-03-21 Thread Praveen C
Dear all In my fortran code, I want to set to TSADAPTNONE by default. But it seems that TSAdaptSetType does not exist in fortran. Is there a solution for this ? Thanks praveen