ceps.cpp

Go to the documentation of this file.

/*~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
   This file is part of CEPS.
 
   CEPS is free software: you can redistribute it and/or modify
   it under the terms of the GNU General Public License as published by
   the Free Software Foundation, either version 3 of the License, or
   (at your option) any later version.
 
   CEPS is distributed in the hope that it will be useful,
   but WITHOUT ANY WARRANTY; without even the implied warranty of
   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
   GNU General Public License for more details.
 
   You should have received a copy of the GNU General Public License
   along with CEPS (see file LICENSE at root of project).
   If not, see <https://www.gnu.org/licenses/>.
 
 
   Copyright 2019-2024 Inria, Universite de Bordeaux
 
   Authors, in alphabetical order:
 
     Pierre-Elliott BECUE, Florian CARO, Yves COUDIERE(*), Andjela DAVIDOVIC, 
     Charlie DOUANLA-LONTSI, Marc FUENTES, Mehdi JUHOOR, Michael LEGUEBE(*), 
     Pauline MIGERDITICHAN, Valentin PANNETIER(*), Nejib ZEMZEMI.
     * : currently active authors
 
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~*/
#include "applications/AllProblems.hpp"
#include "applications/runConvergence.hpp"
#include <unistd.h>
 
Flags* flags;
 
void
doSingleRun(InputParameters* p, Geometry* g)
{
  CEPS_SAYS("Single run detected");
 
  AbstractPdeProblem* pb = nullptr;
 
  CepsString s = ceps::toKey(p->getString("problem type"));
  if (s == ceps::toKey("laplacian"))
    pb = ceps::getNew<LaplacianProblem>(g,p);
  else if (s == ceps::toKey("dirichlet anode cathode"))
    pb = ceps::getNew<DirichletAnodeCathodeProblem>(g,p);
  else if (s == ceps::toKey("flux anode cathode"))
    pb = ceps::getNew<FluxAnodeCathodeProblem>(g,p);
  else if (s == ceps::toKey("heat"))
    pb = ceps::getNew<HeatProblem>(g,p);
  else if (s == ceps::toKey("monodomain"))
    pb = ceps::getNew<MonodomainProblem>(g,p);
  else if (s == ceps::toKey("cl monodomain"))
    pb = ceps::getNew<CLMonodomainProblem>(g,p);
  else if (s == ceps::toKey("bidomain"))
    pb = ceps::getNew<BidomainProblem>(g,p);
  else if (s == ceps::toKey("extended bidomain"))
    pb = ceps::getNew<ExtendedBidomainProblem>(g,p);
  else if (s == ceps::toKey("bilayer monodomain"))
    pb = ceps::getNew<BilayerMonodomainProblem>(g,p);
  else if (s == ceps::toKey("pacemaker bidomain"))
    pb = ceps::getNew<PacemakerBidomainProblem>(g,p);
  else if (s == ceps::toKey("pacemaker poisson"))
    pb = ceps::getNew<PacemakerPoissonProblem>(g,p);
  else 
    CEPS_ABORT(
      "Unknown type of problem " << s << "\n" <<
      "   Valid types are:\n" << validProblems
    );
 
  pb->run();
  ceps::destroyObject(pb);
 
  return; 
 
}
 
int 
main(int argc,CepsChar *argv[])
{
 
  ceps::initialize(argc,argv);
  {
    if (argc == optind) 
    {
      flags->printHelp(argv);
      return 0;
    }
    InputParameters p(argv[optind]);
 
    // Copy parameters file into output directory
    CepsString outDir    = ceps::getDir(p("output file name"));
    CepsString inDir     = ceps::getDir(argv[optind]);
    CepsString oldFile   = inDir  + "/" + ceps::getFilename(argv[optind]);
    CepsString newFile   = outDir + "/" + ceps::getFilename(argv[optind]);
    if (ceps::isMaster() and oldFile != newFile)
    {
      CepsString cpCommand = CepsString("mkdir -p ") + outDir + " 2>/dev/null && " 
                           + "rm "  + newFile + " 2>/dev/null && " 
                           + "cp -f " + argv[optind] + " " +  newFile;
      system(cpCommand.c_str());
    }
 
    CEPS_SEPARATOR;
  
    if (p.hasKey("convergence 3d meshes") 
     or p.hasKey("convergence 2d meshes") 
     or p.hasKey("convergence 1d meshes"))
      doConvergence(&p);
    else
    {
      Geometry g(&p);
      doSingleRun(&p,&g);
    }
  }
 
  return 0;
}