Commit 6de157b0 authored by Michael Turner's avatar Michael Turner

change octree query

parent 1b028cdd
...@@ -106,11 +106,13 @@ NekDouble Octree::Query(Array<OneD, NekDouble> loc) ...@@ -106,11 +106,13 @@ NekDouble Octree::Query(Array<OneD, NekDouble> loc)
// contain the point loc until a leaf is found // contain the point loc until a leaf is found
//first search through sourcepoints //first search through sourcepoints
NekDouble tmp = numeric_limits<double>::max();
for(int i = 0; i < m_lsources.size(); i++) for(int i = 0; i < m_lsources.size(); i++)
{ {
if(m_lsources[i].withinRange(loc)) if(m_lsources[i].withinRange(loc))
{ {
return m_lsources[i].delta; tmp = min(m_lsources[i].delta,tmp);
} }
} }
...@@ -183,7 +185,7 @@ NekDouble Octree::Query(Array<OneD, NekDouble> loc) ...@@ -183,7 +185,7 @@ NekDouble Octree::Query(Array<OneD, NekDouble> loc)
found = true; found = true;
} }
} }
return n->GetDelta(); return min(n->GetDelta(),tmp);
} }
void Octree::WriteOctree(string nm) void Octree::WriteOctree(string nm)
......
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