From 11bf820dd34e763801893da020e57f51229ce713 Mon Sep 17 00:00:00 2001 From: "Roy H. Stogner" Date: Fri, 16 Feb 2018 11:55:23 -0600 Subject: [PATCH] fem_system_ex2 local->distributed solution There doesn't seem to be any way to get this behavior back in reinit() without creating a performance regression: most applications simply never need to move data in that direction. Refs #1593. --- examples/fem_system/fem_system_ex2/fem_system_ex2.C | 8 ++++++-- 1 file changed, 6 insertions(+), 2 deletions(-) diff --git a/examples/fem_system/fem_system_ex2/fem_system_ex2.C b/examples/fem_system/fem_system_ex2/fem_system_ex2.C index 46c12273213..b011034e855 100644 --- a/examples/fem_system/fem_system_ex2/fem_system_ex2.C +++ b/examples/fem_system/fem_system_ex2/fem_system_ex2.C @@ -82,7 +82,9 @@ void setup(EquationSystems & systems, imms.save_initial_mesh(); // Fill global solution vector from local ones - aux_sys.reinit(); + aux_sys.current_local_solution->close(); + *aux_sys.solution = *aux_sys.current_local_solution; + aux_sys.solution->close(); } @@ -119,7 +121,9 @@ void run_timestepping(EquationSystems & systems, GetPot & args) out << "Solving Solid" << std::endl; solid_system.solve(); - aux_system.reinit(); + aux_system.current_local_solution->close(); + *aux_system.solution = *aux_system.current_local_solution; + aux_system.solution->close(); // Carry out the adaptive mesh refinement/coarsening out << "Doing a reinit of the equation systems" << std::endl;