Definition at line 47 of file DiffusionSolverTimeInt.cpp.
Diffusion::Diffusion |
( |
int |
argc, |
|
|
char * |
argv[] |
|
) |
| |
Definition at line 84 of file DiffusionSolverTimeInt.cpp.
87 session = LibUtilities::SessionReader::CreateInstance(argc, argv);
109 fields[0] =
field->UpdatePhys();
112 unsigned int nq =
field->GetNpoints();
114 field->GetCoords(x0,x1,x2);
118 =
session->GetFunction(
"InitialConditions",
"u");
119 icond->Evaluate(x0,x1,x2,0.0,
field->UpdatePhys());
MultiRegions::ContField2DSharedPtr field
General purpose memory allocation routines with the ability to allocate from thread specific memory p...
LibUtilities::FieldIOSharedPtr fld
Array< OneD, Array< OneD, NekDouble > > fields
boost::shared_ptr< Equation > EquationSharedPtr
LibUtilities::SessionReaderSharedPtr session
SpatialDomains::MeshGraphSharedPtr graph
Diffusion::~Diffusion |
( |
| ) |
|
Definition at line 152 of file DiffusionSolverTimeInt.cpp.
References Nektar::StdRegions::eFactorLambda, Nektar::NullFlagList, and Vmath::Smul().
Referenced by TimeIntegrate().
161 for (
int i = 0; i < inarray.num_elements(); ++i)
169 field->HelmSolve(outarray[i],
170 field->UpdateCoeffs(),
174 field->BwdTrans (
field->GetCoeffs(), outarray[i]);
MultiRegions::ContField2DSharedPtr field
std::map< ConstFactorType, NekDouble > ConstFactorMap
void Smul(int n, const T alpha, const T *x, const int incx, T *y, const int incy)
Scalar multiply y = alpha*y.
static FlagList NullFlagList
An empty flag list.
void Diffusion::ExactSolution |
( |
| ) |
|
|
private |
Definition at line 194 of file DiffusionSolverTimeInt.cpp.
196 unsigned int nq =
field->GetNpoints();
198 field->GetCoords(x0,x1,x2);
201 session->GetFunction(
"ExactSolution",0);
210 cout <<
"L inf error: "
211 <<
field->Linf(
field->GetPhys(), exact) << endl;
212 cout <<
"L 2 error: "
213 <<
field->L2(
field->GetPhys(), exact) << endl;
214 cout <<
"H 1 error: "
215 <<
field->H1(
field->GetPhys(), exact) << endl;
MultiRegions::ContField2DSharedPtr field
boost::shared_ptr< Equation > EquationSharedPtr
LibUtilities::SessionReaderSharedPtr session
void Diffusion::TimeIntegrate |
( |
| ) |
|
Definition at line 127 of file DiffusionSolverTimeInt.cpp.
References DoImplicitSolve(), Nektar::LibUtilities::GetTimeIntegrationWrapperFactory(), Vmath::Vcopy(), and Vmath::Zero().
140 for (
int n = 0; n <
nSteps; ++n)
LibUtilities::TimeIntegrationWrapperSharedPtr IntScheme
MultiRegions::ContField2DSharedPtr field
void DefineImplicitSolve(FuncPointerT func, ObjectPointerT obj)
LibUtilities::TimeIntegrationSolutionSharedPtr u
Array< OneD, Array< OneD, NekDouble > > fields
TimeIntegrationWrapperFactory & GetTimeIntegrationWrapperFactory()
void DoImplicitSolve(const Array< OneD, const Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray, const NekDouble time, const NekDouble lambda)
LibUtilities::TimeIntegrationSchemeOperators ode
void Zero(int n, T *x, const int incx)
Zero vector.
void Vcopy(int n, const T *x, const int incx, T *y, const int incy)
void Diffusion::WriteSolution |
( |
| ) |
|
|
private |
Definition at line 178 of file DiffusionSolverTimeInt.cpp.
181 std::vector<LibUtilities::FieldDefinitionsSharedPtr> FieldDef
182 =
field->GetFieldDefinitions();
183 std::vector<std::vector<NekDouble> > FieldData(FieldDef.size());
184 for(
int i = 0; i < FieldDef.size(); ++i)
186 FieldDef[i]->m_fields.push_back(
"u");
187 field->AppendFieldData(FieldDef[i], FieldData[i]);
189 fld->Write(
session->GetSessionName() +
".fld", FieldDef, FieldData);
MultiRegions::ContField2DSharedPtr field
LibUtilities::FieldIOSharedPtr fld
LibUtilities::SessionReaderSharedPtr session
unsigned int Diffusion::nSteps |
|
private |
string Diffusion::sessionName |
|
private |