eos_from_file.C

/*
 * Methods for Eos and file manipulation
 *
 * (see file eos.h for documentation)
 */

/*
 *   Copyright (c) 2000-2001 Eric Gourgoulhon
 *
 *   This file is part of LORENE.
 *
 *   LORENE 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 2 of the License, or
 *   (at your option) any later version.
 *
 *   LORENE 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 LORENE; if not, write to the Free Software
 *   Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
 *
 */


char eos_from_file_C[] = "$Header: /cvsroot/Lorene/C++/Source/Eos/eos_from_file.C,v 1.10 2011/06/16 10:49:18 j_novak Exp $" ;

/*
 * $Id: eos_from_file.C,v 1.10 2011/06/16 10:49:18 j_novak Exp $
 * $Log: eos_from_file.C,v $
 * Revision 1.10  2011/06/16 10:49:18  j_novak
 * New class Eos_mag for EOSs depending on density and magnetic field.
 *
 * Revision 1.9  2010/02/02 13:22:16  j_novak
 * New class Eos_Compstar.
 *
 * Revision 1.8  2005/05/22 20:51:41  k_taniguchi
 * Add a new Eos Eos_fit_AkmalPR.
 *
 * Revision 1.7  2004/09/26 18:53:08  k_taniguchi
 * Introduction of new EOSs: Eos_fit_SLy4 and Eos_fit_FPS
 *
 * Revision 1.6  2004/05/07 08:06:45  k_taniguchi
 * Add the case of Eos_multi_poly.C
 *
 * Revision 1.5  2003/12/08 15:47:03  m_bejger
 * GlendNH3 EOS (Glendenning 1985, case 3) added
 *
 * Revision 1.4  2002/10/16 14:36:35  j_novak
 * Reorganization of #include instructions of standard C++, in order to
 * use experimental version 3 of gcc.
 *
 * Revision 1.3  2002/04/09 14:32:15  e_gourgoulhon
 * 1/ Added extra parameters in EOS computational functions (argument par)
 * 2/ New class MEos for multi-domain EOS
 *
 * Revision 1.2  2001/12/04 21:27:53  e_gourgoulhon
 *
 * All writing/reading to a binary file are now performed according to
 * the big endian convention, whatever the system is big endian or
 * small endian, thanks to the functions fwrite_be and fread_be
 *
 * Revision 1.1.1.1  2001/11/20 15:19:27  e_gourgoulhon
 * LORENE
 *
 * Revision 2.6  2001/09/11  16:23:08  eric
 * Ajout de Eos_AkmalPR, Eos_BBB2 et Eos_BalbN1H1.
 *
 * Revision 2.5  2000/11/23  22:34:10  eric
 * Ajout de Eos_BPAL12.
 *
 * Revision 2.4  2000/11/23  14:46:16  eric
 * Ajout de Eos_strange_cr.
 *
 * Revision 2.3  2000/11/22  19:30:55  eric
 * Ajout des Eos_SLy4 et Eos_FPS
 *
 * Revision 2.2  2000/10/24  15:29:22  eric
 * Ajout de l'EOS matiere etrange (Eos_strange).
 *
 * Revision 2.1  2000/02/14  14:33:41  eric
 * Ajout du constructeur par lecture de fichier formate.
 *
 * Revision 2.0  2000/01/21  15:18:08  eric
 * *** empty log message ***
 *
 *
 * $Header: /cvsroot/Lorene/C++/Source/Eos/eos_from_file.C,v 1.10 2011/06/16 10:49:18 j_novak Exp $
 *
 */
 
// Headers C
#include <stdlib.h>

// Header Lorene
#include "headcpp.h"
#include "eos.h"
#include "eos_multi_poly.h"
#include "eos_fitting.h"
#include "utilitaires.h"

        //--------------------------------------//
        //  Identification virtual functions    //
        //--------------------------------------//


int Eos_poly::identify() const      { return 1; }

int Eos_poly_newt::identify() const { return 2; }

int Eos_incomp::identify() const    { return 3; }

int Eos_incomp_newt::identify() const   { return 4; }

int Eos_strange::identify() const   { return 5; }

int Eos_strange_cr::identify() const    { return 6; }

int Eos_SLy4::identify() const      { return 10; }

int Eos_FPS::identify() const       { return 11; }

int Eos_BPAL12::identify() const    { return 12; }

int Eos_AkmalPR::identify() const   { return 13; }

int Eos_BBB2::identify() const      { return 14; }

int Eos_BalbN1H1::identify() const  { return 15; }

int Eos_GlendNH3::identify() const  { return 16; }

int Eos_Compstar::identify() const  { return 17; }

int Eos_mag::identify() const   { return 18; }

int MEos::identify() const  { return 100; }

int Eos_multi_poly::identify() const    { return 110; }

int Eos_fit_SLy4::identify() const      { return 120; }

int Eos_fit_FPS::identify() const       { return 121; }

int Eos_fit_AkmalPR::identify() const   { return 122; }

        //---------------------------------------------//
        //    EOS construction from a binary file      //
        //---------------------------------------------//

Eos* Eos::eos_from_file(FILE* fich) {
    
    Eos* p_eos ; 
    
    // Type (class) of EOS :
    int identificator ;     
    fread_be(&identificator, sizeof(int), 1, fich) ;        

    switch(identificator) {
    
    case 1 : {
        p_eos = new Eos_poly(fich) ; 
        break ; 
    }
    
    case 2 : {
        p_eos = new Eos_poly_newt(fich) ; 
        break ; 
    }
    
    case 3 : {
        p_eos = new Eos_incomp(fich) ; 
        break ; 
    }
    
    case 4 : {
        p_eos = new Eos_incomp_newt(fich) ; 
        break ; 
    }
    
    case 5 : {
        p_eos = new Eos_strange(fich) ;
        break ;
    }
    
    case 6 : {
        p_eos = new Eos_strange_cr(fich) ;
        break ;
    }
    
    case 10 : {
        p_eos = new Eos_SLy4(fich) ;
        break ;
    }
    
    case 11 : {
        p_eos = new Eos_FPS(fich) ;
        break ;
    }
    
    case 12 : {
        p_eos = new Eos_BPAL12(fich) ;
        break ;
    }
    
    case 13 : {
        p_eos = new Eos_AkmalPR(fich) ;
        break ;
    }
    
    case 14 : {
        p_eos = new Eos_BBB2(fich) ;
        break ;
    }
    
    case 15 : {
        p_eos = new Eos_BalbN1H1(fich) ;
        break ;
    }

    case 16 : {
        p_eos = new Eos_GlendNH3(fich) ;
        break ;
    }

    case 17 : {
        p_eos = new Eos_Compstar(fich) ;
        break ;
    }

    case 18 : {
        p_eos = new Eos_mag(fich) ;
        break ;
    }

    case 100 : {
        p_eos = new MEos(fich) ;
        break ;
    }

    case 110 : {
        p_eos = new Eos_multi_poly(fich) ;
        break ;
    }

    case 120 : {
        p_eos = new Eos_fit_SLy4(fich) ;
        break ;
    }

    case 121 : {
        p_eos = new Eos_fit_FPS(fich) ;
        break ;
    }

    case 122 : {
        p_eos = new Eos_fit_AkmalPR(fich) ;
        break ;
    }

    default : {
        cout << "Eos::eos_from_file : unknown type of EOS !" << endl ; 
        cout << " identificator = " << identificator << endl ; 
        abort() ; 
        break ; 
    }
    
    }
    
    return p_eos ; 
    
}

        //----------------------------------------------//
        //    EOS construction from a formatted file    //
        //----------------------------------------------//

Eos* Eos::eos_from_file(ifstream& fich) {
    
    int identificator ; 
    char blabla[80] ;

    // EOS identificator : 
    fich >> identificator ; fich.getline(blabla, 80) ;

    Eos* p_eos ; 
    
    switch(identificator) {
    
    case 1 : {
        p_eos = new Eos_poly(fich) ; 
        break ; 
    }
    
    case 2 : {
        p_eos = new Eos_poly_newt(fich) ; 
        break ; 
    }
    
    case 3 : {
        p_eos = new Eos_incomp(fich) ; 
        break ; 
    }
    
    case 4 : {
        p_eos = new Eos_incomp_newt(fich) ; 
        break ; 
    }
    
    case 5 : {
        p_eos = new Eos_strange(fich) ;
        break ;
    }
    
    case 6 : {
        p_eos = new Eos_strange_cr(fich) ;
        break ;
    }
    
    case 10 : {
        p_eos = new Eos_SLy4(fich) ;
        break ;
    }
    
    case 11 : {
        p_eos = new Eos_FPS(fich) ;
        break ;
    }
    
    case 12 : {
        p_eos = new Eos_BPAL12(fich) ;
        break ;
    }
    
    case 13 : {
        p_eos = new Eos_AkmalPR(fich) ;
        break ;
    }
    
    case 14 : {
        p_eos = new Eos_BBB2(fich) ;
        break ;
    }
    
    case 15 : {
        p_eos = new Eos_BalbN1H1(fich) ;
        break ;
    }

    case 16 : {
        p_eos = new Eos_GlendNH3(fich) ;
        break ;
    }

    case 17 : {
        p_eos = new Eos_Compstar(fich) ;
        break ;
    }

    case 18 : {
        p_eos = new Eos_mag(fich) ;
        break ;
    }

    case 100 : {
        p_eos = new MEos(fich) ;
        break ;
    }

    case 110 : {
        p_eos = new Eos_multi_poly(fich) ;
        break ;
    }

    case 120 : {
        p_eos = new Eos_fit_SLy4(fich) ;
        break ;
    }

    case 121 : {
        p_eos = new Eos_fit_FPS(fich) ;
        break ;
    }

    case 122 : {
        p_eos = new Eos_fit_AkmalPR(fich) ;
        break ;
    }

    default : {
        cout << "Eos::eos_from_file : unknown type of EOS !" << endl ; 
        cout << " identificator = " << identificator << endl ; 
        abort() ; 
        break ; 
    }
    
    }
    
    return p_eos ; 
    
}

---