Calculation of new time step size using max_retries

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • robertf
    Member
    • Jun 2009
    • 68

    Calculation of new time step size using max_retries

    I'm trying to understand how the automatic time-step controller works. In FETimeStepController.cpp I see
    Code:
    m_ddt = (dt) / (m_maxretries + 1)
    However, the description in Section 3.3.2 of the User's Manual seems to imply that it should be
    Code:
    m_ddt = (dt) / (m_maxretries)
    My C++ is very weak so maybe I'm missing something obvious?

    - Robert
  • maas
    Lead Code Developer
    • Nov 2007
    • 3481

    #2
    Hi Robert,

    It's been a long time since I looked at this part of the code, or its documentation, but I think what's in the code is the correct approach. For instance, imagine the max_retries is set to 1. The documentation's description would suggest that the new time step is set to zero, which would not make a lot of sense. On the other hand, the code's algorithm would set the next time step to 1/2 of the previous time step.

    I'll update the documentation accordingly. Thanks for catching this.

    Best,

    Steve
    Department of Bioengineering, University of Utah
    Scientific Computing and Imaging institute, University of Utah

    Comment

    Working...
    X