dolfin team mailing list archive
-
dolfin team
-
Mailing list archive
-
Message #15240
[Bug 424472] Re: Reaction demo does not converge
** Changed in: dolfin
Milestone: None => 0.93
--
Reaction demo does not converge
https://bugs.launchpad.net/bugs/424472
You received this bug notification because you are subscribed to DOLFIN.
Status in DOLFIN: In Progress
Bug description:
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.
References
