05_ave_num_linked_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;
 
void AveNumLinkedPatches(int argc, char **argv) {
  int num_particles = atof(argv[2]);
 
  filesystem::create_directories("local/polybezier_model/demo_normal_patches");
 
  for (int i = 1; i < 10; i++) {
    STLModel kernel_stl;
    kernel_stl.InitFromSTL(
        "local/polybezier_model/demo_normal_patches/kernel_stl_1.stl");
 
    Polybezier polybezier;
    polybezier.num_cells = i;
    polybezier.InitFromKernelSTL(kernel_stl);
 
    double ave_num{0};
    for (auto cell : polybezier.linked_patches_list) {
      ave_num += cell.size();
    }
    cout << i << ", " << ave_num / polybezier.linked_patches_list.size()
         << ", ";
 
    TriMesh trimesh;
    trimesh.num_cells = i;
    trimesh.use_linked_patches = true;
    trimesh.InitFromSTL(polybezier.GetSTLModel(500));
 
    ave_num = 0;
    for (auto cell : trimesh.linked_vertices) {
      ave_num += cell.size();
    }
    cout << ave_num / trimesh.linked_vertices.size() << endl;
  }
}