Initial commit

BFSSurfaceRecon.xcodeproj/project.pbxproj

+// !$*UTF8*$!
+{
+	archiveVersion = 1;
+	classes = {
+	};
+	objectVersion = 50;
+	objects = {
+
+/* Begin PBXBuildFile section */
+		E1EE8CCB261E2D29000FE218 /* main.cpp in Sources */ = {isa = PBXBuildFile; fileRef = E1EE8CCA261E2D29000FE218 /* main.cpp */; };
+		E1EE8CD4261E2D54000FE218 /* GEL.framework in Frameworks */ = {isa = PBXBuildFile; fileRef = E1EE8CD3261E2D54000FE218 /* GEL.framework */; };
+/* End PBXBuildFile section */
+
+/* Begin PBXCopyFilesBuildPhase section */
+		E1EE8CC5261E2D29000FE218 /* CopyFiles */ = {
+			isa = PBXCopyFilesBuildPhase;
+			buildActionMask = 2147483647;
+			dstPath = /usr/share/man/man1/;
+			dstSubfolderSpec = 0;
+			files = (
+			);
+			runOnlyForDeploymentPostprocessing = 1;
+		};
+/* End PBXCopyFilesBuildPhase section */
+
+/* Begin PBXFileReference section */
+		E1EE8CC7261E2D29000FE218 /* BFSSurfaceRecon */ = {isa = PBXFileReference; explicitFileType = "compiled.mach-o.executable"; includeInIndex = 0; path = BFSSurfaceRecon; sourceTree = BUILT_PRODUCTS_DIR; };
+		E1EE8CCA261E2D29000FE218 /* main.cpp */ = {isa = PBXFileReference; lastKnownFileType = sourcecode.cpp.cpp; path = main.cpp; sourceTree = "<group>"; };
+		E1EE8CD3261E2D54000FE218 /* GEL.framework */ = {isa = PBXFileReference; lastKnownFileType = wrapper.framework; name = GEL.framework; path = ../../../../../../../Library/Frameworks/GEL.framework; sourceTree = "<group>"; };
+/* End PBXFileReference section */
+
+/* Begin PBXFrameworksBuildPhase section */
+		E1EE8CC4261E2D29000FE218 /* Frameworks */ = {
+			isa = PBXFrameworksBuildPhase;
+			buildActionMask = 2147483647;
+			files = (
+				E1EE8CD4261E2D54000FE218 /* GEL.framework in Frameworks */,
+			);
+			runOnlyForDeploymentPostprocessing = 0;
+		};
+/* End PBXFrameworksBuildPhase section */
+
+/* Begin PBXGroup section */
+		E1EE8CBE261E2D29000FE218 = {
+			isa = PBXGroup;
+			children = (
+				E1EE8CC9261E2D29000FE218 /* BFSSurfaceRecon */,
+				E1EE8CC8261E2D29000FE218 /* Products */,
+				E1EE8CD2261E2D54000FE218 /* Frameworks */,
+			);
+			sourceTree = "<group>";
+		};
+		E1EE8CC8261E2D29000FE218 /* Products */ = {
+			isa = PBXGroup;
+			children = (
+				E1EE8CC7261E2D29000FE218 /* BFSSurfaceRecon */,
+			);
+			name = Products;
+			sourceTree = "<group>";
+		};
+		E1EE8CC9261E2D29000FE218 /* BFSSurfaceRecon */ = {
+			isa = PBXGroup;
+			children = (
+				E1EE8CCA261E2D29000FE218 /* main.cpp */,
+			);
+			path = BFSSurfaceRecon;
+			sourceTree = "<group>";
+		};
+		E1EE8CD2261E2D54000FE218 /* Frameworks */ = {
+			isa = PBXGroup;
+			children = (
+				E1EE8CD3261E2D54000FE218 /* GEL.framework */,
+			);
+			name = Frameworks;
+			sourceTree = "<group>";
+		};
+/* End PBXGroup section */
+
+/* Begin PBXNativeTarget section */
+		E1EE8CC6261E2D29000FE218 /* BFSSurfaceRecon */ = {
+			isa = PBXNativeTarget;
+			buildConfigurationList = E1EE8CCE261E2D29000FE218 /* Build configuration list for PBXNativeTarget "BFSSurfaceRecon" */;
+			buildPhases = (
+				E1EE8CC3261E2D29000FE218 /* Sources */,
+				E1EE8CC4261E2D29000FE218 /* Frameworks */,
+				E1EE8CC5261E2D29000FE218 /* CopyFiles */,
+			);
+			buildRules = (
+			);
+			dependencies = (
+			);
+			name = BFSSurfaceRecon;
+			productName = BFSSurfaceRecon;
+			productReference = E1EE8CC7261E2D29000FE218 /* BFSSurfaceRecon */;
+			productType = "com.apple.product-type.tool";
+		};
+/* End PBXNativeTarget section */
+
+/* Begin PBXProject section */
+		E1EE8CBF261E2D29000FE218 /* Project object */ = {
+			isa = PBXProject;
+			attributes = {
+				LastUpgradeCheck = 1240;
+				TargetAttributes = {
+					E1EE8CC6261E2D29000FE218 = {
+						CreatedOnToolsVersion = 12.4;
+					};
+				};
+			};
+			buildConfigurationList = E1EE8CC2261E2D29000FE218 /* Build configuration list for PBXProject "BFSSurfaceRecon" */;
+			compatibilityVersion = "Xcode 9.3";
+			developmentRegion = en;
+			hasScannedForEncodings = 0;
+			knownRegions = (
+				en,
+				Base,
+			);
+			mainGroup = E1EE8CBE261E2D29000FE218;
+			productRefGroup = E1EE8CC8261E2D29000FE218 /* Products */;
+			projectDirPath = "";
+			projectRoot = "";
+			targets = (
+				E1EE8CC6261E2D29000FE218 /* BFSSurfaceRecon */,
+			);
+		};
+/* End PBXProject section */
+
+/* Begin PBXSourcesBuildPhase section */
+		E1EE8CC3261E2D29000FE218 /* Sources */ = {
+			isa = PBXSourcesBuildPhase;
+			buildActionMask = 2147483647;
+			files = (
+				E1EE8CCB261E2D29000FE218 /* main.cpp in Sources */,
+			);
+			runOnlyForDeploymentPostprocessing = 0;
+		};
+/* End PBXSourcesBuildPhase section */
+
+/* Begin XCBuildConfiguration section */
+		E1EE8CCC261E2D29000FE218 /* Debug */ = {
+			isa = XCBuildConfiguration;
+			buildSettings = {
+				ALWAYS_SEARCH_USER_PATHS = NO;
+				CLANG_ANALYZER_NONNULL = YES;
+				CLANG_ANALYZER_NUMBER_OBJECT_CONVERSION = YES_AGGRESSIVE;
+				CLANG_CXX_LANGUAGE_STANDARD = "gnu++14";
+				CLANG_CXX_LIBRARY = "libc++";
+				CLANG_ENABLE_MODULES = YES;
+				CLANG_ENABLE_OBJC_ARC = YES;
+				CLANG_ENABLE_OBJC_WEAK = YES;
+				CLANG_WARN_BLOCK_CAPTURE_AUTORELEASING = YES;
+				CLANG_WARN_BOOL_CONVERSION = YES;
+				CLANG_WARN_COMMA = YES;
+				CLANG_WARN_CONSTANT_CONVERSION = YES;
+				CLANG_WARN_DEPRECATED_OBJC_IMPLEMENTATIONS = YES;
+				CLANG_WARN_DIRECT_OBJC_ISA_USAGE = YES_ERROR;
+				CLANG_WARN_DOCUMENTATION_COMMENTS = YES;
+				CLANG_WARN_EMPTY_BODY = YES;
+				CLANG_WARN_ENUM_CONVERSION = YES;
+				CLANG_WARN_INFINITE_RECURSION = YES;
+				CLANG_WARN_INT_CONVERSION = YES;
+				CLANG_WARN_NON_LITERAL_NULL_CONVERSION = YES;
+				CLANG_WARN_OBJC_IMPLICIT_RETAIN_SELF = YES;
+				CLANG_WARN_OBJC_LITERAL_CONVERSION = YES;
+				CLANG_WARN_OBJC_ROOT_CLASS = YES_ERROR;
+				CLANG_WARN_QUOTED_INCLUDE_IN_FRAMEWORK_HEADER = YES;
+				CLANG_WARN_RANGE_LOOP_ANALYSIS = YES;
+				CLANG_WARN_STRICT_PROTOTYPES = YES;
+				CLANG_WARN_SUSPICIOUS_MOVE = YES;
+				CLANG_WARN_UNGUARDED_AVAILABILITY = YES_AGGRESSIVE;
+				CLANG_WARN_UNREACHABLE_CODE = YES;
+				CLANG_WARN__DUPLICATE_METHOD_MATCH = YES;
+				COPY_PHASE_STRIP = NO;
+				DEBUG_INFORMATION_FORMAT = dwarf;
+				ENABLE_STRICT_OBJC_MSGSEND = YES;
+				ENABLE_TESTABILITY = YES;
+				GCC_C_LANGUAGE_STANDARD = gnu11;
+				GCC_DYNAMIC_NO_PIC = NO;
+				GCC_NO_COMMON_BLOCKS = YES;
+				GCC_OPTIMIZATION_LEVEL = 0;
+				GCC_PREPROCESSOR_DEFINITIONS = (
+					"DEBUG=1",
+					"$(inherited)",
+				);
+				GCC_WARN_64_TO_32_BIT_CONVERSION = YES;
+				GCC_WARN_ABOUT_RETURN_TYPE = YES_ERROR;
+				GCC_WARN_UNDECLARED_SELECTOR = YES;
+				GCC_WARN_UNINITIALIZED_AUTOS = YES_AGGRESSIVE;
+				GCC_WARN_UNUSED_FUNCTION = YES;
+				GCC_WARN_UNUSED_VARIABLE = YES;
+				MACOSX_DEPLOYMENT_TARGET = 11.1;
+				MTL_ENABLE_DEBUG_INFO = INCLUDE_SOURCE;
+				MTL_FAST_MATH = YES;
+				ONLY_ACTIVE_ARCH = YES;
+				SDKROOT = macosx;
+			};
+			name = Debug;
+		};
+		E1EE8CCD261E2D29000FE218 /* Release */ = {
+			isa = XCBuildConfiguration;
+			buildSettings = {
+				ALWAYS_SEARCH_USER_PATHS = NO;
+				CLANG_ANALYZER_NONNULL = YES;
+				CLANG_ANALYZER_NUMBER_OBJECT_CONVERSION = YES_AGGRESSIVE;
+				CLANG_CXX_LANGUAGE_STANDARD = "gnu++14";
+				CLANG_CXX_LIBRARY = "libc++";
+				CLANG_ENABLE_MODULES = YES;
+				CLANG_ENABLE_OBJC_ARC = YES;
+				CLANG_ENABLE_OBJC_WEAK = YES;
+				CLANG_WARN_BLOCK_CAPTURE_AUTORELEASING = YES;
+				CLANG_WARN_BOOL_CONVERSION = YES;
+				CLANG_WARN_COMMA = YES;
+				CLANG_WARN_CONSTANT_CONVERSION = YES;
+				CLANG_WARN_DEPRECATED_OBJC_IMPLEMENTATIONS = YES;
+				CLANG_WARN_DIRECT_OBJC_ISA_USAGE = YES_ERROR;
+				CLANG_WARN_DOCUMENTATION_COMMENTS = YES;
+				CLANG_WARN_EMPTY_BODY = YES;
+				CLANG_WARN_ENUM_CONVERSION = YES;
+				CLANG_WARN_INFINITE_RECURSION = YES;
+				CLANG_WARN_INT_CONVERSION = YES;
+				CLANG_WARN_NON_LITERAL_NULL_CONVERSION = YES;
+				CLANG_WARN_OBJC_IMPLICIT_RETAIN_SELF = YES;
+				CLANG_WARN_OBJC_LITERAL_CONVERSION = YES;
+				CLANG_WARN_OBJC_ROOT_CLASS = YES_ERROR;
+				CLANG_WARN_QUOTED_INCLUDE_IN_FRAMEWORK_HEADER = YES;
+				CLANG_WARN_RANGE_LOOP_ANALYSIS = YES;
+				CLANG_WARN_STRICT_PROTOTYPES = YES;
+				CLANG_WARN_SUSPICIOUS_MOVE = YES;
+				CLANG_WARN_UNGUARDED_AVAILABILITY = YES_AGGRESSIVE;
+				CLANG_WARN_UNREACHABLE_CODE = YES;
+				CLANG_WARN__DUPLICATE_METHOD_MATCH = YES;
+				COPY_PHASE_STRIP = NO;
+				DEBUG_INFORMATION_FORMAT = "dwarf-with-dsym";
+				ENABLE_NS_ASSERTIONS = NO;
+				ENABLE_STRICT_OBJC_MSGSEND = YES;
+				GCC_C_LANGUAGE_STANDARD = gnu11;
+				GCC_NO_COMMON_BLOCKS = YES;
+				GCC_WARN_64_TO_32_BIT_CONVERSION = YES;
+				GCC_WARN_ABOUT_RETURN_TYPE = YES_ERROR;
+				GCC_WARN_UNDECLARED_SELECTOR = YES;
+				GCC_WARN_UNINITIALIZED_AUTOS = YES_AGGRESSIVE;
+				GCC_WARN_UNUSED_FUNCTION = YES;
+				GCC_WARN_UNUSED_VARIABLE = YES;
+				MACOSX_DEPLOYMENT_TARGET = 11.1;
+				MTL_ENABLE_DEBUG_INFO = NO;
+				MTL_FAST_MATH = YES;
+				ONLY_ACTIVE_ARCH = YES;
+				SDKROOT = macosx;
+			};
+			name = Release;
+		};
+		E1EE8CCF261E2D29000FE218 /* Debug */ = {
+			isa = XCBuildConfiguration;
+			buildSettings = {
+				CODE_SIGN_STYLE = Automatic;
+				FRAMEWORK_SEARCH_PATHS = (
+					"$(inherited)",
+					"$(LOCAL_LIBRARY_DIR)/Frameworks",
+				);
+				PRODUCT_NAME = "$(TARGET_NAME)";
+			};
+			name = Debug;
+		};
+		E1EE8CD0261E2D29000FE218 /* Release */ = {
+			isa = XCBuildConfiguration;
+			buildSettings = {
+				CODE_SIGN_STYLE = Automatic;
+				FRAMEWORK_SEARCH_PATHS = (
+					"$(inherited)",
+					"$(LOCAL_LIBRARY_DIR)/Frameworks",
+				);
+				PRODUCT_NAME = "$(TARGET_NAME)";
+			};
+			name = Release;
+		};
+/* End XCBuildConfiguration section */
+
+/* Begin XCConfigurationList section */
+		E1EE8CC2261E2D29000FE218 /* Build configuration list for PBXProject "BFSSurfaceRecon" */ = {
+			isa = XCConfigurationList;
+			buildConfigurations = (
+				E1EE8CCC261E2D29000FE218 /* Debug */,
+				E1EE8CCD261E2D29000FE218 /* Release */,
+			);
+			defaultConfigurationIsVisible = 0;
+			defaultConfigurationName = Release;
+		};
+		E1EE8CCE261E2D29000FE218 /* Build configuration list for PBXNativeTarget "BFSSurfaceRecon" */ = {
+			isa = XCConfigurationList;
+			buildConfigurations = (
+				E1EE8CCF261E2D29000FE218 /* Debug */,
+				E1EE8CD0261E2D29000FE218 /* Release */,
+			);
+			defaultConfigurationIsVisible = 0;
+			defaultConfigurationName = Release;
+		};
+/* End XCConfigurationList section */
+	};
+	rootObject = E1EE8CBF261E2D29000FE218 /* Project object */;
+}

BFSSurfaceRecon.xcodeproj/project.xcworkspace/contents.xcworkspacedata

+<?xml version="1.0" encoding="UTF-8"?>
+<Workspace
+   version = "1.0">
+   <FileRef
+      location = "self:">
+   </FileRef>
+</Workspace>

BFSSurfaceRecon.xcodeproj/project.xcworkspace/xcshareddata/IDEWorkspaceChecks.plist

+<?xml version="1.0" encoding="UTF-8"?>
+<!DOCTYPE plist PUBLIC "-//Apple//DTD PLIST 1.0//EN" "http://www.apple.com/DTDs/PropertyList-1.0.dtd">
+<plist version="1.0">
+<dict>
+	<key>IDEDidComputeMac32BitWarning</key>
+	<true/>
+</dict>
+</plist>

BFSSurfaceRecon.xcodeproj/xcshareddata/xcschemes/BFSSurfaceRecon.xcscheme

+<?xml version="1.0" encoding="UTF-8"?>
+<Scheme
+   LastUpgradeVersion = "1240"
+   version = "1.3">
+   <BuildAction
+      parallelizeBuildables = "YES"
+      buildImplicitDependencies = "YES">
+      <BuildActionEntries>
+         <BuildActionEntry
+            buildForTesting = "YES"
+            buildForRunning = "YES"
+            buildForProfiling = "YES"
+            buildForArchiving = "YES"
+            buildForAnalyzing = "YES">
+            <BuildableReference
+               BuildableIdentifier = "primary"
+               BlueprintIdentifier = "E1EE8CC6261E2D29000FE218"
+               BuildableName = "BFSSurfaceRecon"
+               BlueprintName = "BFSSurfaceRecon"
+               ReferencedContainer = "container:BFSSurfaceRecon.xcodeproj">
+            </BuildableReference>
+         </BuildActionEntry>
+      </BuildActionEntries>
+   </BuildAction>
+   <TestAction
+      buildConfiguration = "Debug"
+      selectedDebuggerIdentifier = "Xcode.DebuggerFoundation.Debugger.LLDB"
+      selectedLauncherIdentifier = "Xcode.DebuggerFoundation.Launcher.LLDB"
+      shouldUseLaunchSchemeArgsEnv = "YES">
+      <Testables>
+      </Testables>
+   </TestAction>
+   <LaunchAction
+      buildConfiguration = "Release"
+      selectedDebuggerIdentifier = "Xcode.DebuggerFoundation.Debugger.LLDB"
+      selectedLauncherIdentifier = "Xcode.DebuggerFoundation.Launcher.LLDB"
+      launchStyle = "0"
+      useCustomWorkingDirectory = "NO"
+      ignoresPersistentStateOnLaunch = "NO"
+      debugDocumentVersioning = "YES"
+      debugServiceExtension = "internal"
+      allowLocationSimulation = "YES">
+      <BuildableProductRunnable
+         runnableDebuggingMode = "0">
+         <BuildableReference
+            BuildableIdentifier = "primary"
+            BlueprintIdentifier = "E1EE8CC6261E2D29000FE218"
+            BuildableName = "BFSSurfaceRecon"
+            BlueprintName = "BFSSurfaceRecon"
+            ReferencedContainer = "container:BFSSurfaceRecon.xcodeproj">
+         </BuildableReference>
+      </BuildableProductRunnable>
+   </LaunchAction>
+   <ProfileAction
+      buildConfiguration = "Release"
+      shouldUseLaunchSchemeArgsEnv = "YES"
+      savedToolIdentifier = ""
+      useCustomWorkingDirectory = "NO"
+      debugDocumentVersioning = "YES">
+      <BuildableProductRunnable
+         runnableDebuggingMode = "0">
+         <BuildableReference
+            BuildableIdentifier = "primary"
+            BlueprintIdentifier = "E1EE8CC6261E2D29000FE218"
+            BuildableName = "BFSSurfaceRecon"
+            BlueprintName = "BFSSurfaceRecon"
+            ReferencedContainer = "container:BFSSurfaceRecon.xcodeproj">
+         </BuildableReference>
+      </BuildableProductRunnable>
+   </ProfileAction>
+   <AnalyzeAction
+      buildConfiguration = "Debug">
+   </AnalyzeAction>
+   <ArchiveAction
+      buildConfiguration = "Release"
+      revealArchiveInOrganizer = "YES">
+   </ArchiveAction>
+</Scheme>

BFSSurfaceRecon.xcodeproj/xcuserdata/andreas.xcuserdatad/xcschemes/xcschememanagement.plist

+<?xml version="1.0" encoding="UTF-8"?>
+<!DOCTYPE plist PUBLIC "-//Apple//DTD PLIST 1.0//EN" "http://www.apple.com/DTDs/PropertyList-1.0.dtd">
+<plist version="1.0">
+<dict>
+	<key>SchemeUserState</key>
+	<dict>
+		<key>BFSSurfaceRecon.xcscheme_^#shared#^_</key>
+		<dict>
+			<key>orderHint</key>
+			<integer>0</integer>
+		</dict>
+	</dict>
+	<key>SuppressBuildableAutocreation</key>
+	<dict>
+		<key>E1EE8CC6261E2D29000FE218</key>
+		<dict>
+			<key>primary</key>
+			<true/>
+		</dict>
+	</dict>
+</dict>
+</plist>

BFSSurfaceRecon/main.cpp

+//
+//  main.cpp
+//  BFSSurfaceRecon
+//
+//  Created by Andreas Bærentzen on 07/04/2021.
+//
+
+#include <future>
+#include <thread>
+#include <iostream>
+#include <fstream>
+#include <sstream>
+#include <GEL/Util/Timer.h>
+#include <GEL/CGLA/CGLA.h>
+#include <GEL/Geometry/rply.h>
+#include <GEL/Geometry/KDTree.h>
+#include <GEL/Geometry/RGrid.h>
+#include <GEL/HMesh/polygonize.h>
+#include <GEL/HMesh/HMesh.h>
+#include <GEL/HMesh/obj_save.h>
+#include <GEL/HMesh/off_save.h>
+#include <GEL/HMesh/x3d_save.h>
+#include <GEL/HMesh/obj_save.h>
+#include <GEL/Util/ArgExtracter.h>
+#include <GEL/Geometry/GridAlgorithm.h>
+
+using namespace std;
+using namespace CGLA;
+using namespace Util;
+using namespace HMesh;
+using namespace Geometry;
+
+
+void save(Manifold& m, const string& file_name) {
+    const string extension = file_name.substr(file_name.length()-4,file_name.length());
+        if(extension==".obj"){
+            obj_save(file_name, m);
+            return;
+        }
+        else if(extension==".off"){
+            off_save(file_name, m);
+            return;
+        }
+        else if(extension==".x3d"){
+            x3d_save(file_name, m);
+            return;
+        }
+        else if(extension==".bhm") {
+            Serialization ser(file_name, std::ios_base::out);
+            m.serialize(ser);
+            return;
+        }
+        cout << "unknown format" << endl;
+}
+
+int main(int argc, char** argv) {
+    int longest_dim = 256;
+    string out_name = "out.bhm";
+    ArgExtracter ae(argc, argv);
+    ae.extract("-d", longest_dim);
+    ae.extract("-o", out_name);
+    string fn = ae.get_last_arg();
+
+    fstream f(fn);
+    string line_buffer;
+    vector<Vec3d> P;
+    vector<Vec3d> N;
+    cout << "Loading file ...." << endl;
+    while(getline(f, line_buffer)) {
+        Vec3d p, n;
+        stringstream sstr(line_buffer);
+        if (sstr >> p) {
+            sstr >> n;
+            if (sstr.good() || sstr.eof()){
+                P.push_back(p);
+                N.push_back(n);
+            }
+        }
+    }
+    
+    Vec3d pmin(1e32), pmax(-1e32);
+    for(auto p: P) {
+        pmin = v_min(pmin,p);
+        pmax = v_max(pmax,p);
+    }
+    
+    Vec3d diagonal = pmax-pmin;
+    double long_side = diagonal.max_coord();
+    diagonal /= long_side;
+    Vec3i dim = Vec3i(double(longest_dim)*diagonal);
+    cout << "Bounding box: " << pmin << " " << pmax << endl;
+    cout << "Volume dimensions: " << dim << endl;
+
+    XForm xform(pmin, pmax, dim);
+    RGrid<float> grid(dim,0);
+    RGrid<float> wgrid(dim,0);
+    cout << "Summing basis functions ..." << endl;
+    for(int i=0;i<P.size();++i) {
+        const Vec3d& p = P[i];
+        const Vec3d& n = N[i];
+        Vec3d p_grid = xform.apply(p);
+        Vec3i pi(p_grid + Vec3d(0.5));
+        Vec3i roi_min = v_max(Vec3i(0),pi-Vec3i(2));
+        Vec3i roi_max = v_min(dim-Vec3i(1), pi+Vec3i(3));
+
+        for(auto pin: Range3D(roi_min, roi_max)) {
+            Vec3d pn(pin);
+            Vec3d v = pn - p_grid;
+            float dist = dot(n, v);
+            float w = exp(-4.0*sqr_length(v));
+            grid[pin] += w * dist;
+            wgrid[pin] += w;
+        }
+    }
+    cout << "Normalizing..." << endl;
+    for(auto pi: Range3D(dim))
+        grid[pi] /= wgrid[pi];
+    
+    cout << "Polygonizing..." << endl;
+    Manifold m;
+    volume_polygonize(xform, grid, m, 0.0, true, false);
+    
+    cout << "Saving to " << out_name << endl;
+    
+    Serialization ser(out_name, std::ios_base::out);
+    m.serialize(ser);
+}