Remove old free surface meshing code

Author: dav-og

demos/sphere/demo_sphere_irreg_waves.cpp

@@ -126,34 +126,6 @@ int main(int argc, char* argv[]) {
     TestHydro hydro_forces(bodies, h5fname);
     hydro_forces.AddWaves(my_hydro_inputs);
 
-    // set up free surface from a mesh
-    auto fse_plane = chrono_types::make_shared<ChBody>();
-    fse_plane->SetPos(ChVector3d(0, 0, 0));
-    fse_plane->SetFixed(true);
-    fse_plane->EnableCollision(false);
-    system.AddBody(fse_plane);
-
-    my_hydro_inputs->SetUpWaveMesh();
-    std::shared_ptr<ChBody> fse_mesh = chrono_types::make_shared<ChBodyEasyMesh>(  //
-        my_hydro_inputs->GetMeshFile(),                                            // file name
-        1000,                                                                      // density
-        false,                                                                     // do not evaluate mass automatically
-        true,                                                                      // create visualization asset
-        false                                                                      // do not collide
-    );
-    fse_mesh->SetMass(1.0);
-    fse_mesh->SetPosDt(my_hydro_inputs->GetWaveMeshVelocity());
-    system.Add(fse_mesh);
-    auto fse_prismatic = chrono_types::make_shared<ChLinkLockPrismatic>();
-    fse_prismatic->Initialize(fse_plane, fse_mesh, ChFramed(ChVector3d(1.0, 0.0, 0.0), QuatFromAngleY(CH_PI_2)));
-    system.AddLink(fse_prismatic);
-
-    // Create a visualization material
-    auto fse_texture = chrono_types::make_shared<ChVisualMaterial>();
-    fse_texture->SetDiffuseColor(ChColor(0.026f, 0.084f, 0.168f));
-    fse_texture->SetOpacity(0.1f);
-    fse_mesh->GetVisualShape(0)->SetMaterial(0, fse_texture);
-
     // for profiling
     auto start = std::chrono::high_resolution_clock::now();
     // main simulation loop

include/hydroc/waves/irregular_wave.h

@@ -44,9 +44,6 @@ class IrregularWaves : public WaveBase {
     WaveMode GetWaveMode() override { return mode_; }
 
     Eigen::VectorXd SetSpectrumFrequencies(double start, double end, int num_steps);
-    void SetUpWaveMesh(std::string filename = "fse_mesh.obj");
-    std::string GetMeshFile();
-    Eigen::Vector3<double> GetWaveMeshVelocity();
 
     void AddH5Data(std::vector<HydroData::IrregularWaveInfo>& irreg_h5_data, HydroData::SimulationParameters& sim_data);
 
@@ -73,7 +70,6 @@ class IrregularWaves : public WaveBase {
     Eigen::VectorXd spectral_widths_;
     Eigen::VectorXd wavenumbers_;
     Eigen::VectorXd wave_phases_;
-    std::string mesh_file_name_;
 
     void InitializeIRFVectors();
     void ReadEtaFromFile();

src/hydro/waves/irregular_wave.cpp

@@ -226,28 +226,6 @@ Eigen::VectorXd IrregularWaves::SetSpectrumFrequencies(double start, double end,
     return result;
 }
 
-void IrregularWaves::SetUpWaveMesh(std::string filename) {
-    mesh_file_name_   = filename;
-    int num_timesteps = static_cast<int>(std::ceil(params_.simulation_duration_ / params_.simulation_dt_));
-    Eigen::VectorXd time_index = Eigen::VectorXd::LinSpaced(num_timesteps + 1, 0, num_timesteps * params_.simulation_dt_);
-    Eigen::VectorXd eta_vec    = Eigen::Map<const Eigen::VectorXd>(free_surface_elevation_sampled_.data(),
-                                                                   static_cast<Eigen::Index>(free_surface_elevation_sampled_.size()));
-
-    std::vector<std::array<double, 3>> free_surface_3d_pts = CreateFreeSurface3DPts(eta_vec, time_index);
-    std::vector<std::array<size_t, 3>> free_surface_triangles =
-        CreateFreeSurfaceTriangles(static_cast<size_t>(time_index.size()));
-
-    WriteFreeSurfaceMeshObj(free_surface_3d_pts, free_surface_triangles, mesh_file_name_);
-}
-
-std::string IrregularWaves::GetMeshFile() {
-    return mesh_file_name_;
-}
-
-Eigen::Vector3<double> IrregularWaves::GetWaveMeshVelocity() {
-    return Eigen::Vector3d(1.0, 0, 0);
-}
-
 void IrregularWaves::CreateSpectrum() {
     int nf;
     if (params_.nfrequencies_ == 0) {

src/hydro/waves/wave_utilities.cpp

@@ -353,67 +353,6 @@ Eigen::VectorXd JONSWAPSpectrumHz(Eigen::VectorXd& f,
     return spectral_densities;
 }
 
-std::vector<std::array<double, 3>> CreateFreeSurface3DPts(const Eigen::VectorXd& eta,
-                                                          const Eigen::VectorXd& t_vec) {
-    std::vector<std::array<double, 3>> surface(t_vec.size() * 2);
-
-    for (int i = 0; i < t_vec.size(); ++i) {
-        double t = -1 * t_vec[i];
-        double z = eta[i];
-
-        surface[2 * i]     = {t, -10.0, z};
-        surface[2 * i + 1] = {t, 10.0, z};
-    }
-
-    return surface;
-}
-
-std::vector<std::array<size_t, 3>> CreateFreeSurfaceTriangles(size_t eta_size) {
-    std::vector<std::array<size_t, 3>> triangles;
-
-    for (size_t i = 0; i < eta_size / 2 - 1; ++i) {
-        triangles.push_back({2 * i, 2 * i + 1, 2 * i + 3});
-        triangles.push_back({2 * i, 2 * i + 3, 2 * i + 2});
-    }
-
-    return triangles;
-}
-
-void WriteFreeSurfaceMeshObj(const std::vector<std::array<double, 3>>& points,
-                             const std::vector<std::array<size_t, 3>>& triangles,
-                             const std::string& file_name) {
-    std::ofstream out(file_name);
-    if (!out) {
-        std::cerr << "Failed to open " << file_name << std::endl;
-        return;
-    }
-
-    auto t  = std::time(nullptr);
-    auto tm = *std::localtime(&t);
-    out << "# Wavefront OBJ file exported by HydroChrono" << std::endl;
-    out << "# File Created: " << std::put_time(&tm, "%Y-%m-%d %H:%M:%S") << std::endl << std::endl;
-
-    out << "# Vertices: " << points.size() << std::endl << std::endl;
-    out << std::fixed << std::setprecision(6);
-    for (const auto& point : points) {
-        out << "v ";
-        out << std::setw(14) << point[0] << ' ';
-        out << std::setw(14) << point[1] << ' ';
-        out << std::setw(14) << point[2] << std::endl;
-    }
-    out << std::endl;
-
-    out << "# Faces: " << triangles.size() << std::endl << std::endl;
-    for (const auto& triangle : triangles) {
-        out << "f ";
-        out << std::setw(9) << triangle[0] + 1;
-        out << std::setw(9) << triangle[1] + 1;
-        out << std::setw(9) << triangle[2] + 1 << std::endl;
-    }
-
-    out.close();
-}
-
 Eigen::VectorXd GetWidthArray(const Eigen::VectorXd& input_array) {
     Eigen::VectorXd width_array(input_array.size());
     for (int ii = 0; ii < width_array.size(); ii++) {

src/hydro/waves/wave_utilities.h

@@ -149,15 +149,6 @@ Eigen::VectorXd JONSWAPSpectrumHz(Eigen::VectorXd& f,
                                   double gamma       = 3.3,
                                   bool is_normalized = false);
 
-std::vector<std::array<double, 3>> CreateFreeSurface3DPts(const Eigen::VectorXd& eta,
-                                                          const Eigen::VectorXd& t_vec);
-
-std::vector<std::array<size_t, 3>> CreateFreeSurfaceTriangles(size_t eta_size);
-
-void WriteFreeSurfaceMeshObj(const std::vector<std::array<double, 3>>& points,
-                             const std::vector<std::array<size_t, 3>>& triangles,
-                             const std::string& file_name);
-
 Eigen::VectorXd GetWidthArray(const Eigen::VectorXd& input_array);
 
 #endif  // HYDRO_WAVES_WAVE_UTILITIES_H