initialise.cpp

/*
 Crown Copyright 2012 AWE.

 This file is part of CloverLeaf.

 CloverLeaf 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.

 CloverLeaf 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
 CloverLeaf. If not, see http://www.gnu.org/licenses/.
 */


//  @brief Top level initialisation routine
//  @author Wayne Gaudin
//  @details Checks for the user input and either invokes the input reader or
//  switches to the internal test problem. It processes the input and strips
//  comments before writing a final input file.
//  It then calls the start routine.

#include "initialise.h"
#include "report.h"
#include "version.h"
#include "read_input.h"
#include "start.h"

#include <fstream>

extern std::ostream g_out;
std::ofstream of;

void initialise(parallel_ &parallel, global_variables& globals) {

  if (parallel.boss) {
    of.open("clover.out");
    if (!of.is_open())
      report_error((char *)"initialise", (char *)"Error opening clover.out file.");

    g_out.rdbuf(of.rdbuf());

  } else {
    g_out.rdbuf(std::cout.rdbuf());
  }


  if (parallel.boss) {
    g_out << "Clover Version " << g_version << std::endl
      << "Kokkos Version" << std::endl
      << "Task Count " << parallel.max_task << std::endl
      << std::endl;

    std::cout << "Output file clover.out opened. All output will go there." << std::endl;
  }

  clover_barrier();

  std::ifstream g_in;
  if (parallel.boss) {
    g_out << "Clover will run from the following input:-" << std::endl
      << std::endl;

    // Try to open clover.in
    g_in.open("clover.in");
    if (!g_in.good()) {
      g_in.close();
      std::ofstream out_unit("clover.in");
      out_unit
        << "*clover" << std::endl
        << " state 1 density=0.2 energy=1.0" << std::endl
        << " state 2 density=1.0 energy=2.5 geometry=rectangle xmin=0.0 xmax=5.0 ymin=0.0 ymax=2.0" << std::endl
        << " x_cells=10" << std::endl
        << " y_cells=2" << std::endl
        << " xmin=0.0" << std::endl
        << " ymin=0.0" << std::endl
        << " xmax=10.0" << std::endl
        << " ymax=2.0" << std::endl
        << " initial_timestep=0.04" << std::endl
        << " timestep_rise=1.5" << std::endl
        << " max_timestep=0.04" << std::endl
        << " end_time=3.0" << std::endl
        << " test_problem 1" << std::endl
        << "*endclover" << std::endl;
      out_unit.close();
      g_in.open("clover.in");
    }
  }

  clover_barrier();

  if (parallel.boss) {
    g_out << std::endl
      << "Initialising and generating" << std::endl
      << std::endl;
  }

  read_input(g_in, parallel, globals);

  clover_barrier();

  globals.step=0;

  start(parallel, globals);

  clover_barrier();

  if (parallel.boss) {
    g_out << "Starting the calculation" << std::endl;
  }

  g_in.close();

}