bin_bhns_extr.h

/*
 *  Definition of Lorene class Bin_bhns_extr
 *
 */

/*
 *   Copyright (c) 2004-2005 Keisuke Taniguchi
 *
 *   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 version 2
 *   as published by the Free Software Foundation.
 *
 *   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
 *
 */

#ifndef __BIN_BHNS_EXTR_H_ 
#define __BIN_BHNS_EXTR_H_ 

/*
 * $Id: bin_bhns_extr.h,v 1.3 2014/10/13 08:52:32 j_novak Exp $
 * $Log: bin_bhns_extr.h,v $
 * Revision 1.3  2014/10/13 08:52:32  j_novak
 * Lorene classes and functions now belong to the namespace Lorene.
 *
 * Revision 1.2  2005/02/28 23:02:51  k_taniguchi
 * Addition of some codes for the conformally flat case
 *
 * Revision 1.1  2004/11/30 20:36:00  k_taniguchi
 * *** empty log message ***
 *
 *
 * $Header: /cvsroot/Lorene/C++/Include/bin_bhns_extr.h,v 1.3 2014/10/13 08:52:32 j_novak Exp $
 *
 */

// Lorene header
#include "et_bin_bhns_extr.h"
#include "etoile.h"

namespace Lorene {
class Bin_bhns_extr {

    // Data : 
    // -----
    private:
        const Base_vect_cart ref_triad ;

    Et_bin_bhns_extr star ;

    double omega ;

    double separ ;

    double mass_bh ;

    // Derived data :
    // ------------
    private :
    mutable double* p_xa_barycenter_extr ;

    mutable double* p_ya_barycenter_extr ;

        mutable double* p_mass_b_extr ;

    // Constructors - Destructor
    // -------------------------
    public:
    Bin_bhns_extr(Map& mp, int nzet, const Eos& eos,
              bool irrot, bool relat, bool kerrs, bool multi) ;

    Bin_bhns_extr(const Bin_bhns_extr& ) ;      

    Bin_bhns_extr(Map& mp, const Eos& eos, FILE* fich) ;

    ~Bin_bhns_extr() ;          
 

    // Memory management
    // -----------------
    private:
    void del_deriv() const ; 
    
    void set_der_0x0() const ; 


    // Mutators / assignment
    // ---------------------
    public:
    void operator=(const Bin_bhns_extr&) ;

    Et_bin_bhns_extr& set_ns()
        { del_deriv() ;
          return star ; } ;

    double& set_omega() { return omega ; } ;

    double& set_separ() { return separ ; } ;

    double& set_mass_bh() { return mass_bh ; } ;

    // Accessors
    // ---------
    public:
    const Et_bin_bhns_extr& get_ns() const
         { return star ; } ;

    double get_omega() const { return omega ; } ;

    double get_separ() const { return separ ; } ;

    double get_mass_bh() const { return mass_bh ; } ;

    // Outputs
    // -------
    public:
    void sauve(FILE* ) const ;      
    
    friend ostream& operator<<(ostream& , const Bin_bhns_extr& ) ;  

    void display_poly(ostream& ) const ;

    private:
    ostream& operator>>(ostream& ) const ;

    // Computational routines
    // ----------------------
    public:
    double xa_barycenter_extr() const ;

    // Absolute coordinate Y of the barycenter of the baryon density
    double ya_barycenter_extr() const ;

    double mass_b_extr() const ;

    void orbit_omega_ks(double fact_omeg_min, double fact_omeg_max) ;

    void orbit_omega_cf(double fact_omeg_min, double fact_omeg_max) ;

    void analytical_omega() ;

    void analytical_shift() ;

};
ostream& operator<<(ostream& , const Bin_bhns_extr& ) ;

}
#endif