main.cpp

Go to the documentation of this file.

/*
 * -------------------------------------------
 * Copyright (c) 2021 - 2026 Prashant K. Jha
 * -------------------------------------------
 * PeriDEM https://github.com/prashjha/PeriDEM
 *
 * Distributed under the Boost Software License, Version 1.0. (See accompanying
 * file LICENSE)
 */
 
#include <PeriDEMConfig.h>
 
// std includes
#include <iostream>
#include <algorithm>
#include <ctime>
#include <filesystem>
#include <stdexcept>
 
// PeriDEM includes
#include "inp/input.h"                          // Input class
#include "periDEMModel.h"                 // Model class
#include "util/io.h"                            // InputParser class
#include "util/json.h"
#include "util/parallelUtil.h"                  // MPI-related functions
#include "util/vecMethods.h"
#include <format>
 
int main(int argc, char *argv[]) {
 
  // init parallel
  util::parallel::initMpi(argc, argv);
  int mpiSize = util::parallel::mpiSize(), mpiRank = util::parallel::mpiRank();
  util::io::print(std::format("Initialized MPI. MPI size = {}, MPI rank = {}\n", mpiSize, mpiRank));
  util::io::print(util::parallel::getMpiStatus()->printStr());
 
  util::io::InputParser input(argc, argv);
 
  if (input.cmdOptionExists("-h") or !input.cmdOptionExists("-i")) {
    // print help
    std::cout << "Syntax to run PeriDEM: PeriDEM -i <input file> -nThreads <number of threads>" << std::endl;
    std::cout << "Example: PeriDEM -i input.json -nThreads 4" << std::endl;
    exit(EXIT_FAILURE);
  }
 
  unsigned int nThreads;
  if (input.cmdOptionExists("-nThreads")) nThreads = std::stoi(input.getCmdOption("-nThreads"));
  else {
    nThreads = std::thread::hardware_concurrency();
    util::io::print(std::format("Running test with default number of threads = {}\n", nThreads));
  }
  // set number of threads
  util::parallel::initNThreads(nThreads);
  util::io::print(std::format("Number of threads = {}\n", util::parallel::getNThreads()));
  
  // print program version
  std::cout << "PeriDEM"
            << " (Version " << MAJOR_VERSION << "." << MINOR_VERSION << "."
            << UPDATE_VERSION << ")" << std::endl;
 
  // current time
  auto begin = steady_clock::now();
 
  // The library reports invalid input by throwing, so that an embedding
  // process is not terminated. Convert that to a message and an exit code.
  try {
    std::string filename = input.getCmdOption("-i");
    if (!std::filesystem::exists(filename)) {
      throw std::runtime_error(
          std::format("Input file {} does not exist.", filename));
    }
    std::ifstream f(filename);
    auto j = json::parse(f);
    auto deck = std::make_shared<inp::Input>(j);
 
    if (deck->isPeriDEM()) {
      PeriDEMModel dem(deck);
      dem.run(deck);
    } else {
      std::cout << "PeriDEM model not found in input file.\n";
    }
  } catch (const std::exception &e) {
    std::cerr << "PeriDEM error: " << e.what() << std::endl;
    util::parallel::finalizeMpi();
    return EXIT_FAILURE;
  }
 
  // get time elapsed
  auto end = steady_clock::now();
 
  std::cout << "Total simulation time (s) = " 
            << util::methods::timeDiff(begin, end, "seconds") 
            << std::endl;
 
  util::parallel::finalizeMpi();
  return 0;
}