fixed bug

This commit is contained in:
Julius Zimmermann 2019-07-15 16:42:13 +02:00
parent ef545c5d86
commit bffdf689ae

View File

@ -4,8 +4,6 @@
#include <mystdlib.h> #include <mystdlib.h>
#include <myadt.hpp> #include <myadt.hpp>
#include <linalg.hpp> #include <linalg.hpp>
#include <csg.hpp> #include <csg.hpp>
@ -242,6 +240,7 @@ namespace netgen
} }
cout << mesh.GetNE() << " elements found" << endl; cout << mesh.GetNE() << " elements found" << endl;
cout << mesh.GetNSE() << " surface elements found" << endl; cout << mesh.GetNSE() << " surface elements found" << endl;
} }
else if(strcmp (reco, "2467") == 0) else if(strcmp (reco, "2467") == 0)
{ {
@ -374,11 +373,11 @@ namespace netgen
if(dim == 2){ if(dim == 2){
// loop through segments to assign default BC to unmarked edges // loop through segments to assign default BC to unmarked edges
int bccounter_tmp = bccounter; int bccounter_tmp = bccounter;
for(int index=0; index < mesh.GetNSeg(); index++){ for(int index=1; index <= mesh.GetNSeg(); index++){
Segment & seg = mesh.LineSegment(get<0>(element_map[index])); Segment & seg = mesh.LineSegment(get<0>(element_map[index]));
if(seg.si == -1){ if(seg.si == -1){
seg.si = bccounter; seg.si = bccounter + 1;
mesh.SetBCName(bccounter, "default"); // could be more efficient if(bccounter_tmp == bccounter) mesh.SetBCName(bccounter, "default"); // could be more efficient
seg.SetBCName(mesh.GetBCNamePtr(bccounter)); seg.SetBCName(mesh.GetBCNamePtr(bccounter));
bccounter_tmp++; bccounter_tmp++;
} }
@ -397,7 +396,6 @@ namespace netgen
for(int i=0; i<bccounter; i++){ for(int i=0; i<bccounter; i++){
cout << mesh.GetBCName(i) << endl; cout << mesh.GetBCName(i) << endl;
} }
} }