temperature_solution = tmp[0];
old_temperature_solution = tmp[1];
+ temperature_constraints.distribute(temperature_solution);
+ temperature_constraints.distribute(old_temperature_solution);
+
stokes_trans.interpolate (x_stokes, stokes_solution);
+ stokes_constraints.distribute(stokes_solution);
+
rebuild_stokes_matrix = true;
rebuild_temperature_matrices = true;
rebuild_stokes_preconditioner = true;
old_saturation_solution = tmp_saturation[1];
saturation_matching_last_computed_darcy_solution = tmp_saturation[2];
+ saturation_constraints.distribute(saturation_solution);
+ saturation_constraints.distribute(old_saturation_solution);
+ saturation_constraints.distribute(saturation_matching_last_computed_darcy_solution);
+
std::vector<TrilinosWrappers::MPI::BlockVector> tmp_darcy (2);
tmp_darcy[0].reinit (darcy_solution);
tmp_darcy[1].reinit (darcy_solution);
last_computed_darcy_solution = tmp_darcy[0];
second_last_computed_darcy_solution = tmp_darcy[1];
+ darcy_constraints.distribute(last_computed_darcy_solution);
+ darcy_constraints.distribute(second_last_computed_darcy_solution);
+
rebuild_saturation_matrix = true;
}
}