28 #include <utilities/aslParametersManager.h> 29 #include <aslDataInc.h> 30 #include <math/aslTemplates.h> 31 #include <aslGeomInc.h> 32 #include <acl/aclGenerators.h> 33 #include <writers/aslVTKFormatWriters.h> 34 #include <num/aslLBGK.h> 35 #include <num/aslLBGKBC.h> 36 #include <utilities/aslTimer.h> 37 #include <acl/aclUtilities.h> 49 int main(
int argc,
char* argv[])
55 appParamsManager.
load(argc, argv);
64 auto gSize(dx.
v()*
AVec<>(size));
67 std::cout <<
"Data initialization... ";
73 .45*
makeAVec(gSize[1],gSize[1],gSize[2])));
77 auto ballMapMem(asl::generateDataContainerACL_SP<FlT>(block, 1, 1u));
81 std::cout <<
"Finished" << endl;
83 std::cout <<
"Numerics initialization... ";
95 std::vector<asl::SPNumMethod> bc;
96 std::vector<asl::SPNumMethod> bcVis;
107 std::cout <<
"Finished" << endl;
108 std::cout <<
"Computing...";
114 writer.
addVector(
"v", *lbgk->getVelocity());
123 for (
unsigned int i(1); i < 1001; ++i)
136 std::cout <<
"Finished" << endl;
138 cout <<
"time=" << timer.
getTime() <<
"; clockTime=" 142 std::cout <<
"Output...";
143 std::cout <<
"Finished" << endl;
144 std::cout <<
"Ok" << endl;
const double getTime() const
const double getProcessorLoad() const
Numerical method for fluid flow.
SPDataWrapperACL generateDataContainer_SP(const Block &b, const acl::VectorOfElements &a)
const double getClockTime() const
void addVector(std::string name, AbstractData &data)
void initAll(std::vector< T * > &v)
SPDistanceFunction generateDFCylinderInf(double r, const AVec< double > &l, const AVec< double > &c)
generates infinite cylinder
const VectorTemplate & d3q15()
Vector template.
int main(int argc, char *argv[])
asl::UValue< double > Param
acl::VectorOfElements dx(const TemplateVE &a)
differential operator
void initData(SPAbstractData d, double a)
void executeAll(std::vector< T * > &v)
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.
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 ...
SPBCond generateBCNoSlip(SPLBGK nm, const std::vector< SlicesNames > &sl)