Commit cbfc4591 authored by Yumnah Mohamied's avatar Yumnah Mohamied

quick (temporary) workaround for boundary conditions from file for 3D

parent a6ad062c
......@@ -2424,22 +2424,23 @@
if(filebcs != "")
{
string var = filebcs.substr(
string varString = filebcs.substr(
0, filebcs.find_last_of("."));
int len = var.length();
var = var.substr(len-1,len);
int len = varString.length();
varString = varString.substr(len-1, len);
int varInt = atoi(varString.c_str());
cout << "Boundary condition from file:"
<< filebcs << endl;
std::vector<LibUtilities::FieldDefinitionsSharedPtr> FieldDef;
std::vector<std::vector<NekDouble> > FieldData;
Import(filebcs,FieldDef, FieldData);
Import(filebcs, FieldDef, FieldData);
// copy FieldData into locExpList
locExpList->ExtractDataToCoeffs(
FieldDef[0], FieldData[0],
FieldDef[0]->m_fields[0], locExpList->UpdateCoeffs());
FieldDef[0]->m_fields[varInt],
locExpList->UpdateCoeffs());
locExpList->BwdTrans_IterPerExp(
locExpList->GetCoeffs(),
locExpList->UpdatePhys());
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment