48 int main(
int argc,
char* argv[])
54 appParamsManager.
load(argc, argv);
60 Param diffCoefNum(diffCoef.
v()*dt.
v()/dx.
v()/dx.
v());
63 auto gSize(dx.
v()*
AVec<>(size));
66 std::cout <<
"Data initialization... ";
76 auto ballMapMem(asl::generateDataContainerACL_SP<FlT>(block, 1, 1u));
80 auto ballBMapMem(asl::generateDataContainerACL_SP<FlT>(block, 1, 1u));
83 auto cField(asl::generateDataContainerACL_SP<FlT>(block, 1, 1u));
87 std::cout <<
"Finished" << endl;
89 std::cout <<
"Numerics initialization... ";
95 std::vector<asl::SPNumMethod> bc;
102 std::cout <<
"Finished" << endl;
103 std::cout <<
"Computing..." <<
flush;
108 writer.addScalars(
"mapE", *ballBMapMem);
109 writer.addScalars(
"c", *cField);
115 for (
unsigned int i(1); i < 201; ++i)
127 std::cout <<
"Finished" << endl;
129 cout <<
"time=" << timer.
getTime() <<
"; clockTime=" 133 std::cout <<
"Output...";
134 std::cout <<
"Finished" << endl;
135 std::cout <<
"Ok" << endl;
SPDistanceFunction generateDFSphere(double r, const AVec< double > &c)
generates sphere
SPNumMethod generateBCConstantGradient2(SPAbstractDataWithGhostNodes d, double v, SPAbstractDataWithGhostNodes map, const VectorTemplate *const t)
Bondary condition that makes fixed gradient, second order accuracy.
const AVec normalize(const AVec< T > &a)
void initAll(std::vector< T *> &v)
asl::UValue< double > Param
const double getProcessorLoad() const
const VectorTemplate & d3q15()
Vector template.
SPFDAdvectionDiffusion generateFDAdvectionDiffusion(SPDataWithGhostNodesACLData c, double diffustionCoeff, SPAbstractDataWithGhostNodes v, const VectorTemplate *vt, bool compressibilityCorrection=false)
acl::VectorOfElements dx(const TemplateVE &a)
differential operator
int main(int argc, char *argv[])
const double getClockTime() const
void executeAll(std::vector< T *> &v)
void initData(SPAbstractData d, double a)
const double getTime() const
void addScalars(std::string name, AbstractData &data)
SPBCond generateBCConstantValue(SPAbstractDataWithGhostNodes d, double v, const std::vector< SlicesNames > &sl)
Bondary condition that puts fixed value in each point.
void load(int argc, char *argv[])