old_temperature_solution = temperature_solution;
if (old_time_step > 0)
{
- TrilinosWrappers::MPI::BlockVector distributed_stokes_solution (stokes_rhs);
- distributed_stokes_solution.block(0).reinit(stokes_solution.block(0),false,true);
- distributed_stokes_solution.block(1).reinit(stokes_solution.block(1),false,true);
- distributed_stokes_solution.sadd (1.+time_step/old_time_step, -time_step/old_time_step,
- old_old_stokes_solution);
- stokes_solution.block(0).reinit(distributed_stokes_solution.block(0), false, true);
- stokes_solution.block(1).reinit(distributed_stokes_solution.block(1), false, true);
-
- TrilinosWrappers::MPI::Vector
- distributed_temperature_solution (temperature_rhs);
- distributed_temperature_solution.reinit(temperature_solution, false, true);
- TrilinosWrappers::MPI::Vector
- distributed_old_old_temperature_solution (temperature_rhs);
- distributed_old_old_temperature_solution.reinit(old_old_temperature_solution, false, true);
-
-
- distributed_temperature_solution.sadd (1.+time_step/old_time_step,
- -time_step/old_time_step,
- distributed_old_old_temperature_solution);
- temperature_solution.reinit(distributed_temperature_solution, false, true);
+ stokes_solution.sadd (1.+time_step/old_time_step, -time_step/old_time_step,
+ old_old_stokes_solution);
+ temperature_solution.sadd (1.+time_step/old_time_step,
+ -time_step/old_time_step,
+ old_old_temperature_solution);
}
}
while (true);