excision_hor.h

/*
 *  Definition of Lorene class Excision_hor, friend class of Spheroid.
 *
 */

/*
 *   Copyright (c) 2009  Jose-Luis Jaramillo & Nicolas Vasset
 *
 *   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 __EXCISIONHOR_H_ 
#define __EXCISIONHOR_H_ 

/*
 * $Header: /cvsroot/Lorene/C++/Include/excision_hor.h,v 1.3 2014/10/13 08:52:34 j_novak Exp $
 *
 */
#include "metric.h"
#include "spheroid.h"

namespace Lorene {
class Excision_hor {


    // Data : 
    // -----
 protected:
  Spheroid sph ;
  
  Scalar conf_fact ; 
  
  Scalar lapse ;
  
  Vector shift ;
  
  Metric gamij ; 
  
  Sym_tensor Kij ;

  double delta_t;

  double no_of_steps;

  Sym_tensor Tij; 
   

    // Derived data : 
    // ------------
    protected:
  mutable Scalar* p_get_BC_conf_fact ; 
  mutable Scalar* p_get_BC_bmN ; 
  mutable Scalar* p_get_BC_bpN ; 
  mutable Vector* p_get_BC_shift ; 


    // Constructors - Destructor
    // -------------------------
    public:

    Excision_hor(const Scalar& h_in, const Metric& gij, const Sym_tensor& Kij2, const Scalar& ppsi, const Scalar& nn, const Vector& beta,const Sym_tensor& Tij2, double timestep, int int_nos = 1) ;
    Excision_hor(const Excision_hor& ) ;        

    Excision_hor(FILE* ) ;          

    virtual ~Excision_hor() ;           
 

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


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

    // Accessors
    // ---------
    public:
    const Spheroid& get_sph() const {return sph; }; 
    
    const Scalar& get_conf_fact() const {return conf_fact; } ;

    const Scalar& get_lapse() const {return lapse ; } ;

    const Vector& get_shift() const {return shift ; } ;

    const Metric& get_gamij() const {return gamij ; } ;

    const Sym_tensor& get_Kij() const {return Kij ; } ;

    double get_delta_t() const {return delta_t ;};

    double get_no_of_steps() const {return no_of_steps ;};
 
        const Sym_tensor& get_Tij() const{return Tij; } ;

    Scalar& set_conf_fact() {del_deriv() ; return conf_fact ; } ;

    Scalar& set_lapse() {del_deriv() ; return lapse ; } ;

    Vector& set_shift() {del_deriv() ; return shift ; } ;

    Metric& set_gamij() {del_deriv() ; return gamij ; } ;

    Sym_tensor& set_Kij() {del_deriv() ; return Kij ; } ;
 
        Sym_tensor& set_Tij() {del_deriv() ; return Tij; } ;
    

    double& set_delta_t() {del_deriv() ; return delta_t ; } ;

    double& set_no_of_steps() {del_deriv() ; return no_of_steps ; } ;

    
    // Computational functions
    // -----------------------
    public:

    const Scalar& get_BC_conf_fact() const ;
    const Scalar& get_BC_bmN(int choice_bmN, double value = 1.) const ;
    const Scalar& get_BC_bpN(int choice_bpN, double c_bpn_lap = 1., double c_bpn_fin = 1., Scalar *bpN_fin = 0x0) const ;
    const Vector& get_BC_shift(double c_V_lap) const ;
  
    // Outputs
    // -------
    public:
    virtual void sauve(FILE *) const ;      
    
    friend ostream& operator<<(ostream& , const Spheroid& ) ;   

};

ostream& operator<<(ostream& , const Spheroid& ) ;


}
#endif