dolfin team mailing list archive

["Benjamin Kehlet" <benjamik@xxxxxxxxxxxxxxxxxxxxx>]

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@xxxxxxxxxx
> http://www.fenics.org/mailman/listinfo/dolfin-dev
>