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 >
signature.asc
Description: Digital signature
_______________________________________________ DOLFIN-dev mailing list [email protected] http://www.fenics.org/mailman/listinfo/dolfin-dev
