#include <math/aslTemplates.h>
#include <aslGeomInc.h>
#include <math/aslPositionFunction.h>
#include <aslDataInc.h>
#include <acl/aclGenerators.h>
#include <writers/aslVTKFormatWriters.h>
#include <num/aslLBGK.h>
#include <num/aslLBGKBC.h>
#include <utilities/aslTimer.h>
#include <acl/aclUtilities.h>
#include <readers/aslVTKFormatReaders.h>
#include <utilities/aslParametersManager.h>
{
}
int main(
int argc,
char* argv[])
{
"1.0");
appParamsManager.
load(argc, argv);
AVec<int> size(
makeAVec(150., 37., 53.2)*(1./dx.
v()));
std::cout << "Data initialization..." << endl;
auto tunnelMap(asl::generateDataContainerACL_SP<FlT>(block, 1, 1u));
auto forceField(asl::generateDataContainerACL_SP<FlT>(block, 3, 1u));
std::cout << "Finished"<<endl;
std::cout << "Numerics initialization..." << endl;
lbgk->init();
std::vector<asl::SPNumMethod> bc;
std::vector<asl::SPNumMethod> bcV;
computeForce->init();
cout << "Finished" << endl;
cout << "Computing..." << endl;
writer.addScalars("tunnel", *tunnelMap);
writer.addScalars("rho", *lbgk->getRho());
writer.addVector("v", *lbgk->getVelocity());
writer.addVector("force", *forceField);
computeForce->execute();
writer.write();
for(unsigned int i(1); i < 2001; ++i)
{
lbgk->execute();
if(!(i%200))
{
cout<<i<<endl;
computeForce->execute();
writer.write();
}
}
cout << "Finished" << endl;
cout << "Computation statistic:" << endl;
cout <<
"Real Time = " << timer.
realTime() <<
"; Processor Time = " return 0;
}