Commit 2b85bb06 authored by Michael Turner's avatar Michael Turner
Browse files

douglas comments

parent 5b755c57
......@@ -538,7 +538,7 @@ To extract a surface use the command:
where the integers are surface IDs to be extracted.
An optional arguemnt of \inltt{detectbnd} can be added to identify the boundary composites as part of the surface extraction.
An optional arguemnt of \inltt{detectbnd} can be added to identify the boundary composites as part of the surface extraction.
\subsection{Linearisation}
......@@ -703,7 +703,7 @@ vertex nodes as boundaries, then the surfaces are meshed in their 2D parameter
plane using the curve meshes as boundaries and finally the 3D volume is meshed
using the surface mesh as the boundary to complete the linear mesh. In \mc, to
achieve greater robustness, the 2D mesh generation library Triangle is used and
the TetGen library for the 3D. Both of which are highly developed Deluanay based
the TetGen library for the 3D. Both of which are highly developed Delaunay based
mesh generators. As with all additional libraies in \nekpp these are
automatically downloaded and installed if needed.
......@@ -726,7 +726,7 @@ the parameter plane of the CAD surface to ensure the high-order triangles have
as little distortion as possible while remaining exactly on the geometric
surface. To do this the system models the high-order edges and triangles as a
network of springs with an associated spring energy which is minimised using a
multidimensional Newton type optimisation procedure with a Guass-Sidel matrix
multidimensional Newton type optimisation procedure with a Gauss-Seidel matrix
solver.
\subsubsection{Mesh Correction}
......
......@@ -171,9 +171,9 @@ Array<OneD, NekDouble> BLMesh::GetNormal(vector<ElementSharedPtr> tris)
Array<OneD, NekDouble> bestN = Np;
NekDouble val = -1.0*numeric_limits<double>::max();
NekDouble dtheta = 3.142/5.0;
NekDouble dphi = 3.142/5.0;
while(dtheta > 3.142/300.0)
NekDouble dtheta = M_PI/5.0;
NekDouble dphi = M_PI/5.0;
while(dtheta > M_PI/300.0)
{
NekDouble theta0 = acos(bestN[2]);
NekDouble phi0 = atan2(bestN[1],bestN[0]);
......@@ -339,16 +339,12 @@ void BLMesh::Mesh()
if(Visability(g->second,bln.N) < 0.0)
{
cout << "failed " << (*it)->m_x << " " << (*it)->m_y << " "
<< (*it)->m_z << " " << Visability(g->second,bln.N) << endl;
cerr << "failed " << (*it)->m_x << " " << (*it)->m_y << " "
<< (*it)->m_z << " "
<< Visability(g->second,bln.N) << endl;
failed++;
}
if(Visability(g->second,bln.N) < 0.1)
{
cout << "Visability low" << endl;
}
Array<OneD, NekDouble> loc = (*it)->GetLoc();
for(int k = 0; k < 3; k++)
{
......
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