grilles.h

/*
 *  Definition of Lorene classes Grille3d
 *               Grille3d_*
 *               Mg3d
 *
 */

/*
 *   Copyright (c) 1999-2000 Jean-Alain Marck
 *   Copyright (c) 1999-2001 Eric Gourgoulhon
 *
 *   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 as published by
 *   the Free Software Foundation; either version 2 of the License, or
 *   (at your option) any later version.
 *
 *   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 __GRILLES_H_ 
#define __GRILLES_H_ 

/*
 * $Id: grilles.h,v 1.19 2012/01/17 10:10:15 j_penner Exp $
 * $Log: grilles.h,v $
 * Revision 1.19  2012/01/17 10:10:15  j_penner
 * added a constructor in which the nucleus and outer domain are both of type FIN
 *
 * Revision 1.18  2008/10/29 08:17:51  jl_cornou
 * Standard spectral bases for pseudo vectors added
 *
 * Revision 1.17  2008/02/18 13:53:37  j_novak
 * Removal of special indentation instructions.
 *
 * Revision 1.16  2007/12/11 15:28:05  jl_cornou
 * Jacobi(0,2) polynomials partially implemented
 *
 * Revision 1.15  2006/05/17 13:17:02  j_novak
 * New member g_angu_1dom, the one-domain angular grid associated with the
 * current grid.
 *
 * Revision 1.14  2005/10/25 08:56:34  p_grandclement
 * addition of std_spectral_base in the case of odd functions near the origin
 *
 * Revision 1.13  2005/10/07 08:47:20  j_novak
 * Addition of the pointer g_non_axi on a grid, with at least 5 points in the
 * theta direction and 4 in the phi one (for tensor rotations).
 *
 * Revision 1.12  2005/03/25 14:54:04  e_gourgoulhon
 * Corrected documentation.
 *
 * Revision 1.11  2004/07/06 13:36:27  j_novak
 * Added methods for desaliased product (operator |) only in r direction.
 *
 * Revision 1.10  2004/06/22 08:49:56  p_grandclement
 * Addition of everything needed for using the logarithmic mapping
 *
 * Revision 1.9  2004/03/22 13:12:41  j_novak
 * Modification of comments to use doxygen instead of doc++
 *
 * Revision 1.8  2003/06/20 14:16:57  f_limousin
 * Add the operator== to compare two Mg3d
 *
 * Revision 1.7  2003/06/18 08:45:26  j_novak
 * In class Mg3d: added the member get_radial, returning only a radial grid
 * For dAlembert solver: the way the coefficients of the operator are defined has been changed.
 *
 * Revision 1.6  2002/10/16 14:36:29  j_novak
 * Reorganization of #include instructions of standard C++, in order to
 * use experimental version 3 of gcc.
 *
 * Revision 1.5  2002/08/13 08:02:45  j_novak
 * Handling of spherical vector/tensor components added in the classes
 * Mg3d and Tenseur. Minor corrections for the class Metconf.
 *
 * Revision 1.4  2002/06/17 14:05:16  j_novak
 * friend functions are now also declared outside the class definition
 *
 * Revision 1.3  2001/12/12 09:23:46  e_gourgoulhon
 * Parameter compact added to the simplified constructor of class Mg3d
 *
 * Revision 1.2  2001/12/11 06:47:42  e_gourgoulhon
 * Simplified constructor for class Mg3d
 *
 * Revision 1.1.1.1  2001/11/20 15:19:27  e_gourgoulhon
 * LORENE
 *
 * Revision 2.10  2001/06/13  14:23:40  eric
 * Les fonctions Mg3d::del_deriv() et Mg3d::set_deriv_0x0() ne sont plus
 * virtuelles puisque Mg3d n'a aucune classe derivee.
 *
 * Revision 2.9  2001/05/26  13:24:49  eric
 * Ajout du membre g_twice (grille double pour le desaliasing)
 * Modif de la declaration de g_angu (pointeur mutable)
 *   g_twice et g_angu ne sont calcules que si necessaire (cad si
 *   on appelle la fonction get_twice() ou get_angu()).
 *
 * Revision 2.8  1999/11/16  14:15:57  eric
 * Ajout de la fonction Mg3d::get_angu().
 *
 * Revision 2.7  1999/10/12  14:54:11  eric
 * Ajout du membre Base_val std_base_scal() const.
 *
 * Revision 2.6  1999/10/01  10:35:42  eric
 * Amelioration des commentaires.
 *
 * Revision 2.5  1999/09/30  14:58:00  eric
 * Operator!= declare const/
 *
 * Revision 2.4  1999/09/30  14:11:43  eric
 * sauve et std_base_vect_cart declarees const.
 *
 * Revision 2.3  1999/09/30  12:52:38  eric
 * Depoussierage.
 * Documentation.
 *
 * Revision 2.2  1999/09/24  14:23:24  eric
 * Declaration de methodes const.
 *
 * Revision 2.1  1999/09/14  15:24:04  phil
 * ajout de std_base_vect_cart
 *
 * Revision 2.0  1999/02/15  10:41:51  hyc
 * *** empty log message ***
 *
 * Revision 2.1  1999/02/15  09:59:50  hyc
 * *** empty log message ***
 *
 * Revision 2.0  1998/12/01  14:28:00  hyc
 * Version 2
 *
 *
 * $Header: /cvsroot/Lorene/C++/Include/grilles.h,v 1.19 2012/01/17 10:10:15 j_penner Exp $
 *
 */

// Classes utilisees

// Fichiers includes
#include <assert.h>
#include <stdio.h>
#include <stdlib.h>
#include "headcpp.h"

#include "type_parite.h"

class Base_val ; 

                //-------------//
                // Mono-grille //
                //-------------//

// Classe de base
class Grille3d {
    protected:
    const int nr ;  
    const int nt ;  
    const int np ;  

    int type_r ;    
    int type_t ;    
    int type_p ;    
    public:
    double* x ; 
    double* tet ;   
    double* phi ;   

    protected: 
    Grille3d(int n_r, int n_t, int n_p) ;
    
    private:
    Grille3d(const Grille3d& ) ;
    
    void operator=(const Grille3d& ) ;
        
    public:
    virtual ~Grille3d() ;       

    public:
        int get_nr() const {return nr ;} ;
        int get_nt() const {return nt ;} ;
        int get_np() const {return np ;} ;

        int get_type_r() const {return type_r ;} ; 
        int get_type_t() const {return type_t ;} ; 
        int get_type_p() const {return type_p ;} ; 

};


            // -------------------- //
                    // Les classes derivees //
            // -------------------- //


// Cas rare + sans symetrie
// ------------------------
class Grille3d_r : public Grille3d {
    public:
    Grille3d_r(int n_r, int n_t, int n_p) ; 
    ~Grille3d_r() ; 
};

// cas fin + sans symetrie
// -----------------------
class Grille3d_f : public Grille3d {
    public:
    Grille3d_f(int n_r, int n_t, int n_p) ; 
    ~Grille3d_f() ; 
};

// cas echantillonnage (fin) en 1/r + sans symetrie
// ------------------------------------------------
class Grille3d_i : public Grille3d {
    public:
    Grille3d_i(int n_r, int n_t, int n_p) ; 
    ~Grille3d_i() ; 
};

// Cas rare + symetrie equatoriale
// -------------------------------
class Grille3d_req : public Grille3d {
    public:
    Grille3d_req(int n_r, int n_t, int n_p) ; 
    ~Grille3d_req() ; 
};

// cas fin + symetrie equatoriale
// ------------------------------
class Grille3d_feq : public Grille3d {
    public:
    Grille3d_feq(int n_r, int n_t, int n_p) ; 
    ~Grille3d_feq() ; 
};

// cas echantillonnage (fin) en 1/r + symetrie equatoriale
// -------------------------------------------------------
class Grille3d_ieq : public Grille3d {
    public:
    Grille3d_ieq(int n_r, int n_t, int n_p) ; 
    ~Grille3d_ieq() ; 
};

// Cas rare supersymetrique
// ------------------------
class Grille3d_rs : public Grille3d {
    public:
    Grille3d_rs(int n_r, int n_t, int n_p) ; 
    ~Grille3d_rs() ; 
};

// cas fin supersymetrique
// -----------------------
class Grille3d_fs : public Grille3d {
    public:
    Grille3d_fs(int n_r, int n_t, int n_p) ; 
    ~Grille3d_fs() ; 
};

// cas echantillonnage (fin) en 1/r supersymetrique
// --------------------------
class Grille3d_is : public Grille3d {
    public:
    Grille3d_is(int n_r, int n_t, int n_p) ; 
    ~Grille3d_is() ; 
};

// Cas rare + 2 phi
// ----------------
class Grille3d_r2p : public Grille3d {
    public:
    Grille3d_r2p(int n_r, int n_t, int n_p) ; 
    ~Grille3d_r2p() ; 
};

// cas fin + 2 phi
// ---------------
class Grille3d_f2p : public Grille3d {
    public:
    Grille3d_f2p(int n_r, int n_t, int n_p) ; 
    ~Grille3d_f2p() ; 
};

// cas echantillonnage (fin) en 1/r + 2 phi
// ----------------------------------------
class Grille3d_i2p : public Grille3d {
    public:
    Grille3d_i2p(int n_r, int n_t, int n_p) ; 
    ~Grille3d_i2p() ; 
};

// cas fin de jacobi + sans symetrie
// ---------------------------------
class Grille3d_fj : public Grille3d {
    public:
    Grille3d_fj(int n_r, int n_t, int n_p) ; 
    ~Grille3d_fj() ; 
};

// cas fin de jacobi + symetrie equatoriale
// ----------------------------------------
class Grille3d_fjeq : public Grille3d {
    public:
    Grille3d_fjeq(int n_r, int n_t, int n_p) ; 
    ~Grille3d_fjeq() ; 
};

// cas fin de jacobi supersymetrique
// ---------------------------------
class Grille3d_fjs : public Grille3d {
    public:
    Grille3d_fjs(int n_r, int n_t, int n_p) ; 
    ~Grille3d_fjs() ; 
};

// cas fin de jacobi + 2 phi
// -------------------------
class Grille3d_fj2p : public Grille3d {
    public:
    Grille3d_fj2p(int n_r, int n_t, int n_p) ; 
    ~Grille3d_fj2p() ; 
};


                //---------------//
                // Multi-grilles //
                //---------------//

class Mg3d {

    // Data  
    // ----
    protected:
    int nzone ; 
    
    int* nr ;   
    int* nt ;   
    int* np ;   
    
    int* type_r ;   
    int type_t ;
    int type_p ;    
    
    Grille3d** g ;  

    mutable Mg3d* g_angu ;  

    mutable Mg3d* g_angu_1dom ; 
    mutable Mg3d* g_radial ; 
    
    mutable Mg3d* g_twice ; 

    mutable Mg3d* g_plus_half ; 

    mutable Mg3d* g_non_axi ; 

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

    Mg3d(int nz, int nbr[], int typr[], int nbt[], int typt, int nbp[],
         int typp) ;

    Mg3d(int nz, int nbr, int nbt, int nbp, int typt, int typp, 
         bool compact) ;

    Mg3d(int nz, int nbr, int nbt, int nbp, int typt, int typp) ;

    Mg3d(FILE* ) ;   
    
   public:
    Mg3d(const Mg3d& ) ;
    
    public:
    
    ~Mg3d() ;   
        
    // Assignement
    // -----------
    private:
    void operator=(const Mg3d& ) ;
        
    // Extraction of information
    // -------------------------
    public:
    int get_nzone() const {      
        return nzone ;
    } ;
    int get_nr(int l) const { 
        assert(l>=0 && l<nzone) ;
        return nr[l] ;
    } ;
    int get_nt(int l) const { 
        assert(l>=0 && l<nzone) ;
        return nt[l] ;
    } ;
    int get_np(int l) const { 
        assert(l>=0 && l<nzone) ;
        return np[l] ;
    } ;
    
    int get_type_r(int l) const {
        assert(l>=0 && l<nzone) ;
        return type_r[l] ;
    } ;

    int get_type_t() const {    
        return type_t ;
    } ;

    int get_type_p() const {
        return type_p ;
    } ;
    
    const Grille3d* get_grille3d(int l) const { 
        assert(l>=0 && l<nzone) ;
        return g[l] ;
    } ;

    const Mg3d* get_angu() const ;
    
    const Mg3d* get_angu_1dom() const ;
    
    const Mg3d* get_radial() const ;
    
    const Mg3d* get_twice() const ;

    const Mg3d* plus_half() const ;

    const Mg3d* get_non_axi() const ;

    bool operator==(const Mg3d& ) const ;  


    
    // Outputs
    // -------
    public: 
    void sauve(FILE* ) const ; 
    
    friend ostream& operator<<(ostream& , const Mg3d & ) ;  
    
    // Management of derived quantities
    // --------------------------------
    protected:
    void del_deriv() const ; 
    
    void set_deriv_0x0() const ; 


    // Miscellaneous
    // -------------
    public:
    bool operator!=(const Mg3d & ) const ;  

    Base_val std_base_scal() const ;    
    
    Base_val std_base_scal_odd() const ;    
    
    Base_val** std_base_vect_cart() const ;

    Base_val** std_base_vect_spher() const ;

    Base_val** pseudo_base_vect_cart() const ;

    Base_val** pseudo_base_vect_spher() const ;

};
ostream& operator<<(ostream& , const Mg3d & ) ;

#endif 

---