wip

WeiqunZhang · Sep 7, 2024 · 477648e · 477648e
1 parent 908c38d
commit 477648e
Show file tree

Hide file tree

Showing 2 changed files with 19 additions and 9 deletions.
diff --git a/Src/EB/AMReX_EB_STL_utils.H b/Src/EB/AMReX_EB_STL_utils.H
@@ -28,6 +28,7 @@ public:
  RealBox boundingbox;
  STLtools::Triangle traiangles[M];
  XDim3 trinorm[M];
+ int ntriangles;
  int children[N];
  };
 

diff --git a/Src/EB/AMReX_EB_STL_utils.cpp b/Src/EB/AMReX_EB_STL_utils.cpp
@@ -458,9 +458,14 @@ STLtools::prepare (Gpu::PinnedVector<Triangle>&& a_tri_pts)
 void
 STLtools::build_bvh (Triangle* begin, Triangle* end)
 {
+ amrex::Abort("xxxxx");
+
  auto ntri = int(end - begin);
 
- if (ntri <= m_bvh_max_size) { return; }
+ if (ntri <= m_bvh_max_size) {
+ // xxxxx update boundingbox .....
+ return;
+ }
 
  RealVect centmin(std::numeric_limits<Real>::max());
  RealVect centmax(std::numeric_limits<Real>::lowest());
@@ -477,16 +482,20 @@ STLtools::build_bvh (Triangle* begin, Triangle* end)
  int tsize = ntri / nsplits;
  int nleft = ntri - tsize*nsplits;
 
- for (int isplit = 0; isplit < nsplits; ++isplit) {
- int tbegin, tend;
- if (isplit < nleft) {
- tbegin = isplit * (tsize+1);
- tend = tbegin + tsize;
+ for (int isplit = 0; isplit < m_bvh_max_splits; ++isplit) {
+ if (isplit < nsplits) {
+ int tbegin, tend;
+ if (isplit < nleft) {
+ tbegin = isplit * (tsize+1);
+ tend = tbegin + tsize;
+ } else {
+ tbegin = isplit * tsize + nleft;
+ tend = tbegin + tsize - 1;
+ }
+ build_bvh(begin+tbegin, begin+tend);
  } else {
- tbegin = isplit * tsize + nleft;
- tend = tbegin + tsize - 1;
+ // xxxxx
  }
- build_bvh(begin+tbegin, begin+tend);
  }
 }