Commit 0cbc9eee authored by Jean-Noël Grad's avatar Jean-Noël Grad
Browse files

Fix source code in tutorials and benchmarks

parent d7b3e0e8
...@@ -1472,7 +1472,7 @@ void Evaluation< LatticeModel_T >::operator()() ...@@ -1472,7 +1472,7 @@ void Evaluation< LatticeModel_T >::operator()()
{ {
WALBERLA_LOG_RESULT_ON_ROOT( "force acting on cylinder (in dimensionless lattice units of the coarsest grid - evaluated in time step " WALBERLA_LOG_RESULT_ON_ROOT( "force acting on cylinder (in dimensionless lattice units of the coarsest grid - evaluated in time step "
<< forceEvaluationExecutionCount_ << "):\n " << force_ << oss.str() << << forceEvaluationExecutionCount_ << "):\n " << force_ << oss.str() <<
"\ndrag and lift coefficients (including extremas of last " << ( coefficients_[0].size() * checkFrequency_ ) << " time steps):" "\ndrag and lift coefficients (including extrema of last " << ( coefficients_[0].size() * checkFrequency_ ) << " time steps):"
"\n \"real\" area:" "\n \"real\" area:"
"\n c_D: " << cDRealArea << " (min = " << coefficientExtremas_[0].first << ", max = " << coefficientExtremas_[0].second << ")" << "\n c_D: " << cDRealArea << " (min = " << coefficientExtremas_[0].first << ", max = " << coefficientExtremas_[0].second << ")" <<
"\n c_L: " << cLRealArea << " (min = " << coefficientExtremas_[1].first << ", max = " << coefficientExtremas_[1].second << ")" << "\n c_L: " << cLRealArea << " (min = " << coefficientExtremas_[1].first << ", max = " << coefficientExtremas_[1].second << ")" <<
...@@ -2920,10 +2920,10 @@ int main( int argc, char **argv ) ...@@ -2920,10 +2920,10 @@ int main( int argc, char **argv )
"// //\n" "// //\n"
"// Schaefer Turek Benchmark //\n" "// Schaefer Turek Benchmark //\n"
"// //\n" "// //\n"
"// Reference: Schaefer, M. and Turek, S. (1996) 'Benchmark computations of laminar flow around a cylinder (with support //\n" "// Reference: Schaefer, M. and Turek, S. (1996) Benchmark computations of laminar flow around a cylinder (with support //\n"
"// by F. Durst, E. Krause and R. Rannacher), in E. Hirschel (Ed.): Flow Simulation with High-Performance //\n" "// by F. Durst, E. Krause and R. Rannacher), in E. Hirschel (Ed.): Flow Simulation with High-Performance //\n"
"// Computers II. DFG Priority Research Program Results 1993-1995, No. 52 in Notes Numer, Fluid Mech., //\n" "// Computers II. DFG Priority Research Program Results 1993-1995, No. 48 in Notes on Numerical Fluid //\n"
"// pp.547-566, Vieweg, Weisbaden. //\n" "// Mechanics, pp.547-566, Vieweg, Weisbaden. //\n"
"// //\n" "// //\n"
"//////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////" ); "//////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////" );
......
...@@ -104,7 +104,7 @@ We will now integrate the generated classes into a waLBerla application. After a ...@@ -104,7 +104,7 @@ We will now integrate the generated classes into a waLBerla application. After a
#include "CumulantMRTNoSlip.h" #include "CumulantMRTNoSlip.h"
#include "CumulantMRTPackInfo.h" #include "CumulantMRTPackInfo.h"
#include "CumulantMRTSweep.h" #include "CumulantMRTSweep.h"
#include "DensityAndVelocityFieldSetter.h" #include "InitialPDFsSetter.h"
\endcode \endcode
We set up typedef aliases for the generated pack info and the D2Q9 stencil. For the PDF and velocity fields, we use instances of the field::GhostLayerField template. The number of entries of the PDF field is specified by the `Stencil_T::Size` parameter. As our domain is two-dimensional, the velocity at each lattice node is a two-dimensional vector. Thus, we set up the velocity field to have two index dimensions passing the stencil's dimension as a template parameter. Finally, we also define a typedef alias for our generated NoSlip boundary. We set up typedef aliases for the generated pack info and the D2Q9 stencil. For the PDF and velocity fields, we use instances of the field::GhostLayerField template. The number of entries of the PDF field is specified by the `Stencil_T::Size` parameter. As our domain is two-dimensional, the velocity at each lattice node is a two-dimensional vector. Thus, we set up the velocity field to have two index dimensions passing the stencil's dimension as a template parameter. Finally, we also define a typedef alias for our generated NoSlip boundary.
......
...@@ -62,7 +62,7 @@ public: ...@@ -62,7 +62,7 @@ public:
WALBERLA_ROOT_SECTION() WALBERLA_ROOT_SECTION()
{ {
std::ofstream fileLocBottom( filename_.c_str()); std::ofstream fileLocBottom( filename_.c_str());
fileLocBottom << "Time,[Locations on the BOTTOM Wall that Reattachment Occures (Normalized with Step Height) ]" << std::endl; fileLocBottom << "Time,[Locations on the BOTTOM Wall that Reattachment Occurs (Normalized with Step Height) ]" << std::endl;
fileLocBottom.close(); fileLocBottom.close();
} }
} }
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment