Re: [petsc-users] FETI-DP implementation and call sequence

2014-09-05 Thread Alp Kalpalp
As I said before, I have checked out master branch and merged with your stefano_zampini/pcbddc-primalfixe branch. configured, compiled successfully. I have used --with-pcbddc option in configure as Barry suggested. However, tests are failed with following reason: akalpalp@a-kalpalp ~/petsc $

Re: [petsc-users] FETI-DP implementation and call sequence

2014-09-05 Thread Matthew Knepley
Please send your configure.log and make.log Thanks, Matt On Fri, Sep 5, 2014 at 4:27 AM, Alp Kalpalp alpkalp...@gmail.com wrote: As I said before, I have checked out master branch and merged with your stefano_zampini/pcbddc-primalfixe branch. configured, compiled successfully. I have

Re: [petsc-users] fieldsplit_0_ monitor in combination with selfp

2014-09-05 Thread Matthew Knepley
On Fri, Sep 5, 2014 at 1:34 AM, Klaij, Christiaan c.kl...@marin.nl wrote: Matt, I think the problem is somehow related to -pc_fieldsplit_schur_precondition selfp. In the example below your are not using that option. Here is the selfp output. It retains the A00 solver. ex62 -run_type

Re: [petsc-users] fieldsplit_0_ monitor in combination with selfp

2014-09-05 Thread Klaij, Christiaan
Thanks! I've spotted another difference: you are setting the fieldsplit_0_ksp_type and I'm not, just relying on the default instead. If I add -fieldsplit_0_ksp_type gmres then is also get the correct answer. Probably, you will get my problem if you remove -fieldsplit_velocity. mpiexec -n 2 ./ex70

[petsc-users] Putting petsc in a namespace

2014-09-05 Thread Florian Lindner
Hello, This may be rather a C/C++ question, but ... I encapsulate some petsc functions into c++ classes. Since I don't want to pull all petsc symbols into the global namespace for anyone using my classes I try to put petsc into it's own namespace: Header petsc.h: namespace petsc { #include

[petsc-users] Using MUMPS and (PT)SCOTCH with PETSC

2014-09-05 Thread Evan Um
Dear PETSC users, I tried to use SCOTCH 5.1.12b in my PETSC codes since MUMPS has compatibility issues with the latest SCOTCH library. I was told that SCOTCH/6.0.0 that comes with PETSC/3.5.0 is automatically downloaded and installed. Is it still possible to use old SCOTCH library 5.1.2b in

Re: [petsc-users] fieldsplit_0_ monitor in combination with selfp

2014-09-05 Thread Matthew Knepley
On Fri, Sep 5, 2014 at 7:31 AM, Klaij, Christiaan c.kl...@marin.nl wrote: Thanks! I've spotted another difference: you are setting the fieldsplit_0_ksp_type and I'm not, just relying on the default instead. If I add -fieldsplit_0_ksp_type gmres then is also get the correct answer. Probably,

Re: [petsc-users] Putting petsc in a namespace

2014-09-05 Thread Jed Brown
Florian Lindner mailingli...@xgm.de writes: Hello, This may be rather a C/C++ question, but ... I encapsulate some petsc functions into c++ classes. Since I don't want to pull all petsc symbols into the global namespace for anyone using my classes I try to put petsc into it's own

Re: [petsc-users] How to run SNES ex62

2014-09-05 Thread Matthew Knepley
On Wed, Aug 13, 2014 at 2:35 PM, Justin Chang jychan...@gmail.com wrote: Hi all, This might seem like a silly question, but whenever I try running ./ex62 i seem to always get a zero solution. I didn't see any script for running ex62 in the makefile so I have tried all the runtime combination

Re: [petsc-users] Using MUMPS and (PT)SCOTCH with PETSC

2014-09-05 Thread Barry Smith
Set a new value for PETSC_ARCH Obtain the Scotch 5.1.12b tar ball. Used ./configure —download-ptscotch=nameoftarball.tar.gz —download-mumps —download-scalapack etc It will use the provided tar ball instead of downloading 6.0.0 Send errors to petsc-ma...@mcs.anl.gov Barry

Re: [petsc-users] Using MUMPS and (PT)SCOTCH with PETSC

2014-09-05 Thread Kai Song
Hi Barry, Thanks for the suggestion. I downloaded the scotch tar ball in the petsc source directory, and set up the configure flag as suggested, but I got the following error: make[1]: Leaving directory

Re: [petsc-users] Using MUMPS and (PT)SCOTCH with PETSC

2014-09-05 Thread Barry Smith
You need the https://gforge.inria.fr/frs/download.php/file/28934/scotch_5.1.12a_esmumps.tar.gz version with mumps in the tar ball name Barry On Sep 5, 2014, at 2:17 PM, Kai Song ks...@lbl.gov wrote: Hi Barry, Thanks for the suggestion. I downloaded the scotch tar ball in the petsc

Re: [petsc-users] Using MUMPS and (PT)SCOTCH with PETSC

2014-09-05 Thread Kai Song
Hi Barry, I got the similar error for scotch_5.1.12b_esmumps.tar.gz: = === Trying to download file://scotch_5.1.12b_esmumps.tar.gz for PTSCOTCH

Re: [petsc-users] Using MUMPS and (PT)SCOTCH with PETSC

2014-09-05 Thread Jed Brown
Barry Smith bsm...@mcs.anl.gov writes: You need the https://gforge.inria.fr/frs/download.php/file/28934/scotch_5.1.12a_esmumps.tar.gz version with mumps in the tar ball name And yes, this is absurd, so don't forget to complain to upstream. One can only hope that overwhelming user

Re: [petsc-users] Using MUMPS and (PT)SCOTCH with PETSC

2014-09-05 Thread Jed Brown
Kai Song ks...@lbl.gov writes: Hi Barry, I got the similar error for scotch_5.1.12b_esmumps.tar.gz: Looks like this commit was only valid with 6.0, so you'll have to revert it to use 5.1.12b. https://bitbucket.org/petsc/petsc/commits/4620623451f619bc2d23a7c3de5bdf5ee5fd0ff2 Is there a

[petsc-users] FW: LaPlacian 3D example 34.c

2014-09-05 Thread Alletto, John M
All, LaPlacian 3D example 34.c - how do you invoke multiGrid ? In general I am also having to try by trial and error which solvers and preconditioners support 3D objects. Is there a list somewhere? Respectfullt John Alletto

Re: [petsc-users] FW: LaPlacian 3D example 34.c

2014-09-05 Thread Jed Brown
Alletto, John M john.m.alle...@lmco.com writes: All, LaPlacian 3D example 34.c - how do you invoke multiGrid ? Did you read my comment on your scicomp post? I recommend using src/ksp/ksp/examples/tutorials/ex45.c since it was written by us (Barry) and is more controllable via run-time

Re: [petsc-users] FETI-DP implementation and call sequence

2014-09-05 Thread Satish Balay
perhaps unrelated - but are you sure you had the latest master? PETSC_VERSION_RELEASE0 PETSC_VERSION_MAJOR 3 PETSC_VERSION_MINOR 5 PETSC_VERSION_SUBMINOR 0 PETSC_VERSION_PATCH 0 PETSC_VERSION_PATCH should be 1 CLINKER

Re: [petsc-users] FETI-DP implementation and call sequence

2014-09-05 Thread Satish Balay
On Wed, 3 Sep 2014, Stefano Zampini wrote: FETIDP is in PETSc as a byproduct of the BDDC preconditioner (one is the dual of the other) and it does not have its own classes so far. That said, you can experiment with FETI-DP only after having setup a BDDC preconditioner with the options and

Re: [petsc-users] FETI-DP implementation and call sequence

2014-09-05 Thread Satish Balay
Ok - looks like this error check in ISLocalToGlobalMappingRestoreBlockInfo() was added a couple of days back. https://bitbucket.org/petsc/petsc/commits/cbc1caf078fb2bf42b82e0b5ac811b1101900405 PetscValidHeaderSpecific(mapping,IS_LTOGM_CLASSID,1); This is breaking PCISDestroy() - which is

Re: [petsc-users] FETI-DP implementation and call sequence

2014-09-05 Thread Jed Brown
Satish Balay ba...@mcs.anl.gov writes: Ok - looks like this error check in ISLocalToGlobalMappingRestoreBlockInfo() was added a couple of days back. https://bitbucket.org/petsc/petsc/commits/cbc1caf078fb2bf42b82e0b5ac811b1101900405 PetscValidHeaderSpecific(mapping,IS_LTOGM_CLASSID,1);

Re: [petsc-users] FETI-DP implementation and call sequence

2014-09-05 Thread Satish Balay
On Fri, 5 Sep 2014, Jed Brown wrote: Satish Balay ba...@mcs.anl.gov writes: Ok - looks like this error check in ISLocalToGlobalMappingRestoreBlockInfo() was added a couple of days back. https://bitbucket.org/petsc/petsc/commits/cbc1caf078fb2bf42b82e0b5ac811b1101900405

Re: [petsc-users] FETI-DP implementation and call sequence

2014-09-05 Thread Jed Brown
Satish Balay ba...@mcs.anl.gov writes: Perhaps the following is the fix [with proper comments, more error checks?]. But someone more familiar with this code should check this.. Satish -- $ git diff |cat diff --git a/src/ksp/pc/impls/is/pcis.c b/src/ksp/pc/impls/is/pcis.c index

[petsc-users] FE discretization in DMPlex

2014-09-05 Thread Justin Chang
Hi all, So I understand how the FEM code works in the DMPlex examples (ex12 and 62). Pardon me if this is a silly question. 1) If I wanted to solve either the poisson or stokes using the discontinuous Galerkin method, is there a way to do this with the built-in DMPlex/FEM functions? Basically