Commit c2c842d2 authored by Michael Turner's avatar Michael Turner

working version, get ho working for quads

parent 486cdf2d
......@@ -286,6 +286,19 @@ void Generator2D::MakeBL(int faceid, vector<EdgeLoop> e)
LibUtilities::eQuadrilateral, conf, qns, tags);
E->m_parentCAD = m_mesh->m_cad->GetSurf(faceid);
for (int j = 0; j < E->GetEdgeCount(); ++j)
{
pair<EdgeSet::iterator,bool> testIns;
EdgeSharedPtr ed = E->GetEdge(j);
// look for edge in m_mesh edgeset from curves
EdgeSet::iterator s = m_mesh->m_edgeSet.find(ed);
if (!(s == m_mesh->m_edgeSet.end()))
{
ed = *s;
E->SetEdge(j, *s);
}
}
m_mesh->m_element[2].push_back(E);
}
}
......
......@@ -886,8 +886,6 @@ void Octree::CompileSourcePointList()
if (C != 0.0)
{
cout << 1.0 / C << endl;
NekDouble del = 2.0 * (1.0 / C) * sqrt(m_eps * (2.0 - m_eps));
if (del > m_maxDelta)
......
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