28 #include <aslDataInc.h>
29 #include <acl/aclGenerators.h>
30 #include <writers/aslVTKFormatWriters.h>
31 #include <num/aslFDPoroElasticity.h>
32 #include <num/aslFDElasticityBC.h>
33 #include <num/aslFDPoroElasticityBC.h>
34 #include <utilities/aslTimer.h>
35 #include <utilities/aslParametersManager.h>
36 #include <math/aslTemplates.h>
37 #include <acl/aclMath/aclVectorOfElements.h>
38 #include <aslGeomInc.h>
39 #include <acl/aclUtilities.h>
47 int main(
int argc,
char* argv[])
63 appParamsManager.
load(argc, argv);
65 Param bulkModulusNum(bulkModulus.
v()/rho.
v()/
dx.v()/
dx.v()*dt.
v()*dt.
v());
66 Param shearModulusNum(shearModulus.
v()/rho.
v()/
dx.v()/
dx.v()*dt.
v()*dt.
v());
70 std::cout <<
"Data initialization... " <<
flush;
73 auto displacement(asl::generateDataContainerACL_SP<FlT>(block, 3, 1u));
75 auto pressureL(asl::generateDataContainerACL_SP<FlT>(block, 1, 1u));
80 auto mapX(asl::generateDataContainerACL_SP<FlT>(block, 1, 1u));
83 std::cout <<
"Finished" << endl;
85 std::cout <<
"Numerics initialization... " <<
flush;
91 hydraulicConductivity.
v(),
96 std::vector<asl::SPNumMethod> bc;
103 writer.
addVector(
"displacement", *displacement);
104 writer.
addScalars(
"pressure", *elasticity->getPressureData());
105 writer.
addScalars(
"pressureL", *elasticity->getLiquidPressureData());
107 std::cout <<
"Finished" << endl;
108 std::cout <<
"Computing..." << endl;
117 for (
unsigned int i(0); i < tsim.
v(); ++i)
120 elasticity->execute();
133 cout <<
"Finished" << endl;
135 cout <<
"Computation statistic:" << endl;
136 cout <<
"Real Time = " << timer.
realTime() <<
"; Processor Time = "
140 cout <<
"timeBulk=" << timerBulk.
realTime() <<
141 "; timeBC=" << timerBC.
realTime() << endl;
void initData(VectorOfElements a, VectorOfElements initializationValue, const KernelConfiguration &kernelConfig)
void load(int argc, char *argv[])
const double realTime() const
const double processorTime() const
const double processorLoad() const
void addVector(std::string name, AbstractData &data)
void addScalars(std::string name, AbstractData &data)
int main(int argc, char *argv[])
SPDataWrapperACL generateDataContainer_SP(const Block &b, const acl::VectorOfElements &a)
acl::VectorOfElements dx(const TemplateVE &a)
differential operator
void addBCZeroStress(std::vector< SPNumMethod > &bcList, SPFDPoroElasticity nm, SPAbstractDataWithGhostNodes map)
Bondary condition set given values to pressure.
SPDistanceFunction generateDFInBlock(const Block &b, unsigned int nG)
generates map corresponding to external (ghost) part of the block
SPDistanceFunction normalize(SPDistanceFunction a, double dx)
const VectorTemplate & d3q15()
Vector template.
VectorOfElements generateVEConstant(T a)
Generates VectorOfElements with 1 Element acl::Constant with value a.
std::shared_ptr< T > map(ElementData m)
void initAll(std::vector< T * > &v)
SPFDPoroElasticity generateFDPoroElasticity(SPDataWithGhostNodesACLData d, SPDataWithGhostNodesACLData pl, double bM, double sM, double k, const VectorTemplate *vT)
void addBCRigidWall(std::vector< SPNumMethod > &bcList, SPFDPoroElasticity nm, const std::vector< SlicesNames > &sl)
void initData(SPAbstractData d, double a)
void executeAll(std::vector< T * > &v)