39 #ifndef OPENVDB_TOOLS_LEVELSETSPHERE_HAS_BEEN_INCLUDED 40 #define OPENVDB_TOOLS_LEVELSETSPHERE_HAS_BEEN_INCLUDED 47 #include <type_traits> 49 #include <tbb/enumerable_thread_specific.h> 50 #include <tbb/parallel_for.h> 51 #include <tbb/parallel_reduce.h> 52 #include <tbb/blocked_range.h> 74 template<
typename Gr
idType,
typename InterruptT>
75 typename GridType::Ptr
78 InterruptT* interrupt =
nullptr,
bool threaded =
true);
93 template<
typename Gr
idType>
94 typename GridType::Ptr
98 return createLevelSetSphere<GridType, util::NullInterrupter>(radius,center,voxelSize,halfWidth,
nullptr,threaded);
111 template<
typename Gr
idT,
typename InterruptT = util::NullInterrupter>
115 using TreeT =
typename GridT::TreeType;
116 using ValueT =
typename GridT::ValueType;
118 static_assert(std::is_floating_point<ValueT>::value,
119 "level set grids must have scalar, floating-point value types");
132 : mRadius(radius), mCenter(center), mInterrupt(interrupt)
144 mGrid = createLevelSet<GridT>(voxelSize, halfWidth);
145 this->rasterSphere(voxelSize, halfWidth, threaded);
157 const ValueT r0 = mRadius/dx, rmax = r0 + w;
160 if (r0 < 1.5f)
return;
163 const Vec3T c(mCenter[0]/dx, mCenter[1]/dx, mCenter[2]/dx);
171 typename GridT::Accessor accessor = mGrid->getAccessor();
173 if (mInterrupt) mInterrupt->start(
"Generating level set of sphere");
175 tbb::enumerable_thread_specific<TreeT> pool(mGrid->tree());
177 auto kernel = [&](
const tbb::blocked_range<int>& r) {
179 int &i = ijk[0], &j = ijk[1], &k = ijk[2], m=1;
180 TreeT &tree = pool.local();
181 typename GridT::Accessor acc(tree);
183 for (i = r.begin(); i <= r.end(); ++i) {
186 for (j = jmin; j <= jmax; ++j) {
188 for (k = kmin; k <= kmax; k += m) {
194 acc.setValue(ijk, dx*v);
208 tbb::parallel_for(tbb::blocked_range<int>(imin, imax, 128), kernel);
209 using RangeT = tbb::blocked_range<typename tbb::enumerable_thread_specific<TreeT>::iterator>;
213 Op(
TreeT &tree) : mDelete(
false), mTree(&tree) {}
214 Op(
const Op& other, tbb::split) : mDelete(
true), mTree(
new TreeT(other.mTree->background())) {}
215 ~Op() {
if (mDelete)
delete mTree; }
216 void operator()(RangeT &r) {
for (
auto i=r.begin(); i!=r.end(); ++i) this->merge(*i);}
217 void join(Op &other) { this->merge(*(other.mTree)); }
219 } op( mGrid->tree() );
220 tbb::parallel_reduce(RangeT(pool.begin(), pool.end(), 4), op);
222 kernel(tbb::blocked_range<int>(imin, imax));
229 if (mInterrupt) mInterrupt->end();
234 InterruptT* mInterrupt;
235 typename GridT::Ptr mGrid;
242 template<
typename Gr
idType,
typename InterruptT>
243 typename GridType::Ptr
245 float halfWidth, InterruptT* interrupt,
bool threaded)
248 static_assert(std::is_floating_point<typename GridType::ValueType>::value,
249 "level set grids must have scalar, floating-point value types");
251 using ValueT =
typename GridType::ValueType;
253 return factory.
getLevelSet(ValueT(voxelSize), ValueT(halfWidth), threaded);
260 #endif // OPENVDB_TOOLS_LEVELSETSPHERE_HAS_BEEN_INCLUDED
General-purpose arithmetic and comparison routines, most of which accept arbitrary value types (or at...
Coord Abs(const Coord &xyz)
Definition: Coord.h:542
#define OPENVDB_THROW(exception, message)
Definition: Exceptions.h:109
Signed (x, y, z) 32-bit integer coordinates.
Definition: Coord.h:52
static const Real LEVEL_SET_HALF_WIDTH
Definition: Types.h:510
float Sqrt(float x)
Return the square root of a floating-point value.
Definition: Math.h:735
#define OPENVDB_VERSION_NAME
The version namespace name for this library version.
Definition: version.h:128
Propagate the signs of distance values from the active voxels in the narrow band to the inactive valu...
Type Pow2(Type x)
Return x2.
Definition: Math.h:522
Definition: Exceptions.h:40
Definition: Exceptions.h:92
#define OPENVDB_USE_VERSION_NAMESPACE
Definition: version.h:180
int Ceil(float x)
Return the ceiling of x.
Definition: Math.h:830
int Floor(float x)
Return the floor of x.
Definition: Math.h:822
bool wasInterrupted(T *i, int percent=-1)
Definition: NullInterrupter.h:76