Skip to content

Commit 31589f6

Browse files
authored
Merge pull request #432 from mmorphew/master
GLM: Add 0th iteration output 0.par and rei0 files
2 parents b8871dd + d83ab43 commit 31589f6

1 file changed

Lines changed: 19 additions & 1 deletion

File tree

src/libs/pestpp_common/SVDSolver.cpp

Lines changed: 19 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -249,7 +249,25 @@ ModelRun SVDSolver::solve(RunManagerAbstract &run_manager, TerminationController
249249
return best_upgrade_run;
250250
if (restart_runs) restart_controller.get_restart_option() = RestartController::RestartOption::NONE;
251251
}
252-
252+
if (global_iter_num == 1)
253+
{
254+
stringstream filename;
255+
string complete_filename;
256+
257+
// Write rei0
258+
filename << "rei" << 0;
259+
output_file_writer.write_rei(file_manager.open_ofile_ext(filename.str()), global_iter_num-1,
260+
*(best_upgrade_run.get_obj_func_ptr()->get_obs_ptr()),
261+
best_upgrade_run.get_obs(), *(best_upgrade_run.get_obj_func_ptr()),
262+
best_upgrade_run.get_ctl_pars());
263+
file_manager.close_file(filename.str());
264+
filename.str(""); // reset the stringstream
265+
// Write 0.par
266+
filename << 0 << ".par";
267+
output_file_writer.write_par(file_manager.open_ofile_ext(filename.str()), best_upgrade_run.get_ctl_pars(), *(par_transform.get_offset_ptr()),
268+
*(par_transform.get_scale_ptr()));
269+
file_manager.close_file(filename.str());
270+
}
253271
// Update Regularization weights if REG_FRAC is used
254272
ModelRun prev_run(best_upgrade_run);
255273
os << endl;

0 commit comments

Comments
 (0)