52 ProcessCreateExp::create,
53 "dummy module used to create m_exp.");
70 if (
m_f->m_comm->TreatAsRankZero())
77 bool fldfilegiven = (
m_f->m_fielddef.size() != 0);
78 bool expFromFld = fldfilegiven && !vm.count(
"useSessionExpansion");
82 m_f->m_numHomogeneousDir = 0;
85 m_f->m_numHomogeneousDir =
m_f->m_fielddef[0]->m_numHomogeneousDir;
88 m_f->m_graph->SetExpansionInfo(
m_f->m_fielddef);
92 if (
m_f->m_session->DefinesSolverInfo(
"HOMOGENEOUS"))
95 m_f->m_session->GetSolverInfo(
"HOMOGENEOUS");
97 if ((HomoStr ==
"HOMOGENEOUS1D") ||
98 (HomoStr ==
"Homogeneous1D") ||
99 (HomoStr ==
"1D") || (HomoStr ==
"Homo1D"))
101 m_f->m_numHomogeneousDir = 1;
103 if ((HomoStr ==
"HOMOGENEOUS2D") ||
104 (HomoStr ==
"Homogeneous2D") ||
105 (HomoStr ==
"2D") || (HomoStr ==
"Homo2D"))
107 m_f->m_numHomogeneousDir = 2;
112 m_f->m_exp.resize(1);
115 auto domain =
m_f->m_graph->GetDomain();
116 for(
int d = 0; d < domain.size(); ++d)
118 for (
auto &compIter : domain[d])
120 for (
auto &x : compIter.second->m_geomVec)
122 IDs.push_back(x->GetGlobalID());
137 if (vm.count(
"output-points"))
139 int nPointsNew = vm[
"output-points"].as<
int>();
140 m_f->m_graph->SetExpansionInfoToPointOrder(nPointsNew);
145 if (
m_f->m_comm->TreatAsRankZero())
151 ss << cpuTime <<
"s";
152 cout <<
"\t ProcessCreateExp setexpansion CPU Time: "
159 if (
m_f->m_numHomogeneousDir == 1 && vm.count(
"output-points-hom-z"))
161 int expdim =
m_f->m_graph->GetMeshDimension();
162 m_f->m_fielddef[0]->m_numModes[expdim] =
163 vm[
"output-points-hom-z"].as<
int>();
165 m_f->m_exp[0] =
m_f->SetUpFirstExpList(
m_f->m_numHomogeneousDir,
169 if (
m_f->m_comm->TreatAsRankZero())
176 ss1 << cpuTime <<
"s";
177 cout <<
"\t ProcessCreateExp set first exp CPU Time: "
179 << ss1.str() << endl;
194 int nfields, nstrips;
195 m_f->m_session->LoadParameter(
"Strip_Z", nstrips, 1);
196 vector<string> vars =
m_f->m_session->GetVariables();
199 if (useSessionVariables)
201 m_f->m_variables = vars;
203 nfields =
m_f->m_variables.size();
205 m_f->m_exp.resize(nfields * nstrips);
207 for (
int s = 0; s < nstrips; ++s)
209 for (i = 0; i < nfields; ++i)
214 if (!
m_f->m_exp[s * nfields + i])
216 m_f->m_exp[s * nfields + i] =
m_f->AppendExpList(
217 m_f->m_numHomogeneousDir, vars[i]);
224 m_f->m_exp[s * nfields + i] =
m_f->AppendExpList(
225 m_f->m_numHomogeneousDir, vars[0]);
229 m_f->m_exp[s * nfields + i] =
m_f->AppendExpList(
230 m_f->m_numHomogeneousDir);
237 for (
int s = 0; s < nstrips; ++s)
239 for (j = 0; j < nfields; ++j)
241 for (i = 0; i <
m_f->m_data.size() / nstrips; ++i)
243 int n = i * nstrips + s;
246 auto it =
find (
m_f->m_fielddef[n]->m_fields.begin(),
247 m_f->m_fielddef[n]->m_fields.end(),
248 m_f->m_variables[j]);
249 if(it !=
m_f->m_fielddef[n]->m_fields.end())
251 m_f->m_exp[s * nfields + j]->ExtractDataToCoeffs(
255 m_f->m_exp[s * nfields + j]->UpdateCoeffs());
258 m_f->m_exp[s * nfields + j]->BwdTrans(
259 m_f->m_exp[s * nfields + j]->GetCoeffs(),
260 m_f->m_exp[s * nfields + j]->UpdatePhys());
265 m_f->m_fielddef = vector<LibUtilities::FieldDefinitionsSharedPtr>();
266 m_f->m_data = vector<std::vector<NekDouble> >();
FieldSharedPtr m_f
Field object.
virtual void Process(po::variables_map &vm)
Write mesh to output file.
virtual ~ProcessCreateExp()
void LoadFieldData(bool useSessionVariables=false)
Abstract base class for processing modules.
tKey RegisterCreatorFunction(tKey idKey, CreatorFunction classCreator, std::string pDesc="")
Register a class with the factory.
NekDouble TimePerTest(unsigned int n)
Returns amount of seconds per iteration in a test with n iterations.
static std::shared_ptr< DataType > AllocateSharedPtr(const Args &...args)
Allocate a shared pointer from the memory pool.
std::shared_ptr< Field > FieldSharedPtr
std::pair< ModuleType, std::string > ModuleKey
ModuleFactory & GetModuleFactory()
InputIterator find(InputIterator first, InputIterator last, InputIterator startingpoint, const EqualityComparable &value)
The above copyright notice and this permission notice shall be included.