From bb6adceb6bee0f5e40feb75810dbf072cc1fe2d9 Mon Sep 17 00:00:00 2001 From: Wolfgang Bangerth Date: Fri, 15 Jan 2010 18:38:50 +0000 Subject: [PATCH] Add missing include file. git-svn-id: https://svn.dealii.org/trunk@20374 0785d39b-7218-0410-832d-ea1e28bc413d --- .../numerics/theta_timestepping.templates.h | 22 ++++++++++--------- 1 file changed, 12 insertions(+), 10 deletions(-) diff --git a/deal.II/deal.II/include/numerics/theta_timestepping.templates.h b/deal.II/deal.II/include/numerics/theta_timestepping.templates.h index 2d324b80b8..485fca03ae 100644 --- a/deal.II/deal.II/include/numerics/theta_timestepping.templates.h +++ b/deal.II/deal.II/include/numerics/theta_timestepping.templates.h @@ -11,6 +11,8 @@ //--------------------------------------------------------------------------- #include + +#include #include DEAL_II_NAMESPACE_OPEN @@ -22,7 +24,7 @@ namespace Algorithms : op_explicit(&e), op_implicit(&i) {} - + template void ThetaTimestepping::notify(const Event& e) @@ -59,14 +61,14 @@ namespace Algorithms ThetaTimestepping::operator() (NamedData& out, const NamedData& in) { Assert(!adaptive, ExcNotImplemented()); - + deallog.push ("Theta"); GrowingVectorMemory mem; typename VectorMemory::Pointer aux(mem); aux->reinit(*out(0)); - + control.restart(); - + bool first = true; d_explicit.time = control.now(); @@ -85,12 +87,12 @@ namespace Algorithms NamedData src2; src2.add(p, "Previous time data"); src2.merge(in); - + if (output != 0) (*output) << 0U << out; - + for (unsigned int count = 1; d_explicit.time < control.final(); ++count) - { + { const bool step_change = control.advance(); d_implicit.time = control.now(); d_explicit.step = (1.-vtheta)*control.step(); @@ -104,15 +106,15 @@ namespace Algorithms op_explicit->notify(Events::new_timestep_size); op_implicit->notify(Events::new_timestep_size); } - + // Compute // (I + (1-theta)dt A) u (*op_explicit)(out1, src1); (*op_implicit)(out, src2); - + if (output != 0 && control.print()) (*output) << count << out; - + first = false; d_explicit.time = control.now(); } -- 2.39.5