boundary layer

This commit is contained in:
Joachim Schoeberl 2015-01-07 10:58:09 +00:00
parent 9a403ce44b
commit 6b0079cca3
2 changed files with 5 additions and 5 deletions

View File

@ -482,13 +482,13 @@ namespace netgen
mesh.AddVolumeElement(el); mesh.AddVolumeElement(el);
numprisms++; numprisms++;
*/ */
cout << "add element: " << endl; // cout << "add element: " << endl;
int classify = 0; int classify = 0;
for (int j = 0; j < 3; j++) for (int j = 0; j < 3; j++)
if (mapto[sel[j]]) if (mapto[sel[j]])
classify += (1 << j); classify += (1 << j);
cout << "classify = " << classify << endl; // cout << "classify = " << classify << endl;
ELEMENT_TYPE types[] = { PRISM, TET, TET, PYRAMID, ELEMENT_TYPE types[] = { PRISM, TET, TET, PYRAMID,
TET, PYRAMID, PYRAMID, PRISM }; TET, PYRAMID, PYRAMID, PRISM };
@ -588,15 +588,14 @@ namespace netgen
volel.Delete(); volel.Delete();
} }
} }
}
mesh.Compress();
}
else else
{ {
mesh.Point(i).SetPoint(pointtomove + layerht * growthvectors.Elem(i)); mesh.Point(i).SetPoint(pointtomove + layerht * growthvectors.Elem(i));
} }
} }
} }
mesh.Compress();
} }
// Optimise the tet part of the volume mesh after all the modifications // Optimise the tet part of the volume mesh after all the modifications

View File

@ -17,6 +17,7 @@ public:
int matnr = 1; int matnr = 1;
double hfirst = 0.01; double hfirst = 0.01;
double growthfactor = 1; double growthfactor = 1;
bool optimize = true;
}; };
extern void GenerateBoundaryLayer (Mesh & mesh, BoundaryLayerParameters & blp); extern void GenerateBoundaryLayer (Mesh & mesh, BoundaryLayerParameters & blp);