level_set_function.hpp

Go to the documentation of this file.

#pragma once
 
#include "sdf_calculator.hpp"
 
namespace netdem {
 
class LevelSetFunction {
protected:
  Vec3d corner{-0.5, -0.5, -0.5};
 
  double spacing{0.05};
 
  Vec3i dim{21, 21, 21};
 
  VecXT<VecXT<VecXT<double>>> signed_distance_table;
 
public:
  LevelSetFunction();
 
  // return the lower corner of the box
  Vec3d const &GetCorner() const;
 
  // return the spacing of the level set grid
  double GetSpacing() const;
 
  // return the dimesions of the level set grid
  Vec3i const &GetDimensions() const;
 
  double GetSignedDistance(int id_x, int id_y, int id_z);
 
  // return SDF of the nodes of the level set grid
  VecXT<VecXT<VecXT<double>>> const &GetSignedDistanceTable() const;
 
  void SetCorner(double corner_x, double corner_y, double corner_z);
 
  void SetSpacing(double sp);
 
  void SetDimension(double dim_x, double dim_y, double dim_z);
 
  void SetSignedDistance(int id_x, int id_y, int id_z, double sdf);
 
  void InitFromSDFCalculator(const SDFCalculator &sdf_calculator);
 
  double SignedDistance(Vec3d const &pos) const;
 
  Vec3d GradientInterpolate(Vec3d const &pos);
 
  Vec3d GradientMinus(int i, int j, int k);
 
  Vec3d GradientPlus(int i, int j, int k);
 
  void Reinitialization(int iter, double dt);
 
  void Reinitialization();
};
 
} // namespace netdem