diff options
Adapt refinement apps to meta-descriptor
Diffstat (limited to 'apps/adrian/poiseuille2d')
-rw-r--r-- | apps/adrian/poiseuille2d/poiseuille2d.cpp | 12 |
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; } |