THM restart with thermal strains test required.
The following discussion from !3742 (merged) should be addressed: discussion:
Since
eps
is considered in the IP output and input now, for an exact restarting computationeps
,eps_prev
andeps_m_prev
have to be considered in IP output and input () as well because of the thermal strain being computed in the increment manner as:eps_m.noalias() = eps_m_prev + eps - eps_prev - dthermal_strain;
If these IP data are considered in IP I/O, the setting of
eps_prev
andeps_m_prev
have to be avoided in callingIntegrationPointData:: pushBackState()
.
I'd think storing eps_m is sufficient but this is to be checked. Best way would be an example with thermal strains requiring the restart.