#include <utilities/aslParametersManager.h>
#include <math/aslTemplates.h>
#include <aslGeomInc.h>
#include <aslDataInc.h>
#include <math/aslPositionFunction.h>
#include <acl/aclGenerators.h>
#include <writers/aslVTKFormatWriters.h>
#include <num/aslLBGK.h>
#include <num/aslLBGKBC.h>
#include <utilities/aslTimer.h>
#include <readers/aslVTKFormatReaders.h>
{
}
int main(
int argc,
char* argv[])
{
"1.0");
appParamsManager.
load(argc, argv);
AVec<int> size(
makeAVec(33., 33., 90.)*(1./dx.
v()));
std::cout <<
"Data initialization..." <<
flush;
asl::Block block(compressorMap->getInternalBlock());
auto comprCaseMap(asl::generateDataContainerACL_SP<FlT>(block, 1, 1u));
std::cout << "Finished" << endl;
std::cout <<
"Numerics initialization..." <<
flush;
lbgk->init();
std::vector<asl::SPNumMethod> bc;
std::vector<asl::SPNumMethod> bcVis;
.5*(block.getBPosition() + block.position)));
std::cout << "Finished" << endl;
std::cout << "Computing..." << endl;
for (unsigned int i(1); i < 10001; ++i)
{
lbgk->execute();
if (!(i%2000))
{
cout << i << endl;
}
}
cout << "Finished" << endl;
cout << "Computation statistic:" << endl;
cout <<
"Real Time = " << timer.
realTime() <<
"; Processor Time = " return 0;
}