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

* Improved the Extended STL export code to hopefully speed it up substantially

Browse Source
This commit is contained in:
Philippose Rajan 2010-08-18 19:54:12 +00:00
parent 196caf7b67
commit 4fbca00961
1 changed files with 21 additions and 10 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

31
libsrc/interface/writeuser.cpp
View File

@ -349,7 +349,6 @@ void WriteSTLFormat (const Mesh & mesh,
* STL into named boundary faces * STL into named boundary faces
* when using a third-party mesher * when using a third-party mesher
*/ */
void WriteSTLExtFormat (const Mesh & mesh, void WriteSTLExtFormat (const Mesh & mesh,
const string & filename) const string & filename)
{ {
@ -359,33 +358,45 @@ void WriteSTLExtFormat (const Mesh & mesh,
outfile.precision(10); outfile.precision(10);
int numBCs = 0;
Array<int> faceBCs; Array<int> faceBCs;
TABLE<int> faceBCMapping;
faceBCs.SetSize(mesh.GetNFD());
faceBCMapping.SetSize(mesh.GetNFD());
faceBCs = -1;
// Collect the BC numbers used in the mesh // Collect the BC numbers used in the mesh
for(int faceNr = 1; faceNr <= mesh.GetNFD(); faceNr++) for(int faceNr = 1; faceNr <= mesh.GetNFD(); faceNr++)
{ {
int bcNum = mesh.GetFaceDescriptor(faceNr).BCProperty(); int bcNum = mesh.GetFaceDescriptor(faceNr).BCProperty();
if(!faceBCs.Contains(bcNum)) if(faceBCs.Pos(bcNum) < 0)
{ {
faceBCs.Append(bcNum); numBCs++;
faceBCs.Set(numBCs,bcNum);
faceBCMapping.Add1(numBCs,faceNr);
}
else
{
faceBCMapping.Add1(faceBCs.Pos(bcNum)+1,faceNr);
} }
} }
faceBCs.SetSize(numBCs);
faceBCMapping.ChangeSize(numBCs);
// Now actually write the data to file // Now actually write the data to file
for(int bcInd = 1; bcInd <= faceBCs.Size(); bcInd++) for(int bcInd = 1; bcInd <= faceBCs.Size(); bcInd++)
{ {
outfile << "solid Boundary_" << faceBCs.Elem(bcInd) << "\n"; outfile << "solid Boundary_" << faceBCs.Elem(bcInd) << "\n";
for(int faceNr = 1;faceNr <= mesh.GetNFD(); faceNr++) for(int faceNr = 1;faceNr <= faceBCMapping.EntrySize(bcInd); faceNr++)
{ {
if(mesh.GetFaceDescriptor(faceNr).BCProperty() != faceBCs.Elem(bcInd))
{
continue;
}
Array<SurfaceElementIndex> faceSei; Array<SurfaceElementIndex> faceSei;
mesh.GetSurfaceElementsOfFace(faceNr,faceSei); mesh.GetSurfaceElementsOfFace(faceBCMapping.Get(bcInd,faceNr),faceSei);
for (int i = 0; i < faceSei.Size(); i++) for (int i = 0; i < faceSei.Size(); i++)
{ {