ForkMaster/netgen
1
0
Fork 0
mirror of https://github.com/NGSolve/netgen.git synced 2024-12-24 21:10:33 +05:00
Code Issues Packages Projects Releases Wiki Activity

boundarylayer - calculate edge tangents on the fly

Browse Source
This commit is contained in:
mhochsteger@cerbsim.com 2022-03-07 17:24:36 +01:00
parent 97b9dae812
commit 8dd4c0c9e3
1 changed files with 12 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

15
libsrc/meshing/boundarylayer.cpp
View File

@ -16,6 +16,17 @@ namespace netgen
return Cross(mesh[el[1]]-v0, mesh[el[2]]-v0).Normalize(); return Cross(mesh[el[1]]-v0, mesh[el[2]]-v0).Normalize();
}; };
Vec<3> getEdgeTangent(const Mesh & mesh, PointIndex pi)
{
Vec<3> tangent = 0.0;
for(auto segi : mesh.GetTopology().GetVertexSegments(pi))
{
auto & seg = mesh[segi];
tangent += (mesh[seg[1]] - mesh[seg[0]]);
}
return tangent.Normalize();
}
void InsertVirtualBoundaryLayer (Mesh & mesh) void InsertVirtualBoundaryLayer (Mesh & mesh)
{ {
@ -324,7 +335,6 @@ namespace netgen
} }
const auto & p2sel = mesh.CreatePoint2SurfaceElementTable(); const auto & p2sel = mesh.CreatePoint2SurfaceElementTable();
Array<Vec<3>, PointIndex> edge_tangents(np);
for(auto pi : mesh.Points().Range()) for(auto pi : mesh.Points().Range())
{ {
@ -365,7 +375,6 @@ namespace netgen
auto npnp = np * np; auto npnp = np * np;
auto nn = n * n; auto nn = n * n;
if(nn-npn*npn/npnp == 0) { np = n; continue; } if(nn-npn*npn/npnp == 0) { np = n; continue; }
edge_tangents[pi] = Cross(np, n).Normalize();
np += (nn - npn)/(nn - npn*npn/npnp) * (n - npn/npnp * np); np += (nn - npn)/(nn - npn*npn/npnp) * (n - npn/npnp * np);
} }
} }
@ -546,7 +555,7 @@ namespace netgen
{ {
auto pi = points[i]; auto pi = points[i];
len += (mesh[pi] - mesh[points[i-1]]).Length(); len += (mesh[pi] - mesh[points[i-1]]).Length();
auto t = edge_tangents[pi]; auto t = getEdgeTangent(mesh, pi);
auto lam = len/edge_len; auto lam = len/edge_len;
auto interpol = (1-lam) * (gt1 * t) * t + lam * (gt2 * t) * t; auto interpol = (1-lam) * (gt1 * t) * t + lam * (gt2 * t) * t;
growthvectors[pi] += interpol; growthvectors[pi] += interpol;