connection.h

/*
 *  Definition of Lorene class Connection
 *
 */

/*
 *   Copyright (c) 2003-2004 Eric Gourgoulhon & Jerome Novak
 *
 *   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 __CONNECTION_H_ 
#define __CONNECTION_H_ 

/*
 * $Id: connection.h,v 1.13 2004/03/22 13:12:40 j_novak Exp $
 * $Log: connection.h,v $
 * Revision 1.13  2004/03/22 13:12:40  j_novak
 * Modification of comments to use doxygen instead of doc++
 *
 * Revision 1.12  2004/01/04 20:50:24  e_gourgoulhon
 * Class Connection: data member delta is now of type Tensor_sym (and no
 *  longer of type Tensor_delta).
 *
 * Revision 1.11  2003/12/30 22:56:40  e_gourgoulhon
 * Replaced member flat_conn (flat connection) by flat_met (flat metric)
 * Added argument flat_met to the constructors of Connection.
 * Suppressed method fait_ricci() (the computation of the Ricci is
 * now devoted to the virtual method ricci()).
 *
 * Revision 1.10  2003/12/27 14:56:20  e_gourgoulhon
 * -- Method derive_cov() suppressed.
 * -- Change of the position of the derivation index from the first one
 *    to the last one in methods p_derive_cov() and p_divergence().
 *
 * Revision 1.9  2003/10/16 14:21:33  j_novak
 * The calculation of the divergence of a Tensor is now possible.
 *
 * Revision 1.8  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.7  2003/10/06 06:52:26  e_gourgoulhon
 * Corrected documentation.
 *
 * Revision 1.6  2003/10/05 21:04:25  e_gourgoulhon
 * Improved comments
 *
 * Revision 1.5  2003/10/03 14:07:23  e_gourgoulhon
 * Added derived class Connection_fcart.
 *
 * Revision 1.4  2003/10/02 21:31:11  e_gourgoulhon
 * Added methods fait_delta and update
 * flat_conn is now a modifiable pointer.
 *
 * Revision 1.3  2003/10/02 15:44:23  j_novak
 * The destructor is now public...
 *
 * Revision 1.2  2003/10/01 15:41:49  e_gourgoulhon
 * Added mapping
 *
 * Revision 1.1  2003/09/29 21:14:10  e_gourgoulhon
 * First version --- not ready yet.
 *
 *
 *
 * $Header: /cvsroot/Lorene/C++/Include/connection.h,v 1.13 2004/03/22 13:12:40 j_novak Exp $
 *
 */


// Lorene headers
#include "tensor.h"

class Metric ; 
class Metric_flat ; 

                //--------------------------//
                //     class Connection     // 
                //--------------------------//

class Connection {

    // Data : 
    // -----
    protected:

    const Map* const mp ;   

    const Base_vect* const triad ; 
    
    Tensor_sym delta ; 

    bool assoc_metric ;


    private:

    const Metric_flat* flat_met ;


    // Derived data : 
    // ------------
    protected:

    mutable Tensor* p_ricci ;   

    // Constructors - Destructor
    // -------------------------
    public:
    
    Connection(const Tensor_sym& delta_i, const Metric_flat& flat_met_i) ;      
    
    Connection(const Metric& met, const Metric_flat& flat_met_i) ;      
    
    Connection(const Connection& ) ;        
    
    protected:

    Connection(const Map&, const Base_vect& ) ;         

        public:
    virtual ~Connection() ;         
 

    // Memory management
    // -----------------
    protected:      

    void del_deriv() const ; 
    
    void set_der_0x0() const ; 


    // Mutators / assignment
    // ---------------------
    public:

    void operator=(const Connection&) ; 
    
    void update(const Tensor_sym& delta_i) ;        
    
    void update(const Metric& met) ;
    

    // Accessors
    // ---------
    public:
    const Map& get_mp() const {return *mp; } ;  


    const Tensor_sym& get_delta() const {return delta; } ; 

    // Computational methods
    // ---------------------
    
    public: 

    virtual Tensor* p_derive_cov(const Tensor& tens) const ; 

    virtual Tensor* p_divergence(const Tensor& tens) const ; 

    virtual const Tensor& ricci() const ; 
    
    private:
    void fait_delta(const Metric& ) ;  

};


                //-------------------------------//
                //     class Connection_flat     // 
                //-------------------------------//

class Connection_flat : public Connection {

  // Constructors - Destructor
  // -------------------------
 protected:
  
  Connection_flat(const Map&, const Base_vect&) ; 

 public:

    Connection_flat(const Connection_flat & ) ; 

  virtual ~Connection_flat() ; 


  // Mutators / assignment
  // ---------------------
 public:

  void operator=(const Connection_flat&) ;  
  

  // Computational methods
  // ---------------------
  
 public: 

        virtual Tensor* p_derive_cov(const Tensor& tens) const = 0 ; 

        virtual Tensor* p_divergence(const Tensor& tens) const = 0 ; 

    virtual const Tensor& ricci() const ; 
    
};


                //-------------------------------//
                //     class Connection_fspher   // 
                //-------------------------------//

class Connection_fspher : public Connection_flat {

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

    Connection_fspher(const Map&, const Base_vect_spher&) ; 

    Connection_fspher(const Connection_fspher& ) ;      
    
 public:

  virtual ~Connection_fspher() ; 


  // Mutators / assignment
  // ---------------------
 public:

  void operator=(const Connection_fspher&) ;    
  

  // Computational methods
  // ---------------------
  
 public: 
    virtual Tensor* p_derive_cov(const Tensor& tens) const ; 

    virtual Tensor* p_divergence(const Tensor& tens) const ; 
  
};



                //-------------------------------//
                //     class Connection_fcart    // 
                //-------------------------------//

class Connection_fcart : public Connection_flat {

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

    Connection_fcart(const Map&, const Base_vect_cart&) ; 

    Connection_fcart(const Connection_fcart& ) ;        
    
 public:

  virtual ~Connection_fcart() ; 


  // Mutators / assignment
  // ---------------------
 public:

  void operator=(const Connection_fcart&) ; 
  

  // Computational methods
  // ---------------------
  
 public: 
    virtual Tensor* p_derive_cov(const Tensor& tens) const ; 

    virtual Tensor* p_divergence(const Tensor& tens) const ; 

  
};




#endif

---