#include <aslDataInc.h>
#include <acl/aclGenerators.h>
#include <writers/aslVTKFormatWriters.h>
#include <num/aslFDElasticity.h>
#include <num/aslFDElasticityBC.h>
#include <utilities/aslTimer.h>
#include <utilities/aslParametersManager.h>
#include <math/aslTemplates.h>
#include <acl/aclMath/aclVectorOfElements.h>
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();
}
std::cout << "Finished" << endl;
cout <<
"time=" << timer.
getTime() <<
"; clockTime=" std::cout << "Output...";
std::cout << "Finished" << endl;
std::cout << "Ok" << endl;
return 0;
}