int main(
int argc,
char* argv[])
{
appParamsManager.
load(argc, argv);
Param bulkModulusNum(bulkModulus.
v()/rho.
v()/
dx.v()/
dx.v());
Param shearModulusNum(shearModulus.
v()/rho.
v()/
dx.v()/
dx.v());
std::cout << "Data initialization... ";
auto displacement(asl::generateDataContainerACL_SP<FlT>(block, 3, 1u));
writer.addVector("displacement", *displacement);
writer.write();
std::cout << "Finished" << endl;
std::cout << "Numerics initialization... ";
elasticity->init();
bcRigidWall->init();
std::cout << "Finished" << endl;
std::cout << "Computing...";
bcRigidWall->execute();
for (
unsigned int i(0); i < tsim.
v(); ++i)
{
elasticity->execute();
bcRigidWall->execute();
writer.write();
}
cout << "Finished" << endl;
cout << "Computation statistic:" << endl;
cout <<
"Real Time = " << timer.
realTime() <<
"; Processor Time = " return 0;
}