Open shishiousan opened 7 months ago
now you can set SDAlgoParam_
for HHT alpha and Collocation methods by
TYPE(SDAlgoParam_) :: obj
call obj%HHTAlpha(alpha=alpha, beta=beta, gamma=gamma)
call obj%Collocation(beta=beta, gamma=gamma, theta=theta)
Since the Houbolt method needs a somewhat different structure, its implementation is future work.
Implementing some well-known semi discrete time integration schemes