29 #include <utilities/aslParametersManager.h> 30 #include <math/aslTemplates.h> 31 #include <aslGeomInc.h> 32 #include <math/aslPositionFunction.h> 33 #include <aslDataInc.h> 34 #include <acl/aclGenerators.h> 35 #include <writers/aslVTKFormatWriters.h> 36 #include <num/aslLBGK.h> 37 #include <num/aslLBGKBC.h> 38 #include <utilities/aslTimer.h> 39 #include <readers/aslVTKFormatReaders.h> 60 center[1] = bl.
position[1] + .25 * rTunnel;
79 int main(
int argc,
char* argv[])
109 appParamsManager.
load(argc, argv);
120 FlT nuNum(nu.
v() * dt.
v() / dx.
v() / dx.
v());
122 cout <<
"Data initialization... " <<
flush;
128 asl::Block block(locomotive->getInternalBlock());
131 auto tunnelMap(asl::generateDataContainerACL_SP<FlT>(block, 1, 1u));
136 auto forceField(asl::generateDataContainerACL_SP<FlT>(block, 3, 1u));
140 cout <<
"Finished" << endl;
142 cout <<
"Numerics initialization... " <<
flush;
158 std::vector<asl::SPNumMethod> bc;
159 std::vector<asl::SPNumMethod> bcV;
184 computeForce->init();
186 cout <<
"Finished" << endl;
187 cout <<
"Computing..." << endl;
194 writer.addScalars(
"tunnel", *tunnelMap);
195 writer.addScalars(
"rho", *lbgk->getRho());
196 writer.addVector(
"v", *lbgk->getVelocity());
197 writer.addVector(
"force", *forceField);
202 computeForce->execute();
209 for (
unsigned int i(1); i < iterations.
v(); ++i)
222 computeForce->execute();
229 cout <<
"Finished" << endl;
231 cout <<
"Computation statistic:" << endl;
232 cout <<
"time = " << timer.
getTime() <<
"; clockTime = " int main(int argc, char *argv[])
const double getTime() const
const double getProcessorLoad() const
const double getClockTime() const
const AVec normalize(const AVec< T > &a)
SPNumMethod generateBCNoSlipRho(SPLBGK nmU, SPAbstractDataWithGhostNodes map)
void initAll(std::vector< T * > &v)
SPDataWithGhostNodesACLData readSurface(const string &fileName, double dx, acl::CommandQueue queue=acl::hardware.defaultQueue)
SPBCond generateBCConstantPressureVelocity(SPLBGK nm, double p, AVec<> v, const std::vector< SlicesNames > &sl)
const V getBPosition() const
SPDistanceFunction generateDFPlane(const AVec< double > &n, const AVec< double > &p0)
const VectorTemplate & d3q15()
Vector template.
acl::VectorOfElements dx(const TemplateVE &a)
differential operator
asl::SPDistanceFunction generateTunnel(asl::Block &bl)
SPNumMethod generateBCVelocityVel(SPLBGK nm, SPPositionFunction v, SPAbstractDataWithGhostNodes map)
void initData(SPAbstractData d, double a)
void executeAll(std::vector< T * > &v)
void addScalars(std::string name, AbstractData &data)
VectorOfElements generateVEConstant(T a)
Generates VectorOfElements with 1 Element acl::Constant with value a.
SPNumMethod generateComputeSurfaceForce(SPLBGK nm, SPDataWithGhostNodesACLData fF, SPAbstractDataWithGhostNodes map)
SPDistanceFunction generateDFCylinder(double r, const AVec< double > &l, const AVec< double > &c)
generates cylinder
void load(int argc, char *argv[])
SPNumMethod generateBCNoSlipVel(SPLBGK nmU, SPAbstractDataWithGhostNodes map)
for velocity field
SPPositionFunction generatePFConstant(const AVec< double > &a)
contains different kernels for preprocessing and posprocessing of data used by LBGK ...
SPNumMethod generateBCVelocity(SPLBGK nm, SPPositionFunction v, SPAbstractDataWithGhostNodes map)
SPBCond generateBCNoSlip(SPLBGK nm, const std::vector< SlicesNames > &sl)