[dune-pdelab] some problems in dune
Rebecca de Cuveland
rebecca.cuveland at iwr.uni-heidelberg.de
Thu Oct 30 15:00:24 CET 2014
Hey Robert,
perhaps you can show us your program and output? (verbosity 3 or 4 for AMG and 1 for linear solver)
AMG is quite sensible for the used options. I only use Q1, but as far as I know Q2 should also work...
Bye,
Rebecca
--
Rebecca de Cuveland
Universität Heidelberg
Interdisziplinäres Zentrum für Wissenschaftliches Rechnen
Im Neuenheimer Feld 368, 69120 Heidelberg
Tel: 06221/548890 Fax: 06221/548884
e-mail: <rebecca.cuveland at iwr.uni-heidelberg.de>
> On 23.10.2014 10:16, Robert Jendrny wrote:
> Moin moin,
>
> thank you for these options ....
>
> I have one last question:
> If I calculate Poisson with Q1-elements, I get mesh-independent convergence with AMG, i.e. I needed the same iterations for meshsize h, 1/2*h, 1/4*h,....
> But I have problems with Q2-elements: The iteration number "explodes"
>
> How can I fix it? Is it a consequence from mesh-transfer?
>
> It would be nice if you can help me,
>
> Robert
>
>
>> On 10/15/2014 09:22 PM, Rebecca de Cuveland wrote:
>> Hey Robert,
>> there should be another possibility, Markus implemented this some time ago.
>> You can use the parameters class to change the variable, I used it in my code like this:
>>
>> typedef Dune::PDELab::ISTLBackend_BCGS_AMG_ILU0<IGO> LS;
>> LS ls(tpgfs,250,1);
>> Dune::Amg::Parameters params;
>> params.setAccumulate(Dune::Amg::noAccu);
>> ls.setparams(params);
>>
>> Bye,
>> Rebecca
More information about the dune-pdelab
mailing list