Re: [petsc-users] monitoring the convergence of fieldsplit 0 and 1

2017-01-18 Thread Lawrence Mitchell
On 18/01/17 08:40, Klaij, Christiaan wrote: > Barry, > > I've managed to replicate the problem with 3.7.4 > snes/examples/tutorials/ex70.c. Basically I've added > KSPGetTotalIterations to main (file is attached): PCFieldSplitGetSubKSP returns, in the Schur case: MatSchurComplementGet(pc->schur,

Re: [petsc-users] monitoring the convergence of fieldsplit 0 and 1

2017-01-18 Thread Klaij, Christiaan
Barry, I've managed to replicate the problem with 3.7.4 snes/examples/tutorials/ex70.c. Basically I've added KSPGetTotalIterations to main (file is attached): $ diff -u ex70.c.bak ex70.c --- ex70.c.bak2017-01-18 09:25:46.286174830 +0100 +++ ex70.c2017-01-18 09:03:40.904483434 +0100 @@ -669,6

Re: [petsc-users] monitoring the convergence of fieldsplit 0 and 1

2017-01-18 Thread Klaij, Christiaan
Thanks Lawrence, that nicely explains the unexpected behaviour! I guess in general there ought to be getters for the four ksp(A00)'s that occur in the full factorization. Chris dr. ir. Christiaan Klaij | CFD Researcher | Research & Development MARIN | T +31 317 49 33 44 |

Re: [petsc-users] monitoring the convergence of fieldsplit 0 and 1

2017-01-18 Thread Matthew Knepley
On Wed, Jan 18, 2017 at 4:42 AM, Klaij, Christiaan wrote: > Thanks Lawrence, that nicely explains the unexpected behaviour! > > I guess in general there ought to be getters for the four > ksp(A00)'s that occur in the full factorization. Yes, we will fix it. I think that the

[petsc-users] About GMRES Solver in PETSc

2017-01-18 Thread Chih-Hao Chen
Hello, Sorry for this bother. I"ve checked GMRES solver in your libary supports parallel computing in your website: http://www.mcs.anl.gov/petsc/documentation/linearsolvertable.html. Summary of Sparse Linear Solvers

Re: [petsc-users] About GMRES Solver in PETSc

2017-01-18 Thread Hong
Yes. Hong On Wed, Jan 18, 2017 at 2:22 PM, Chih-Hao Chen < chih-hao.ch...@mail.mcgill.ca> wrote: > Hello, > > > Sorry for this bother. > > I"ve checked GMRES solver in your libary supports parallel computing in > your website: > >

Re: [petsc-users] About GMRES Solver in PETSc

2017-01-18 Thread Chih-Hao Chen
Hello, Thanks for the quick reply. So if I specify matrix A and RHS vector b in parallel forms, ie. having them distributed across several cores, when I using the functions like MatMult, VecAYPX and etc., all the operations would be executed in parallel too? Thanks very much. Best,

Re: [petsc-users] About GMRES Solver in PETSc

2017-01-18 Thread Matthew Knepley
On Wed, Jan 18, 2017 at 2:32 PM, Chih-Hao Chen < chih-hao.ch...@mail.mcgill.ca> wrote: > Hello, > > > Thanks for the quick reply. > > So if I specify matrix A and RHS vector b in parallel forms, > > ie. having them distributed across several cores, > > when I using the functions like MatMult,