Commit 9775a4f5 authored by Julian Marcon's avatar Julian Marcon

Transfer BL end info to CurveMesh

parent 11b1632a
......@@ -102,6 +102,34 @@ void Generator2D::Process()
{
m_curvemeshes[i] =
MemoryManager<CurveMesh>::AllocateSharedPtr(i, m_mesh);
if (m_blends.count(i))
{
vector<CADVertSharedPtr> vertices =
m_mesh->m_cad->GetCurve(i)->GetVertex();
Array<OneD, NekDouble> loc;
NekDouble t;
if (m_blends[i] == 2)
{
loc = vertices[0]->GetLoc();
t = m_thickness.Evaluate(m_thickness_ID, loc[0], loc[1],
loc[2], 0.0);
m_curvemeshes[i]->SetOffset(0, t);
loc = vertices[1]->GetLoc();
t = m_thickness.Evaluate(m_thickness_ID, loc[0], loc[1],
loc[2], 0.0);
m_curvemeshes[i]->SetOffset(1, t);
}
else
{
loc = vertices[m_blends[i]]->GetLoc();
t = m_thickness.Evaluate(m_thickness_ID, loc[0], loc[1],
loc[2], 0.0);
m_curvemeshes[i]->SetOffset(m_blends[i], t);
}
}
}
else
{
......
......@@ -136,6 +136,11 @@ public:
return m_id;
}
void SetOffset(unsigned i, NekDouble offset)
{
m_endoffset[i] = offset;
}
private:
/**
* @brief get node spacing sampling function
......@@ -188,6 +193,8 @@ private:
std::vector<NodeSharedPtr> m_meshpoints;
LibUtilities::AnalyticExpressionEvaluator m_bl;
int m_blID;
/// offset of second point at each end
std::map<unsigned, NekDouble> m_endoffset;
};
}
......
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