et_bin_nsbh.h

/*
 *  Definition of Lorene class Etoile_bin_nsbh
 *
 */

/*
 *   Copyright (c) 2003 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 __ET_BIN_NSBH_H_ 
#define __ET_BIN_NSBH_H_ 

/*
 * $Id: et_bin_nsbh.h,v 1.9 2014/10/13 08:52:34 j_novak Exp $
 * $Log: et_bin_nsbh.h,v $
 * Revision 1.9  2014/10/13 08:52:34  j_novak
 * Lorene classes and functions now belong to the namespace Lorene.
 *
 * Revision 1.8  2006/09/25 10:01:45  p_grandclement
 * Addition of N-dimensional Tbl
 *
 * Revision 1.7  2006/06/01 12:47:50  p_grandclement
 * update of the Bin_ns_bh project
 *
 * Revision 1.6  2006/04/25 07:21:54  p_grandclement
 * Various changes for the NS_BH project
 *
 * Revision 1.5  2005/10/18 13:12:31  p_grandclement
 * update of the mixted binary codes
 *
 * Revision 1.4  2005/08/29 15:10:12  p_grandclement
 * Addition of things needed :
 *   1) For BBH with different masses
 *   2) Provisory files for the mixted binaries (Bh and NS) : THIS IS NOT
 *   WORKING YET !!!
 *
 * Revision 1.3  2004/06/09 06:32:51  k_taniguchi
 * Introduce set_n_auto() and set_confpsi_auto().
 *
 * Revision 1.2  2003/10/24 12:25:19  k_taniguchi
 * Introduce the method of update metric for NS-BH
 *
 * Revision 1.1  2003/10/21 11:46:13  k_taniguchi
 * Definition of class Et_bin_nsbh
 *
 *
 * $Header: /cvsroot/Lorene/C++/Include/et_bin_nsbh.h,v 1.9 2014/10/13 08:52:34 j_novak Exp $
 *
 */

// Lorene headers
#include "etoile.h"
#include "bhole.h"

namespace Lorene {
class Et_bin_nsbh : public Etoile_bin {

    // Data : 
    // -----
    protected:
    Tenseur n_auto ;

    Tenseur n_comp ;

    Tenseur d_n_auto ;

    Tenseur d_n_comp ;

    Tenseur confpsi ;

    Tenseur confpsi_auto ;

    Tenseur confpsi_comp ;

    Tenseur d_confpsi_auto ;

    Tenseur d_confpsi_comp ;

    Tenseur_sym taij_auto ;

    Tenseur_sym taij_comp ;

    Tenseur_sym taij_tot ;

    Tenseur_sym tkij_auto ;

    Tenseur_sym tkij_tot ;

    Cmp ssjm1_lapse ;

    Cmp ssjm1_confpsi ;

    // Constructors - Destructor
    // -------------------------
    public:
    Et_bin_nsbh(Map& mp_i, int nzet_i, bool relat, const Eos& eos_i,
        bool irrot, const Base_vect& ref_triad_i) ;

    Et_bin_nsbh(const Et_bin_nsbh& ) ;      

    Et_bin_nsbh(Map& mp_i, const Eos& eos_i, const Base_vect& ref_triad_i,
        FILE* fich, bool old = false) ;         

    virtual ~Et_bin_nsbh() ;            
 

    // Mutators / assignment
    // ---------------------
    public:
    void operator=(const Et_bin_nsbh&) ;    
    
    Tenseur& set_n_auto() ;

    Tenseur& set_n_comp() ;

    Tenseur& set_confpsi_auto() ;

    Tenseur& set_confpsi_comp() ;

    // Accessors
    // ---------
    public:
    const Tenseur& get_n_auto() const {return n_auto;} ;

    const Tenseur& get_n_comp() const {return n_comp;} ;

    const Tenseur& get_d_n_auto() const {return d_n_auto;} ;

    const Tenseur& get_d_n_comp() const {return d_n_comp;} ;

    const Tenseur& get_confpsi() const {return confpsi;} ;

    const Tenseur& get_confpsi_auto() const {return confpsi_auto;} ;

    const Tenseur& get_confpsi_comp() const {return confpsi_comp;} ;

    const Tenseur& get_d_confpsi_auto() const {return d_confpsi_auto;} ;

    const Tenseur& get_d_confpsi_comp() const {return d_confpsi_comp;} ;

    const Tenseur_sym& get_taij_auto() const {return taij_auto;} ;

    const Tenseur_sym& get_taij_comp() const {return tkij_comp;} ;

    const Tenseur_sym& get_taij_tot() const {return taij_tot;} ;

    const Tenseur_sym& get_tkij_auto() const {return tkij_auto;} ;

    const Tenseur_sym& get_tkij_tot() const {return tkij_tot;} ;


    // Outputs
    // -------
    public:
    virtual void sauve(FILE *) const ;      

    protected:
    virtual ostream& operator>>(ostream& ) const ;


    // Global quantities
    // -----------------

    // Computational routines
    // ----------------------
    public:
        
    void fait_taij_auto() ;
    void update_metric(const Bhole& comp) ;

    void update_metric_der_comp(const Bhole& comp) ;

    virtual void equilibrium_nsbh(double ent_c, int mermax,
                  int mermax_poisson, double relax_poisson,
                  int mermax_potvit, double relax_potvit,
                  double thres_adapt,
                  const Tbl& fact, Tbl& diff) ;

    void equilibrium_nsbh (bool, double, int&, int, int, double, int, double, Tbl&) ;
    
    virtual void kinematics(double omega, double x_axe) ; 
    
    double compute_angul() const ;
    double compute_axe(double) const ;
    
    friend class Bin_ns_bh ;

};

}
#endif