Example: Flow around rotating cylinders
int main(
int argc,
char* argv[])
{
"1.0");
appParamsManager.
load(argc, argv);
std::cout << "Data initialization... ";
auto exCylMap(asl::generateDataContainerACL_SP<FlT>(block, 1, 1u));
.5*gSize),
auto inCylMap(asl::generateDataContainerACL_SP<FlT>(block, 1, 1u));
.5*gSize),
auto computationalDomainMap(asl::generateDataContainerACL_SP<FlT>(block, 1, 1u));
.5*gSize) |
.5*gSize) |
std::cout << "Finished" << endl;
std::cout << "Numerics initialization... ";
lbgk->init();
std::vector<asl::SPNumMethod> bc;
std::vector<asl::SPNumMethod> bcV;
std::cout << "Finished" << endl;
std::cout << "Computing...";
writer.
addScalars(
"map", *computationalDomainMap);
for (unsigned int i(0); i < 10001 ; ++i)
{
lbgk->execute();
if (!(i%1000))
{
cout << i << endl;
}
}
cout << "Finished" << endl;
cout << "Computation statistic:" << endl;
cout <<
"Real Time = " << timer.
realTime() <<
"; Processor Time = " return 0;
}