04_demo_normal_patches.cpp

This is an example of how to use the netdem library.

#include "data_dumper.hpp"
#include "gen_pack.hpp"
#include "gen_wall_box_plane.hpp"
#include "gen_wall_box_plate.hpp"
#include "gravity.hpp"
#include "model_linear_spring.hpp"
#include "shape_polybezier.hpp"
#include "shape_sphere.hpp"
#include "shape_trimesh.hpp"
#include "simulation.hpp"
#include <filesystem>
#include <iostream>
#include <unordered_map>
 
using namespace netdem;
using namespace std;
 
STLModel GetKernelPolyhedron();
 
void DemoNormalPatches(int argc, char **argv) {
  int num_particles = atof(argv[2]);
 
  filesystem::create_directories("local/polybezier_model/demo_normal_patches");
 
  for (int i = 0; i < num_particles; i++) {
    auto kernel_stl = GetKernelPolyhedron();
    kernel_stl.SaveAsSTL(
        "local/polybezier_model/demo_normal_patches/kernel_stl_" +
        to_string(i) + ".stl");
 
    Polybezier polybezier;
    polybezier.InitFromKernelSTL(kernel_stl);
    polybezier.SaveNormalPatchesCubic(
        "local/polybezier_model/demo_normal_patches/polybezier_normal_" +
        to_string(i) + ".stl");
 
    TriMesh trimesh;
    trimesh.use_linked_patches = true;
    trimesh.InitFromSTL(polybezier.GetSTLModel(500));
    trimesh.SaveAsSTL("local/polybezier_model/demo_normal_patches/mesh_stl_" +
                      to_string(i) + ".stl");
    trimesh.SaveNormalPatchesCubic(
        "local/polybezier_model/demo_normal_patches/mesh_normal_" +
        to_string(i) + ".stl");
  }
}