Re: [petsc-users] explanations on DM_BOUNDARY_PERIODIC

2017-04-28 Thread neok m4700
2:00 Barry Smith <bsm...@mcs.anl.gov>: > > > On Apr 27, 2017, at 12:43 PM, neok m4700 <neok.m4...@gmail.com> wrote: > > > > Hi Matthew, > > > > Thank you for the clarification, however, it is unclear why there is an > additional unknown in the case

Re: [petsc-users] explanations on DM_BOUNDARY_PERIODIC

2017-04-27 Thread neok m4700
. And in the case of dirichlet or neumann bcs, the extremum cell add information to the RHS, they do not appear in the matrix formulation. Hope I was clear enough, thanks 2017-04-27 16:15 GMT+02:00 Matthew Knepley <knep...@gmail.com>: > On Thu, Apr 27, 2017 at 3:46 AM, neok m4700

[petsc-users] explanations on DM_BOUNDARY_PERIODIC

2017-04-27 Thread neok m4700
Hi, I am trying to change my problem to using periodic boundary conditions. However, when I use DMDASetUniformCoordinates on the DA, the spacing changes. This is due to an additional point e.g. in dm/impls/da/gr1.c else if (dim == 2) { if (bx == DM_BOUNDARY_PERIODIC) hx = (xmax-xmin)/(M);

Re: [petsc-users] (no subject)

2016-06-02 Thread neok m4700
Re, Makes sense to read the documentation, I will try with another preconditioners. Thanks for the support. 2016-06-02 18:15 GMT+02:00 Matthew Knepley <knep...@gmail.com>: > On Thu, Jun 2, 2016 at 11:10 AM, neok m4700 <neok.m4...@gmail.com> wrote: > >>

Re: [petsc-users] (no subject)

2016-06-02 Thread neok m4700
etsc-eab7b92 with petsc4py-6e8e093 > > Satish > > > On Thu, 2 Jun 2016, neok m4700 wrote: > > > Hi Matthew, > > > > I've rebuilt petsc // petsc4py with following versions: > > > > 3.7.0 // 3.7.0 => same runtime error > > 00c67f3 // 3.7.1 => fails to bu

Re: [petsc-users] (no subject)

2016-06-02 Thread neok m4700
__pyx_v_cdim, __pyx_v_interp, (&__pyx_v_newdm))); if (unlikely(__pyx_t_4 == -1)) __PYX_ERR(42, 49, __pyx_L1_error) using - numpy 1.11.0 - openblas 0.2.18 - openmpi 1.10.2 Thanks 2016-06-02 16:39 GMT+02:00 Matthew Knepley <knep...@gmail.com>: > On Thu, Jun 2, 2016 at 9:12 AM, neok m4700 <

[petsc-users] (no subject)

2016-06-02 Thread neok m4700
Hi, I built petsc 3.7.1 and petsc4py 3.7.0 (with openmpi 1.10.2) and ran the examples in the demo directory. $ python test_mat_ksp.py => runs as expected (serial) $ mpiexec -np 2 python test_mat_ksp.py => fails with the following output: Traceback (most recent call last): File