From: Wolfgang Bangerth Date: Fri, 15 Jan 2010 18:38:50 +0000 (+0000) Subject: Add missing include file. X-Git-Tag: v8.0.0~6621 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=bb6adceb6bee0f5e40feb75810dbf072cc1fe2d9;p=dealii.git Add missing include file. git-svn-id: https://svn.dealii.org/trunk@20374 0785d39b-7218-0410-832d-ea1e28bc413d --- 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(); }