29#include <utilities/aslParametersManager.h>
30#include <math/aslTemplates.h>
31#include <aslGeomInc.h>
32#include <aslDataInc.h>
33#include <acl/aclGenerators.h>
34#include <writers/aslVTKFormatWriters.h>
35#include <num/aslFDStefanMaxwell.h>
36#include <num/aslBasicBC.h>
37#include <utilities/aslTimer.h>
38#include <acl/aclMath/aclVectorOfElements.h>
48int main(
int argc,
char* argv[])
54 appParamsManager.
load(argc, argv);
60 Param diffCoefNum(diffCoef.
v()*dt.
v()/dx.v()/dx.v());
64 auto gSize(dx.v()*
AVec<>(size));
67 std::cout <<
"Flow: Data initialization... ";
79 std::cout <<
"Finished" << endl;
81 std::cout <<
"Flow: Numerics initialization... ";
84 auto nm(generateFDStefanMaxwell(c1Field, c2Field, diffCoefNum.
v(), templ));
92 std::vector<asl::SPNumMethod> bc;
105 std::cout <<
"Finished" << endl;
106 std::cout <<
"Computing..." << flush;
118 for(
unsigned int i(1); i < 401; ++i)
130 cout <<
"Finished" << endl;
132 cout <<
"Computation statistic:" << endl;
133 cout <<
"Real Time = " << timer.
realTime() <<
"; Processor Time = "
void load(int argc, char *argv[])
const double realTime() const
const double processorTime() const
const double processorLoad() const
Updatable value. This class stores value and its TimeStamp.
void addScalars(std::string name, AbstractData &data)
SPDataWrapperACLData generateDataContainerACL_SP(const Block &b, unsigned int n=1)
generates pointer to ACL Data field with n components
SPBCond generateBCConstantGradient(SPAbstractDataWithGhostNodes d, double v, const VectorTemplate *const t, const std::vector< SlicesNames > &sl)
Bondary condition that makes fixed gradient <>
SPBCond generateBCConstantValue(SPAbstractDataWithGhostNodes d, double v, const std::vector< SlicesNames > &sl)
Bondary condition that puts fixed value in each point.
const VectorTemplate & d3q7()
Vector template.
VectorOfElements generateVEConstant(T a)
Generates VectorOfElements with 1 Element acl::Constant with value a.
void initData(SPAbstractData d, double a)
asl::UValue< double > Param