main.cpp File Reference

#include "geom/geomIncludes.h"
#include "inp/deckIncludes.h"
#include "inp/input.h"
#include "material/materialUtil.h"
#include "particle/baseParticle.h"
#include "particle/particleMpi.h"
#include "periDEMModel.h"
#include "postprocess/postprocess.h"
#include "util/function.h"
#include "util/io.h"
#include "util/parallelUtil.h"
#include "util/point.h"
#include <mpi.h>
#include <algorithm>
#include <cmath>
#include <cstdint>
#include <filesystem>
#include <format>
#include <fstream>
#include <stdexcept>
#include <string>
#include <thread>
#include <vector>

Include dependency graph for main.cpp:

Go to the source code of this file.

Data Structures
class	anonymous_namespace{main.cpp}::ContactWitness

Namespaces
namespace	anonymous_namespace{main.cpp}

Functions
std::string	anonymous_namespace{main.cpp}::directoryPathWithTrailingSep (const std::filesystem::path &dir)
bool	anonymous_namespace{main.cpp}::ownsNode (const data::ModelData &data, size_t i)
void	anonymous_namespace{main.cpp}::dumpFinalNodal (const data::ModelData &data, const std::filesystem::path &out_dir)
json	anonymous_namespace{main.cpp}::buildInputJson (const std::string &output_path, const std::filesystem::path &mesh_cir, const std::filesystem::path &mesh_wall, const std::string &mpi_strategy, double final_time, size_t num_steps)
int	main (int argc, char *argv[])

Function Documentation

main()

int main	(	int	argc,
		char *	argv[]
	)

Definition at line 345 of file main.cpp.

                                 {
  util::parallel::initMpi(argc, argv);
  util::io::InputParser input(argc, argv);
 
  unsigned n_threads = 1;
  if (input.cmdOptionExists("-nThreads"))
    n_threads = static_cast<unsigned>(std::stoi(input.getCmdOption("-nThreads")));
  util::parallel::initNThreads(n_threads);
 
  std::string mpi_strategy = "none";
  if (input.cmdOptionExists("-mpiStrategy"))
    mpi_strategy = input.getCmdOption("-mpiStrategy");
 
  namespace fs = std::filesystem;
  fs::path out_dir = fs::current_path() / "out";
  fs::path inp_dir = fs::current_path() / "inp";
  if (input.cmdOptionExists("-outputDir"))
    out_dir = input.getCmdOption("-outputDir");
  if (input.cmdOptionExists("-inputDir"))
    inp_dir = input.getCmdOption("-inputDir");
  fs::create_directories(out_dir);
  fs::create_directories(inp_dir);
 
  // Short stable window with soft contact.
  const double final_time = 2.0e-4;
  const size_t num_steps = 1000;
 
  auto input_json =
      buildInputJson(directoryPathWithTrailingSep(out_dir),
                     inp_dir / "mesh_cir.msh", inp_dir / "mesh_wall.msh",
                     mpi_strategy, final_time, num_steps);
  {
    std::ofstream os(inp_dir / "input.json");
    os << input_json.dump(2);
  }
 
  auto deck = std::make_shared<inp::Input>(input_json);
  PeriDEMModel dem(deck);
  auto witness = std::make_unique<ContactWitness>();
  ContactWitness *witness_p = witness.get();
  dem.setPostprocess(std::move(witness));
  dem.init();
  dem.integrate();
  witness_p->assertOk();
  dumpFinalNodal(dem, out_dir);
  dem.close();
 
  if (util::parallel::mpiRank() == 0)
    util::io::print(std::format("mpi_identity_twop_wall done (strategy={})\n",
                                mpi_strategy));
  util::parallel::finalizeMpi();
  return 0;
}

References anonymous_namespace{main.cpp}::buildInputJson(), PeriDEMModel::close(), util::io::InputParser::cmdOptionExists(), anonymous_namespace{main.cpp}::dumpFinalNodal(), util::parallel::finalizeMpi(), util::io::InputParser::getCmdOption(), PeriDEMModel::init(), util::parallel::initMpi(), util::parallel::initNThreads(), PeriDEMModel::integrate(), util::parallel::mpiRank(), util::io::print(), and PeriDEMModel::setPostprocess().

Here is the call graph for this function: