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

Matt, Just a reminder. With petsc-3.8.4 the issue is still there. Chris dr. ir. Christiaan Klaij | Senior Researcher | Research & Development MARIN | T +31 317 49 33 44 | c.kl...@marin.nl | www.marin.nl

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

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

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

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

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

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

Well, that's it, all the rest was hard coded. Here's the relevant part of the code: CALL PCSetType(pc_system,PCFIELDSPLIT,ierr); CHKERRQ(ierr) CALL PCFieldSplitSetType(pc_system,PC_COMPOSITE_SCHUR,ierr); CHKERRQ(ierr) CALL PCFieldSplitSetIS(pc_system,"0",isgs(1),ierr); CHKERRQ(ierr)

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

Please send all the command line options you use. > On Jan 16, 2017, at 1:47 AM, Klaij, Christiaan wrote: > > Barry, > > Sure, here's the output with: > > -sys_ksp_view -sys_ksp_converged_reason > -sys_fieldsplit_0_ksp_converged_reason

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

Barry, Sure, here's the output with: -sys_ksp_view -sys_ksp_converged_reason -sys_fieldsplit_0_ksp_converged_reason -sys_fieldsplit_1_ksp_converged_reason (In my previous email, I rearranged 0 & 1 for easy summing.) Chris Linear sys_fieldsplit_0_ solve converged due to CONVERGED_RTOL

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

Yes, I would have expected this to work. Could you send the output from -ksp_view in this case? > On Jan 13, 2017, at 3:46 AM, Klaij, Christiaan wrote: > > Barry, > > It's been a while but I'm finally using this function in > 3.7.4. Is it supposed to work with

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

Barry, It's been a while but I'm finally using this function in 3.7.4. Is it supposed to work with fieldsplit? Here's why. I'm solving a Navier-Stokes system with fieldsplit (pc has one velocity solve and one pressure solve) and trying to retrieve the totals like this: CALL

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

Barry, Matt Thanks, the total iteration count is already very useful. Looking forward to the custom monitors as well. Chris From: Barry Smith bsm...@mcs.anl.gov Sent: Saturday, April 11, 2015 12:27 AM To: Klaij, Christiaan Cc: petsc-users@mcs.anl.gov

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

Chris, I have added KSPGetTotalIterations() to the branch barry/add-ksp-total-iterations/master and next. After tests it will go into master Barry On Apr 10, 2015, at 8:07 AM, Klaij, Christiaan c.kl...@marin.nl wrote: Barry, Sure, I can call PCFieldSplitGetSubKSP() to get the