summaryrefslogtreecommitdiff
path: root/apps/adrian/poiseuille2d/poiseuille2d.cpp
diff options
context:
space:
mode:
Diffstat (limited to 'apps/adrian/poiseuille2d/poiseuille2d.cpp')
-rw-r--r--apps/adrian/poiseuille2d/poiseuille2d.cpp12
1 files changed, 5 insertions, 7 deletions
diff --git a/apps/adrian/poiseuille2d/poiseuille2d.cpp b/apps/adrian/poiseuille2d/poiseuille2d.cpp
index d4d8e0f..81870a8 100644
--- a/apps/adrian/poiseuille2d/poiseuille2d.cpp
+++ b/apps/adrian/poiseuille2d/poiseuille2d.cpp
@@ -24,9 +24,7 @@
*/
#include "olb2D.h"
-#ifndef OLB_PRECOMPILED
#include "olb2D.hh"
-#endif
#include <vector>
@@ -34,7 +32,7 @@ using namespace olb;
typedef double T;
-#define DESCRIPTOR descriptors::D2Q9Descriptor
+#define DESCRIPTOR descriptors::D2Q9<>
const T lx = 4.0; // length of the channel
const T ly = 1.0; // height of the channel
@@ -123,7 +121,7 @@ void prepareLattice(Grid2D<T,DESCRIPTOR>& grid)
const T Lx = converter.getLatticeLength(lx);
const T Ly = converter.getLatticeLength(ly);
const T p0 = 8.*converter.getLatticeViscosity()*converter.getCharLatticeVelocity()*Lx/(Ly*Ly);
- AnalyticalLinear2D<T,T> rho(-p0/lx*DESCRIPTOR<T>::invCs2, 0, p0*DESCRIPTOR<T>::invCs2+1);
+ AnalyticalLinear2D<T,T> rho(-p0/lx*descriptors::invCs2<T,DESCRIPTOR>(), 0, p0*descriptors::invCs2<T,DESCRIPTOR>()+1);
const T maxVelocity = converter.getCharLatticeVelocity();
const T radius = ly/2;
@@ -248,8 +246,8 @@ int main(int argc, char* argv[])
prepareLattice(fineGrid);
- clout << "Total number of active cells: " << coarseGrid.getActiveVoxelN() << endl;
- clout << "starting simulation..." << endl;
+ clout << "Total number of active cells: " << coarseGrid.getActiveVoxelN() << std::endl;
+ clout << "starting simulation..." << std::endl;
Timer<T> timer(
coarseGrid.getConverter().getLatticeTime(maxPhysT),
coarseGrid.getSuperGeometry().getStatistics().getNvoxel());
@@ -260,7 +258,7 @@ int main(int argc, char* argv[])
for (int iT = 0; iT < coarseGrid.getConverter().getLatticeTime(maxPhysT); ++iT) {
if (converge.hasConverged()) {
- clout << "Simulation converged." << endl;
+ clout << "Simulation converged." << std::endl;
break;
}