ForkMaster/netgen
1
0
Fork 0
mirror of https://github.com/NGSolve/netgen.git synced 2024-11-11 16:49:16 +05:00
Code Issues Packages Projects Releases Wiki Activity

Merge branch 'surface_vectors_on_quad' into 'master'

Browse Source 
Draw surface vectors on quads

See merge request ngsolve/netgen!558
This commit is contained in:
Hochsteger, Matthias 2023-03-08 10:04:27 +01:00
commit 3e52c44aeb
2 changed files with 16 additions and 11 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

25
libsrc/visualization/vssolution.cpp
View File

@ -2202,7 +2202,7 @@ namespace netgen
void VisualSceneSolution :: DrawTrigSurfaceVectors(const NgArray< Point<3> > & lp,
const Point<3> & pmin, const Point<3> & pmax,
const int sei, const SolData * vsol)
const int sei, const SolData * vsol, bool swap_lam)
{
shared_ptr<Mesh> mesh = GetMesh();
@ -2267,12 +2267,17 @@ namespace netgen
if (lam1 >= 0 && lam2 >= 0 && lam1+lam2 <= 1)
{
if(swap_lam)
{
Swap(lam1, lam2);
lam1 = 1.0-lam1;
lam2 = 1.0-lam2;
}
Point<3> cp;
for (k = 0; k < 3; k++)
cp(k) = lp[0](k) +
lam1 * (lp[1](k)-lp[0](k)) +
lam2 * (lp[2](k)-lp[0](k));
Point<2> xref(lam1, lam2);
if (mesh->GetCurvedElements().IsHighOrder())
mesh->GetCurvedElements().
@ -2460,22 +2465,21 @@ namespace netgen
}
else if (el.GetType() == QUAD)
{
/*
NgArray < Point<3> > lp(3);
lp[0] = mesh->Point(el[0]);
lp[1] = mesh->Point(el[1]);
lp[2] = mesh->Point(el[2]);
DrawTrigSurfaceVectors(lp,pmin,pmax,sei,vsol);
lp[0] = mesh->Point(el[0]);
lp[1] = mesh->Point(el[2]);
lp[2] = mesh->Point(el[3]);
DrawTrigSurfaceVectors(lp,pmin,pmax,sei,vsol);
*/
lp[0] = mesh->Point(el[2]);
lp[1] = mesh->Point(el[1]);
lp[2] = mesh->Point(el[3]);
DrawTrigSurfaceVectors(lp,pmin,pmax,sei,vsol, true);
/*
Point<3> lp[4];
Point<2> p2d[4];
@ -2574,6 +2578,7 @@ namespace netgen
}
}
*/
}
}
}

2
libsrc/visualization/vssolution.hpp
View File

@ -328,7 +328,7 @@ public:
void DrawSurfaceVectors ();
void DrawTrigSurfaceVectors(const NgArray< Point<3> > & lp, const Point<3> & pmin, const Point<3> & pmax,
const int sei, const SolData * vsol);
const int sei, const SolData * vsol, bool swap_lam=false);
void DrawIsoSurface(const SolData * sol, const SolData * grad, int comp);
void DrawIsoLines (const Point<3> & p1,