GCC Code Coverage Report


Directory: ./
File: openvdb/openvdb/tree/LeafNode.h
Date: 2022-07-25 17:40:05
Exec Total Coverage
Lines: 391 415 94.2%
Functions: 559 1220 45.8%
Branches: 600 1568 38.3%

Line Branch Exec Source
1 // Copyright Contributors to the OpenVDB Project
2 // SPDX-License-Identifier: MPL-2.0
3
4 #ifndef OPENVDB_TREE_LEAFNODE_HAS_BEEN_INCLUDED
5 #define OPENVDB_TREE_LEAFNODE_HAS_BEEN_INCLUDED
6
7 #include <openvdb/Types.h>
8 #include <openvdb/util/NodeMasks.h>
9 #include <openvdb/io/Compression.h> // for io::readData(), etc.
10 #include "Iterator.h"
11 #include "LeafBuffer.h"
12 #include <algorithm> // for std::nth_element()
13 #include <iostream>
14 #include <memory>
15 #include <sstream>
16 #include <string>
17 #include <type_traits>
18 #include <vector>
19
20
21 class TestLeaf;
22 template<typename> class TestLeafIO;
23
24 namespace openvdb {
25 OPENVDB_USE_VERSION_NAMESPACE
26 namespace OPENVDB_VERSION_NAME {
27 namespace tree {
28
29 template<Index, typename> struct SameLeafConfig; // forward declaration
30
31
32 /// @brief Templated block class to hold specific data types and a fixed
33 /// number of values determined by Log2Dim. The actual coordinate
34 /// dimension of the block is 2^Log2Dim, i.e. Log2Dim=3 corresponds to
35 /// a LeafNode that spans a 8^3 block.
36 template<typename T, Index Log2Dim>
37 16 class LeafNode
38 {
39 public:
40 using BuildType = T;
41 using ValueType = T;
42 using Buffer = LeafBuffer<ValueType, Log2Dim>;
43 using LeafNodeType = LeafNode<ValueType, Log2Dim>;
44 using NodeMaskType = util::NodeMask<Log2Dim>;
45 using Ptr = SharedPtr<LeafNode>;
46
47 static const Index
48 LOG2DIM = Log2Dim, // needed by parent nodes
49 TOTAL = Log2Dim, // needed by parent nodes
50 DIM = 1 << TOTAL, // dimension along one coordinate direction
51 NUM_VALUES = 1 << 3 * Log2Dim,
52 NUM_VOXELS = NUM_VALUES, // total number of voxels represented by this node
53 SIZE = NUM_VALUES,
54 LEVEL = 0; // level 0 = leaf
55
56 /// @brief ValueConverter<T>::Type is the type of a LeafNode having the same
57 /// dimensions as this node but a different value type, T.
58 template<typename OtherValueType>
59 struct ValueConverter { using Type = LeafNode<OtherValueType, Log2Dim>; };
60
61 /// @brief SameConfiguration<OtherNodeType>::value is @c true if and only if
62 /// OtherNodeType is the type of a LeafNode with the same dimensions as this node.
63 template<typename OtherNodeType>
64 struct SameConfiguration {
65 static const bool value = SameLeafConfig<LOG2DIM, OtherNodeType>::value;
66 };
67
68
69 /// Default constructor
70 LeafNode();
71
72 /// @brief Constructor
73 /// @param coords the grid index coordinates of a voxel
74 /// @param value a value with which to fill the buffer
75 /// @param active the active state to which to initialize all voxels
76 explicit LeafNode(const Coord& coords,
77
1/2
✓ Branch 1 taken 1 times.
✗ Branch 2 not taken.
1 const ValueType& value = zeroVal<ValueType>(),
78 bool active = false);
79
80 /// @brief "Partial creation" constructor used during file input
81 /// @param coords the grid index coordinates of a voxel
82 /// @param value a value with which to fill the buffer
83 /// @param active the active state to which to initialize all voxels
84 /// @details This constructor does not allocate memory for voxel values.
85 LeafNode(PartialCreate,
86 const Coord& coords,
87 const ValueType& value = zeroVal<ValueType>(),
88 bool active = false);
89
90 /// Deep copy constructor
91 LeafNode(const LeafNode&);
92
93 /// Deep assignment operator
94 LeafNode& operator=(const LeafNode&) = default;
95
96 /// Value conversion copy constructor
97 template<typename OtherValueType>
98 explicit LeafNode(const LeafNode<OtherValueType, Log2Dim>& other);
99
100 /// Topology copy constructor
101 template<typename OtherValueType>
102 LeafNode(const LeafNode<OtherValueType, Log2Dim>& other,
103 const ValueType& offValue, const ValueType& onValue, TopologyCopy);
104
105 /// Topology copy constructor
106 template<typename OtherValueType>
107 LeafNode(const LeafNode<OtherValueType, Log2Dim>& other,
108 const ValueType& background, TopologyCopy);
109
110 /// Destructor.
111 ~LeafNode();
112
113 //
114 // Statistics
115 //
116 /// Return log2 of the dimension of this LeafNode, e.g. 3 if dimensions are 8^3
117 static Index log2dim() { return Log2Dim; }
118 /// Return the number of voxels in each coordinate dimension.
119 static Index dim() { return DIM; }
120 /// Return the total number of voxels represented by this LeafNode
121 static Index size() { return SIZE; }
122 /// Return the total number of voxels represented by this LeafNode
123 static Index numValues() { return SIZE; }
124 /// Return the level of this node, which by definition is zero for LeafNodes
125 static Index getLevel() { return LEVEL; }
126 /// Append the Log2Dim of this LeafNode to the specified vector
127
2/4
✓ Branch 1 taken 1 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 1 times.
✗ Branch 5 not taken.
18792 static void getNodeLog2Dims(std::vector<Index>& dims) { dims.push_back(Log2Dim); }
128 /// Return the dimension of child nodes of this LeafNode, which is one for voxels.
129 static Index getChildDim() { return 1; }
130 /// Return the leaf count for this node, which is one.
131 static Index32 leafCount() { return 1; }
132 /// no-op
133 void nodeCount(std::vector<Index32> &) const {}
134 /// Return the non-leaf count for this node, which is zero.
135 static Index32 nonLeafCount() { return 0; }
136 /// Return the child count for this node, which is zero.
137 static Index32 childCount() { return 0; }
138
139 /// Return the number of voxels marked On.
140
33/60
✓ Branch 0 taken 11711 times.
✓ Branch 1 taken 52642 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 1 times.
✓ Branch 4 taken 21370 times.
✓ Branch 5 taken 1 times.
✗ Branch 6 not taken.
✓ Branch 7 taken 3 times.
✗ Branch 8 not taken.
✓ Branch 9 taken 1 times.
✓ Branch 10 taken 2 times.
✓ Branch 11 taken 1 times.
✓ Branch 13 taken 2 times.
✗ Branch 14 not taken.
✓ Branch 16 taken 2 times.
✗ Branch 17 not taken.
✓ Branch 19 taken 2701 times.
✗ Branch 20 not taken.
✓ Branch 22 taken 1 times.
✗ Branch 23 not taken.
✓ Branch 25 taken 1 times.
✗ Branch 26 not taken.
✓ Branch 28 taken 1 times.
✗ Branch 29 not taken.
✓ Branch 31 taken 1 times.
✗ Branch 32 not taken.
✓ Branch 34 taken 1 times.
✗ Branch 35 not taken.
✓ Branch 37 taken 1 times.
✗ Branch 38 not taken.
✓ Branch 40 taken 1 times.
✗ Branch 41 not taken.
✓ Branch 43 taken 1 times.
✗ Branch 44 not taken.
✓ Branch 46 taken 1 times.
✗ Branch 47 not taken.
✓ Branch 49 taken 1 times.
✗ Branch 50 not taken.
✓ Branch 52 taken 1 times.
✗ Branch 53 not taken.
✓ Branch 55 taken 1 times.
✗ Branch 56 not taken.
✓ Branch 58 taken 1 times.
✗ Branch 59 not taken.
✓ Branch 61 taken 1 times.
✗ Branch 62 not taken.
✓ Branch 64 taken 1 times.
✗ Branch 65 not taken.
✓ Branch 67 taken 1 times.
✗ Branch 68 not taken.
✓ Branch 70 taken 1 times.
✗ Branch 71 not taken.
✓ Branch 73 taken 1 times.
✗ Branch 74 not taken.
✓ Branch 76 taken 1 times.
✗ Branch 77 not taken.
✓ Branch 79 taken 1 times.
✗ Branch 80 not taken.
✓ Branch 82 taken 1 times.
✗ Branch 83 not taken.
112566 Index64 onVoxelCount() const { return mValueMask.countOn(); }
141 /// Return the number of voxels marked Off.
142
1/2
✓ Branch 1 taken 20520 times.
✗ Branch 2 not taken.
41100 Index64 offVoxelCount() const { return mValueMask.countOff(); }
143 Index64 onLeafVoxelCount() const { return onVoxelCount(); }
144 Index64 offLeafVoxelCount() const { return offVoxelCount(); }
145 static Index64 onTileCount() { return 0; }
146 static Index64 offTileCount() { return 0; }
147 /// Return @c true if this node has no active voxels.
148 bool isEmpty() const { return mValueMask.isOff(); }
149 /// Return @c true if this node contains only active voxels.
150 bool isDense() const { return mValueMask.isOn(); }
151 /// Return @c true if memory for this node's buffer has been allocated.
152
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 4 times.
8 bool isAllocated() const { return !mBuffer.isOutOfCore() && !mBuffer.empty(); }
153 /// Allocate memory for this node's buffer if it has not already been allocated.
154 106 bool allocate() { return mBuffer.allocate(); }
155
156 /// Return the memory in bytes occupied by this node.
157 Index64 memUsage() const;
158 Index64 memUsageIfLoaded() const;
159
160 /// Expand the given bounding box so that it includes this leaf node's active voxels.
161 /// If visitVoxels is false this LeafNode will be approximated as dense, i.e. with all
162 /// voxels active. Else the individual active voxels are visited to produce a tight bbox.
163 void evalActiveBoundingBox(CoordBBox& bbox, bool visitVoxels = true) const;
164
165 /// @brief Return the bounding box of this node, i.e., the full index space
166 /// spanned by this leaf node.
167 1314416 CoordBBox getNodeBoundingBox() const { return CoordBBox::createCube(mOrigin, DIM); }
168
169 /// Set the grid index coordinates of this node's local origin.
170
1/2
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
130719136 void setOrigin(const Coord& origin) { mOrigin = origin; }
171 //@{
172 /// Return the grid index coordinates of this node's local origin.
173
19/79
✓ Branch 1 taken 1419 times.
✓ Branch 2 taken 1 times.
✓ Branch 3 taken 11139 times.
✗ Branch 4 not taken.
✓ Branch 5 taken 1895 times.
✓ Branch 6 taken 62 times.
✓ Branch 7 taken 2 times.
✓ Branch 8 taken 4 times.
✓ Branch 9 taken 62 times.
✓ Branch 10 taken 4 times.
✓ Branch 11 taken 60 times.
✗ Branch 12 not taken.
✗ Branch 13 not taken.
✓ Branch 14 taken 1430 times.
✓ Branch 15 taken 1 times.
✓ Branch 16 taken 1 times.
✗ Branch 17 not taken.
✗ Branch 18 not taken.
✗ Branch 19 not taken.
✓ Branch 20 taken 9 times.
✗ Branch 21 not taken.
✗ Branch 22 not taken.
✗ Branch 23 not taken.
✗ Branch 24 not taken.
✗ Branch 25 not taken.
✓ Branch 26 taken 1 times.
✗ Branch 27 not taken.
✓ Branch 28 taken 1 times.
✓ Branch 29 taken 5 times.
✓ Branch 30 taken 1 times.
✗ Branch 31 not taken.
✗ Branch 32 not taken.
✗ Branch 33 not taken.
✗ Branch 34 not taken.
✗ Branch 35 not taken.
✗ Branch 36 not taken.
✗ Branch 37 not taken.
✗ Branch 38 not taken.
✗ Branch 39 not taken.
✗ Branch 40 not taken.
✗ Branch 41 not taken.
✗ Branch 42 not taken.
✓ Branch 43 taken 45 times.
✗ Branch 44 not taken.
✗ Branch 45 not taken.
✗ Branch 46 not taken.
✗ Branch 47 not taken.
✗ Branch 48 not taken.
✗ Branch 49 not taken.
✗ Branch 50 not taken.
✗ Branch 51 not taken.
✗ Branch 52 not taken.
✗ Branch 53 not taken.
✗ Branch 54 not taken.
✗ Branch 55 not taken.
✗ Branch 56 not taken.
✗ Branch 57 not taken.
✗ Branch 58 not taken.
✗ Branch 59 not taken.
✗ Branch 60 not taken.
✗ Branch 61 not taken.
✗ Branch 62 not taken.
✗ Branch 63 not taken.
✗ Branch 64 not taken.
✗ Branch 65 not taken.
✗ Branch 66 not taken.
✗ Branch 67 not taken.
✗ Branch 68 not taken.
✗ Branch 69 not taken.
✗ Branch 70 not taken.
✗ Branch 71 not taken.
✗ Branch 72 not taken.
✗ Branch 73 not taken.
✗ Branch 74 not taken.
✗ Branch 75 not taken.
✗ Branch 76 not taken.
✗ Branch 77 not taken.
✗ Branch 78 not taken.
✗ Branch 79 not taken.
273381 const Coord& origin() const { return mOrigin; }
174 void getOrigin(Coord& origin) const { origin = mOrigin; }
175 void getOrigin(Int32& x, Int32& y, Int32& z) const { mOrigin.asXYZ(x, y, z); }
176 //@}
177
178 /// Return the linear table offset of the given global or local coordinates.
179 static Index coordToOffset(const Coord& xyz);
180 /// @brief Return the local coordinates for a linear table offset,
181 /// where offset 0 has coordinates (0, 0, 0).
182 static Coord offsetToLocalCoord(Index n);
183 /// Return the global coordinates for a linear table offset.
184 Coord offsetToGlobalCoord(Index n) const;
185
186 #if OPENVDB_ABI_VERSION_NUMBER >= 9
187 /// Return the transient data value.
188
3/6
✓ Branch 1 taken 1 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 1 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 1 times.
✗ Branch 8 not taken.
3 Index32 transientData() const { return mTransientData; }
189 /// Set the transient data value.
190
1/2
✓ Branch 1 taken 1 times.
✗ Branch 2 not taken.
1 void setTransientData(Index32 transientData) { mTransientData = transientData; }
191 #endif
192
193 /// Return a string representation of this node.
194 std::string str() const;
195
196 /// @brief Return @c true if the given node (which may have a different @c ValueType
197 /// than this node) has the same active value topology as this node.
198 template<typename OtherType, Index OtherLog2Dim>
199 bool hasSameTopology(const LeafNode<OtherType, OtherLog2Dim>* other) const;
200
201 /// Check for buffer, state and origin equivalence.
202 bool operator==(const LeafNode& other) const;
203
4/8
✓ Branch 1 taken 1 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✓ Branch 4 taken 1 times.
✓ Branch 6 taken 1 times.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 1 times.
2 bool operator!=(const LeafNode& other) const { return !(other == *this); }
204
205 protected:
206 using MaskOnIterator = typename NodeMaskType::OnIterator;
207 using MaskOffIterator = typename NodeMaskType::OffIterator;
208 using MaskDenseIterator = typename NodeMaskType::DenseIterator;
209
210 // Type tags to disambiguate template instantiations
211 struct ValueOn {}; struct ValueOff {}; struct ValueAll {};
212 struct ChildOn {}; struct ChildOff {}; struct ChildAll {};
213
214 template<typename MaskIterT, typename NodeT, typename ValueT, typename TagT>
215 struct ValueIter:
216 // Derives from SparseIteratorBase, but can also be used as a dense iterator,
217 // if MaskIterT is a dense mask iterator type.
218 public SparseIteratorBase<
219 MaskIterT, ValueIter<MaskIterT, NodeT, ValueT, TagT>, NodeT, ValueT>
220 {
221 using BaseT = SparseIteratorBase<MaskIterT, ValueIter, NodeT, ValueT>;
222
223 ValueIter() {}
224 ValueIter(const MaskIterT& iter, NodeT* parent): BaseT(iter, parent) {}
225
226 130791909 ValueT& getItem(Index pos) const { return this->parent().getValue(pos); }
227 422721684 ValueT& getValue() const { return this->parent().getValue(this->pos()); }
228
229 // Note: setItem() can't be called on const iterators.
230 void setItem(Index pos, const ValueT& value) const
231 {
232 this->parent().setValueOnly(pos, value);
233 }
234 // Note: setValue() can't be called on const iterators.
235 50653358 void setValue(const ValueT& value) const
236 {
237 50653358 this->parent().setValueOnly(this->pos(), value);
238 50653358 }
239
240 // Note: modifyItem() can't be called on const iterators.
241 template<typename ModifyOp>
242 void modifyItem(Index n, const ModifyOp& op) const { this->parent().modifyValue(n, op); }
243 // Note: modifyValue() can't be called on const iterators.
244 template<typename ModifyOp>
245 2754 void modifyValue(const ModifyOp& op) const { this->parent().modifyValue(this->pos(), op); }
246 };
247
248 /// Leaf nodes have no children, so their child iterators have no get/set accessors.
249 template<typename MaskIterT, typename NodeT, typename TagT>
250 struct ChildIter:
251 public SparseIteratorBase<MaskIterT, ChildIter<MaskIterT, NodeT, TagT>, NodeT, ValueType>
252 {
253 ChildIter() {}
254 ChildIter(const MaskIterT& iter, NodeT* parent): SparseIteratorBase<
255 MaskIterT, ChildIter<MaskIterT, NodeT, TagT>, NodeT, ValueType>(iter, parent) {}
256 };
257
258 template<typename NodeT, typename ValueT, typename TagT>
259 struct DenseIter: public DenseIteratorBase<
260 MaskDenseIterator, DenseIter<NodeT, ValueT, TagT>, NodeT, /*ChildT=*/void, ValueT>
261 {
262 using BaseT = DenseIteratorBase<MaskDenseIterator, DenseIter, NodeT, void, ValueT>;
263 using NonConstValueT = typename BaseT::NonConstValueType;
264
265 DenseIter() {}
266 DenseIter(const MaskDenseIterator& iter, NodeT* parent): BaseT(iter, parent) {}
267
268 32768 bool getItem(Index pos, void*& child, NonConstValueT& value) const
269 {
270 49152 value = this->parent().getValue(pos);
271 32768 child = nullptr;
272 32768 return false; // no child
273 }
274
275 // Note: setItem() can't be called on const iterators.
276 //void setItem(Index pos, void* child) const {}
277
278 // Note: unsetItem() can't be called on const iterators.
279 void unsetItem(Index pos, const ValueT& value) const
280 {
281 this->parent().setValueOnly(pos, value);
282 }
283 };
284
285 public:
286 using ValueOnIter = ValueIter<MaskOnIterator, LeafNode, const ValueType, ValueOn>;
287 using ValueOnCIter = ValueIter<MaskOnIterator, const LeafNode, const ValueType, ValueOn>;
288 using ValueOffIter = ValueIter<MaskOffIterator, LeafNode, const ValueType, ValueOff>;
289 using ValueOffCIter = ValueIter<MaskOffIterator,const LeafNode,const ValueType,ValueOff>;
290 using ValueAllIter = ValueIter<MaskDenseIterator, LeafNode, const ValueType, ValueAll>;
291 using ValueAllCIter = ValueIter<MaskDenseIterator,const LeafNode,const ValueType,ValueAll>;
292 using ChildOnIter = ChildIter<MaskOnIterator, LeafNode, ChildOn>;
293 using ChildOnCIter = ChildIter<MaskOnIterator, const LeafNode, ChildOn>;
294 using ChildOffIter = ChildIter<MaskOffIterator, LeafNode, ChildOff>;
295 using ChildOffCIter = ChildIter<MaskOffIterator, const LeafNode, ChildOff>;
296 using ChildAllIter = DenseIter<LeafNode, ValueType, ChildAll>;
297 using ChildAllCIter = DenseIter<const LeafNode, const ValueType, ChildAll>;
298
299
45/139
✗ Branch 1 not taken.
✓ Branch 2 taken 26 times.
✗ Branch 4 not taken.
✓ Branch 5 taken 6391 times.
✗ Branch 7 not taken.
✓ Branch 8 taken 334 times.
✗ Branch 10 not taken.
✓ Branch 11 taken 11877 times.
✗ Branch 13 not taken.
✓ Branch 14 taken 2688 times.
✗ Branch 16 not taken.
✓ Branch 17 taken 564 times.
✗ Branch 19 not taken.
✓ Branch 20 taken 542 times.
✗ Branch 22 not taken.
✓ Branch 23 taken 676 times.
✗ Branch 25 not taken.
✓ Branch 26 taken 592 times.
✗ Branch 28 not taken.
✓ Branch 29 taken 27 times.
✗ Branch 31 not taken.
✗ Branch 32 not taken.
✓ Branch 33 taken 566 times.
✗ Branch 34 not taken.
✓ Branch 35 taken 5 times.
✗ Branch 36 not taken.
✗ Branch 37 not taken.
✗ Branch 38 not taken.
✗ Branch 40 not taken.
✗ Branch 41 not taken.
✗ Branch 43 not taken.
✗ Branch 44 not taken.
✗ Branch 46 not taken.
✗ Branch 47 not taken.
✗ Branch 48 not taken.
✗ Branch 49 not taken.
✗ Branch 50 not taken.
✗ Branch 51 not taken.
✗ Branch 52 not taken.
✗ Branch 53 not taken.
✗ Branch 54 not taken.
✗ Branch 55 not taken.
✗ Branch 56 not taken.
✗ Branch 57 not taken.
✗ Branch 58 not taken.
✗ Branch 59 not taken.
✗ Branch 60 not taken.
✗ Branch 61 not taken.
✗ Branch 62 not taken.
✗ Branch 63 not taken.
✗ Branch 64 not taken.
✗ Branch 65 not taken.
✗ Branch 66 not taken.
✗ Branch 67 not taken.
✗ Branch 68 not taken.
✗ Branch 70 not taken.
✗ Branch 71 not taken.
✗ Branch 73 not taken.
✗ Branch 74 not taken.
✓ Branch 75 taken 1 times.
✗ Branch 76 not taken.
✗ Branch 77 not taken.
✓ Branch 79 taken 1 times.
✗ Branch 80 not taken.
✗ Branch 82 not taken.
✓ Branch 83 taken 1 times.
✗ Branch 84 not taken.
✗ Branch 85 not taken.
✗ Branch 86 not taken.
✓ Branch 87 taken 1 times.
✗ Branch 88 not taken.
✗ Branch 89 not taken.
✓ Branch 91 taken 1 times.
✗ Branch 92 not taken.
✗ Branch 94 not taken.
✓ Branch 95 taken 1 times.
✗ Branch 96 not taken.
✗ Branch 97 not taken.
✗ Branch 98 not taken.
✓ Branch 99 taken 1 times.
✗ Branch 100 not taken.
✗ Branch 101 not taken.
✓ Branch 103 taken 1 times.
✗ Branch 104 not taken.
✗ Branch 106 not taken.
✓ Branch 107 taken 1 times.
✗ Branch 108 not taken.
✗ Branch 109 not taken.
✗ Branch 110 not taken.
✓ Branch 111 taken 1 times.
✗ Branch 112 not taken.
✗ Branch 113 not taken.
✓ Branch 115 taken 1 times.
✗ Branch 116 not taken.
✗ Branch 118 not taken.
✓ Branch 119 taken 1 times.
✗ Branch 120 not taken.
✓ Branch 123 taken 1 times.
✗ Branch 124 not taken.
✓ Branch 127 taken 1 times.
✗ Branch 128 not taken.
✓ Branch 131 taken 1 times.
✗ Branch 132 not taken.
✓ Branch 135 taken 1 times.
✗ Branch 136 not taken.
✓ Branch 139 taken 1 times.
✗ Branch 140 not taken.
✓ Branch 143 taken 1 times.
✗ Branch 144 not taken.
✓ Branch 147 taken 1 times.
✗ Branch 148 not taken.
✓ Branch 151 taken 1 times.
✗ Branch 152 not taken.
✓ Branch 155 taken 1 times.
✗ Branch 156 not taken.
✓ Branch 159 taken 1 times.
✗ Branch 160 not taken.
✓ Branch 163 taken 1 times.
✗ Branch 164 not taken.
✓ Branch 167 taken 1 times.
✗ Branch 168 not taken.
✓ Branch 171 taken 1 times.
✗ Branch 172 not taken.
✓ Branch 175 taken 1 times.
✗ Branch 176 not taken.
✓ Branch 179 taken 1 times.
✗ Branch 180 not taken.
✓ Branch 183 taken 1 times.
✗ Branch 184 not taken.
✓ Branch 187 taken 1 times.
✗ Branch 188 not taken.
✓ Branch 191 taken 1 times.
✗ Branch 192 not taken.
✓ Branch 195 taken 1 times.
✗ Branch 196 not taken.
✓ Branch 199 taken 1 times.
✗ Branch 200 not taken.
✓ Branch 203 taken 1 times.
✗ Branch 204 not taken.
2779551 ValueOnCIter cbeginValueOn() const { return ValueOnCIter(mValueMask.beginOn(), this); }
300 20124 ValueOnCIter beginValueOn() const { return ValueOnCIter(mValueMask.beginOn(), this); }
301
6/14
✗ Branch 1 not taken.
✓ Branch 2 taken 58 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 56 times.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 10 taken 1 times.
✗ Branch 11 not taken.
✓ Branch 14 taken 1 times.
✗ Branch 15 not taken.
✓ Branch 18 taken 1 times.
✗ Branch 19 not taken.
604089 ValueOnIter beginValueOn() { return ValueOnIter(mValueMask.beginOn(), this); }
302 4764 ValueOffCIter cbeginValueOff() const { return ValueOffCIter(mValueMask.beginOff(), this); }
303 ValueOffCIter beginValueOff() const { return ValueOffCIter(mValueMask.beginOff(), this); }
304 40409 ValueOffIter beginValueOff() { return ValueOffIter(mValueMask.beginOff(), this); }
305 107698 ValueAllCIter cbeginValueAll() const { return ValueAllCIter(mValueMask.beginDense(), this); }
306 ValueAllCIter beginValueAll() const { return ValueAllCIter(mValueMask.beginDense(), this); }
307 48092 ValueAllIter beginValueAll() { return ValueAllIter(mValueMask.beginDense(), this); }
308
309 ValueOnCIter cendValueOn() const { return ValueOnCIter(mValueMask.endOn(), this); }
310 ValueOnCIter endValueOn() const { return ValueOnCIter(mValueMask.endOn(), this); }
311 ValueOnIter endValueOn() { return ValueOnIter(mValueMask.endOn(), this); }
312 ValueOffCIter cendValueOff() const { return ValueOffCIter(mValueMask.endOff(), this); }
313 ValueOffCIter endValueOff() const { return ValueOffCIter(mValueMask.endOff(), this); }
314 ValueOffIter endValueOff() { return ValueOffIter(mValueMask.endOff(), this); }
315 ValueAllCIter cendValueAll() const { return ValueAllCIter(mValueMask.endDense(), this); }
316 ValueAllCIter endValueAll() const { return ValueAllCIter(mValueMask.endDense(), this); }
317 ValueAllIter endValueAll() { return ValueAllIter(mValueMask.endDense(), this); }
318
319 // Note that [c]beginChildOn() and [c]beginChildOff() actually return end iterators,
320 // because leaf nodes have no children.
321 4877008 ChildOnCIter cbeginChildOn() const { return ChildOnCIter(mValueMask.endOn(), this); }
322 ChildOnCIter beginChildOn() const { return ChildOnCIter(mValueMask.endOn(), this); }
323 8426 ChildOnIter beginChildOn() { return ChildOnIter(mValueMask.endOn(), this); }
324 ChildOffCIter cbeginChildOff() const { return ChildOffCIter(mValueMask.endOff(), this); }
325 ChildOffCIter beginChildOff() const { return ChildOffCIter(mValueMask.endOff(), this); }
326 ChildOffIter beginChildOff() { return ChildOffIter(mValueMask.endOff(), this); }
327 ChildAllCIter cbeginChildAll() const { return ChildAllCIter(mValueMask.beginDense(), this); }
328 64 ChildAllCIter beginChildAll() const { return ChildAllCIter(mValueMask.beginDense(), this); }
329 198 ChildAllIter beginChildAll() { return ChildAllIter(mValueMask.beginDense(), this); }
330
331 ChildOnCIter cendChildOn() const { return ChildOnCIter(mValueMask.endOn(), this); }
332 ChildOnCIter endChildOn() const { return ChildOnCIter(mValueMask.endOn(), this); }
333 ChildOnIter endChildOn() { return ChildOnIter(mValueMask.endOn(), this); }
334 ChildOffCIter cendChildOff() const { return ChildOffCIter(mValueMask.endOff(), this); }
335 ChildOffCIter endChildOff() const { return ChildOffCIter(mValueMask.endOff(), this); }
336 ChildOffIter endChildOff() { return ChildOffIter(mValueMask.endOff(), this); }
337 ChildAllCIter cendChildAll() const { return ChildAllCIter(mValueMask.endDense(), this); }
338 ChildAllCIter endChildAll() const { return ChildAllCIter(mValueMask.endDense(), this); }
339 ChildAllIter endChildAll() { return ChildAllIter(mValueMask.endDense(), this); }
340
341 //
342 // Buffer management
343 //
344 /// @brief Exchange this node's data buffer with the given data buffer
345 /// without changing the active states of the values.
346 41621 void swap(Buffer& other) { mBuffer.swap(other); }
347
9/135
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
✓ Branch 6 taken 126 times.
✓ Branch 7 taken 3959 times.
✗ Branch 8 not taken.
✗ Branch 10 not taken.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 252 times.
✗ Branch 14 not taken.
✗ Branch 15 not taken.
✗ Branch 16 not taken.
✗ Branch 17 not taken.
✗ Branch 18 not taken.
✓ Branch 19 taken 11327 times.
✓ Branch 20 taken 1 times.
✓ Branch 21 taken 1395 times.
✗ Branch 22 not taken.
✗ Branch 23 not taken.
✗ Branch 24 not taken.
✗ Branch 25 not taken.
✗ Branch 26 not taken.
✗ Branch 27 not taken.
✗ Branch 28 not taken.
✓ Branch 29 taken 8 times.
✗ Branch 30 not taken.
✓ Branch 31 taken 688 times.
✗ Branch 32 not taken.
✗ Branch 33 not taken.
✓ Branch 34 taken 824 times.
✗ Branch 35 not taken.
✗ Branch 36 not taken.
✗ Branch 37 not taken.
✗ Branch 38 not taken.
✗ Branch 39 not taken.
✗ Branch 40 not taken.
✗ Branch 41 not taken.
✗ Branch 42 not taken.
✗ Branch 43 not taken.
✗ Branch 44 not taken.
✗ Branch 45 not taken.
✗ Branch 46 not taken.
✗ Branch 47 not taken.
✗ Branch 48 not taken.
✗ Branch 49 not taken.
✗ Branch 50 not taken.
✗ Branch 51 not taken.
✗ Branch 52 not taken.
✗ Branch 53 not taken.
✗ Branch 54 not taken.
✗ Branch 55 not taken.
✗ Branch 56 not taken.
✗ Branch 57 not taken.
✗ Branch 58 not taken.
✗ Branch 59 not taken.
✗ Branch 60 not taken.
✗ Branch 61 not taken.
✗ Branch 62 not taken.
✗ Branch 63 not taken.
✗ Branch 64 not taken.
✗ Branch 65 not taken.
✗ Branch 66 not taken.
✗ Branch 67 not taken.
✗ Branch 68 not taken.
✗ Branch 69 not taken.
✗ Branch 70 not taken.
✗ Branch 71 not taken.
✗ Branch 72 not taken.
✗ Branch 73 not taken.
✗ Branch 74 not taken.
✗ Branch 75 not taken.
✗ Branch 76 not taken.
✗ Branch 77 not taken.
✗ Branch 78 not taken.
✗ Branch 79 not taken.
✗ Branch 80 not taken.
✗ Branch 81 not taken.
✗ Branch 82 not taken.
✗ Branch 83 not taken.
✗ Branch 84 not taken.
✗ Branch 85 not taken.
✗ Branch 86 not taken.
✗ Branch 87 not taken.
✗ Branch 88 not taken.
✗ Branch 89 not taken.
✗ Branch 90 not taken.
✗ Branch 91 not taken.
✗ Branch 92 not taken.
✗ Branch 93 not taken.
✗ Branch 94 not taken.
✗ Branch 95 not taken.
✗ Branch 96 not taken.
✗ Branch 97 not taken.
✗ Branch 98 not taken.
✗ Branch 99 not taken.
✗ Branch 100 not taken.
✗ Branch 101 not taken.
✗ Branch 102 not taken.
✗ Branch 103 not taken.
✗ Branch 104 not taken.
✗ Branch 105 not taken.
✗ Branch 106 not taken.
✗ Branch 107 not taken.
✗ Branch 108 not taken.
✗ Branch 109 not taken.
✗ Branch 110 not taken.
✗ Branch 112 not taken.
✗ Branch 113 not taken.
✗ Branch 115 not taken.
✗ Branch 116 not taken.
✗ Branch 118 not taken.
✗ Branch 119 not taken.
✗ Branch 121 not taken.
✗ Branch 122 not taken.
✗ Branch 124 not taken.
✗ Branch 125 not taken.
✗ Branch 127 not taken.
✗ Branch 128 not taken.
✗ Branch 130 not taken.
✗ Branch 131 not taken.
✗ Branch 133 not taken.
✗ Branch 134 not taken.
✗ Branch 136 not taken.
✗ Branch 137 not taken.
✗ Branch 139 not taken.
✗ Branch 140 not taken.
✗ Branch 142 not taken.
✗ Branch 143 not taken.
✗ Branch 145 not taken.
✗ Branch 146 not taken.
✗ Branch 148 not taken.
✗ Branch 149 not taken.
1404380 const Buffer& buffer() const { return mBuffer; }
348
79/253
✓ Branch 1 taken 2191 times.
✓ Branch 2 taken 1395 times.
✓ Branch 3 taken 1024 times.
✓ Branch 4 taken 135343 times.
✓ Branch 5 taken 8 times.
✓ Branch 6 taken 4096 times.
✓ Branch 7 taken 4025 times.
✓ Branch 8 taken 8 times.
✓ Branch 9 taken 4096 times.
✓ Branch 10 taken 1342 times.
✗ Branch 11 not taken.
✓ Branch 12 taken 4608 times.
✓ Branch 13 taken 115200 times.
✗ Branch 14 not taken.
✓ Branch 15 taken 4119 times.
✓ Branch 16 taken 41 times.
✓ Branch 17 taken 252 times.
✓ Branch 18 taken 5396 times.
✗ Branch 19 not taken.
✓ Branch 20 taken 11579 times.
✓ Branch 21 taken 7168 times.
✗ Branch 22 not taken.
✓ Branch 23 taken 11579 times.
✓ Branch 24 taken 1024 times.
✗ Branch 25 not taken.
✓ Branch 26 taken 1300 times.
✓ Branch 27 taken 4096 times.
✗ Branch 28 not taken.
✗ Branch 29 not taken.
✓ Branch 30 taken 1842688 times.
✗ Branch 31 not taken.
✗ Branch 32 not taken.
✓ Branch 33 taken 4096 times.
✓ Branch 34 taken 2080 times.
✓ Branch 35 taken 11327 times.
✓ Branch 36 taken 4096 times.
✗ Branch 37 not taken.
✗ Branch 38 not taken.
✗ Branch 39 not taken.
✗ Branch 40 not taken.
✗ Branch 41 not taken.
✓ Branch 42 taken 4096 times.
✓ Branch 43 taken 2671 times.
✗ Branch 44 not taken.
✓ Branch 45 taken 5266 times.
✓ Branch 46 taken 2671 times.
✗ Branch 47 not taken.
✗ Branch 48 not taken.
✓ Branch 49 taken 2671 times.
✗ Branch 50 not taken.
✓ Branch 51 taken 4096 times.
✓ Branch 52 taken 2671 times.
✗ Branch 53 not taken.
✓ Branch 54 taken 19968 times.
✗ Branch 55 not taken.
✗ Branch 56 not taken.
✓ Branch 57 taken 574 times.
✗ Branch 58 not taken.
✗ Branch 59 not taken.
✗ Branch 60 not taken.
✓ Branch 61 taken 4096 times.
✓ Branch 62 taken 1170 times.
✗ Branch 63 not taken.
✗ Branch 64 not taken.
✓ Branch 65 taken 4096 times.
✗ Branch 66 not taken.
✗ Branch 67 not taken.
✗ Branch 68 not taken.
✓ Branch 69 taken 4096 times.
✗ Branch 70 not taken.
✗ Branch 71 not taken.
✗ Branch 72 not taken.
✓ Branch 73 taken 4096 times.
✓ Branch 74 taken 67 times.
✗ Branch 75 not taken.
✗ Branch 76 not taken.
✓ Branch 77 taken 4096 times.
✗ Branch 78 not taken.
✓ Branch 79 taken 1170 times.
✗ Branch 80 not taken.
✓ Branch 81 taken 4096 times.
✗ Branch 82 not taken.
✗ Branch 83 not taken.
✗ Branch 84 not taken.
✓ Branch 85 taken 4096 times.
✗ Branch 86 not taken.
✗ Branch 87 not taken.
✗ Branch 88 not taken.
✓ Branch 89 taken 4096 times.
✗ Branch 90 not taken.
✓ Branch 91 taken 60 times.
✗ Branch 92 not taken.
✓ Branch 93 taken 4627 times.
✗ Branch 94 not taken.
✗ Branch 95 not taken.
✓ Branch 96 taken 1170 times.
✓ Branch 97 taken 6144 times.
✗ Branch 98 not taken.
✗ Branch 99 not taken.
✗ Branch 100 not taken.
✗ Branch 101 not taken.
✗ Branch 102 not taken.
✗ Branch 103 not taken.
✗ Branch 104 not taken.
✗ Branch 105 not taken.
✗ Branch 106 not taken.
✗ Branch 107 not taken.
✓ Branch 108 taken 67 times.
✗ Branch 109 not taken.
✗ Branch 110 not taken.
✗ Branch 112 not taken.
✓ Branch 113 taken 1170 times.
✗ Branch 114 not taken.
✗ Branch 115 not taken.
✗ Branch 116 not taken.
✗ Branch 117 not taken.
✗ Branch 118 not taken.
✗ Branch 119 not taken.
✗ Branch 120 not taken.
✗ Branch 121 not taken.
✗ Branch 122 not taken.
✗ Branch 123 not taken.
✗ Branch 124 not taken.
✓ Branch 125 taken 63 times.
✗ Branch 126 not taken.
✗ Branch 127 not taken.
✗ Branch 128 not taken.
✓ Branch 130 taken 1170 times.
✗ Branch 131 not taken.
✗ Branch 133 not taken.
✗ Branch 134 not taken.
✗ Branch 136 not taken.
✗ Branch 137 not taken.
✗ Branch 139 not taken.
✗ Branch 140 not taken.
✓ Branch 142 taken 64 times.
✗ Branch 143 not taken.
✗ Branch 145 not taken.
✗ Branch 146 not taken.
✓ Branch 147 taken 1560 times.
✗ Branch 148 not taken.
✗ Branch 149 not taken.
✗ Branch 150 not taken.
✗ Branch 151 not taken.
✗ Branch 152 not taken.
✗ Branch 153 not taken.
✗ Branch 154 not taken.
✗ Branch 155 not taken.
✗ Branch 156 not taken.
✗ Branch 157 not taken.
✗ Branch 158 not taken.
✓ Branch 159 taken 67 times.
✗ Branch 160 not taken.
✗ Branch 161 not taken.
✗ Branch 163 not taken.
✓ Branch 164 taken 1560 times.
✗ Branch 165 not taken.
✗ Branch 166 not taken.
✗ Branch 167 not taken.
✗ Branch 168 not taken.
✗ Branch 169 not taken.
✗ Branch 170 not taken.
✗ Branch 171 not taken.
✗ Branch 172 not taken.
✗ Branch 173 not taken.
✗ Branch 174 not taken.
✗ Branch 175 not taken.
✓ Branch 176 taken 71 times.
✗ Branch 177 not taken.
✗ Branch 178 not taken.
✗ Branch 179 not taken.
✓ Branch 181 taken 1560 times.
✗ Branch 182 not taken.
✗ Branch 184 not taken.
✗ Branch 185 not taken.
✗ Branch 187 not taken.
✗ Branch 188 not taken.
✓ Branch 190 taken 112 times.
✗ Branch 191 not taken.
✗ Branch 195 not taken.
✗ Branch 196 not taken.
✗ Branch 198 not taken.
✗ Branch 199 not taken.
✗ Branch 201 not taken.
✗ Branch 202 not taken.
✓ Branch 204 taken 128 times.
✗ Branch 205 not taken.
✓ Branch 209 taken 1040 times.
✗ Branch 210 not taken.
✗ Branch 212 not taken.
✗ Branch 213 not taken.
✗ Branch 215 not taken.
✗ Branch 216 not taken.
✗ Branch 218 not taken.
✗ Branch 219 not taken.
✓ Branch 221 taken 47 times.
✗ Branch 222 not taken.
✓ Branch 224 taken 2080 times.
✗ Branch 225 not taken.
✗ Branch 227 not taken.
✗ Branch 228 not taken.
✗ Branch 230 not taken.
✗ Branch 231 not taken.
✗ Branch 233 not taken.
✗ Branch 234 not taken.
✓ Branch 236 taken 47 times.
✗ Branch 237 not taken.
✓ Branch 241 taken 2080 times.
✗ Branch 242 not taken.
✗ Branch 244 not taken.
✗ Branch 245 not taken.
✓ Branch 247 taken 4097 times.
✗ Branch 248 not taken.
✗ Branch 250 not taken.
✗ Branch 251 not taken.
✓ Branch 253 taken 47 times.
✗ Branch 254 not taken.
✓ Branch 258 taken 4455205 times.
✗ Branch 259 not taken.
✗ Branch 261 not taken.
✗ Branch 262 not taken.
✓ Branch 266 taken 3640 times.
✗ Branch 267 not taken.
✗ Branch 269 not taken.
✗ Branch 270 not taken.
✗ Branch 272 not taken.
✗ Branch 273 not taken.
✓ Branch 275 taken 34 times.
✗ Branch 276 not taken.
✗ Branch 280 not taken.
✗ Branch 281 not taken.
✗ Branch 283 not taken.
✗ Branch 284 not taken.
✓ Branch 286 taken 280 times.
✗ Branch 287 not taken.
✗ Branch 291 not taken.
✗ Branch 292 not taken.
✗ Branch 294 not taken.
✗ Branch 295 not taken.
✓ Branch 297 taken 93 times.
✗ Branch 298 not taken.
✓ Branch 300 taken 1 times.
✗ Branch 301 not taken.
✓ Branch 303 taken 1 times.
✗ Branch 304 not taken.
✓ Branch 306 taken 280 times.
✗ Branch 307 not taken.
✗ Branch 309 not taken.
✗ Branch 310 not taken.
✗ Branch 312 not taken.
✗ Branch 313 not taken.
✓ Branch 315 taken 392 times.
✗ Branch 316 not taken.
15415992 Buffer& buffer() { return mBuffer; }
349
350 //
351 // I/O methods
352 //
353 /// @brief Read in just the topology.
354 /// @param is the stream from which to read
355 /// @param fromHalf if true, floating-point input values are assumed to be 16-bit
356 void readTopology(std::istream& is, bool fromHalf = false);
357 /// @brief Write out just the topology.
358 /// @param os the stream to which to write
359 /// @param toHalf if true, output floating-point values as 16-bit half floats
360 void writeTopology(std::ostream& os, bool toHalf = false) const;
361
362 /// @brief Read buffers from a stream.
363 /// @param is the stream from which to read
364 /// @param fromHalf if true, floating-point input values are assumed to be 16-bit
365 void readBuffers(std::istream& is, bool fromHalf = false);
366 /// @brief Read buffers that intersect the given bounding box.
367 /// @param is the stream from which to read
368 /// @param bbox an index-space bounding box
369 /// @param fromHalf if true, floating-point input values are assumed to be 16-bit
370 void readBuffers(std::istream& is, const CoordBBox& bbox, bool fromHalf = false);
371 /// @brief Write buffers to a stream.
372 /// @param os the stream to which to write
373 /// @param toHalf if true, output floating-point values as 16-bit half floats
374 void writeBuffers(std::ostream& os, bool toHalf = false) const;
375
376 size_t streamingSize(bool toHalf = false) const;
377
378 //
379 // Accessor methods
380 //
381 /// Return the value of the voxel at the given coordinates.
382 const ValueType& getValue(const Coord& xyz) const;
383 /// Return the value of the voxel at the given linear offset.
384 const ValueType& getValue(Index offset) const;
385
386 /// @brief Return @c true if the voxel at the given coordinates is active.
387 /// @param xyz the coordinates of the voxel to be probed
388 /// @param[out] val the value of the voxel at the given coordinates
389 bool probeValue(const Coord& xyz, ValueType& val) const;
390 /// @brief Return @c true if the voxel at the given offset is active.
391 /// @param offset the linear offset of the voxel to be probed
392 /// @param[out] val the value of the voxel at the given coordinates
393 bool probeValue(Index offset, ValueType& val) const;
394
395 /// Return the level (i.e., 0) at which leaf node values reside.
396 static Index getValueLevel(const Coord&) { return LEVEL; }
397
398 /// Set the active state of the voxel at the given coordinates but don't change its value.
399 void setActiveState(const Coord& xyz, bool on);
400 /// Set the active state of the voxel at the given offset but don't change its value.
401
3/4
✗ Branch 0 not taken.
✓ Branch 1 taken 438757 times.
✓ Branch 3 taken 241316 times.
✓ Branch 4 taken 197441 times.
877514 void setActiveState(Index offset, bool on) { assert(offset<SIZE); mValueMask.set(offset, on); }
402
403 /// Set the value of the voxel at the given coordinates but don't change its active state.
404 void setValueOnly(const Coord& xyz, const ValueType& val);
405 /// Set the value of the voxel at the given offset but don't change its active state.
406 void setValueOnly(Index offset, const ValueType& val);
407
408 /// Mark the voxel at the given coordinates as inactive but don't change its value.
409 6240 void setValueOff(const Coord& xyz) { mValueMask.setOff(LeafNode::coordToOffset(xyz)); }
410 /// Mark the voxel at the given offset as inactive but don't change its value.
411
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 438345 times.
438603 void setValueOff(Index offset) { assert(offset < SIZE); mValueMask.setOff(offset); }
412
413 /// Set the value of the voxel at the given coordinates and mark the voxel as inactive.
414 void setValueOff(const Coord& xyz, const ValueType& val);
415 /// Set the value of the voxel at the given offset and mark the voxel as inactive.
416 void setValueOff(Index offset, const ValueType& val);
417
418 /// Mark the voxel at the given coordinates as active but don't change its value.
419 void setValueOn(const Coord& xyz) { mValueMask.setOn(LeafNode::coordToOffset(xyz)); }
420 /// Mark the voxel at the given offset as active but don't change its value.
421
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 9461142 times.
18922294 void setValueOn(Index offset) { assert(offset < SIZE); mValueMask.setOn(offset); }
422 /// Set the value of the voxel at the given coordinates and mark the voxel as active.
423 155689312 void setValueOn(const Coord& xyz, const ValueType& val) {
424 this->setValueOn(LeafNode::coordToOffset(xyz), val);
425 155689312 }
426 /// Set the value of the voxel at the given coordinates and mark the voxel as active.
427
9/18
✓ Branch 1 taken 1 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 1 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 1 times.
✗ Branch 8 not taken.
✓ Branch 10 taken 1 times.
✗ Branch 11 not taken.
✓ Branch 13 taken 1 times.
✗ Branch 14 not taken.
✓ Branch 16 taken 1 times.
✗ Branch 17 not taken.
✓ Branch 19 taken 1 times.
✗ Branch 20 not taken.
✓ Branch 22 taken 1 times.
✗ Branch 23 not taken.
✓ Branch 25 taken 1 times.
✗ Branch 26 not taken.
9 void setValue(const Coord& xyz, const ValueType& val) { this->setValueOn(xyz, val); }
428 /// Set the value of the voxel at the given offset and mark the voxel as active.
429 void setValueOn(Index offset, const ValueType& val) {
430
16/32
✓ Branch 1 taken 89414 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 4 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 397 times.
✗ Branch 8 not taken.
✓ Branch 9 taken 169518 times.
✓ Branch 10 taken 7 times.
✗ Branch 11 not taken.
✓ Branch 12 taken 169518 times.
✓ Branch 13 taken 81359 times.
✗ Branch 14 not taken.
✓ Branch 15 taken 57263 times.
✓ Branch 16 taken 6 times.
✗ Branch 17 not taken.
✗ Branch 18 not taken.
✓ Branch 19 taken 6 times.
✗ Branch 20 not taken.
✗ Branch 21 not taken.
✓ Branch 22 taken 6 times.
✗ Branch 23 not taken.
✗ Branch 24 not taken.
✓ Branch 25 taken 1 times.
✗ Branch 26 not taken.
✓ Branch 28 taken 1 times.
✗ Branch 29 not taken.
✓ Branch 31 taken 1 times.
✗ Branch 32 not taken.
✓ Branch 34 taken 1 times.
✗ Branch 35 not taken.
✓ Branch 37 taken 1 times.
✗ Branch 38 not taken.
92443662 mBuffer.setValue(offset, val);
431 96474851 mValueMask.setOn(offset);
432 1761058 }
433
434 /// @brief Apply a functor to the value of the voxel at the given offset
435 /// and mark the voxel as active.
436 template<typename ModifyOp>
437 1079684540 void modifyValue(Index offset, const ModifyOp& op)
438 {
439 1079684540 mBuffer.loadValues();
440 if (!mBuffer.empty()) {
441 // in-place modify value
442 ValueType& val = const_cast<ValueType&>(mBuffer[offset]);
443 1026 op(val);
444 1079684540 mValueMask.setOn(offset);
445 }
446 1079684540 }
447
448 /// @brief Apply a functor to the value of the voxel at the given coordinates
449 /// and mark the voxel as active.
450 template<typename ModifyOp>
451 1078979841 void modifyValue(const Coord& xyz, const ModifyOp& op)
452 {
453 1078979841 this->modifyValue(this->coordToOffset(xyz), op);
454 1078979841 }
455
456 /// Apply a functor to the voxel at the given coordinates.
457 template<typename ModifyOp>
458 16 void modifyValueAndActiveState(const Coord& xyz, const ModifyOp& op)
459 {
460 16 mBuffer.loadValues();
461 if (!mBuffer.empty()) {
462 const Index offset = this->coordToOffset(xyz);
463 16 bool state = mValueMask.isOn(offset);
464 // in-place modify value
465 ValueType& val = const_cast<ValueType&>(mBuffer[offset]);
466 6 op(val, state);
467
1/2
✓ Branch 0 taken 3 times.
✗ Branch 1 not taken.
6 mValueMask.set(offset, state);
468 }
469 16 }
470
471 /// Mark all voxels as active but don't change their values.
472 void setValuesOn() { mValueMask.setOn(); }
473 /// Mark all voxels as inactive but don't change their values.
474 void setValuesOff() { mValueMask.setOff(); }
475
476 /// Return @c true if the voxel at the given coordinates is active.
477 44324056 bool isValueOn(const Coord& xyz) const {return this->isValueOn(LeafNode::coordToOffset(xyz));}
478 /// Return @c true if the voxel at the given offset is active.
479 40614225 bool isValueOn(Index offset) const { return mValueMask.isOn(offset); }
480
481 /// Return @c false since leaf nodes never contain tiles.
482 static bool hasActiveTiles() { return false; }
483
484 /// Set all voxels that lie outside the given axis-aligned box to the background.
485 void clip(const CoordBBox&, const ValueType& background);
486
487 /// Set all voxels within an axis-aligned box to the specified value and active state.
488 void fill(const CoordBBox& bbox, const ValueType&, bool active = true);
489 /// Set all voxels within an axis-aligned box to the specified value and active state.
490 void denseFill(const CoordBBox& bbox, const ValueType& value, bool active = true)
491 {
492 11920 this->fill(bbox, value, active);
493 }
494
495 /// Set all voxels to the specified value but don't change their active states.
496 void fill(const ValueType& value);
497 /// Set all voxels to the specified value and active state.
498 void fill(const ValueType& value, bool active);
499
500 /// @brief Copy into a dense grid the values of the voxels that lie within
501 /// a given bounding box.
502 ///
503 /// @param bbox inclusive bounding box of the voxels to be copied into the dense grid
504 /// @param dense dense grid with a stride in @e z of one (see tools::Dense
505 /// in tools/Dense.h for the required API)
506 ///
507 /// @note @a bbox is assumed to be identical to or contained in the coordinate domains
508 /// of both the dense grid and this node, i.e., no bounds checking is performed.
509 /// @note Consider using tools::CopyToDense in tools/Dense.h
510 /// instead of calling this method directly.
511 template<typename DenseT>
512 void copyToDense(const CoordBBox& bbox, DenseT& dense) const;
513
514 /// @brief Copy from a dense grid into this node the values of the voxels
515 /// that lie within a given bounding box.
516 /// @details Only values that are different (by more than the given tolerance)
517 /// from the background value will be active. Other values are inactive
518 /// and truncated to the background value.
519 ///
520 /// @param bbox inclusive bounding box of the voxels to be copied into this node
521 /// @param dense dense grid with a stride in @e z of one (see tools::Dense
522 /// in tools/Dense.h for the required API)
523 /// @param background background value of the tree that this node belongs to
524 /// @param tolerance tolerance within which a value equals the background value
525 ///
526 /// @note @a bbox is assumed to be identical to or contained in the coordinate domains
527 /// of both the dense grid and this node, i.e., no bounds checking is performed.
528 /// @note Consider using tools::CopyFromDense in tools/Dense.h
529 /// instead of calling this method directly.
530 template<typename DenseT>
531 void copyFromDense(const CoordBBox& bbox, const DenseT& dense,
532 const ValueType& background, const ValueType& tolerance);
533
534 /// @brief Return the value of the voxel at the given coordinates.
535 /// @note Used internally by ValueAccessor.
536 template<typename AccessorT>
537 const ValueType& getValueAndCache(const Coord& xyz, AccessorT&) const
538 {
539 844314473 return this->getValue(xyz);
540 }
541
542 /// @brief Return @c true if the voxel at the given coordinates is active.
543 /// @note Used internally by ValueAccessor.
544 template<typename AccessorT>
545 9005668 bool isValueOnAndCache(const Coord& xyz, AccessorT&) const { return this->isValueOn(xyz); }
546
547 /// @brief Change the value of the voxel at the given coordinates and mark it as active.
548 /// @note Used internally by ValueAccessor.
549 template<typename AccessorT>
550 void setValueAndCache(const Coord& xyz, const ValueType& val, AccessorT&)
551 {
552 93772111 this->setValueOn(xyz, val);
553 93772111 }
554
555 /// @brief Change the value of the voxel at the given coordinates
556 /// but preserve its state.
557 /// @note Used internally by ValueAccessor.
558 template<typename AccessorT>
559 void setValueOnlyAndCache(const Coord& xyz, const ValueType& val, AccessorT&)
560 {
561 8136073 this->setValueOnly(xyz, val);
562 8136073 }
563
564 /// @brief Apply a functor to the value of the voxel at the given coordinates
565 /// and mark the voxel as active.
566 /// @note Used internally by ValueAccessor.
567 template<typename ModifyOp, typename AccessorT>
568 void modifyValueAndCache(const Coord& xyz, const ModifyOp& op, AccessorT&)
569 {
570 539489919 this->modifyValue(xyz, op);
571 539489919 }
572
573 /// Apply a functor to the voxel at the given coordinates.
574 /// @note Used internally by ValueAccessor.
575 template<typename ModifyOp, typename AccessorT>
576 void modifyValueAndActiveStateAndCache(const Coord& xyz, const ModifyOp& op, AccessorT&)
577 {
578 8 this->modifyValueAndActiveState(xyz, op);
579 8 }
580
581 /// @brief Change the value of the voxel at the given coordinates and mark it as inactive.
582 /// @note Used internally by ValueAccessor.
583 template<typename AccessorT>
584 void setValueOffAndCache(const Coord& xyz, const ValueType& value, AccessorT&)
585 {
586 5685953 this->setValueOff(xyz, value);
587 5685953 }
588
589 /// @brief Set the active state of the voxel at the given coordinates
590 /// without changing its value.
591 /// @note Used internally by ValueAccessor.
592 template<typename AccessorT>
593 void setActiveStateAndCache(const Coord& xyz, bool on, AccessorT&)
594 {
595 177034 this->setActiveState(xyz, on);
596 177034 }
597
598 /// @brief Return @c true if the voxel at the given coordinates is active
599 /// and return the voxel value in @a val.
600 /// @note Used internally by ValueAccessor.
601 template<typename AccessorT>
602 bool probeValueAndCache(const Coord& xyz, ValueType& val, AccessorT&) const
603 {
604 999273928 return this->probeValue(xyz, val);
605 }
606
607 /// @brief Return the value of the voxel at the given coordinates and return
608 /// its active state and level (i.e., 0) in @a state and @a level.
609 /// @note Used internally by ValueAccessor.
610 template<typename AccessorT>
611 const ValueType& getValue(const Coord& xyz, bool& state, int& level, AccessorT&) const
612 {
613 const Index offset = this->coordToOffset(xyz);
614 state = mValueMask.isOn(offset);
615 level = LEVEL;
616 return mBuffer[offset];
617 }
618
619 /// @brief Return the LEVEL (=0) at which leaf node values reside.
620 /// @note Used internally by ValueAccessor (note last argument is a dummy).
621 template<typename AccessorT>
622 static Index getValueLevelAndCache(const Coord&, AccessorT&) { return LEVEL; }
623
624 /// @brief Return a const reference to the first value in the buffer.
625 /// @note Though it is potentially risky you can convert this
626 /// to a non-const pointer by means of const_case<ValueType*>&.
627
2/4
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✓ Branch 5 taken 1 times.
✗ Branch 6 not taken.
2441 const ValueType& getFirstValue() const { return mBuffer[0]; }
628 /// Return a const reference to the last value in the buffer.
629
5/8
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 4 times.
✓ Branch 5 taken 1 times.
✗ Branch 6 not taken.
✓ Branch 7 taken 1 times.
✗ Branch 8 not taken.
385583 const ValueType& getLastValue() const { return mBuffer[SIZE - 1]; }
630
631 /// @brief Replace inactive occurrences of @a oldBackground with @a newBackground,
632 /// and inactive occurrences of @a -oldBackground with @a -newBackground.
633 void resetBackground(const ValueType& oldBackground, const ValueType& newBackground);
634
635 void negate();
636
637 /// @brief No-op
638 /// @details This function exists only to enable template instantiation.
639 void voxelizeActiveTiles(bool = true) {}
640
641 template<MergePolicy Policy> void merge(const LeafNode&);
642 template<MergePolicy Policy> void merge(const ValueType& tileValue, bool tileActive);
643 template<MergePolicy Policy>
644 void merge(const LeafNode& other, const ValueType& /*bg*/, const ValueType& /*otherBG*/);
645
646 /// @brief Union this node's set of active values with the active values
647 /// of the other node, whose @c ValueType may be different. So a
648 /// resulting voxel will be active if either of the original voxels
649 /// were active.
650 ///
651 /// @note This operation modifies only active states, not values.
652 template<typename OtherType>
653 void topologyUnion(const LeafNode<OtherType, Log2Dim>& other, const bool preserveTiles = false);
654
655 /// @brief Intersect this node's set of active values with the active values
656 /// of the other node, whose @c ValueType may be different. So a
657 /// resulting voxel will be active only if both of the original voxels
658 /// were active.
659 ///
660 /// @details The last dummy argument is required to match the signature
661 /// for InternalNode::topologyIntersection.
662 ///
663 /// @note This operation modifies only active states, not
664 /// values. Also note that this operation can result in all voxels
665 /// being inactive so consider subsequently calling prune.
666 template<typename OtherType>
667 void topologyIntersection(const LeafNode<OtherType, Log2Dim>& other, const ValueType&);
668
669 /// @brief Difference this node's set of active values with the active values
670 /// of the other node, whose @c ValueType may be different. So a
671 /// resulting voxel will be active only if the original voxel is
672 /// active in this LeafNode and inactive in the other LeafNode.
673 ///
674 /// @details The last dummy argument is required to match the signature
675 /// for InternalNode::topologyDifference.
676 ///
677 /// @note This operation modifies only active states, not values.
678 /// Also, because it can deactivate all of this node's voxels,
679 /// consider subsequently calling prune.
680 template<typename OtherType>
681 void topologyDifference(const LeafNode<OtherType, Log2Dim>& other, const ValueType&);
682
683 template<typename CombineOp>
684 void combine(const LeafNode& other, CombineOp& op);
685 template<typename CombineOp>
686 void combine(const ValueType& value, bool valueIsActive, CombineOp& op);
687
688 template<typename CombineOp, typename OtherType /*= ValueType*/>
689 void combine2(const LeafNode& other, const OtherType&, bool valueIsActive, CombineOp&);
690 template<typename CombineOp, typename OtherNodeT /*= LeafNode*/>
691 void combine2(const ValueType&, const OtherNodeT& other, bool valueIsActive, CombineOp&);
692 template<typename CombineOp, typename OtherNodeT /*= LeafNode*/>
693 void combine2(const LeafNode& b0, const OtherNodeT& b1, CombineOp&);
694
695 /// @brief Calls the templated functor BBoxOp with bounding box
696 /// information. An additional level argument is provided to the
697 /// callback.
698 ///
699 /// @note The bounding boxes are guaranteed to be non-overlapping.
700 template<typename BBoxOp> void visitActiveBBox(BBoxOp&) const;
701
702 template<typename VisitorOp> void visit(VisitorOp&);
703 template<typename VisitorOp> void visit(VisitorOp&) const;
704
705 template<typename OtherLeafNodeType, typename VisitorOp>
706 void visit2Node(OtherLeafNodeType& other, VisitorOp&);
707 template<typename OtherLeafNodeType, typename VisitorOp>
708 void visit2Node(OtherLeafNodeType& other, VisitorOp&) const;
709 template<typename IterT, typename VisitorOp>
710 void visit2(IterT& otherIter, VisitorOp&, bool otherIsLHS = false);
711 template<typename IterT, typename VisitorOp>
712 void visit2(IterT& otherIter, VisitorOp&, bool otherIsLHS = false) const;
713
714 //@{
715 /// This function exists only to enable template instantiation.
716 void prune(const ValueType& /*tolerance*/ = zeroVal<ValueType>()) {}
717 void addLeaf(LeafNode*) {}
718 template<typename AccessorT>
719 void addLeafAndCache(LeafNode*, AccessorT&) {}
720 template<typename NodeT>
721 NodeT* stealNode(const Coord&, const ValueType&, bool) { return nullptr; }
722 template<typename NodeT>
723 NodeT* probeNode(const Coord&) { return nullptr; }
724 template<typename NodeT>
725 const NodeT* probeConstNode(const Coord&) const { return nullptr; }
726 template<typename ArrayT> void getNodes(ArrayT&) const {}
727 template<typename ArrayT> void stealNodes(ArrayT&, const ValueType&, bool) {}
728 //@}
729
730 void addTile(Index level, const Coord&, const ValueType&, bool);
731 void addTile(Index offset, const ValueType&, bool);
732 template<typename AccessorT>
733 void addTileAndCache(Index, const Coord&, const ValueType&, bool, AccessorT&);
734
735 //@{
736 /// @brief Return a pointer to this node.
737 LeafNode* touchLeaf(const Coord&) { return this; }
738 template<typename AccessorT>
739 LeafNode* touchLeafAndCache(const Coord&, AccessorT&) { return this; }
740 template<typename NodeT, typename AccessorT>
741 NodeT* probeNodeAndCache(const Coord&, AccessorT&)
742 {
743 OPENVDB_NO_UNREACHABLE_CODE_WARNING_BEGIN
744 if (!(std::is_same<NodeT, LeafNode>::value)) return nullptr;
745 return reinterpret_cast<NodeT*>(this);
746 OPENVDB_NO_UNREACHABLE_CODE_WARNING_END
747 }
748 LeafNode* probeLeaf(const Coord&) { return this; }
749 template<typename AccessorT>
750 LeafNode* probeLeafAndCache(const Coord&, AccessorT&) { return this; }
751 //@}
752 //@{
753 /// @brief Return a @const pointer to this node.
754 const LeafNode* probeConstLeaf(const Coord&) const { return this; }
755 template<typename AccessorT>
756 const LeafNode* probeConstLeafAndCache(const Coord&, AccessorT&) const { return this; }
757 template<typename AccessorT>
758 const LeafNode* probeLeafAndCache(const Coord&, AccessorT&) const { return this; }
759 const LeafNode* probeLeaf(const Coord&) const { return this; }
760 template<typename NodeT, typename AccessorT>
761 const NodeT* probeConstNodeAndCache(const Coord&, AccessorT&) const
762 {
763 OPENVDB_NO_UNREACHABLE_CODE_WARNING_BEGIN
764 if (!(std::is_same<NodeT, LeafNode>::value)) return nullptr;
765 return reinterpret_cast<const NodeT*>(this);
766 OPENVDB_NO_UNREACHABLE_CODE_WARNING_END
767 }
768 //@}
769
770 /// Return @c true if all of this node's values have the same active state
771 /// and are in the range this->getFirstValue() +/- @a tolerance.
772 ///
773 ///
774 /// @param firstValue Is updated with the first value of this leaf node.
775 /// @param state Is updated with the state of all values IF method
776 /// returns @c true. Else the value is undefined!
777 /// @param tolerance The tolerance used to determine if values are
778 /// approximately equal to the for value.
779 bool isConstant(ValueType& firstValue, bool& state,
780 const ValueType& tolerance = zeroVal<ValueType>()) const;
781
782 /// Return @c true if all of this node's values have the same active state
783 /// and the range (@a maxValue - @a minValue) < @a tolerance.
784 ///
785 /// @param minValue Is updated with the minimum of all values IF method
786 /// returns @c true. Else the value is undefined!
787 /// @param maxValue Is updated with the maximum of all values IF method
788 /// returns @c true. Else the value is undefined!
789 /// @param state Is updated with the state of all values IF method
790 /// returns @c true. Else the value is undefined!
791 /// @param tolerance The tolerance used to determine if values are
792 /// approximately constant.
793 bool isConstant(ValueType& minValue, ValueType& maxValue,
794 bool& state, const ValueType& tolerance = zeroVal<ValueType>()) const;
795
796
797 /// @brief Computes the median value of all the active AND inactive voxels in this node.
798 /// @return The median value of all values in this node.
799 ///
800 /// @param tmp Optional temporary storage that can hold at least NUM_VALUES values
801 /// Use of this temporary storage can improve performance
802 /// when this method is called multiple times.
803 ///
804 /// @note If tmp = this->buffer().data() then the median
805 /// value is computed very efficiently (in place) but
806 /// the voxel values in this node are re-shuffled!
807 ///
808 /// @warning If tmp != nullptr then it is the responsibility of
809 /// the client code that it points to enough memory to
810 /// hold NUM_VALUES elements of type ValueType.
811 ValueType medianAll(ValueType *tmp = nullptr) const;
812
813 /// @brief Computes the median value of all the active voxels in this node.
814 /// @return The number of active voxels.
815 ///
816 /// @param value If the return value is non zero @a value is updated
817 /// with the median value.
818 ///
819 /// @param tmp Optional temporary storage that can hold at least
820 /// as many values as there are active voxels in this node.
821 /// Use of this temporary storage can improve performance
822 /// when this method is called multiple times.
823 ///
824 /// @warning If tmp != nullptr then it is the responsibility of
825 /// the client code that it points to enough memory to
826 /// hold the number of active voxels of type ValueType.
827 Index medianOn(ValueType &value, ValueType *tmp = nullptr) const;
828
829 /// @brief Computes the median value of all the inactive voxels in this node.
830 /// @return The number of inactive voxels.
831 ///
832 /// @param value If the return value is non zero @a value is updated
833 /// with the median value.
834 ///
835 /// @param tmp Optional temporary storage that can hold at least
836 /// as many values as there are inactive voxels in this node.
837 /// Use of this temporary storage can improve performance
838 /// when this method is called multiple times.
839 ///
840 /// @warning If tmp != nullptr then it is the responsibility of
841 /// the client code that it points to enough memory to
842 /// hold the number of inactive voxels of type ValueType.
843 Index medianOff(ValueType &value, ValueType *tmp = nullptr) const;
844
845 /// Return @c true if all of this node's values are inactive.
846 bool isInactive() const { return mValueMask.isOff(); }
847
848 protected:
849 friend class ::TestLeaf;
850 template<typename> friend class ::TestLeafIO;
851
852 // During topology-only construction, access is needed
853 // to protected/private members of other template instances.
854 template<typename, Index> friend class LeafNode;
855
856 friend struct ValueIter<MaskOnIterator, LeafNode, ValueType, ValueOn>;
857 friend struct ValueIter<MaskOffIterator, LeafNode, ValueType, ValueOff>;
858 friend struct ValueIter<MaskDenseIterator, LeafNode, ValueType, ValueAll>;
859 friend struct ValueIter<MaskOnIterator, const LeafNode, ValueType, ValueOn>;
860 friend struct ValueIter<MaskOffIterator, const LeafNode, ValueType, ValueOff>;
861 friend struct ValueIter<MaskDenseIterator, const LeafNode, ValueType, ValueAll>;
862
863 // Allow iterators to call mask accessor methods (see below).
864 /// @todo Make mask accessors public?
865 friend class IteratorBase<MaskOnIterator, LeafNode>;
866 friend class IteratorBase<MaskOffIterator, LeafNode>;
867 friend class IteratorBase<MaskDenseIterator, LeafNode>;
868
869 // Mask accessors
870 public:
871 7540057 bool isValueMaskOn(Index n) const { return mValueMask.isOn(n); }
872 bool isValueMaskOn() const { return mValueMask.isOn(); }
873 bool isValueMaskOff(Index n) const { return mValueMask.isOff(n); }
874 bool isValueMaskOff() const { return mValueMask.isOff(); }
875
4/32
✗ Branch 0 not taken.
✗ Branch 1 not taken.
✓ Branch 2 taken 1894 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 41 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✗ Branch 9 not taken.
✓ Branch 10 taken 41 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✗ Branch 13 not taken.
✗ Branch 14 not taken.
✗ Branch 15 not taken.
✓ Branch 16 taken 41 times.
✗ Branch 17 not taken.
✗ Branch 18 not taken.
✗ Branch 19 not taken.
✗ Branch 20 not taken.
✗ Branch 21 not taken.
✗ Branch 22 not taken.
✗ Branch 23 not taken.
✗ Branch 24 not taken.
✗ Branch 25 not taken.
✗ Branch 26 not taken.
✗ Branch 27 not taken.
✗ Branch 28 not taken.
✗ Branch 29 not taken.
✗ Branch 30 not taken.
✗ Branch 31 not taken.
329759 const NodeMaskType& getValueMask() const { return mValueMask; }
876
3/32
✗ Branch 0 not taken.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✓ Branch 4 taken 41 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✗ Branch 9 not taken.
✓ Branch 10 taken 41 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✗ Branch 13 not taken.
✗ Branch 14 not taken.
✗ Branch 15 not taken.
✓ Branch 16 taken 41 times.
✗ Branch 17 not taken.
✗ Branch 18 not taken.
✗ Branch 19 not taken.
✗ Branch 20 not taken.
✗ Branch 21 not taken.
✗ Branch 22 not taken.
✗ Branch 23 not taken.
✗ Branch 24 not taken.
✗ Branch 25 not taken.
✗ Branch 26 not taken.
✗ Branch 27 not taken.
✗ Branch 28 not taken.
✗ Branch 29 not taken.
✗ Branch 30 not taken.
✗ Branch 31 not taken.
394727 NodeMaskType& getValueMask() { return mValueMask; }
877 50973 const NodeMaskType& valueMask() const { return mValueMask; }
878 void setValueMask(const NodeMaskType& mask) { mValueMask = mask; }
879 bool isChildMaskOn(Index) const { return false; } // leaf nodes have no children
880 bool isChildMaskOff(Index) const { return true; }
881 bool isChildMaskOff() const { return true; }
882 protected:
883
10/32
✓ Branch 0 taken 54 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 3326941 times.
✗ Branch 6 not taken.
✓ Branch 7 taken 3 times.
✓ Branch 8 taken 278 times.
✓ Branch 9 taken 7918 times.
✓ Branch 10 taken 3894 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 2 times.
✗ Branch 14 not taken.
✗ Branch 15 not taken.
✗ Branch 16 not taken.
✗ Branch 17 not taken.
✗ Branch 18 not taken.
✗ Branch 19 not taken.
✓ Branch 20 taken 3326941 times.
✗ Branch 21 not taken.
✗ Branch 22 not taken.
✗ Branch 23 not taken.
✓ Branch 24 taken 4 times.
✗ Branch 25 not taken.
✗ Branch 26 not taken.
✗ Branch 27 not taken.
✓ Branch 28 taken 2 times.
✗ Branch 29 not taken.
✗ Branch 30 not taken.
✗ Branch 31 not taken.
6666037 void setValueMask(Index n, bool on) { mValueMask.set(n, on); }
884 5498740 void setValueMaskOn(Index n) { mValueMask.setOn(n); }
885 void setValueMaskOff(Index n) { mValueMask.setOff(n); }
886
887 inline void skipCompressedValues(bool seekable, std::istream&, bool fromHalf);
888
889 /// Compute the origin of the leaf node that contains the voxel with the given coordinates.
890 static void evalNodeOrigin(Coord& xyz) { xyz &= ~(DIM - 1); }
891
892 template<typename NodeT, typename VisitorOp, typename ChildAllIterT>
893 static inline void doVisit(NodeT&, VisitorOp&);
894
895 template<typename NodeT, typename OtherNodeT, typename VisitorOp,
896 typename ChildAllIterT, typename OtherChildAllIterT>
897 static inline void doVisit2Node(NodeT& self, OtherNodeT& other, VisitorOp&);
898
899 template<typename NodeT, typename VisitorOp,
900 typename ChildAllIterT, typename OtherChildAllIterT>
901 static inline void doVisit2(NodeT& self, OtherChildAllIterT&, VisitorOp&, bool otherIsLHS);
902
903 private:
904 /// Buffer containing the actual data values
905 Buffer mBuffer;
906 /// Bitmask that determines which voxels are active
907 NodeMaskType mValueMask;
908 /// Global grid index coordinates (x,y,z) of the local origin of this node
909 Coord mOrigin;
910 #if OPENVDB_ABI_VERSION_NUMBER >= 9
911 /// Transient data (not serialized)
912 Index32 mTransientData = 0;
913 #endif
914 }; // end of LeafNode class
915
916
917 ////////////////////////////////////////
918
919
920 //@{
921 /// Helper metafunction used to implement LeafNode::SameConfiguration
922 /// (which, as an inner class, can't be independently specialized)
923 template<Index Dim1, typename NodeT2>
924 struct SameLeafConfig { static const bool value = false; };
925
926 template<Index Dim1, typename T2>
927 struct SameLeafConfig<Dim1, LeafNode<T2, Dim1> > { static const bool value = true; };
928 //@}
929
930
931 ////////////////////////////////////////
932
933
934 template<typename T, Index Log2Dim>
935 inline
936 44883 LeafNode<T, Log2Dim>::LeafNode():
937 mValueMask(),//default is off!
938
5/13
✓ Branch 1 taken 333 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 12 times.
✓ Branch 4 taken 400 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 333 times.
✗ Branch 8 not taken.
✗ Branch 9 not taken.
✓ Branch 10 taken 400 times.
✗ Branch 11 not taken.
✗ Branch 13 not taken.
✗ Branch 14 not taken.
46362 mOrigin(0, 0, 0)
939 {
940 }
941
942
943 template<typename T, Index Log2Dim>
944 inline
945 11780157 LeafNode<T, Log2Dim>::LeafNode(const Coord& xyz, const ValueType& val, bool active):
946 mBuffer(val),
947 mValueMask(active),
948 11780157 mOrigin(xyz & (~(DIM - 1)))
949 {
950 11780157 }
951
952
953 template<typename T, Index Log2Dim>
954 inline
955
1/2
✓ Branch 0 taken 23074 times.
✗ Branch 1 not taken.
46148 LeafNode<T, Log2Dim>::LeafNode(PartialCreate, const Coord& xyz, const ValueType& val, bool active):
956 mBuffer(PartialCreate(), val),
957 mValueMask(active),
958 46148 mOrigin(xyz & (~(DIM - 1)))
959 {
960 46148 }
961
962
963 template<typename T, Index Log2Dim>
964 inline
965 87898 LeafNode<T, Log2Dim>::LeafNode(const LeafNode& other)
966 87898 : mBuffer(other.mBuffer)
967 , mValueMask(other.valueMask())
968 , mOrigin(other.mOrigin)
969 #if OPENVDB_ABI_VERSION_NUMBER >= 9
970
0/2
✗ Branch 1 not taken.
✗ Branch 2 not taken.
87898 , mTransientData(other.mTransientData)
971 #endif
972 {
973 87898 }
974
975
976 // Copy-construct from a leaf node with the same configuration but a different ValueType.
977 template<typename T, Index Log2Dim>
978 template<typename OtherValueType>
979 inline
980 33625 LeafNode<T, Log2Dim>::LeafNode(const LeafNode<OtherValueType, Log2Dim>& other)
981 : mValueMask(other.valueMask())
982 , mOrigin(other.mOrigin)
983 #if OPENVDB_ABI_VERSION_NUMBER >= 9
984 67250 , mTransientData(other.mTransientData)
985 #endif
986 {
987 struct Local {
988 /// @todo Consider using a value conversion functor passed as an argument instead.
989 17219584 static inline ValueType convertValue(const OtherValueType& val) { return ValueType(val); }
990 };
991
992
2/2
✓ Branch 0 taken 17217536 times.
✓ Branch 1 taken 33621 times.
17253209 for (Index i = 0; i < SIZE; ++i) {
993
1/2
✓ Branch 1 taken 17217536 times.
✗ Branch 2 not taken.
34439168 mBuffer[i] = Local::convertValue(other.mBuffer[i]);
994 }
995 33625 }
996
997
998 template<typename T, Index Log2Dim>
999 template<typename OtherValueType>
1000 inline
1001 85826 LeafNode<T, Log2Dim>::LeafNode(const LeafNode<OtherValueType, Log2Dim>& other,
1002 const ValueType& background, TopologyCopy)
1003 : mBuffer(background)
1004 , mValueMask(other.valueMask())
1005 , mOrigin(other.mOrigin)
1006 #if OPENVDB_ABI_VERSION_NUMBER >= 9
1007
23/100
✓ Branch 1 taken 5669 times.
✓ Branch 2 taken 3848 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 9515 times.
✓ Branch 5 taken 2526 times.
✓ Branch 6 taken 32 times.
✓ Branch 7 taken 751 times.
✓ Branch 8 taken 136 times.
✗ Branch 9 not taken.
✓ Branch 10 taken 4078 times.
✗ Branch 11 not taken.
✓ Branch 12 taken 10519 times.
✓ Branch 13 taken 16510 times.
✓ Branch 14 taken 16 times.
✓ Branch 15 taken 3254 times.
✗ Branch 16 not taken.
✗ Branch 17 not taken.
✗ Branch 18 not taken.
✓ Branch 19 taken 16049 times.
✗ Branch 20 not taken.
✗ Branch 21 not taken.
✓ Branch 22 taken 1833 times.
✗ Branch 23 not taken.
✗ Branch 24 not taken.
✗ Branch 25 not taken.
✓ Branch 26 taken 793 times.
✗ Branch 27 not taken.
✓ Branch 28 taken 3517 times.
✓ Branch 29 taken 6 times.
✗ Branch 30 not taken.
✗ Branch 31 not taken.
✗ Branch 32 not taken.
✓ Branch 33 taken 4 times.
✓ Branch 34 taken 24 times.
✗ Branch 35 not taken.
✗ Branch 36 not taken.
✗ Branch 37 not taken.
✓ Branch 38 taken 72 times.
✗ Branch 39 not taken.
✗ Branch 40 not taken.
✗ Branch 41 not taken.
✗ Branch 42 not taken.
✗ Branch 43 not taken.
✗ Branch 44 not taken.
✗ Branch 45 not taken.
✗ Branch 46 not taken.
✗ Branch 47 not taken.
✗ Branch 48 not taken.
✗ Branch 49 not taken.
✗ Branch 50 not taken.
✓ Branch 51 taken 18 times.
✗ Branch 52 not taken.
✗ Branch 53 not taken.
✗ Branch 55 not taken.
✗ Branch 56 not taken.
✗ Branch 58 not taken.
✗ Branch 59 not taken.
✗ Branch 60 not taken.
✗ Branch 61 not taken.
✗ Branch 62 not taken.
✗ Branch 63 not taken.
✗ Branch 64 not taken.
✗ Branch 65 not taken.
✗ Branch 67 not taken.
✗ Branch 68 not taken.
✗ Branch 70 not taken.
✗ Branch 71 not taken.
✗ Branch 72 not taken.
✗ Branch 73 not taken.
✗ Branch 74 not taken.
✗ Branch 76 not taken.
✗ Branch 77 not taken.
✗ Branch 79 not taken.
✗ Branch 80 not taken.
✗ Branch 82 not taken.
✗ Branch 83 not taken.
✗ Branch 85 not taken.
✗ Branch 86 not taken.
✗ Branch 88 not taken.
✗ Branch 89 not taken.
✗ Branch 91 not taken.
✗ Branch 92 not taken.
✗ Branch 94 not taken.
✗ Branch 95 not taken.
✗ Branch 97 not taken.
✗ Branch 98 not taken.
✗ Branch 100 not taken.
✗ Branch 101 not taken.
✗ Branch 103 not taken.
✗ Branch 104 not taken.
✗ Branch 106 not taken.
✗ Branch 107 not taken.
✗ Branch 109 not taken.
✗ Branch 110 not taken.
✓ Branch 112 taken 1408 times.
✗ Branch 113 not taken.
✓ Branch 115 taken 5248 times.
✗ Branch 116 not taken.
✗ Branch 118 not taken.
✗ Branch 119 not taken.
171652 , mTransientData(other.mTransientData)
1008 #endif
1009 {
1010 }
1011
1012
1013 template<typename T, Index Log2Dim>
1014 template<typename OtherValueType>
1015 inline
1016 9 LeafNode<T, Log2Dim>::LeafNode(const LeafNode<OtherValueType, Log2Dim>& other,
1017 const ValueType& offValue, const ValueType& onValue, TopologyCopy)
1018 : mValueMask(other.valueMask())
1019 , mOrigin(other.mOrigin)
1020 #if OPENVDB_ABI_VERSION_NUMBER >= 9
1021 18 , mTransientData(other.mTransientData)
1022 #endif
1023 {
1024
2/2
✓ Branch 0 taken 2560 times.
✓ Branch 1 taken 5 times.
4617 for (Index i = 0; i < SIZE; ++i) {
1025
2/2
✓ Branch 1 taken 512 times.
✓ Branch 2 taken 2048 times.
9216 mBuffer[i] = (mValueMask.isOn(i) ? onValue : offValue);
1026 }
1027 9 }
1028
1029
1030 template<typename T, Index Log2Dim>
1031 inline
1032 LeafNode<T, Log2Dim>::~LeafNode()
1033 {
1034
1/2
✓ Branch 11 taken 1 times.
✗ Branch 12 not taken.
6079330 }
1035
1036
1037 template<typename T, Index Log2Dim>
1038 inline std::string
1039 LeafNode<T, Log2Dim>::str() const
1040 {
1041 std::ostringstream ostr;
1042 ostr << "LeafNode @" << mOrigin << ": " << mBuffer;
1043 return ostr.str();
1044 }
1045
1046
1047 ////////////////////////////////////////
1048
1049
1050 template<typename T, Index Log2Dim>
1051 inline Index
1052 LeafNode<T, Log2Dim>::coordToOffset(const Coord& xyz)
1053 {
1054 2545225891 assert ((xyz[0] & (DIM-1u)) < DIM && (xyz[1] & (DIM-1u)) < DIM && (xyz[2] & (DIM-1u)) < DIM);
1055 2545225891 return ((xyz[0] & (DIM-1u)) << 2*Log2Dim)
1056 2545225891 + ((xyz[1] & (DIM-1u)) << Log2Dim)
1057
19/82
✗ Branch 0 not taken.
✓ Branch 1 taken 5640 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 264 times.
✓ Branch 4 taken 17 times.
✓ Branch 5 taken 4 times.
✓ Branch 6 taken 257026 times.
✗ Branch 7 not taken.
✓ Branch 8 taken 17 times.
✗ Branch 9 not taken.
✓ Branch 10 taken 41344 times.
✓ Branch 11 taken 4 times.
✗ Branch 12 not taken.
✗ Branch 13 not taken.
✓ Branch 14 taken 3 times.
✗ Branch 15 not taken.
✗ Branch 16 not taken.
✓ Branch 17 taken 15 times.
✗ Branch 18 not taken.
✗ Branch 19 not taken.
✓ Branch 20 taken 1 times.
✗ Branch 21 not taken.
✗ Branch 22 not taken.
✓ Branch 23 taken 8 times.
✓ Branch 24 taken 2 times.
✗ Branch 25 not taken.
✓ Branch 26 taken 6 times.
✗ Branch 27 not taken.
✗ Branch 28 not taken.
✓ Branch 29 taken 15 times.
✗ Branch 30 not taken.
✗ Branch 31 not taken.
✓ Branch 32 taken 1654104 times.
✗ Branch 33 not taken.
✗ Branch 34 not taken.
✗ Branch 35 not taken.
✗ Branch 37 not taken.
✗ Branch 38 not taken.
✗ Branch 40 not taken.
✗ Branch 41 not taken.
✗ Branch 43 not taken.
✗ Branch 44 not taken.
✗ Branch 46 not taken.
✗ Branch 47 not taken.
✗ Branch 49 not taken.
✗ Branch 50 not taken.
✗ Branch 52 not taken.
✗ Branch 53 not taken.
✗ Branch 55 not taken.
✗ Branch 56 not taken.
✗ Branch 58 not taken.
✗ Branch 59 not taken.
✗ Branch 61 not taken.
✗ Branch 62 not taken.
✗ Branch 64 not taken.
✗ Branch 65 not taken.
✗ Branch 67 not taken.
✗ Branch 68 not taken.
✗ Branch 70 not taken.
✗ Branch 71 not taken.
✗ Branch 73 not taken.
✗ Branch 74 not taken.
✓ Branch 76 taken 4 times.
✗ Branch 77 not taken.
✓ Branch 79 taken 4 times.
✗ Branch 80 not taken.
✗ Branch 82 not taken.
✗ Branch 83 not taken.
✗ Branch 85 not taken.
✗ Branch 86 not taken.
✗ Branch 88 not taken.
✗ Branch 89 not taken.
✗ Branch 91 not taken.
✗ Branch 92 not taken.
✗ Branch 94 not taken.
✗ Branch 95 not taken.
✗ Branch 97 not taken.
✗ Branch 98 not taken.
✗ Branch 100 not taken.
✗ Branch 101 not taken.
✓ Branch 103 taken 3837100 times.
✗ Branch 104 not taken.
704680159 + (xyz[2] & (DIM-1u));
1058 }
1059
1060 template<typename T, Index Log2Dim>
1061 inline Coord
1062 8868062148 LeafNode<T, Log2Dim>::offsetToLocalCoord(Index n)
1063 {
1064
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 4444854948 times.
8868062148 assert(n<(1<< 3*Log2Dim));
1065 Coord xyz;
1066 8868062148 xyz.setX(n >> 2*Log2Dim);
1067 8868062148 n &= ((1<<2*Log2Dim)-1);
1068 8868062148 xyz.setY(n >> Log2Dim);
1069 8868062148 xyz.setZ(n & ((1<<Log2Dim)-1));
1070 8868062148 return xyz;
1071 }
1072
1073
1074 template<typename T, Index Log2Dim>
1075 inline Coord
1076 508882843 LeafNode<T, Log2Dim>::offsetToGlobalCoord(Index n) const
1077 {
1078 508882843 return (this->offsetToLocalCoord(n) + this->origin());
1079 }
1080
1081
1082 ////////////////////////////////////////
1083
1084
1085 template<typename ValueT, Index Log2Dim>
1086 inline const ValueT&
1087 1446046736 LeafNode<ValueT, Log2Dim>::getValue(const Coord& xyz) const
1088 {
1089 1446046736 return this->getValue(LeafNode::coordToOffset(xyz));
1090 }
1091
1092 template<typename ValueT, Index Log2Dim>
1093 inline const ValueT&
1094 2395982389 LeafNode<ValueT, Log2Dim>::getValue(Index offset) const
1095 {
1096
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 1200523297 times.
2395982389 assert(offset < SIZE);
1097
27/54
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 4 taken 1 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 1 times.
✗ Branch 8 not taken.
✓ Branch 10 taken 1 times.
✗ Branch 11 not taken.
✓ Branch 12 taken 1 times.
✓ Branch 13 taken 1 times.
✗ Branch 14 not taken.
✓ Branch 15 taken 1 times.
✓ Branch 16 taken 1 times.
✗ Branch 17 not taken.
✓ Branch 19 taken 1 times.
✗ Branch 20 not taken.
✓ Branch 22 taken 1 times.
✗ Branch 23 not taken.
✓ Branch 25 taken 1 times.
✗ Branch 26 not taken.
✓ Branch 28 taken 1 times.
✗ Branch 29 not taken.
✓ Branch 31 taken 1 times.
✗ Branch 32 not taken.
✓ Branch 34 taken 1 times.
✗ Branch 35 not taken.
✓ Branch 37 taken 1 times.
✗ Branch 38 not taken.
✓ Branch 40 taken 1 times.
✗ Branch 41 not taken.
✓ Branch 43 taken 1 times.
✗ Branch 44 not taken.
✓ Branch 46 taken 1 times.
✗ Branch 47 not taken.
✓ Branch 49 taken 1 times.
✗ Branch 50 not taken.
✓ Branch 52 taken 1 times.
✗ Branch 53 not taken.
✓ Branch 55 taken 1 times.
✗ Branch 56 not taken.
✓ Branch 58 taken 1 times.
✗ Branch 59 not taken.
✓ Branch 61 taken 1 times.
✗ Branch 62 not taken.
✓ Branch 64 taken 1 times.
✗ Branch 65 not taken.
✓ Branch 67 taken 1 times.
✗ Branch 68 not taken.
✓ Branch 70 taken 1 times.
✗ Branch 71 not taken.
✓ Branch 73 taken 1 times.
✗ Branch 74 not taken.
✓ Branch 76 taken 1 times.
✗ Branch 77 not taken.
2395993176 return mBuffer[offset];
1098 }
1099
1100
1101 template<typename T, Index Log2Dim>
1102 inline bool
1103 1994932138 LeafNode<T, Log2Dim>::probeValue(const Coord& xyz, ValueType& val) const
1104 {
1105 1994932138 return this->probeValue(LeafNode::coordToOffset(xyz), val);
1106 }
1107
1108 template<typename T, Index Log2Dim>
1109 inline bool
1110 1994932138 LeafNode<T, Log2Dim>::probeValue(Index offset, ValueType& val) const
1111 {
1112
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 999284969 times.
1994932138 assert(offset < SIZE);
1113 1994932138 val = mBuffer[offset];
1114 1994932138 return mValueMask.isOn(offset);
1115 }
1116
1117
1118 template<typename T, Index Log2Dim>
1119 inline void
1120 8801555 LeafNode<T, Log2Dim>::setValueOff(const Coord& xyz, const ValueType& val)
1121 {
1122 8801555 this->setValueOff(LeafNode::coordToOffset(xyz), val);
1123 8801555 }
1124
1125 template<typename T, Index Log2Dim>
1126 inline void
1127 11557955 LeafNode<T, Log2Dim>::setValueOff(Index offset, const ValueType& val)
1128 {
1129
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 5779029 times.
11557955 assert(offset < SIZE);
1130
13/26
✓ Branch 1 taken 1 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 1 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 1 times.
✗ Branch 8 not taken.
✓ Branch 10 taken 1 times.
✗ Branch 11 not taken.
✓ Branch 13 taken 1 times.
✗ Branch 14 not taken.
✓ Branch 16 taken 1 times.
✗ Branch 17 not taken.
✓ Branch 19 taken 1 times.
✗ Branch 20 not taken.
✓ Branch 22 taken 1 times.
✗ Branch 23 not taken.
✓ Branch 25 taken 1 times.
✗ Branch 26 not taken.
✓ Branch 28 taken 1 times.
✗ Branch 29 not taken.
✓ Branch 31 taken 1 times.
✗ Branch 32 not taken.
✓ Branch 34 taken 1 times.
✗ Branch 35 not taken.
✓ Branch 37 taken 1 times.
✗ Branch 38 not taken.
11557968 mBuffer.setValue(offset, val);
1131 11557968 mValueMask.setOff(offset);
1132 11557955 }
1133
1134
1135 template<typename T, Index Log2Dim>
1136 inline void
1137 354102 LeafNode<T, Log2Dim>::setActiveState(const Coord& xyz, bool on)
1138 {
1139
2/2
✓ Branch 0 taken 177059 times.
✓ Branch 1 taken 1 times.
354102 mValueMask.set(this->coordToOffset(xyz), on);
1140 354102 }
1141
1142
1143 template<typename T, Index Log2Dim>
1144 inline void
1145 16272169 LeafNode<T, Log2Dim>::setValueOnly(const Coord& xyz, const ValueType& val)
1146 {
1147 16272169 this->setValueOnly(LeafNode::coordToOffset(xyz), val);
1148 16272169 }
1149
1150 template<typename T, Index Log2Dim>
1151 inline void
1152 104273471 LeafNode<T, Log2Dim>::setValueOnly(Index offset, const ValueType& val)
1153 {
1154
56/113
✗ Branch 0 not taken.
✓ Branch 1 taken 52441495 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 1 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 1 times.
✗ Branch 8 not taken.
✓ Branch 10 taken 1 times.
✗ Branch 11 not taken.
✓ Branch 13 taken 1 times.
✗ Branch 14 not taken.
✓ Branch 16 taken 1 times.
✗ Branch 17 not taken.
✓ Branch 19 taken 1 times.
✗ Branch 20 not taken.
✓ Branch 22 taken 1 times.
✗ Branch 23 not taken.
✓ Branch 25 taken 1 times.
✗ Branch 26 not taken.
✓ Branch 28 taken 1 times.
✗ Branch 29 not taken.
✓ Branch 31 taken 1 times.
✗ Branch 32 not taken.
✓ Branch 34 taken 1 times.
✗ Branch 35 not taken.
✓ Branch 37 taken 1 times.
✗ Branch 38 not taken.
✓ Branch 40 taken 1 times.
✗ Branch 41 not taken.
✓ Branch 43 taken 1 times.
✗ Branch 44 not taken.
✓ Branch 46 taken 1 times.
✗ Branch 47 not taken.
✓ Branch 49 taken 1 times.
✗ Branch 50 not taken.
✓ Branch 52 taken 1 times.
✗ Branch 53 not taken.
✓ Branch 55 taken 1 times.
✗ Branch 56 not taken.
✓ Branch 58 taken 1 times.
✗ Branch 59 not taken.
✓ Branch 61 taken 1 times.
✗ Branch 62 not taken.
✓ Branch 64 taken 1 times.
✗ Branch 65 not taken.
✓ Branch 67 taken 1 times.
✗ Branch 68 not taken.
✓ Branch 70 taken 1 times.
✗ Branch 71 not taken.
✓ Branch 73 taken 1 times.
✗ Branch 74 not taken.
✓ Branch 76 taken 1 times.
✗ Branch 77 not taken.
✓ Branch 79 taken 1 times.
✗ Branch 80 not taken.
✓ Branch 82 taken 1 times.
✗ Branch 83 not taken.
✓ Branch 85 taken 1 times.
✗ Branch 86 not taken.
✓ Branch 88 taken 1 times.
✗ Branch 89 not taken.
✓ Branch 91 taken 1 times.
✗ Branch 92 not taken.
✓ Branch 94 taken 1 times.
✗ Branch 95 not taken.
✓ Branch 97 taken 1 times.
✗ Branch 98 not taken.
✓ Branch 100 taken 1 times.
✗ Branch 101 not taken.
✓ Branch 103 taken 1 times.
✗ Branch 104 not taken.
✓ Branch 106 taken 1 times.
✗ Branch 107 not taken.
✓ Branch 109 taken 1 times.
✗ Branch 110 not taken.
✓ Branch 112 taken 1 times.
✗ Branch 113 not taken.
✓ Branch 115 taken 1 times.
✗ Branch 116 not taken.
✓ Branch 118 taken 1 times.
✗ Branch 119 not taken.
✓ Branch 121 taken 1 times.
✗ Branch 122 not taken.
✓ Branch 124 taken 1 times.
✗ Branch 125 not taken.
✓ Branch 127 taken 1 times.
✗ Branch 128 not taken.
✓ Branch 130 taken 1 times.
✗ Branch 131 not taken.
✓ Branch 133 taken 1 times.
✗ Branch 134 not taken.
✓ Branch 136 taken 1 times.
✗ Branch 137 not taken.
✓ Branch 139 taken 1 times.
✗ Branch 140 not taken.
✓ Branch 142 taken 1 times.
✗ Branch 143 not taken.
✓ Branch 145 taken 1 times.
✗ Branch 146 not taken.
✓ Branch 148 taken 1 times.
✗ Branch 149 not taken.
✓ Branch 151 taken 1 times.
✗ Branch 152 not taken.
✓ Branch 154 taken 1 times.
✗ Branch 155 not taken.
✓ Branch 157 taken 1 times.
✗ Branch 158 not taken.
✓ Branch 160 taken 1 times.
✗ Branch 161 not taken.
✓ Branch 163 taken 1 times.
✗ Branch 164 not taken.
✓ Branch 166 taken 1 times.
✗ Branch 167 not taken.
104273527 assert(offset<SIZE); mBuffer.setValue(offset, val);
1155 104273527 }
1156
1157
1158 ////////////////////////////////////////
1159
1160
1161 template<typename T, Index Log2Dim>
1162 inline void
1163 8268 LeafNode<T, Log2Dim>::clip(const CoordBBox& clipBBox, const T& background)
1164 {
1165 8268 CoordBBox nodeBBox = this->getNodeBoundingBox();
1166 if (!clipBBox.hasOverlap(nodeBBox)) {
1167 // This node lies completely outside the clipping region. Fill it with the background.
1168 this->fill(background, /*active=*/false);
1169 } else if (clipBBox.isInside(nodeBBox)) {
1170 // This node lies completely inside the clipping region. Leave it intact.
1171 8222 return;
1172 }
1173
1174 // This node isn't completely contained inside the clipping region.
1175 // Set any voxels that lie outside the region to the background value.
1176
1177 // Construct a boolean mask that is on inside the clipping region and off outside it.
1178 NodeMaskType mask;
1179 46 nodeBBox.intersect(clipBBox);
1180 Coord xyz;
1181 int &x = xyz.x(), &y = xyz.y(), &z = xyz.z();
1182
2/2
✓ Branch 0 taken 82 times.
✓ Branch 1 taken 23 times.
210 for (x = nodeBBox.min().x(); x <= nodeBBox.max().x(); ++x) {
1183
2/2
✓ Branch 0 taken 368 times.
✓ Branch 1 taken 82 times.
900 for (y = nodeBBox.min().y(); y <= nodeBBox.max().y(); ++y) {
1184
2/2
✓ Branch 0 taken 1698 times.
✓ Branch 1 taken 368 times.
4132 for (z = nodeBBox.min().z(); z <= nodeBBox.max().z(); ++z) {
1185 3396 mask.setOn(static_cast<Index32>(this->coordToOffset(xyz)));
1186 }
1187 }
1188 }
1189
1190 // Set voxels that lie in the inactive region of the mask (i.e., outside
1191 // the clipping region) to the background value.
1192
2/2
✓ Branch 1 taken 10078 times.
✓ Branch 2 taken 23 times.
20248 for (MaskOffIterator maskIter = mask.beginOff(); maskIter; ++maskIter) {
1193 20156 this->setValueOff(maskIter.pos(), background);
1194 }
1195 }
1196
1197
1198 ////////////////////////////////////////
1199
1200
1201 template<typename T, Index Log2Dim>
1202 inline void
1203
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 169903 times.
339600 LeafNode<T, Log2Dim>::fill(const CoordBBox& bbox, const ValueType& value, bool active)
1204 {
1205
1/2
✓ Branch 0 taken 90 times.
✗ Branch 1 not taken.
180 if (!this->allocate()) return;
1206
1207 339600 auto clippedBBox = this->getNodeBoundingBox();
1208 339600 clippedBBox.intersect(bbox);
1209 if (!clippedBBox) return;
1210
1211
2/2
✓ Branch 0 taken 909285 times.
✓ Branch 1 taken 169993 times.
2157424 for (Int32 x = clippedBBox.min().x(); x <= clippedBBox.max().x(); ++x) {
1212 1817824 const Index offsetX = (x & (DIM-1u)) << 2*Log2Dim;
1213
2/2
✓ Branch 0 taken 4732938 times.
✓ Branch 1 taken 909285 times.
11280490 for (Int32 y = clippedBBox.min().y(); y <= clippedBBox.max().y(); ++y) {
1214 9462666 const Index offsetXY = offsetX + ((y & (DIM-1u)) << Log2Dim);
1215
2/2
✓ Branch 0 taken 20825678 times.
✓ Branch 1 taken 4732938 times.
51094840 for (Int32 z = clippedBBox.min().z(); z <= clippedBBox.max().z(); ++z) {
1216 41632174 const Index offset = offsetXY + (z & (DIM-1u));
1217 41631994 mBuffer[offset] = value;
1218
2/2
✓ Branch 0 taken 15176741 times.
✓ Branch 1 taken 5648937 times.
41632174 mValueMask.set(offset, active);
1219 }
1220 }
1221 }
1222 }
1223
1224 template<typename T, Index Log2Dim>
1225 inline void
1226 LeafNode<T, Log2Dim>::fill(const ValueType& value)
1227 {
1228 8693 mBuffer.fill(value);
1229 }
1230
1231 template<typename T, Index Log2Dim>
1232 inline void
1233 1119 LeafNode<T, Log2Dim>::fill(const ValueType& value, bool active)
1234 {
1235 1119 mBuffer.fill(value);
1236 mValueMask.set(active);
1237 1119 }
1238
1239
1240 ////////////////////////////////////////
1241
1242
1243 template<typename T, Index Log2Dim>
1244 template<typename DenseT>
1245 inline void
1246 5202 LeafNode<T, Log2Dim>::copyToDense(const CoordBBox& bbox, DenseT& dense) const
1247 {
1248 5202 mBuffer.loadValues();
1249
1250 using DenseValueType = typename DenseT::ValueType;
1251
1252 const size_t xStride = dense.xStride(), yStride = dense.yStride(), zStride = dense.zStride();
1253 const Coord& min = dense.bbox().min();
1254 5202 DenseValueType* t0 = dense.data() + zStride * (bbox.min()[2] - min[2]); // target array
1255 5202 const T* s0 = &mBuffer[bbox.min()[2] & (DIM-1u)]; // source array
1256
2/2
✓ Branch 0 taken 11323 times.
✓ Branch 1 taken 2601 times.
27848 for (Int32 x = bbox.min()[0], ex = bbox.max()[0] + 1; x < ex; ++x) {
1257 22646 DenseValueType* t1 = t0 + xStride * (x - min[0]);
1258 22646 const T* s1 = s0 + ((x & (DIM-1u)) << 2*Log2Dim);
1259
2/2
✓ Branch 0 taken 51262 times.
✓ Branch 1 taken 11323 times.
125170 for (Int32 y = bbox.min()[1], ey = bbox.max()[1] + 1; y < ey; ++y) {
1260 102524 DenseValueType* t2 = t1 + yStride * (y - min[1]);
1261 102524 const T* s2 = s1 + ((y & (DIM-1u)) << Log2Dim);
1262
2/2
✓ Branch 0 taken 203768 times.
✓ Branch 1 taken 51262 times.
510060 for (Int32 z = bbox.min()[2], ez = bbox.max()[2] + 1; z < ez; ++z, t2 += zStride) {
1263 407536 *t2 = DenseValueType(*s2++);
1264 }
1265 }
1266 }
1267 5202 }
1268
1269
1270 template<typename T, Index Log2Dim>
1271 template<typename DenseT>
1272 inline void
1273
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 1130 times.
2260 LeafNode<T, Log2Dim>::copyFromDense(const CoordBBox& bbox, const DenseT& dense,
1274 const ValueType& background, const ValueType& tolerance)
1275 {
1276 if (!this->allocate()) return;
1277
1278 using DenseValueType = typename DenseT::ValueType;
1279
1280 const size_t xStride = dense.xStride(), yStride = dense.yStride(), zStride = dense.zStride();
1281 const Coord& min = dense.bbox().min();
1282
1283 2260 const DenseValueType* s0 = dense.data() + zStride * (bbox.min()[2] - min[2]); // source
1284 2260 const Int32 n0 = bbox.min()[2] & (DIM-1u);
1285
2/2
✓ Branch 0 taken 7314 times.
✓ Branch 1 taken 1130 times.
16888 for (Int32 x = bbox.min()[0], ex = bbox.max()[0]+1; x < ex; ++x) {
1286 14628 const DenseValueType* s1 = s0 + xStride * (x - min[0]);
1287 14628 const Int32 n1 = n0 + ((x & (DIM-1u)) << 2*LOG2DIM);
1288
2/2
✓ Branch 0 taken 47458 times.
✓ Branch 1 taken 7314 times.
109544 for (Int32 y = bbox.min()[1], ey = bbox.max()[1]+1; y < ey; ++y) {
1289 94916 const DenseValueType* s2 = s1 + yStride * (y - min[1]);
1290 94916 Int32 n2 = n1 + ((y & (DIM-1u)) << LOG2DIM);
1291
2/2
✓ Branch 0 taken 313626 times.
✓ Branch 1 taken 47458 times.
722168 for (Int32 z = bbox.min()[2], ez = bbox.max()[2]+1; z < ez; ++z, ++n2, s2 += zStride) {
1292
2/2
✓ Branch 0 taken 105064 times.
✓ Branch 1 taken 208562 times.
627252 if (math::isApproxEqual(background, ValueType(*s2), tolerance)) {
1293 210128 mValueMask.setOff(n2);
1294 210128 mBuffer[n2] = background;
1295 } else {
1296 417124 mValueMask.setOn(n2);
1297 417124 mBuffer[n2] = ValueType(*s2);
1298 }
1299 }
1300 }
1301 }
1302 }
1303
1304
1305 ////////////////////////////////////////
1306
1307
1308 template<typename T, Index Log2Dim>
1309 inline void
1310 LeafNode<T, Log2Dim>::readTopology(std::istream& is, bool /*fromHalf*/)
1311 {
1312 mValueMask.load(is);
1313 }
1314
1315
1316 template<typename T, Index Log2Dim>
1317 inline void
1318 LeafNode<T, Log2Dim>::writeTopology(std::ostream& os, bool /*toHalf*/) const
1319 {
1320 mValueMask.save(os);
1321 }
1322
1323
1324 ////////////////////////////////////////
1325
1326
1327
1328 template<typename T, Index Log2Dim>
1329 inline void
1330 46112 LeafNode<T,Log2Dim>::skipCompressedValues(bool seekable, std::istream& is, bool fromHalf)
1331 {
1332
2/2
✓ Branch 0 taken 22873 times.
✓ Branch 1 taken 183 times.
46112 if (seekable) {
1333 // Seek over voxel values.
1334 45746 io::readCompressedValues<ValueType, NodeMaskType>(
1335 45746 is, nullptr, SIZE, mValueMask, fromHalf);
1336 } else {
1337 // Read and discard voxel values.
1338 732 Buffer temp;
1339
1/2
✓ Branch 1 taken 183 times.
✗ Branch 2 not taken.
366 io::readCompressedValues(is, temp.mData, SIZE, mValueMask, fromHalf);
1340 }
1341 46112 }
1342
1343
1344 template<typename T, Index Log2Dim>
1345 inline void
1346 LeafNode<T,Log2Dim>::readBuffers(std::istream& is, bool fromHalf)
1347 {
1348
12/25
✓ Branch 2 taken 10221 times.
✗ Branch 3 not taken.
✓ Branch 5 taken 10221 times.
✗ Branch 6 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 1 times.
✗ Branch 10 not taken.
✓ Branch 12 taken 1 times.
✗ Branch 13 not taken.
✓ Branch 16 taken 1 times.
✗ Branch 17 not taken.
✓ Branch 19 taken 1 times.
✗ Branch 20 not taken.
✓ Branch 23 taken 1 times.
✗ Branch 24 not taken.
✓ Branch 26 taken 1 times.
✗ Branch 27 not taken.
✓ Branch 30 taken 1 times.
✗ Branch 31 not taken.
✓ Branch 33 taken 1 times.
✗ Branch 34 not taken.
✓ Branch 37 taken 1 times.
✗ Branch 38 not taken.
✓ Branch 40 taken 1 times.
✗ Branch 41 not taken.
26872 this->readBuffers(is, CoordBBox::inf(), fromHalf);
1349 }
1350
1351
1352 template<typename T, Index Log2Dim>
1353 inline void
1354 54334 LeafNode<T,Log2Dim>::readBuffers(std::istream& is, const CoordBBox& clipBBox, bool fromHalf)
1355 {
1356 54334 SharedPtr<io::StreamMetadata> meta = io::getStreamMetadataPtr(is);
1357
5/6
✓ Branch 0 taken 23057 times.
✓ Branch 1 taken 4110 times.
✓ Branch 3 taken 23057 times.
✗ Branch 4 not taken.
✓ Branch 5 taken 184 times.
✓ Branch 6 taken 22873 times.
54334 const bool seekable = meta && meta->seekable();
1358
1359
1/2
✓ Branch 1 taken 27167 times.
✗ Branch 2 not taken.
54334 std::streamoff maskpos = is.tellg();
1360
1361
2/2
✓ Branch 0 taken 22873 times.
✓ Branch 1 taken 4294 times.
54334 if (seekable) {
1362 // Seek over the value mask.
1363 mValueMask.seek(is);
1364 } else {
1365 // Read in the value mask.
1366 mValueMask.load(is);
1367 }
1368
1369 54334 int8_t numBuffers = 1;
1370
2/4
✓ Branch 1 taken 27167 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✓ Branch 4 taken 27167 times.
54334 if (io::getFormatVersion(is) < OPENVDB_FILE_VERSION_NODE_MASK_COMPRESSION) {
1371 // Read in the origin.
1372 is.read(reinterpret_cast<char*>(&mOrigin), sizeof(Coord::ValueType) * 3);
1373
1374 // Read in the number of buffers, which should now always be one.
1375 is.read(reinterpret_cast<char*>(&numBuffers), sizeof(int8_t));
1376 }
1377
1378 54334 CoordBBox nodeBBox = this->getNodeBoundingBox();
1379 if (!clipBBox.hasOverlap(nodeBBox)) {
1380 // This node lies completely outside the clipping region.
1381
1/2
✓ Branch 1 taken 295 times.
✗ Branch 2 not taken.
590 skipCompressedValues(seekable, is, fromHalf);
1382 mValueMask.setOff();
1383 mBuffer.setOutOfCore(false);
1384 } else {
1385 // If this node lies completely inside the clipping region and it is being read
1386 // from a memory-mapped file, delay loading of its buffer until the buffer
1387 // is actually accessed. (If this node requires clipping, its buffer
1388 // must be accessed and therefore must be loaded.)
1389
1/2
✓ Branch 1 taken 26872 times.
✗ Branch 2 not taken.
53744 io::MappedFile::Ptr mappedFile = io::getMappedFilePtr(is);
1390
2/2
✓ Branch 0 taken 22761 times.
✓ Branch 1 taken 4111 times.
53744 const bool delayLoad = ((mappedFile.get() != nullptr) && clipBBox.isInside(nodeBBox));
1391
1392 if (delayLoad) {
1393 mBuffer.setOutOfCore(true);
1394
1/2
✓ Branch 1 taken 22761 times.
✗ Branch 2 not taken.
45522 mBuffer.mFileInfo = new typename Buffer::FileInfo;
1395 mBuffer.mFileInfo->meta = meta;
1396
1/2
✓ Branch 1 taken 22761 times.
✗ Branch 2 not taken.
45522 mBuffer.mFileInfo->bufpos = is.tellg();
1397 45522 mBuffer.mFileInfo->mapping = mappedFile;
1398 // Save the offset to the value mask, because the in-memory copy
1399 // might change before the value buffer gets read.
1400 45522 mBuffer.mFileInfo->maskpos = maskpos;
1401 // Skip over voxel values.
1402
1/2
✓ Branch 1 taken 22761 times.
✗ Branch 2 not taken.
45522 skipCompressedValues(seekable, is, fromHalf);
1403 } else {
1404
1/2
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
4 mBuffer.allocate();
1405
1/2
✓ Branch 1 taken 4111 times.
✗ Branch 2 not taken.
8222 io::readCompressedValues(is, mBuffer.mData, SIZE, mValueMask, fromHalf);
1406 mBuffer.setOutOfCore(false);
1407
1408 // Get this tree's background value.
1409 7194 T background = zeroVal<T>();
1410
3/4
✓ Branch 1 taken 4111 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 8 times.
✓ Branch 4 taken 4103 times.
8222 if (const void* bgPtr = io::getGridBackgroundValuePtr(is)) {
1411 16 background = *static_cast<const T*>(bgPtr);
1412 }
1413
1/2
✓ Branch 1 taken 4111 times.
✗ Branch 2 not taken.
8222 this->clip(clipBBox, background);
1414 }
1415 }
1416
1417
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 27167 times.
54334 if (numBuffers > 1) {
1418 // Read in and discard auxiliary buffers that were created with earlier
1419 // versions of the library. (Auxiliary buffers are not mask compressed.)
1420 const bool zipped = io::getDataCompression(is) & io::COMPRESS_ZIP;
1421 Buffer temp;
1422 for (int i = 1; i < numBuffers; ++i) {
1423 if (fromHalf) {
1424 io::HalfReader<io::RealToHalf<T>::isReal, T>::read(is, temp.mData, SIZE, zipped);
1425 } else {
1426 io::readData<T>(is, temp.mData, SIZE, zipped);
1427 }
1428 }
1429 }
1430
1431 // increment the leaf number
1432
4/6
✓ Branch 0 taken 23057 times.
✓ Branch 1 taken 4110 times.
✓ Branch 3 taken 23057 times.
✗ Branch 4 not taken.
✓ Branch 6 taken 23057 times.
✗ Branch 7 not taken.
54334 if (meta) meta->setLeaf(meta->leaf() + 1);
1433 54334 }
1434
1435
1436 template<typename T, Index Log2Dim>
1437 inline void
1438 80124 LeafNode<T, Log2Dim>::writeBuffers(std::ostream& os, bool toHalf) const
1439 {
1440 // Write out the value mask.
1441 mValueMask.save(os);
1442
1443 80124 mBuffer.loadValues();
1444
1445 80124 io::writeCompressedValues(os, mBuffer.mData, SIZE,
1446 80124 mValueMask, /*childMask=*/NodeMaskType(), toHalf);
1447 80124 }
1448
1449
1450 ////////////////////////////////////////
1451
1452
1453 template<typename T, Index Log2Dim>
1454 inline bool
1455
1/2
✓ Branch 0 taken 628 times.
✗ Branch 1 not taken.
628 LeafNode<T, Log2Dim>::operator==(const LeafNode& other) const
1456 {
1457
2/2
✓ Branch 0 taken 625 times.
✓ Branch 1 taken 2 times.
627 return mOrigin == other.mOrigin &&
1458
2/2
✓ Branch 0 taken 1 times.
✓ Branch 1 taken 624 times.
625 mValueMask == other.valueMask() &&
1459 625 mBuffer == other.mBuffer;
1460 }
1461
1462
1463 template<typename T, Index Log2Dim>
1464 inline Index64
1465
2/2
✓ Branch 0 taken 56914 times.
✓ Branch 1 taken 11607 times.
137042 LeafNode<T, Log2Dim>::memUsage() const
1466 {
1467 // Use sizeof(*this) to capture alignment-related padding
1468 // (but note that sizeof(*this) includes sizeof(mBuffer)).
1469 137042 return sizeof(*this) + mBuffer.memUsage() - sizeof(mBuffer);
1470 }
1471
1472
1473 template<typename T, Index Log2Dim>
1474 inline Index64
1475 LeafNode<T, Log2Dim>::memUsageIfLoaded() const
1476 {
1477 // Use sizeof(*this) to capture alignment-related padding
1478 // (but note that sizeof(*this) includes sizeof(mBuffer)).
1479 return sizeof(*this) + mBuffer.memUsageIfLoaded() - sizeof(mBuffer);
1480 }
1481
1482
1483 template<typename T, Index Log2Dim>
1484 inline void
1485 236725 LeafNode<T, Log2Dim>::evalActiveBoundingBox(CoordBBox& bbox, bool visitVoxels) const
1486 {
1487 236725 CoordBBox this_bbox = this->getNodeBoundingBox();
1488 194842 if (bbox.isInside(this_bbox)) return;//this LeafNode is already enclosed in the bbox
1489
2/2
✓ Branch 0 taken 21779 times.
✓ Branch 1 taken 156 times.
41883 if (ValueOnCIter iter = this->cbeginValueOn()) {//any active values?
1490
2/2
✓ Branch 0 taken 20774 times.
✓ Branch 1 taken 1005 times.
41571 if (visitVoxels) {//use voxel granularity?
1491 39561 this_bbox.reset();
1492
2/2
✓ Branch 0 taken 1739965 times.
✓ Branch 1 taken 20774 times.
3220402 for(; iter; ++iter) this_bbox.expand(this->offsetToLocalCoord(iter.pos()));
1493 this_bbox.translate(this->origin());
1494 }
1495 41571 bbox.expand(this_bbox);
1496 }
1497 }
1498
1499
1500 template<typename T, Index Log2Dim>
1501 template<typename OtherType, Index OtherLog2Dim>
1502 inline bool
1503 9474 LeafNode<T, Log2Dim>::hasSameTopology(const LeafNode<OtherType, OtherLog2Dim>* other) const
1504 {
1505
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 4373 times.
9474 assert(other);
1506 9474 return (Log2Dim == OtherLog2Dim && mValueMask == other->getValueMask());
1507 }
1508
1509 template<typename T, Index Log2Dim>
1510 inline bool
1511
2/2
✓ Branch 0 taken 6889 times.
✓ Branch 1 taken 16169 times.
44986 LeafNode<T, Log2Dim>::isConstant(ValueType& firstValue,
1512 bool& state,
1513 const ValueType& tolerance) const
1514 {
1515
2/2
✓ Branch 0 taken 15853 times.
✓ Branch 1 taken 2436 times.
35916 if (!mValueMask.isConstant(state)) return false;// early termination
1516 31280 firstValue = mBuffer[0];
1517
2/2
✓ Branch 0 taken 286242 times.
✓ Branch 1 taken 525 times.
402284 for (Index i = 1; i < SIZE; ++i) {
1518
2/3
✓ Branch 0 taken 249963 times.
✓ Branch 1 taken 5946 times.
✗ Branch 2 not taken.
340898 if ( !math::isApproxEqual(mBuffer[i], firstValue, tolerance) ) return false;// early termination
1519 }
1520 return true;
1521 }
1522
1523 template<typename T, Index Log2Dim>
1524 inline bool
1525
2/2
✓ Branch 0 taken 22370 times.
✓ Branch 1 taken 12632 times.
65298 LeafNode<T, Log2Dim>::isConstant(ValueType& minValue,
1526 ValueType& maxValue,
1527 bool& state,
1528 const ValueType& tolerance) const
1529 {
1530
2/2
✓ Branch 0 taken 16512 times.
✓ Branch 1 taken 5412 times.
39142 if (!mValueMask.isConstant(state)) return false;// early termination
1531 28318 minValue = maxValue = mBuffer[0];
1532
2/2
✓ Branch 0 taken 8002731 times.
✓ Branch 1 taken 15448 times.
13626886 for (Index i = 1; i < SIZE; ++i) {
1533 const T& v = mBuffer[i];
1534
3/3
✓ Branch 0 taken 472 times.
✓ Branch 1 taken 7218729 times.
✓ Branch 2 taken 783530 times.
13600696 if (v < minValue) {
1535
3/3
✓ Branch 0 taken 14 times.
✓ Branch 1 taken 458 times.
✓ Branch 2 taken 73 times.
1090 if ((maxValue - v) > tolerance) return false;// early termination
1536 28 minValue = v;
1537
3/3
✓ Branch 0 taken 473 times.
✓ Branch 1 taken 7218255 times.
✓ Branch 2 taken 783458 times.
13599606 } else if (v > maxValue) {
1538
3/3
✓ Branch 0 taken 12 times.
✓ Branch 1 taken 461 times.
✓ Branch 2 taken 72 times.
1090 if ((v - minValue) > tolerance) return false;// early termination
1539 24 maxValue = v;
1540 }
1541 }
1542 return true;
1543 }
1544
1545 template<typename T, Index Log2Dim>
1546 inline T
1547 15488 LeafNode<T, Log2Dim>::medianAll(T *tmp) const
1548 {
1549 15488 std::unique_ptr<T[]> data(nullptr);
1550
2/2
✓ Branch 0 taken 13 times.
✓ Branch 1 taken 4740 times.
15488 if (tmp == nullptr) {//allocate temporary storage
1551
1/4
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
13 data.reset(new T[NUM_VALUES]);
1552 tmp = data.get();
1553 }
1554
3/4
✓ Branch 1 taken 4753 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 13 times.
✓ Branch 4 taken 4740 times.
15488 if (tmp != mBuffer.data()) {
1555
1/2
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
13 const T* src = mBuffer.data();
1556
2/2
✓ Branch 0 taken 6656 times.
✓ Branch 1 taken 13 times.
6669 for (T* dst = tmp; dst-tmp < NUM_VALUES;) *dst++ = *src++;
1557 }
1558 static const size_t midpoint = (NUM_VALUES - 1) >> 1;
1559 15488 std::nth_element(tmp, tmp + midpoint, tmp + NUM_VALUES);
1560
2/2
✓ Branch 0 taken 13 times.
✓ Branch 1 taken 4740 times.
15501 return tmp[midpoint];
1561 }
1562
1563 template<typename T, Index Log2Dim>
1564 inline Index
1565 11 LeafNode<T, Log2Dim>::medianOn(T &value, T *tmp) const
1566 {
1567 11 const Index count = mValueMask.countOn();
1568
2/2
✓ Branch 0 taken 1 times.
✓ Branch 1 taken 10 times.
11 if (count == NUM_VALUES) {//special case: all voxels are active
1569 1 value = this->medianAll(tmp);
1570 1 return NUM_VALUES;
1571
2/2
✓ Branch 0 taken 9 times.
✓ Branch 1 taken 1 times.
10 } else if (count == 0) {
1572 return 0;
1573 }
1574 9 std::unique_ptr<T[]> data(nullptr);
1575
1/2
✓ Branch 0 taken 9 times.
✗ Branch 1 not taken.
9 if (tmp == nullptr) {//allocate temporary storage
1576
1/2
✓ Branch 1 taken 9 times.
✗ Branch 2 not taken.
9 data.reset(new T[count]);// 0 < count < NUM_VALUES
1577 tmp = data.get();
1578 }
1579
2/2
✓ Branch 0 taken 45 times.
✓ Branch 1 taken 9 times.
99 for (auto iter=this->cbeginValueOn(); iter; ++iter) *tmp++ = *iter;
1580 9 T *begin = tmp - count;
1581 9 const size_t midpoint = (count - 1) >> 1;
1582 9 std::nth_element(begin, begin + midpoint, tmp);
1583
1/2
✓ Branch 0 taken 9 times.
✗ Branch 1 not taken.
9 value = begin[midpoint];
1584 return count;
1585 }
1586
1587 template<typename T, Index Log2Dim>
1588 inline Index
1589 11 LeafNode<T, Log2Dim>::medianOff(T &value, T *tmp) const
1590 {
1591
2/2
✓ Branch 0 taken 1 times.
✓ Branch 1 taken 10 times.
11 const Index count = mValueMask.countOff();
1592
2/2
✓ Branch 0 taken 1 times.
✓ Branch 1 taken 10 times.
11 if (count == NUM_VALUES) {//special case: all voxels are inactive
1593 1 value = this->medianAll(tmp);
1594 1 return NUM_VALUES;
1595
2/2
✓ Branch 0 taken 9 times.
✓ Branch 1 taken 1 times.
10 } else if (count == 0) {
1596 return 0;
1597 }
1598 9 std::unique_ptr<T[]> data(nullptr);
1599
1/2
✓ Branch 0 taken 9 times.
✗ Branch 1 not taken.
9 if (tmp == nullptr) {//allocate temporary storage
1600
1/2
✓ Branch 1 taken 9 times.
✗ Branch 2 not taken.
9 data.reset(new T[count]);// 0 < count < NUM_VALUES
1601 tmp = data.get();
1602 }
1603
2/2
✓ Branch 0 taken 4563 times.
✓ Branch 1 taken 9 times.
9135 for (auto iter=this->cbeginValueOff(); iter; ++iter) *tmp++ = *iter;
1604 9 T *begin = tmp - count;
1605 9 const size_t midpoint = (count - 1) >> 1;
1606 9 std::nth_element(begin, begin + midpoint, tmp);
1607
1/2
✓ Branch 0 taken 9 times.
✗ Branch 1 not taken.
9 value = begin[midpoint];
1608 return count;
1609 }
1610
1611 ////////////////////////////////////////
1612
1613
1614 template<typename T, Index Log2Dim>
1615 inline void
1616 4 LeafNode<T, Log2Dim>::addTile(Index /*level*/, const Coord& xyz, const ValueType& val, bool active)
1617 {
1618 4 this->addTile(this->coordToOffset(xyz), val, active);
1619 }
1620
1621 template<typename T, Index Log2Dim>
1622 inline void
1623 4 LeafNode<T, Log2Dim>::addTile(Index offset, const ValueType& val, bool active)
1624 {
1625
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
4 assert(offset < SIZE);
1626 4 setValueOnly(offset, val);
1627 4 setActiveState(offset, active);
1628 }
1629
1630 template<typename T, Index Log2Dim>
1631 template<typename AccessorT>
1632 inline void
1633 LeafNode<T, Log2Dim>::addTileAndCache(Index level, const Coord& xyz,
1634 const ValueType& val, bool active, AccessorT&)
1635 {
1636 this->addTile(level, xyz, val, active);
1637 }
1638
1639
1640 ////////////////////////////////////////
1641
1642
1643 template<typename T, Index Log2Dim>
1644 inline void
1645
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 20866 times.
21710 LeafNode<T, Log2Dim>::resetBackground(const ValueType& oldBackground,
1646 const ValueType& newBackground)
1647 {
1648 if (!this->allocate()) return;
1649
1650 typename NodeMaskType::OffIterator iter;
1651 // For all inactive values...
1652
2/2
✓ Branch 1 taken 6829447 times.
✓ Branch 2 taken 20866 times.
7255439 for (iter = this->mValueMask.beginOff(); iter; ++iter) {
1653 7212019 ValueType &inactiveValue = mBuffer[iter.pos()];
1654
2/2
✓ Branch 0 taken 5754789 times.
✓ Branch 1 taken 1074658 times.
7212019 if (math::isApproxEqual(inactiveValue, oldBackground)) {
1655 6137361 inactiveValue = newBackground;
1656
1/2
✓ Branch 0 taken 1074658 times.
✗ Branch 1 not taken.
1074658 } else if (math::isApproxEqual(inactiveValue, math::negative(oldBackground))) {
1657 1074658 inactiveValue = math::negative(newBackground);
1658 }
1659 }
1660 }
1661
1662
1663 template<typename T, Index Log2Dim>
1664 template<MergePolicy Policy>
1665 inline void
1666
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 553 times.
1106 LeafNode<T, Log2Dim>::merge(const LeafNode& other)
1667 {
1668 if (!this->allocate()) return;
1669
1670 OPENVDB_NO_UNREACHABLE_CODE_WARNING_BEGIN
1671 if (Policy == MERGE_NODES) return;
1672 1106 typename NodeMaskType::OnIterator iter = other.valueMask().beginOn();
1673
2/2
✓ Branch 0 taken 74022 times.
✓ Branch 1 taken 553 times.
149150 for (; iter; ++iter) {
1674 const Index n = iter.pos();
1675
2/2
✓ Branch 1 taken 68581 times.
✓ Branch 2 taken 5441 times.
148044 if (mValueMask.isOff(n)) {
1676 137162 mBuffer[n] = other.mBuffer[n];
1677 137162 mValueMask.setOn(n);
1678 }
1679 }
1680 OPENVDB_NO_UNREACHABLE_CODE_WARNING_END
1681 }
1682
1683 template<typename T, Index Log2Dim>
1684 template<MergePolicy Policy>
1685 inline void
1686 LeafNode<T, Log2Dim>::merge(const LeafNode& other,
1687 const ValueType& /*bg*/, const ValueType& /*otherBG*/)
1688 {
1689 553 this->template merge<Policy>(other);
1690 553 }
1691
1692 template<typename T, Index Log2Dim>
1693 template<MergePolicy Policy>
1694 inline void
1695
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 1 times.
1 LeafNode<T, Log2Dim>::merge(const ValueType& tileValue, bool tileActive)
1696 {
1697 if (!this->allocate()) return;
1698
1699 OPENVDB_NO_UNREACHABLE_CODE_WARNING_BEGIN
1700 if (Policy != MERGE_ACTIVE_STATES_AND_NODES) return;
1701
1/2
✓ Branch 0 taken 1 times.
✗ Branch 1 not taken.
1 if (!tileActive) return;
1702 // Replace all inactive values with the active tile value.
1703
2/2
✓ Branch 1 taken 512 times.
✓ Branch 2 taken 1 times.
514 for (typename NodeMaskType::OffIterator iter = mValueMask.beginOff(); iter; ++iter) {
1704 const Index n = iter.pos();
1705 512 mBuffer[n] = tileValue;
1706 512 mValueMask.setOn(n);
1707 }
1708 OPENVDB_NO_UNREACHABLE_CODE_WARNING_END
1709 }
1710
1711
1712 template<typename T, Index Log2Dim>
1713 template<typename OtherType>
1714 inline void
1715 LeafNode<T, Log2Dim>::topologyUnion(const LeafNode<OtherType, Log2Dim>& other, bool)
1716 {
1717 mValueMask |= other.valueMask();
1718 }
1719
1720 template<typename T, Index Log2Dim>
1721 template<typename OtherType>
1722 inline void
1723 LeafNode<T, Log2Dim>::topologyIntersection(const LeafNode<OtherType, Log2Dim>& other,
1724 const ValueType&)
1725 {
1726 mValueMask &= other.valueMask();
1727 }
1728
1729 template<typename T, Index Log2Dim>
1730 template<typename OtherType>
1731 inline void
1732 11 LeafNode<T, Log2Dim>::topologyDifference(const LeafNode<OtherType, Log2Dim>& other,
1733 const ValueType&)
1734 {
1735 11 mValueMask &= !other.valueMask();
1736 11 }
1737
1738 template<typename T, Index Log2Dim>
1739 inline void
1740
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 334 times.
334 LeafNode<T, Log2Dim>::negate()
1741 {
1742 if (!this->allocate()) return;
1743
1744
2/2
✓ Branch 0 taken 171008 times.
✓ Branch 1 taken 334 times.
171342 for (Index i = 0; i < SIZE; ++i) {
1745 171008 mBuffer[i] = -mBuffer[i];
1746 }
1747 }
1748
1749
1750 ////////////////////////////////////////
1751
1752
1753 template<typename T, Index Log2Dim>
1754 template<typename CombineOp>
1755 inline void
1756
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 42 times.
73 LeafNode<T, Log2Dim>::combine(const LeafNode& other, CombineOp& op)
1757 {
1758 if (!this->allocate()) return;
1759
1760 CombineArgs<T> args;
1761
2/2
✓ Branch 0 taken 21504 times.
✓ Branch 1 taken 42 times.
37449 for (Index i = 0; i < SIZE; ++i) {
1762 70656 op(args.setARef(mBuffer[i])
1763 .setAIsActive(mValueMask.isOn(i))
1764 74752 .setBRef(other.mBuffer[i])
1765 .setBIsActive(other.valueMask().isOn(i))
1766 37376 .setResultRef(mBuffer[i]));
1767 mValueMask.set(i, args.resultIsActive());
1768 }
1769 }
1770
1771
1772 template<typename T, Index Log2Dim>
1773 template<typename CombineOp>
1774 inline void
1775
0/2
✗ Branch 0 not taken.
✗ Branch 1 not taken.
673 LeafNode<T, Log2Dim>::combine(const ValueType& value, bool valueIsActive, CombineOp& op)
1776 {
1777 if (!this->allocate()) return;
1778
1779 CombineArgs<T> args;
1780 args.setBRef(value).setBIsActive(valueIsActive);
1781
0/2
✗ Branch 0 not taken.
✗ Branch 1 not taken.
345249 for (Index i = 0; i < SIZE; ++i) {
1782 689152 op(args.setARef(mBuffer[i])
1783 .setAIsActive(mValueMask.isOn(i))
1784 344576 .setResultRef(mBuffer[i]));
1785 mValueMask.set(i, args.resultIsActive());
1786 }
1787 }
1788
1789
1790 ////////////////////////////////////////
1791
1792
1793 template<typename T, Index Log2Dim>
1794 template<typename CombineOp, typename OtherType>
1795 inline void
1796 LeafNode<T, Log2Dim>::combine2(const LeafNode& other, const OtherType& value,
1797 bool valueIsActive, CombineOp& op)
1798 {
1799 if (!this->allocate()) return;
1800
1801 CombineArgs<T, OtherType> args;
1802 args.setBRef(value).setBIsActive(valueIsActive);
1803 for (Index i = 0; i < SIZE; ++i) {
1804 op(args.setARef(other.mBuffer[i])
1805 .setAIsActive(other.valueMask().isOn(i))
1806 .setResultRef(mBuffer[i]));
1807 mValueMask.set(i, args.resultIsActive());
1808 }
1809 }
1810
1811
1812 template<typename T, Index Log2Dim>
1813 template<typename CombineOp, typename OtherNodeT>
1814 inline void
1815 LeafNode<T, Log2Dim>::combine2(const ValueType& value, const OtherNodeT& other,
1816 bool valueIsActive, CombineOp& op)
1817 {
1818 if (!this->allocate()) return;
1819
1820 CombineArgs<T, typename OtherNodeT::ValueType> args;
1821 args.setARef(value).setAIsActive(valueIsActive);
1822 for (Index i = 0; i < SIZE; ++i) {
1823 op(args.setBRef(other.mBuffer[i])
1824 .setBIsActive(other.valueMask().isOn(i))
1825 .setResultRef(mBuffer[i]));
1826 mValueMask.set(i, args.resultIsActive());
1827 }
1828 }
1829
1830
1831 template<typename T, Index Log2Dim>
1832 template<typename CombineOp, typename OtherNodeT>
1833 inline void
1834 4 LeafNode<T, Log2Dim>::combine2(const LeafNode& b0, const OtherNodeT& b1, CombineOp& op)
1835 {
1836 if (!this->allocate()) return;
1837
1838 CombineArgs<T, typename OtherNodeT::ValueType> args;
1839 2052 for (Index i = 0; i < SIZE; ++i) {
1840 2048 mValueMask.set(i, b0.valueMask().isOn(i) || b1.valueMask().isOn(i));
1841 2048 op(args.setARef(b0.mBuffer[i])
1842 .setAIsActive(b0.valueMask().isOn(i))
1843 3584 .setBRef(b1.mBuffer[i])
1844 512 .setBIsActive(b1.valueMask().isOn(i))
1845 4096 .setResultRef(mBuffer[i]));
1846 mValueMask.set(i, args.resultIsActive());
1847 }
1848 }
1849
1850
1851 ////////////////////////////////////////
1852
1853
1854 template<typename T, Index Log2Dim>
1855 template<typename BBoxOp>
1856 inline void
1857 2 LeafNode<T, Log2Dim>::visitActiveBBox(BBoxOp& op) const
1858 {
1859 if (op.template descent<LEVEL>()) {
1860 for (ValueOnCIter i=this->cbeginValueOn(); i; ++i) {
1861 op.template operator()<LEVEL>(CoordBBox::createCube(i.getCoord(), 1));
1862 }
1863 } else {
1864 2 op.template operator()<LEVEL>(this->getNodeBoundingBox());
1865 }
1866 2 }
1867
1868
1869 template<typename T, Index Log2Dim>
1870 template<typename VisitorOp>
1871 inline void
1872 LeafNode<T, Log2Dim>::visit(VisitorOp& op)
1873 {
1874
4/8
✓ Branch 1 taken 8 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 8 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 8 times.
✗ Branch 8 not taken.
✓ Branch 10 taken 8 times.
✗ Branch 11 not taken.
32 doVisit<LeafNode, VisitorOp, ChildAllIter>(*this, op);
1875 32 }
1876
1877
1878 template<typename T, Index Log2Dim>
1879 template<typename VisitorOp>
1880 inline void
1881 LeafNode<T, Log2Dim>::visit(VisitorOp& op) const
1882 {
1883
4/8
✓ Branch 1 taken 8 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 8 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 8 times.
✗ Branch 8 not taken.
✓ Branch 10 taken 8 times.
✗ Branch 11 not taken.
32 doVisit<const LeafNode, VisitorOp, ChildAllCIter>(*this, op);
1884 32 }
1885
1886
1887 template<typename T, Index Log2Dim>
1888 template<typename NodeT, typename VisitorOp, typename ChildAllIterT>
1889 inline void
1890 128 LeafNode<T, Log2Dim>::doVisit(NodeT& self, VisitorOp& op)
1891 {
1892
2/2
✓ Branch 1 taken 32768 times.
✓ Branch 2 taken 64 times.
65792 for (ChildAllIterT iter = self.beginChildAll(); iter; ++iter) {
1893 65536 op(iter);
1894 }
1895 128 }
1896
1897
1898 ////////////////////////////////////////
1899
1900
1901 template<typename T, Index Log2Dim>
1902 template<typename OtherLeafNodeType, typename VisitorOp>
1903 inline void
1904 LeafNode<T, Log2Dim>::visit2Node(OtherLeafNodeType& other, VisitorOp& op)
1905 {
1906 doVisit2Node<LeafNode, OtherLeafNodeType, VisitorOp, ChildAllIter,
1907 24 typename OtherLeafNodeType::ChildAllIter>(*this, other, op);
1908 24 }
1909
1910
1911 template<typename T, Index Log2Dim>
1912 template<typename OtherLeafNodeType, typename VisitorOp>
1913 inline void
1914 LeafNode<T, Log2Dim>::visit2Node(OtherLeafNodeType& other, VisitorOp& op) const
1915 {
1916 doVisit2Node<const LeafNode, OtherLeafNodeType, VisitorOp, ChildAllCIter,
1917 typename OtherLeafNodeType::ChildAllCIter>(*this, other, op);
1918 }
1919
1920
1921 template<typename T, Index Log2Dim>
1922 template<
1923 typename NodeT,
1924 typename OtherNodeT,
1925 typename VisitorOp,
1926 typename ChildAllIterT,
1927 typename OtherChildAllIterT>
1928 inline void
1929 48 LeafNode<T, Log2Dim>::doVisit2Node(NodeT& self, OtherNodeT& other, VisitorOp& op)
1930 {
1931 // Allow the two nodes to have different ValueTypes, but not different dimensions.
1932 static_assert(OtherNodeT::SIZE == NodeT::SIZE,
1933 "can't visit nodes of different sizes simultaneously");
1934 static_assert(OtherNodeT::LEVEL == NodeT::LEVEL,
1935 "can't visit nodes at different tree levels simultaneously");
1936
1937 48 ChildAllIterT iter = self.beginChildAll();
1938 48 OtherChildAllIterT otherIter = other.beginChildAll();
1939
1940
3/4
✓ Branch 0 taken 12288 times.
✓ Branch 1 taken 24 times.
✓ Branch 2 taken 12288 times.
✗ Branch 3 not taken.
49200 for ( ; iter && otherIter; ++iter, ++otherIter) {
1941 24576 op(iter, otherIter);
1942 }
1943 48 }
1944
1945
1946 ////////////////////////////////////////
1947
1948
1949 template<typename T, Index Log2Dim>
1950 template<typename IterT, typename VisitorOp>
1951 inline void
1952 LeafNode<T, Log2Dim>::visit2(IterT& otherIter, VisitorOp& op, bool otherIsLHS)
1953 {
1954 19 doVisit2<LeafNode, VisitorOp, ChildAllIter, IterT>(
1955 *this, otherIter, op, otherIsLHS);
1956 19 }
1957
1958
1959 template<typename T, Index Log2Dim>
1960 template<typename IterT, typename VisitorOp>
1961 inline void
1962 LeafNode<T, Log2Dim>::visit2(IterT& otherIter, VisitorOp& op, bool otherIsLHS) const
1963 {
1964 doVisit2<const LeafNode, VisitorOp, ChildAllCIter, IterT>(
1965 *this, otherIter, op, otherIsLHS);
1966 }
1967
1968
1969 template<typename T, Index Log2Dim>
1970 template<
1971 typename NodeT,
1972 typename VisitorOp,
1973 typename ChildAllIterT,
1974 typename OtherChildAllIterT>
1975 inline void
1976 19 LeafNode<T, Log2Dim>::doVisit2(NodeT& self, OtherChildAllIterT& otherIter,
1977 VisitorOp& op, bool otherIsLHS)
1978 {
1979 19 if (!otherIter) return;
1980
1981 19 if (otherIsLHS) {
1982 5140 for (ChildAllIterT iter = self.beginChildAll(); iter; ++iter) {
1983 5120 op(otherIter, iter);
1984 }
1985 } else {
1986 4626 for (ChildAllIterT iter = self.beginChildAll(); iter; ++iter) {
1987 4608 op(iter, otherIter);
1988 }
1989 }
1990 }
1991
1992
1993 ////////////////////////////////////////
1994
1995
1996 template<typename T, Index Log2Dim>
1997 inline std::ostream&
1998 operator<<(std::ostream& os, const typename LeafNode<T, Log2Dim>::Buffer& buf)
1999 {
2000 for (Index32 i = 0, N = buf.size(); i < N; ++i) os << buf.mData[i] << ", ";
2001 return os;
2002 }
2003
2004 } // namespace tree
2005 } // namespace OPENVDB_VERSION_NAME
2006 } // namespace openvdb
2007
2008
2009 ////////////////////////////////////////
2010
2011
2012 // Specialization for LeafNodes of type bool
2013 #include "LeafNodeBool.h"
2014
2015 // Specialization for LeafNodes with mask information only
2016 #include "LeafNodeMask.h"
2017
2018 #endif // OPENVDB_TREE_LEAFNODE_HAS_BEEN_INCLUDED
2019