From 3b50a2df3303711e8a69d958feb3d631da76b7cc Mon Sep 17 00:00:00 2001 From: Guido Kanschat Date: Thu, 8 Jul 2010 04:50:45 +0000 Subject: [PATCH] make theta_timestepping.cc work and make OutputOperator gnuplot compatible git-svn-id: https://svn.dealii.org/trunk@21462 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/deal.II/include/numerics/operator.templates.h | 5 +++-- deal.II/examples/doxygen/theta_timestepping.cc | 8 +++----- 2 files changed, 6 insertions(+), 7 deletions(-) diff --git a/deal.II/deal.II/include/numerics/operator.templates.h b/deal.II/deal.II/include/numerics/operator.templates.h index 79d4f5980d..cf27cc2a91 100644 --- a/deal.II/deal.II/include/numerics/operator.templates.h +++ b/deal.II/deal.II/include/numerics/operator.templates.h @@ -68,9 +68,10 @@ namespace Algorithms } else { - (*os) << ' ' << step << ' '; + (*os) << ' ' << step; for (unsigned int i=0;iprint(*os); + for (unsigned int j=0;jsize();++j) + (*os) << ' ' << (*vectors(i))(j); (*os) << std::endl; } return *this; diff --git a/deal.II/examples/doxygen/theta_timestepping.cc b/deal.II/examples/doxygen/theta_timestepping.cc index ecf62e2d3f..f4daa804ae 100644 --- a/deal.II/examples/doxygen/theta_timestepping.cc +++ b/deal.II/examples/doxygen/theta_timestepping.cc @@ -61,8 +61,8 @@ class Implicit int main() { FullMatrix matrix(2); - matrix(0,0) = 0.; - matrix(1,1) = 0.; + matrix(0,0) = 1.; + matrix(1,1) = 1.; matrix(0,1) = 31.4; matrix(1,0) = -31.4; @@ -109,10 +109,9 @@ Explicit::operator() (NamedData*>& out, const NamedDatanotifications.test(Events::initial) || this->notifications.test(Events::new_timestep_size)) { - m.equ(timestep_data->step, *matrix); + m.equ(-timestep_data->step, *matrix); for (unsigned int i=0;inotifications.clear(); unsigned int i = in.find("Previous time"); @@ -144,7 +143,6 @@ Implicit::operator() (NamedData*>& out, const NamedDatanotifications.clear(); -- 2.39.5