30 const std::vector<std::unique_ptr<Process>>& processes,
31 const std::map<std::string, std::unique_ptr<NumLib::NonlinearSolverBase>>&
33 std::vector<std::unique_ptr<MeshLib::Mesh>>& meshes,
34 bool const compensate_non_equilibrium_initial_residuum)
38 if (!output_config_tree)
40 INFO(
"No output section found.");
44 ?
createOutput(*output_config_tree, output_directory, meshes)
47 output_directory, meshes);
48 auto const fixed_times_for_output =
51 if (
auto const submesh_residuum_output_config_tree =
54 submesh_residuum_output_config_tree)
57 *submesh_residuum_output_config_tree, output_directory, meshes);
59 for (
auto& process : processes)
61 auto const& residuum_vector_names =
62 process->initializeAssemblyOnSubmeshes(smroc.meshes);
64 for (
auto const& name : residuum_vector_names)
66 smroc.output.doNotProjectFromBulkMeshToSubmeshes(
71 outputs.push_back(std::move(smroc.output));
76 for (
auto& process : processes)
78 process->initializeAssemblyOnSubmeshes({});
85 compensate_non_equilibrium_initial_residuum, fixed_times_for_output);
87 const bool use_staggered_scheme =
88 ranges::any_of(processes.begin(), processes.end(),
89 [](
auto const& process)
90 { return !(process->isMonolithicSchemeUsed()); });
92 std::unique_ptr<NumLib::StaggeredCoupling> staggered_coupling =
nullptr;
93 if (use_staggered_scheme)
95 staggered_coupling = NumLib::createStaggeredCoupling<ProcessData>(
96 config, per_process_data);
100 if (per_process_data.size() > 1)
103 "The monolithic scheme is used. However more than one "
104 "process data tags (by name \"process\") inside tag "
105 "\"time_loop\" are defined for the staggered scheme. If you "
106 "want to use staggered scheme, please set the element of tag "
107 "\"<coupling_scheme>\" to \"staggered\".");
111 const auto minmax_iter =
112 std::minmax_element(per_process_data.begin(),
113 per_process_data.end(),
114 [](std::unique_ptr<ProcessData>
const& a,
115 std::unique_ptr<ProcessData>
const& b) {
116 return (a->timestep_algorithm->end() <
117 b->timestep_algorithm->end());
119 const double start_time =
120 per_process_data[minmax_iter.first - per_process_data.begin()]
121 ->timestep_algorithm->begin();
122 const double end_time =
123 per_process_data[minmax_iter.second - per_process_data.begin()]
124 ->timestep_algorithm->end();
126 return std::make_unique<TimeLoop>(
127 std::move(outputs), std::move(per_process_data),
128 std::move(staggered_coupling), start_time, end_time);
std::vector< std::unique_ptr< ProcessData > > createPerProcessData(BaseLib::ConfigTree const &config, std::vector< std::unique_ptr< Process > > const &processes, std::map< std::string, std::unique_ptr< NumLib::NonlinearSolverBase > > const &nonlinear_solvers, bool const compensate_non_equilibrium_initial_residuum, std::vector< double > const &fixed_times_for_output)
std::unique_ptr< TimeLoop > createTimeLoop(BaseLib::ConfigTree const &config, std::string const &output_directory, const std::vector< std::unique_ptr< Process > > &processes, const std::map< std::string, std::unique_ptr< NumLib::NonlinearSolverBase > > &nonlinear_solvers, std::vector< std::unique_ptr< MeshLib::Mesh > > &meshes, bool const compensate_non_equilibrium_initial_residuum)
Builds a TimeLoop from the given configuration.