main_CahnHilliard.C

// $Id$
//==============================================================================
//!
//! \file main_CahnHilliard.C
//!
//! \date Jul 13 2015
//!
//! \author Arne Morten Kvarving
//!
//! \brief Main program for an isogeometric solver for Cahn-Hilliard.
//!
//==============================================================================

#include "IFEM.h"
#include "Profiler.h"
#include "SIM1D.h"
#include "SIM2D.h"
#include "SIM3D.h"
#include "SIMPhaseField.h"
#include "SIMSolver.h"


/*!
  \brief Creates the Cahn-Hilliard simulator and launches the simulation.
  \param[in] infile The input file to parse
*/

template<class Dim> int runSimulator (char* infile)
{
  typedef SIMPhaseField<Dim> PhaseFieldDriver;

  utl::profiler->start("Model input");
  IFEM::cout <<"\n\n0. Parsing input file(s)."
             <<"\n========================="<< std::endl;

  PhaseFieldDriver phaseSim;
  if (!phaseSim.read(infile))
    return 1;

  phaseSim.opt.print(IFEM::cout) << std::endl;

  SIMSolverStat<PhaseFieldDriver> solver(phaseSim);
  if (!solver.read(infile))
    return 1;

  utl::profiler->stop("Model input");
  IFEM::cout <<"\n\n10. Preprocessing the finite element model:"
             <<"\n==========================================="<< std::endl;

  // Preprocess the model and establish data structures for the algebraic system
  if (!phaseSim.preprocess())
    return 2;

  // Initialize the linear solvers
  if (!phaseSim.initSystem(phaseSim.opt.solver))
    return 2;

  // Initialize the solution field
  phaseSim.init(TimeStep());

  if (phaseSim.opt.dumpHDF5(infile))
    solver.handleDataOutput(phaseSim.opt.hdf5,phaseSim.getProcessAdm(),
                            phaseSim.opt.saveInc);

  return solver.solveProblem(infile,"100. Starting the simulation");
}


/*!
  \brief Main program for the isogeometric Cahn-Hilliard solver.
*/

int main (int argc, char** argv)
{
  Profiler prof(argv[0]);

  int ndim = 3;
  char* infile = nullptr;

  IFEM::Init(argc,argv,"Cahn-Hilliard equation solver");

  for (int i = 1; i < argc; i++)
    if (SIMoptions::ignoreOldOptions(argc,argv,i))
      ; // ignore the obsolete option
    else if (!strcmp(argv[i],"-2D"))
      ndim = 2;
    else if (!strcmp(argv[i],"-1D"))
      ndim = 1;
    else if (!infile)
      infile = argv[i];
    else
      std::cerr <<"  ** Unknown option ignored: "<< argv[i] << std::endl;

  if (!infile)
  {
    std::cout <<"usage: "<< argv[0]
              <<" <inputfile> [-dense|-spr|-superlu[<nt>]|-samg|-petsc]\n"
              <<"       [-lag|-spec|-LR] [-1D|-2D] [-nGauss <n>]\n"
              <<"       [-vtf <format> [-nviz <nviz>]"
              <<" [-nu <nu>] [-nv <nv>] [-nw <nw>]] [-hdf5]"<< std::endl;
    return 0;
  }

  IFEM::cout <<"\nInput file: "<< infile;
  IFEM::getOptions().print(IFEM::cout) << std::endl;

  if (ndim == 3)
    return runSimulator<SIM3D>(infile);
  else if (ndim == 2)
    return runSimulator<SIM2D>(infile);
  else
    return runSimulator<SIM1D>(infile);
}