From 57de7d6157dbf7e8981dcca19965f8c5dc15a1b9 Mon Sep 17 00:00:00 2001
From: Tomas Oberhuber <tomas.oberhuber@fjfi.cvut.cz>
Date: Wed, 13 Feb 2013 09:39:41 +0100
Subject: [PATCH] Fixing PDE solver.
---
src/implementation/mesh/tnlGrid2D_impl.h | 2 +-
src/implementation/solvers/pde/tnlPDESolver_impl.h | 7 ++++++-
2 files changed, 7 insertions(+), 2 deletions(-)
diff --git a/src/implementation/mesh/tnlGrid2D_impl.h b/src/implementation/mesh/tnlGrid2D_impl.h
index 4977931659..f626a05f1e 100644
--- a/src/implementation/mesh/tnlGrid2D_impl.h
+++ b/src/implementation/mesh/tnlGrid2D_impl.h
@@ -250,7 +250,7 @@ bool tnlGrid< 2, Real, Device, Index> :: write( const MeshFunction& function,
for( IndexType i = 0; i < getDimensions(). x(); i++ )
{
const RealType x = this -> getLowerCorner(). x() + i * hx;
- const RealType y = this -> getLowerCorner(). y() + i * hy;
+ const RealType y = this -> getLowerCorner(). y() + j * hy;
file << x << " " << " " << y << " " << function[ this -> getNodeIndex( j, i ) ] << endl;
}
file << endl;
diff --git a/src/implementation/solvers/pde/tnlPDESolver_impl.h b/src/implementation/solvers/pde/tnlPDESolver_impl.h
index 7e5917d3a7..54d11bf02d 100644
--- a/src/implementation/solvers/pde/tnlPDESolver_impl.h
+++ b/src/implementation/solvers/pde/tnlPDESolver_impl.h
@@ -90,11 +90,16 @@ bool tnlPDESolver< Problem, TimeStepper > :: solve()
IndexType step( 0 );
IndexType allSteps = ceil( this -> finalTime / this -> snapshotTau );
this -> timeStepper -> setProblem( * ( this -> problem ) );
+ if( ! this -> problem -> makeSnapshot( t, step ) )
+ {
+ cerr << "Making the snapshot failed." << endl;
+ return false;
+ }
while( step < allSteps )
{
RealType tau = Min( this -> snapshotTau,
this -> finalTime - t );
- this -> timeStepper -> setTau( tau );
+ //this -> timeStepper -> setTau( tau );
if( ! this -> timeStepper -> solve( t, t + tau ) )
return false;
step ++;
--
GitLab