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

visualize deformation of 1d elements

Browse Source
This commit is contained in:
Joachim Schöberl 2018-04-30 18:41:19 +02:00
parent 819fca6940
commit 1629d6bf69
1 changed files with 15 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

16
libsrc/visualization/vssolution.cpp
View File

@ -1164,8 +1164,12 @@ namespace netgen
const SolData * sol = NULL; const SolData * sol = NULL;
if (scalfunction != -1) sol = soldata[scalfunction]; if (scalfunction != -1) sol = soldata[scalfunction];
const SolData * vsol = NULL;
if (deform && vecfunction != -1) vsol = soldata[vecfunction];
int ncomp = 0; int ncomp = 0;
if (sol) ncomp = sol->components; if (sol) ncomp = sol->components;
if (vsol) ncomp = vsol->components;
Array<double> mvalues(ncomp); Array<double> mvalues(ncomp);
@ -1180,7 +1184,17 @@ namespace netgen
for (int j = 0; j < npt; j++) for (int j = 0; j < npt; j++)
mesh->GetCurvedElements(). mesh->GetCurvedElements().
CalcSegmentTransformation (pref[j], i, points[j]); CalcSegmentTransformation (pref[j], i, points[j]);
if (sol) if (vsol)
{
for (int j = 0; j < npt; j++)
{
vsol->solclass->GetSegmentValue (i, pref[j], &mvalues[0]);
// values[j] = ExtractValue (sol, scalcomp, &mvalues[0]);
points[j](0) += scaledeform * mvalues[0];
points[j](1) += scaledeform * mvalues[1];
}
}
else if (sol)
{ {
for (int j = 0; j < npt; j++) for (int j = 0; j < npt; j++)
{ {