input.h

Go to the documentation of this file.

/* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *\
 *                                                                           *
 *                       / /   / /    __    \ \  / /                         *
 *                      / /__ / /   / _ \    \ \/ /                          *
 *                     /  ___  /   | |/_/    / /\ \                          *
 *                    / /   / /    \_\      / /  \ \                         *
 *                                                                           *
 *                         Jakub Benda (c) 2014                              *
 *                     Charles University in Prague                          *
 *                                                                           *
\* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * */

#ifndef HEX_INPUT
#define HEX_INPUT

#include <cstdio>
#include <fstream>
#include <vector>
#include <string>

#include "arrays.h"

class CommandLine
{
    public:
        
        typedef enum {
            StgNone    = 0x00,
            StgRadial  = 0x01,
            StgSolve   = 0x02,
            StgExtract = 0x04
        } HexEcsStg;
        
        // constructor
        CommandLine (int argc, char* argv[])
            : zipcount(0), zipmax(-1), parallel(false), preconditioner(0),
              droptol(1e-15), itinerary(StgNone), outofcore(false), itertol(1e-10),
              parallel_dot(false), parallel_block(true)
        {
            // get command line options
            parse(argc, argv);
            
            // run the whole sequence if nothing specified
            if (itinerary == StgNone)
                itinerary = StgRadial | StgSolve | StgExtract;
        }
        
        void parse (int argc, char* argv[]);
        
        //
        // public attributes
        //
        
        std::ifstream inputfile;
        
        std::string zipfile;
        
        int  zipcount;
        
        double zipmax;
        
        bool parallel;
        
        int preconditioner;
        
        double droptol;
        
        int itinerary;
        
        bool outofcore;
        
        double itertol;
        
        bool parallel_dot;
        
        bool parallel_block;
};

class InputFile
{
    public:
        
        // constructor
        InputFile (std::ifstream & inputfile)
        {
            // read inputfile
            read(inputfile);
            
            // compute angular momentum limit
            maxell = levels + L;
            
            // projectile momenta
            ki.resize(Ei.size());
            for (size_t i = 0; i < Ei.size(); i++)
                ki[i] = sqrt(Ei[i]);
        }
        
        // read data from file
        void read (std::ifstream & inputfile);
        
        //
        // public attributes
        //
        
        int order, ni, L, Spin, Pi, levels, maxell;
        double ecstheta, B;
        rArray rknots, cknots, Ei, ki;
        std::vector<std::tuple<int,int,int>> instates, outstates;
};

#endif