metric.h

 /*
 *  Definition of Lorene class Metric
 *
 */

/*
 *   Copyright (c) 2003 Eric Gourgoulhon & Jerome Novak
 *
 *   Copyright (c) 1999-2001 Philippe Grandclement (for previous class Metrique)
 *
 *   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 __METRIC_H_ 
#define __METRIC_H_ 

/*
 * $Id: metric.h,v 1.9 2014/10/13 08:52:35 j_novak Exp $
 * $Log: metric.h,v $
 * Revision 1.9  2014/10/13 08:52:35  j_novak
 * Lorene classes and functions now belong to the namespace Lorene.
 *
 * Revision 1.8  2004/11/18 12:23:42  jl_jaramillo
 * radial vector normal to a spherical slicing and pointing towards
 * spatial infinity
 *
 * Revision 1.7  2004/03/22 13:12:42  j_novak
 * Modification of comments to use doxygen instead of doc++
 *
 * Revision 1.6  2003/12/30 23:04:58  e_gourgoulhon
 * Important reorganization of class Metric:
 *  -- suppression of virtual methods fait_* : the actual computations
 *     are now performed via the virtual methods con(), cov(), connect(),
 *     ricci(), ricci_scal(), determinant()
 *  -- the member p_connect is now treated as an ordinary derived data
 *     member
 *  -- the construction of the associated connection (member p_connect)
 *     is performed thanks to the new methods Map::flat_met_spher() and
 *     Map::flat_met_cart().
 *
 * Revision 1.5  2003/11/06 14:43:37  e_gourgoulhon
 * Gave a name to const arguments in certain method prototypes (e.g.
 * constructors) to correct a bug of DOC++.
 *
 * Revision 1.4  2003/10/06 15:30:32  j_novak
 * Defined methods for flat metric.
 *
 * Revision 1.3  2003/10/06 13:58:45  j_novak
 * The memory management has been improved.
 * Implementation of the covariant derivative with respect to the exact Tensor
 * type.
 *
 * Revision 1.2  2003/10/03 11:21:45  j_novak
 * More methods for the class Metric
 *
 * Revision 1.1  2003/10/02 15:45:48  j_novak
 * New class Metric
 *
 *
 *
 *
 * $Header: /cvsroot/Lorene/C++/Include/metric.h,v 1.9 2014/10/13 08:52:35 j_novak Exp $
 *
 */

// Lorene headers
#include "connection.h"

#define N_TENSOR_DEPEND 200

namespace Lorene {
class Metric {

    // Data : 
    // -----
    protected:
    const Map* const mp ;   

    mutable Sym_tensor* p_met_cov ;

    mutable Sym_tensor* p_met_con ;


    // Derived data : 
    // ------------
    protected:
    
    mutable Connection* p_connect ; 

    mutable Scalar* p_ricci_scal ;

    mutable Vector* p_radial_vect ;


    
    mutable Scalar* p_determinant ;
    
    mutable const Tensor* tensor_depend[N_TENSOR_DEPEND] ;
    
    // Constructors - Destructor
    // -------------------------
    public:
    explicit Metric(const Sym_tensor& tens) ;  
         
    Metric(const Metric& met) ;     

    Metric(const Map&, FILE* ) ;            

    protected:
    explicit Metric(const Map& mpi) ;

    public:
    virtual ~Metric() ;         
 

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

    void del_tensor_depend() const ;
        
    void set_tensor_depend_0x0() const ;
        

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

    virtual void operator=(const Sym_tensor& tens) ;
    
    // Accessors
    // ---------
    public:
    const Map& get_mp() const {return *mp ; } ;

    virtual const Sym_tensor& cov() const ;

    virtual const Sym_tensor& con() const ;

    virtual const Connection& connect() const ;

    const Sym_tensor& ricci() const ;
    
    virtual const Scalar& ricci_scal() const ;

    virtual const Vector& radial_vect() const ;


    virtual const Scalar& determinant() const ;





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

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


    friend class Tensor ;

};

class Metric_flat: public Metric {

    // Data : 
    // -----
    protected:
    
  const Base_vect* triad ; 

    // Constructors - Destructor
    // -------------------------
    public:
    Metric_flat(const Map&, const Base_vect& ) ;   

    Metric_flat(const Metric_flat& ) ;      

    Metric_flat(const Map&, FILE* ) ;           

   public:
    virtual ~Metric_flat() ;            
 

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

    virtual void operator=(const Sym_tensor& tens) ;


    // Accessors
    // ---------
    public:
    const Base_vect* get_triad() const {return triad;} ; 
    
    virtual const Sym_tensor& cov() const ;

    virtual const Sym_tensor& con() const ;

    virtual const Connection& connect() const ;

    virtual const Scalar& ricci_scal() const ;

    virtual const Scalar& determinant() const ;


    // Outputs
    // -------
    public:
    virtual void sauve(FILE *) const ;      
    
    protected:
    virtual ostream& operator>>(ostream& ) const ;    


};



}
#endif