dolfin team mailing list archive
-
dolfin team
-
Mailing list archive
-
Message #15126
Re: [Bug 424472] [NEW] Reaction demo does not converge
Benjamin, can you register on launchpad so I can assign this bug to
you? ;-)
--
Anders
On Fri, Sep 04, 2009 at 05:22:00PM -0000, logg wrote:
> Public bug reported:
>
> The reaction demo does not converge after merge with Benjamin (changeset
> 7a94d1c7bf41).
>
> Before the merge, it converged nicely, but now it seems to diverge and
> runs much slower (additional damping needed). It might just be something
> with the choice of tolerance or the number of iterations.
>
> Earlier output:
> ------------------
>
> Creating ODE of size 5000.
> Solving ODE over the time interval [0.0, 0.01]
> Solving primal problem
> Initializing continuous Galerkin method cG(1).
> Using multi-adaptive fixed-point solver.
> Using discrete tolerance tol = 1e-06.
> Computing transpose (inverse) of dependency pattern.
> Time-stepping over the time interval [0, 0.01]
> *** Warning: Time slab system seems to be diverging.
> Direct fixed-point iteration does not converge.
> Trying diagonally damped fixed-point iteration.
> Rejecting time slab K = 1.000e-05, trying again.
> Rejecting time slab K = 6.969e-06, trying again.
> Rejecting time slab K = 4.679e-06, trying again.
> Rejecting time slab K = 4.417e-06, trying again.
> | Time-stepping |
> |=======|---------------------------------------------------------| 11.1%
> | Time-stepping |
> |==============================|----------------------------------| 47.5%
> | Time-stepping |
> |=================================================================| 100.0%
> Multi-adaptive efficiency index: 24.511
> Average number of global iterations per step: 35.106
> Average number of local iterations per global iteration: 1.957
> Total number of (macro) time steps: 47
> Number of rejected time steps: 4
> Not solving the dual problem as requested.
> ODE solution computed in 2.820 seconds.
>
> New output:
> ----------------
>
> lving ODE over the time interval [0.0, 0.01]
> Solving primal problem
> Initializing continuous Galerkin method cG(1).
> Using multi-adaptive fixed-point solver.
> Using discrete tolerance tol = 1e-06....
> Computing transpose (inverse) of dependency pattern.
> Time-stepping over the time interval [0, 0.01]
> *** Warning: Time slab system seems to be diverging.
> Direct fixed-point iteration does not converge.
> Trying diagonally damped fixed-point iteration.
> Rejecting time slab K = 1.000e-05, trying again.
> Rejecting time slab K = 6.969e-06, trying again.
> Rejecting time slab K = 4.679e-06, trying again.
> Rejecting time slab K = 4.417e-06, trying again.
> | Time-stepping |
> |-----------------------------------------------------------------| 0.0%
> | Time-stepping |
> |-----------------------------------------------------------------| 0.2%
> | Time-stepping |
> |-----------------------------------------------------------------| 0.4%
> | Time-stepping |
> |-----------------------------------------------------------------| 0.6%
> | Time-stepping |
> |-----------------------------------------------------------------| 1.0%
> | Time-stepping |
> |-----------------------------------------------------------------| 1.4%
> | Time-stepping |
> |=|---------------------------------------------------------------| 1.8%
> | Time-stepping |
> |=|---------------------------------------------------------------| 2.1%
> | Time-stepping |
> |=|---------------------------------------------------------------| 2.6%
> | Time-stepping |
> |==|--------------------------------------------------------------| 3.1%
> | Time-stepping |
> |==|--------------------------------------------------------------| 3.7%
> | Time-stepping |
> |==|--------------------------------------------------------------| 4.5%
> | Time-stepping |
> |===|-------------------------------------------------------------| 5.4%
> | Time-stepping |
> |====|------------------------------------------------------------| 6.4%
> | Time-stepping |
> |=====|-----------------------------------------------------------| 7.7%
> | Time-stepping |
> |======|----------------------------------------------------------| 9.3%
> | Time-stepping |
> |=======|---------------------------------------------------------| 11.1%
> | Time-stepping |
> |========|--------------------------------------------------------| 13.3%
> | Time-stepping |
> |==========|------------------------------------------------------| 15.9%
> | Time-stepping |
> |============|----------------------------------------------------| 19.1%
> | Time-stepping |
> |==============|--------------------------------------------------| 22.9%
> | Time-stepping |
> |=================|-----------------------------------------------| 27.4%
> | Time-stepping |
> |=====================|-------------------------------------------| 32.5%
> | Time-stepping |
> |========================|----------------------------------------| 37.5%
> | Time-stepping |
> |===========================|-------------------------------------| 42.5%
> | Time-stepping |
> |==============================|----------------------------------| 47.5%
>
> | Time-stepping |
> |==================================|------------------------------| 52.6%
> | Time-stepping |
> |=====================================|---------------------------| 57.6%
> | Time-stepping |
> |========================================|------------------------| 62.6%
> | Time-stepping |
> |============================================|--------------------| 67.7%
> | Time-stepping |
> |===============================================|-----------------| 73.0%
> | Time-stepping |
> |==================================================|--------------| 78.3%
> | Time-stepping |
> |======================================================|----------| 83.4%
> | Time-stepping |
> |============================================================|----| 93.7%
> | Time-stepping |
> |=================================================================| 100.0%
> Multi-adaptive efficiency index: 24.511
> Average number of global iterations per step: 35.106
> Average number of local iterations per global iteration: 1.957
> Total number of (macro) time steps: 47
> Number of rejected time steps: 4
> Not solving the dual problem as requested.
> ODE solution computed in 113.380 seconds.
>
> ** Affects: dolfin
> Importance: Undecided
> Status: New
>
Attachment:
signature.asc
Description: Digital signature
Follow ups
References
