Added visualization for contains

samples/04-object-properties/main.cc

@@ -22,9 +22,18 @@ void test_obj(Obj const& obj)
             size = 100.f + length(obj);
 
         auto uniformPts = std::vector<tg::pos3>();
-        for (tg::u64 i = 0; i < size; ++i)
+        for (auto i = 0; i < size; ++i)
             uniformPts.emplace_back(uniform(rng, obj));
 
+        auto const aabb = aabb_of(obj);
+        tg::aabb3 aabb3;
+        if constexpr (domainD == 3)
+            aabb3 = aabb;
+        else
+            aabb3 = tg::aabb3(tg::pos3(aabb.min), tg::pos3(aabb.max));
+
+        auto g = gv::grid();
+        {
             auto v = gv::view();
 
             // uniform
@@ -41,16 +50,33 @@ void test_obj(Obj const& obj)
                 gv::view(gv::points(tg::pos3(centroid)).point_size_px(12.f), tg::color4::green);
 
             // aabb_of
-        auto const aabb = aabb_of(obj);
-        if constexpr (domainD == 3)
-            gv::view(gv::lines(aabb));
-        else
-            gv::view(gv::lines(tg::aabb3(tg::pos3(aabb.min), tg::pos3(aabb.max))));
+            gv::view(gv::lines(aabb3));
 
             // normal_of
             if constexpr (domainD == 3 && ((objectD == 2 && !tg::object_traits<Obj>::is_boundary) || (objectD == 1 && tg::object_traits<Obj>::is_boundary)))
                 gv::view(tg::segment3(centroid, centroid + normal_of(obj)), tg::color4::green);
         }
+        {
+            auto const aabbBigger = tg::aabb(aabb.min - 0.25f, aabb.max + 0.25f);
+            auto pointsInside = std::vector<tg::pos3>();
+            auto pointsOutside = std::vector<tg::pos3>();
+            auto const numPts = tg::pow(20.f, domainD);
+            auto const tolerance = objectD < domainD ? 0.1f : 0.f;
+            for (auto i = 0; i < numPts; ++i)
+            {
+                auto const p = uniform(rng, aabbBigger);
+                if (contains(obj, p, tolerance))
+                    pointsInside.emplace_back(p);
+                else
+                    pointsOutside.emplace_back(p);
+            }
+
+            // contains
+            auto v = gv::view();
+            view(pointsInside, gv::maybe_empty, tg::color4::green);
+            gv::view(gv::points(pointsOutside).point_size_px(5.f));
+        }
+    }
 }
 
 int main()