Currently the build_adiabatic_slfm_library method will actually solve the extinction problem before identifying that extinction has occurred. This is wasteful as the extinction solve can be difficult. We should make the compute_steady_state method smarter. A caveat here is that we could "predict" extinction simply due to a nonlinear solve being too aggressive. We would need to be careful here, perhaps rerunning the solve more conservatively to ensure that extinction is indeed expected.
Currently the
build_adiabatic_slfm_librarymethod will actually solve the extinction problem before identifying that extinction has occurred. This is wasteful as the extinction solve can be difficult. We should make thecompute_steady_statemethod smarter. A caveat here is that we could "predict" extinction simply due to a nonlinear solve being too aggressive. We would need to be careful here, perhaps rerunning the solve more conservatively to ensure that extinction is indeed expected.