[CIG-SHORT] Adaptive time stepping in PyLith

[Brad Aagaard]

Hi all:

Tabrez did indeed find a bug in the adaptive time stepping. I am testing 
a fix right now. The bug is related to the Jacobian of the matrix being 
formed at the wrong time when the time step is not uniform.

Brad



Tabrez Ali wrote:
> - I just ran an example having velocity BC's (total_time=1.0*year)  
> with a fixed time step of 0.1*year and it runs fine.
> 
> - With adaptive time stepping (max_dt=1.0*year; adapt_skip=10) it  
> also runs fine and seems to take a time step of 0.088*year.
> 
> - But if I use max_dt=5.0*year (adapt_skip=10) then things go wrong  
> as the solution fails to converge. Even though it appears as if the  
> time step it takes is 0.088*year (from looking at the "Finishing  
> advancing solution from ..." info).