@@ -806,7 +806,12 @@ def _get_bos_data_from_restart(self, source_rate, bos_conc):
806806 if res .source_rate != 0.0 :
807807 # Scale reaction rates by ratio of source rates
808808 rates *= source_rate / res .source_rate
809- return bos_conc , OperatorResult (k , rates )
809+
810+ if self ._keff_search_control is not None and source_rate != 0.0 :
811+ keff_search_root = self ._restore_keff_search_control (res )
812+ else :
813+ keff_search_root = None
814+ return bos_conc , OperatorResult (k , rates ), keff_search_root
810815
811816 def _get_start_data (self ) -> tuple [float , int ]:
812817 """
@@ -861,12 +866,8 @@ def _get_bos_data(self, step_index, source_rate, bos_conc):
861866 bos_conc , res = self ._get_bos_data_from_operator (
862867 step_index , source_rate , bos_conc )
863868 else :
864- bos_conc , res = self ._get_bos_data_from_restart (
869+ bos_conc , res , keff_search_root = self ._get_bos_data_from_restart (
865870 source_rate , bos_conc )
866- if self ._keff_search_control is not None and source_rate != 0.0 :
867- keff_search_root = self ._restore_keff_search_control (res )
868- else :
869- keff_search_root = None
870871
871872 return bos_conc , res , keff_search_root
872873
0 commit comments