opal.cpp

Go to the documentation of this file.

#include "opal.h"
 
extern Ippl *ippl;
extern Inform *gmsg;
 
#include "AbstractObjects/OpalData.h"
#include "OpalConfigure/Configure.h"
#include "OpalParser/OpalParser.h"
#include "Parser/FileStream.h"
#include "Utilities/OpalException.h"
#include "Fields/Fieldmap.h"
#include "Structure/IpplInfoWrapper.h"
#include "Utilities/Options.h"
 
#include "OPALconfig.h"
 
#ifdef ENABLE_AMR
    #include <AMReX.H>
#endif
 
#include "Message/Communicate.h"
 
#include "GSLErrorHandling.h"
 
#include <gsl/gsl_errno.h>
 
#include <fstream>
#include <iostream>
#include <limits>
#include <string>
 
int run_opal(char */*args*/[], std::string inputfile, int restartStep,
             int infoLevel, int warnLevel, MPI_Comm comm)
{
    std::string::size_type startExtension    = inputfile.find_last_of('.');
    std::string outputFileName = inputfile.substr(0,startExtension) + ".out";
    std::ofstream output(outputFileName.c_str());
 
    MPI_Barrier(comm);
 
    IpplInfoWrapper *newippl = new IpplInfoWrapper(inputfile, infoLevel, warnLevel, comm);
    gmsg = new Inform("OPAL ", output);
    IpplInfo::Info->setDestination(output);
    IpplInfo::Error->setDestination(output);
    IpplInfo::Warn->setDestination(output);
 
#ifdef ENABLE_AMR
    if (Options::amr)
        amrex::Initialize(comm);
#endif
 
    gsl_set_error_handler(&handleGSLErrors);
 
    OpalData *opal = OpalData::getInstance();
 
    Configure::configure();
    opal->storeInputFn(inputfile);
 
    // FileStream is a RCObject
    FileStream *is = 0;
    try {
        is = new FileStream(inputfile);
    } catch(...) {
        is = 0;
        throw new OpalException("run_opal", "Could not open inputfile: " + inputfile);
    }
 
    // run simulation
    OpalParser *parser = new OpalParser();
 
    if (restartStep > std::numeric_limits<int>::min()) {
        opal->setRestartRun();
        opal->setRestartStep(restartStep);
        opal->setRestartFileName(inputfile.substr(0,startExtension) + ".h5");
    }
 
    if(is) parser->run(is);
 
    Ippl::Comm->barrier();
 
    IpplInfo::Info->setDestination(std::cout);
    IpplInfo::Error->setDestination(std::cout);
    IpplInfo::Warn->setDestination(std::cout);
 
    // cleanup
    //OPAL->reset();
    OpalData::deleteInstance();
    _Fieldmap::clearDictionary();
    delete parser;
    delete gmsg;
 
#ifdef ENABLE_AMR
    if (Options::amr)
        amrex::Finalize(true);
#endif
 
    //FIXME: strange side effects
    //ippl = 0;
    //delete aippl;
 
    //XXX: seems like Ippl is always returning the same instance after the
    //     initial instantiation.
    delete newippl;
 
    output.close();
    return 0;
}