61 int main(
int argc,
char* argv[])
68 appParamsManager.
load(argc, argv);
74 Param w(6.*3.14*2./60.);
84 std::cout <<
"Data initialization..." <<
flush;
89 asl::Block block(compressorMap->getInternalBlock());
91 auto comprCaseMap(asl::generateDataContainerACL_SP<FlT>(block, 1, 1u));
95 std::cout <<
"Finished" << endl;
97 std::cout <<
"Numerics initialization..." <<
flush;
110 std::vector<asl::SPNumMethod> bc;
111 std::vector<asl::SPNumMethod> bcVis;
115 .5*(block.getBPosition() + block.position)));
128 std::cout <<
"Finished" << endl;
129 std::cout <<
"Computing..." << endl;
133 writer.
addScalars(
"compressor", *compressorMap);
136 writer.
addVector(
"v", *lbgk->getVelocity());
144 for (
unsigned int i(1); i < 10001; ++i)
157 std::cout <<
"Finished" << endl;
159 cout <<
"time=" << timer.
getTime() <<
"; clockTime=" 163 std::cout <<
"Output...";
164 std::cout <<
"Finished" << endl;
165 std::cout <<
"Ok" << endl;
Numerical method for fluid flow.
const AVec normalize(const AVec< T > &a)
int main(int argc, char *argv[])
void addVector(std::string name, AbstractData &data)
asl::SPDistanceFunction generateCase(asl::Block &bl)
void initAll(std::vector< T *> &v)
SPDataWithGhostNodesACLData readSurface(const string &fileName, double dx, acl::CommandQueue queue=acl::hardware.defaultQueue)
asl::UValue< double > Param
SPDistanceFunction generateDFCylinderInf(double r, const AVec< double > &l, const AVec< double > &c)
generates infinite cylinder
const double getProcessorLoad() const
const VectorTemplate & d3q15()
Vector template.
acl::VectorOfElements dx(const TemplateVE &a)
differential operator
const double getClockTime() const
void executeAll(std::vector< T *> &v)
SPNumMethod generateBCVelocityVel(SPLBGK nm, SPPositionFunction v, SPAbstractDataWithGhostNodes map)
void initData(SPAbstractData d, double a)
const double getTime() const
void addScalars(std::string name, AbstractData &data)
SPBCond generateBCConstantPressure(SPLBGK nm, double p, const std::vector< SlicesNames > &sl)
VectorOfElements generateVEConstant(T a)
Generates VectorOfElements with 1 Element acl::Constant with value a.
SPPositionFunction generatePFRotationField(const AVec< double > &axis, const AVec< double > &c)
const V getBPosition() const
void load(int argc, char *argv[])
SPNumMethod generateBCNoSlipVel(SPLBGK nmU, SPAbstractDataWithGhostNodes map)
for velocity field
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)