Dear Prof. Arndt,
Thank you very much!
> when the solver fails to converge within 1000 iterations.
This is just because my solver code is :
template <int dim>
void StokesProblem<dim>::solve_preconditioner ()
{
SolverControl solver_control (1000, 1e-8*system_rhs.l2_norm());
SolverGMRES<BlockVector<double> > gmres (solver_control);
gmres.solve (system_matrix, solution, system_rhs,
PreconditionIdentity());//, preconditioner_H
}
when I change "1000" above into larger number, like "100000", the original
code can run. And I rewrite the whole code, and the code can run now.
Thankyou very much!
Best,
Chucui
在 2019年1月16日星期三 UTC+8下午7:42:28，Daniel Arndt写道：
>
> Chucui,
>
>
> Thanks for your quick answer! In my deal.ii (version 8.5.1), the code also
>> fails in the 9th time step. Shall I need to reinstall my deal.ii into the
>> newest version?
>>
>
> Normally, we only fix bugs in the library in the developer version which
> is also very stable unless you are trying to use the newest features.
> However, I would expect that there is a way to write your code in such a
> way that it also works with version 8.5 correctly.
>
>
>
>
>> For the shortest test code, please give me some time. Because I have
>> deleted about 2000 lines. If I delete some other lines, the error
>> information will change.
>>
>
> The faster you do that, the faster we can help you. It is totally up to
> you. ;-)
>
> Best,
> Daniel
>
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
For more options, visit https://groups.google.com/d/optout.
