have you tried 'fieldsplit_1_ksp_convergence_test':'skip' ? I think then it only does a V-cycle to precondition the Schur-complement. Sorry, didn't read your email carefully enough, I've only tried those options for HDiv. Eike ________________________________ From: firedrake-bounces@imperial.ac.uk <firedrake-bounces@imperial.ac.uk> on behalf of Justin Chang <jychang48@gmail.com> Sent: 15 October 2015 14:57 To: Firedrake Project Subject: Re: [firedrake] Preconditioner for Least-squares FEM Btw, I am not working in HDiv space. The LSFEM uses equal order CG1 elements (VectorFunctionSpace(...) for v and FunctionSpace(...) for p). I tried every combination of this: -ksp_type gmres/cg -pc_type fieldsplit -pc_fieldsplit_type schur -pc_fieldsplit_schur_precondition selfp -fieldsplit_0_ksp_type preonly/cg -fieldsplit_0_pc_type gamg -fieldsplir_1_ksp_type preonly/cg -fieldsplit_1_pc_type gamg But get PETSc errors (Segmentation Violation). Perhaps I am not doing this right? Eike's, the parameters you suggested still give me the DIVERGED_MAX_ITER error. Colin, what do you mean by hybridizing the equations? Thanks all for your input, Justin On Thu, Oct 15, 2015 at 7:32 AM, Colin Cotter <colin.cotter@imperial.ac.uk<mailto:colin.cotter@imperial.ac.uk>> wrote: This is a tricky system to precondition in general. Approaches that are known to work well are: 1) Use a block H(div) preconditioner, i.e. precondition by the operator ( 1 - grad(div .) 0 ) ( 0 1 ) This converges fast with direct solvers to obtain the inverse of this operator, but classical iterative methods do not converge well for the top block, and so you need to use Schwartz preconditioners to get a fast parallel implementation. We don't currently have Firedrake support for that. 2) Hybridise the equation. The resulting reduced system behaves well with classical iterative methods or AMG. We have been able to construct the hybridised system using petsc4py. all the best --cjc On 15 October 2015 at 14:18, Andrew McRae <A.T.T.McRae@bath.ac.uk<mailto:A.T.T.McRae@bath.ac.uk>> wrote: As usual, the Firedrake person in the best position to answer your question is Lawrence, but he's currently on holiday. I won't try to answer your question, but in an old email<http://mailman.ic.ac.uk/pipermail/firedrake/2015-July/000953.html>, Lawrence suggested you try using multigrid on the blocks; e.g. 'fieldsplit_0_pc_type': 'gamg', 'fieldsplit_1_pc_type': 'gamg', My linear solver knowledge is pretty dire, so I'll stop talking here :) Andrew On 15 October 2015 at 13:58, Justin Chang <jychang48@gmail.com<mailto:jychang48@gmail.com>> wrote: I did ask them, but it went unanswered for the last 1.5 months... Let me clarify my question a bit though, I tried using -ksp_type cg -pc_type hypre -pc_hypre_type boomeramg, but it is not supported for matnest. I added nest=False but it seems to break my solver somehow. The only thing that works for me is -ksp_type cg -pc_type bjacobi plus MixedVectorSpaceBasis(W, ...) but the number of solver iterations increases with problem size... Thanks, Justin On Thu, Oct 15, 2015 at 6:50 AM, Andrew McRae <A.T.T.McRae@bath.ac.uk<mailto:A.T.T.McRae@bath.ac.uk>> wrote: I think the PETSc mailing list is a much better place to ask questions like this? On 15 October 2015 at 13:47, Justin Chang <jychang48@gmail.com<mailto:jychang48@gmail.com>> wrote: Hi all, I am attempting to solve Darcy's equation: u + grad[p] = 0 div[u] = f The weak form under the least-squares finite element method (LSFEM) looks like this: (u + grad[p]; v + grad[q]) + div[u]*div[v] = (f; div[v]) For H(div) elements like RT0, these options worked nicely: -ksp_type gmres -pc_type fieldsplit -pc_fieldsplit_type schur -pc_fieldsplit_schur_precondition selfp -fieldsplit_0_ksp_type preonly -fieldsplit_0_pc_type bjacobi -fieldsplit_0_sub_pc_type ilu -fieldsplir_1_ksp_type preonly -fieldsplit_1_pc_type hypre but for the above LSFEM, it does not work (get a DIVERGED_MAX_IT error). I heard multigrid methods are good for these types of problems, so how do I tweak the above parameters? Thanks, Justin _______________________________________________ firedrake mailing list firedrake@imperial.ac.uk<mailto:firedrake@imperial.ac.uk> https://mailman.ic.ac.uk/mailman/listinfo/firedrake -- http://www.imperial.ac.uk/people/colin.cotter www.cambridge.org/9781107663916<http://www.cambridge.org/9781107663916> [http://assets.cambridge.org/97811076/63916/cover/9781107663916.jpg] _______________________________________________ firedrake mailing list firedrake@imperial.ac.uk<mailto:firedrake@imperial.ac.uk> https://mailman.ic.ac.uk/mailman/listinfo/firedrake