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

2014-09-03 Thread Alp Kalpalp
Hi, I need to investigate the performance of FETI-DP on a heteregenous problem with a proposed preconditioner/scaling and different parameters. I do not want to reinvent the known. I need to investigate the outcomes of my research on top of FETI-DP. I have seen several papers in literature

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

2014-09-03 Thread Matthew Knepley
On Wed, Sep 3, 2014 at 3:10 AM, Alp Kalpalp alpkalp...@gmail.com wrote: Hi, I need to investigate the performance of FETI-DP on a heteregenous problem with a proposed preconditioner/scaling and different parameters. I do not want to reinvent the known. I need to investigate the outcomes of

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

2014-09-03 Thread KarinNiKo
You could perhaps have a look at the PCBDDC : http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/PC/PCBDDC.html contributed by Stefano Zampini. Best, Nicolas 2014-09-03 13:09 GMT+02:00 Matthew Knepley knep...@gmail.com: On Wed, Sep 3, 2014 at 3:10 AM, Alp Kalpalp

Re: [petsc-users] PCFieldSplitSetSchurPre in fortran

2014-09-03 Thread Klaij, Christiaan
Matt, Thanks for the fix. If I understand correctly, in an existing install of petsc-3.5.1, I would only need to replace the file finclude/petscpc.h by the new file for the fix to work? (instead of downloading dev, configuring, installing on various machines). Chris MARIN news: Bas Buchner

Re: [petsc-users] PCFieldSplitSetSchurPre in fortran

2014-09-03 Thread Matthew Knepley
On Wed, Sep 3, 2014 at 7:00 AM, Klaij, Christiaan c.kl...@marin.nl wrote: Matt, Thanks for the fix. If I understand correctly, in an existing install of petsc-3.5.1, I would only need to replace the file finclude/petscpc.h by the new file for the fix to work? (instead of downloading dev,

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

2014-09-03 Thread Alp Kalpalp
Dear Matt, While searcing inside of the petsc bundle, I have seen PCBDDCMatFETIDPGetRHS PCBDDCMatFETIDPGetSolution PCBDDCCreateFETIDPOperators also there is an example in petsc-3.5.1\src\ksp\ksp\examples\tutorials\ex59.c however I could not see it on webbrowser interface. Is this an

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

2014-09-03 Thread Matthew Knepley
On Wed, Sep 3, 2014 at 8:51 AM, Alp Kalpalp alpkalp...@gmail.com wrote: Dear Matt, While searcing inside of the petsc bundle, I have seen PCBDDCMatFETIDPGetRHS PCBDDCMatFETIDPGetSolution PCBDDCCreateFETIDPOperators also there is an example in

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

2014-09-03 Thread Jed Brown
Matthew Knepley knep...@gmail.com writes: 1- Is it possible to complete a FETI-DP solution with the provided functions in current PetSc release? There is no FETI-DP in PETSc. Wrong. There is PCBDDC, which has the same eigenvalues as FETI-DP. You can enable it by configuring --with-pcbddc.

Re: [petsc-users] PCFieldSplitSetSchurPre in fortran

2014-09-03 Thread Klaij, Christiaan
Matt, Thanks, after applying the fix to my petsc-3.5.1 install, the small Fortran program works as expected. Now, I would like to change the fortran strategy to the option 3) Using Fortran modules. So, in the small fortran program I replace these seven lines #include finclude/petscsys.h

Re: [petsc-users] PCFieldSplitSetSchurPre in fortran

2014-09-03 Thread Matthew Knepley
On Wed, Sep 3, 2014 at 9:23 AM, Klaij, Christiaan c.kl...@marin.nl wrote: Matt, Thanks, after applying the fix to my petsc-3.5.1 install, the small Fortran program works as expected. Now, I would like to change the fortran strategy to the option 3) Using Fortran modules. So, in the small

Re: [petsc-users] PCFieldSplitSetSchurPre in fortran

2014-09-03 Thread Klaij, Christiaan
I'm sorry, how do I do that? Chris MARIN news: Applied Hydrodynamics of Floating Offshore Structures course, Oct 8 - 10, Houstonhttp://www.marin.nl/web/News/News-items/Applied-Hydrodynamics-of-Floating-Offshore-Structures-course-Oct-8-10-Houston.htm This e-mail may be confidential,

Re: [petsc-users] PCFieldSplitSetSchurPre in fortran

2014-09-03 Thread Matthew Knepley
On Wed, Sep 3, 2014 at 9:43 AM, Klaij, Christiaan c.kl...@marin.nl wrote: I'm sorry, how do I do that? print it Matt Chris MARIN news: Applied Hydrodynamics of Floating Offshore Structures course, Oct 8 - 10, Houston

Re: [petsc-users] PCFieldSplitSetSchurPre in fortran

2014-09-03 Thread Klaij, Christiaan
print *, PC_FIELDSPLIT_SCHUR_PRE_USER gives: 2. (changing to *_SELFP gives: This name does not have a type, and must have an explicit type. [PC_FIELDSPLIT_SCHUR_PRE_SELFP]) MARIN news: MARIN Report 112:

Re: [petsc-users] PCFieldSplitSetSchurPre in fortran

2014-09-03 Thread Matthew Knepley
On Wed, Sep 3, 2014 at 9:46 AM, Klaij, Christiaan c.kl...@marin.nl wrote: print *, PC_FIELDSPLIT_SCHUR_PRE_USER gives: 2. (changing to *_SELFP gives: This name does not have a type, and must have an explicit type. [PC_FIELDSPLIT_SCHUR_PRE_SELFP]) I have

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

2014-09-03 Thread Stefano Zampini
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 customization you prefer. Use

[petsc-users] Loosing a flux when using the asm preconditioned for 2 cores

2014-09-03 Thread Welland, Michael J.
Hi all, I'm simulating a problem with small fluxes, using the asm preconditioner and lu as the sub preconditioner. The simulation runs fine using 2 cores, but when I use more the fluxes disappear and the desired effect goes with them. Does anyone have an idea of a suitable tolerance or

Re: [petsc-users] Loosing a flux when using the asm preconditioned for 2 cores

2014-09-03 Thread Barry Smith
Start by running with -snes_converged_reason -snes_monitor -ksp_converged_reason -ksp_monitor_true_residual -snes_linesearch_monitor on two and then three cores. By default, “not converging” of an iterative solver does not generate an error so the most likely cause is that the iterative

[petsc-users] Loosing a flux when using the asm preconditioned for 2 cores

2014-09-03 Thread Welland, Michael J.
Thanks Barry, I attach the terminal outputs. From my reading, it seems like everything converged alright, no? Usually I get an error message when the linear or snes solver doesn't converge. By the way, I don't think it should matter but for your reference I am doing an operator split method in

Re: [petsc-users] Loosing a flux when using the asm preconditioned for 2 cores

2014-09-03 Thread Matthew Knepley
On Wed, Sep 3, 2014 at 2:05 PM, Welland, Michael J. mwell...@anl.gov wrote: Thanks Barry, I attach the terminal outputs. From my reading, it seems like everything converged alright, no? Usually I get an error message when the linear or snes solver doesn't converge. By the way, I don't

Re: [petsc-users] Loosing a flux when using the asm preconditioned for 2 cores

2014-09-03 Thread Barry Smith
Force really tight convergence on the nonlinear solver -snes_rtol 1.e-12 Note also that initial function norm 0 SNES Function norm 2.418625311837e+01 0 SNES Function norm 3.156023513435e+04 is very different in the two cases. This means either you are solving a different nonlinear

[petsc-users] Postdoctoral position at ANL/MCS (CHiMaD)

2014-09-03 Thread Dmitry Karpeyev
*Postdoctoral position at the Center for Hierarchical Materials Design* We have an opening for a postdoctoral position in the Center for Hierarchical Materials Design (CHiMaD), a joint center between Northwestern University, the University of Chicago, and Argonne National Laboratory, sponsored