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

drawing solution at the boundary

Browse Source
This commit is contained in:
Joachim Schöberl 2019-03-16 06:30:10 +01:00
parent f57ad036b5
commit bf59c1d490
1 changed files with 7 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

11
libsrc/visualization/vssolution.cpp
View File

@ -1175,7 +1175,7 @@ namespace netgen
for (int i = 0; i < npt; i++) for (int i = 0; i < npt; i++)
pref[i] = double(i) / (npt-1); pref[i] = double(i) / (npt-1);
int meshdim = mesh->GetDimension();
for (SegmentIndex i = 0; i < mesh -> GetNSeg(); i++) for (SegmentIndex i = 0; i < mesh -> GetNSeg(); i++)
{ {
// mesh->GetCurvedElements(). // mesh->GetCurvedElements().
@ -1190,8 +1190,10 @@ namespace netgen
{ {
vsol->solclass->GetSegmentValue (i, pref[j], &mvalues[0]); vsol->solclass->GetSegmentValue (i, pref[j], &mvalues[0]);
// values[j] = ExtractValue (sol, scalcomp, &mvalues[0]); // values[j] = ExtractValue (sol, scalcomp, &mvalues[0]);
points[j](0) += scaledeform * mvalues[0]; for (int k = 0; k < min(ncomp, 3); k++)
points[j](1) += scaledeform * mvalues[1]; points[j](k) += scaledeform * mvalues[0];
// points[j](0) += scaledeform * mvalues[0];
// points[j](1) += scaledeform * mvalues[1];
} }
} }
else if (sol) else if (sol)
@ -1200,7 +1202,8 @@ namespace netgen
{ {
sol->solclass->GetSegmentValue (i, pref[j], &mvalues[0]); sol->solclass->GetSegmentValue (i, pref[j], &mvalues[0]);
values[j] = ExtractValue (sol, scalcomp, &mvalues[0]); values[j] = ExtractValue (sol, scalcomp, &mvalues[0]);
points[j](1) += scaledeform * values[j]; if (meshdim <= 2)
points[j](meshdim) += scaledeform * values[j];
} }
} }