et_bin_bhns_extr.h

/*
 *  Definition of Lorene class Et_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 __ET_BIN_BHNS_EXTR_H_ 
#define __ET_BIN_BHNS_EXTR_H_ 

/*
 * $Id: et_bin_bhns_extr.h,v 1.5 2014/10/13 08:52:34 j_novak Exp $
 * $Log: et_bin_bhns_extr.h,v $
 * Revision 1.5  2014/10/13 08:52:34  j_novak
 * Lorene classes and functions now belong to the namespace Lorene.
 *
 * Revision 1.4  2005/02/28 23:04:56  k_taniguchi
 * Addition of two indicators for the backfround metric and the boundary
 * condition, and some codes for the conformally flat case
 *
 * Revision 1.3  2005/01/31 20:25:22  k_taniguchi
 * Change the argument of equil_bhns_extr_ylm.
 *
 * Revision 1.2  2004/12/29 16:26:02  k_taniguchi
 * Addition of a comutational method to calculate equilibrium figures
 * with a multipole falloff condition at the outer boundary.
 *
 * Revision 1.1  2004/11/30 20:36:31  k_taniguchi
 * *** empty log message ***
 *
 *
 * $Header: /cvsroot/Lorene/C++/Include/et_bin_bhns_extr.h,v 1.5 2014/10/13 08:52:34 j_novak Exp $
 *
 */

// Lorene headers
#include "etoile.h"

namespace Lorene {
class Et_bin_bhns_extr : public Etoile_bin {

    // Data
    // ----
    protected:

        bool kerrschild ;

    bool multipole ;

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

    Et_bin_bhns_extr(const Et_bin_bhns_extr& ) ;    

    Et_bin_bhns_extr(Map& mp_i, const Eos& eos_i,
             const Base_vect& ref_triad_i, FILE* fich) ;

    virtual ~Et_bin_bhns_extr() ;           
 

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

    // Accessors
    // ---------
    public:

    bool in_kerrschild() const {return kerrschild ;} ;

    bool with_multipole() const {return multipole ;} ;

    // Outputs
    // -------
    public:
    virtual void sauve(FILE *) const ;      
    
    // Computational routines
    // ----------------------
    public:
    void hydro_euler_extr(const double& mass, const double& sepa) ;

    void update_metric_extr(const double& mass, const double& sepa) ;

     void update_metric_der_comp_extr(const double& mass,
                      const double& sepa) ;

    void kinematics_extr(double omega, const double& mass,
                 const double& sepa) ;

    void extrinsic_curv_extr(const double& mass, const double& sepa) ;

    void equil_bhns_extr_ks(double ent_c, const double& mass,
                const double& sepa, int mermax,
                int mermax_poisson, 
                double relax_poisson, int mermax_potvit, 
                double relax_potvit, int np_filter,
                double thres_adapt, Tbl& diff) ;

    void equil_bhns_extr_cf(double ent_c, const double& mass,
                const double& sepa, int mermax,
                int mermax_poisson, 
                double relax_poisson, int mermax_potvit, 
                double relax_potvit, int np_filter,
                double thres_adapt, Tbl& diff) ;

    void equil_bhns_extr_ylm_ks(double ent_c, const double& mass,
                    const double& sepa, double* nu_int,
                    double* beta_int, double* shift_int,
                    int mermax, int mermax_poisson, 
                    double relax_poisson, double relax_ylm,
                    int mermax_potvit, double relax_potvit,
                    int np_filter,
                    double thres_adapt, Tbl& diff) ;

    void equil_bhns_extr_ylm_cf(double ent_c, const double& mass,
                    const double& sepa, double* nu_int,
                    double* beta_int, double* shift_int,
                    int mermax, int mermax_poisson, 
                    double relax_poisson, double relax_ylm,
                    int mermax_potvit, double relax_potvit,
                    int np_filter,
                    double thres_adapt, Tbl& diff) ;

    void test_bhns_extr(const double& mass,
                const double& sepa, int mermax_poisson,
                double relax_poisson, int mermax_potvit,
                double relax_potvit, Tbl& diff) ;

    double velocity_pot_extr(const double& mass, const double& sepa,
                 int mermax, double precis, double relax) ;

    void ent_max_search(double& xx, double& yy) const ;

    double phi_longest_rad(double x_max, double y_max) const ;

    void get_ylm(int nylm, Cmp** ylmvec) const ;

    void get_integrals(int nylm, double* intvec, Cmp** ylmvec,
               Cmp source) const ;

    friend class Bin_bhns_extr ;

};

}
#endif