Done. Benjamin
> 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 >> > _______________________________________________ > DOLFIN-dev mailing list > DOLFIN-dev@fenics.org > http://www.fenics.org/mailman/listinfo/dolfin-dev > _______________________________________________ DOLFIN-dev mailing list DOLFIN-dev@fenics.org http://www.fenics.org/mailman/listinfo/dolfin-dev