tenseur.h
/*
* Definition of Lorene classes Tenseur
* Tenseur_sym
*
*/
/*
* Copyright (c) 1999-2001 Philippe Grandclement
* Copyright (c) 2000-2001 Eric Gourgoulhon
* Copyright (c) 2002 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 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 __TENSEUR_H_
#define __TENSEUR_H_
/*
* $Id: tenseur.h,v 1.19 2016/09/19 15:26:22 j_novak Exp $
* $Log: tenseur.h,v $
* Revision 1.19 2016/09/19 15:26:22 j_novak
* Correction of several bugs preventing the shared library compilation.
*
* Revision 1.18 2014/10/13 08:52:37 j_novak
* Lorene classes and functions now belong to the namespace Lorene.
*
* Revision 1.17 2010/02/02 13:34:12 e_gourgoulhon
* Marked DEPRECATED (in the documentation).
*
* Revision 1.16 2005/08/30 08:35:10 p_grandclement
* Addition of the Tau version of the vectorial Poisson equation for the Tensors
*
* Revision 1.15 2004/12/29 16:24:03 k_taniguchi
* Addition of the method for vectorial Poisson equations with a multipole
* falloff condition at the outer boundary.
*
* Revision 1.14 2004/12/22 18:25:12 k_taniguchi
* Change an argument of poisson_vect_falloff
*
* Revision 1.13 2004/11/30 20:43:32 k_taniguchi
* Addition of the method for vectorial Poisson equations with falloff
* condition at the outer boundary.
*
* Revision 1.12 2004/03/22 13:12:43 j_novak
* Modification of comments to use doxygen instead of doc++
*
* Revision 1.11 2003/11/06 12:17:31 r_prix
* fixed mini-bug in documentation: without explicit argument in function-prototype,
* doc++ seemed to merge docu of Tensor::operator=(Cmp&) and operator=(Tenseur&)
*
* Revision 1.10 2003/06/20 14:23:38 f_limousin
* Add the functions compare().
*
* Revision 1.9 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.8 2002/09/19 14:12:37 e_gourgoulhon
* Modif documentation for LaTeX compliance.
*
* Revision 1.7 2002/09/10 13:44:17 j_novak
* The method "manipule" of one indice has been removed for Tenseur_sym objects
* (the result cannot be a Tenseur_sym).
* The method "sans_trace" now computes the traceless part of a Tenseur (or
* Tenseur_sym) of valence 2.
*
* Revision 1.6 2002/09/06 14:49:25 j_novak
* Added method lie_derive for Tenseur and Tenseur_sym.
* Corrected various errors for derive_cov and arithmetic.
*
* Revision 1.5 2002/08/14 13:46:14 j_novak
* Derived quantities of a Tenseur can now depend on several Metrique's
*
* Revision 1.4 2002/08/08 15:10:44 j_novak
* The flag "plat" has been added to the class Metrique to show flat metrics.
*
* Revision 1.3 2002/08/07 16:14:11 j_novak
* class Tenseur can now also handle tensor densities, this should be transparent to older codes
*
* Revision 1.2 2002/06/17 14:05:17 j_novak
* friend functions are now also declared outside the class definition
*
* Revision 1.1.1.1 2001/11/20 15:19:27 e_gourgoulhon
* LORENE
*
* Revision 2.46 2001/08/27 10:03:56 eric
* Ajout de l'operator% (produit tensoriel avec desaliasing)
*
* Revision 2.45 2001/06/19 15:38:00 eric
* Modif commentaires: mise en conformite Doc++ 3.4.8.
*
* Revision 2.44 2001/06/18 13:56:07 novak
* Ajout de la fonction abs()
*
* Revision 2.43 2001/05/29 16:10:43 eric
* Modif commentaires (mise en conformite Doc++ 3.4.7).
*
* Revision 2.42 2001/05/26 15:48:17 eric
* *** empty log message ***
*
* Revision 2.41 2001/05/26 15:42:54 eric
* Ajout de la fonction flat_scalar_prod_desal (desaliasage)
*
* Revision 2.40 2000/10/19 10:37:51 phil
* *** empty log message ***
*
* Revision 2.39 2000/10/19 09:47:15 phil
* *** empty log message ***
*
* Revision 2.38 2000/10/19 09:41:48 phil
* ajout de inverse_poisson_vect
*
* Revision 2.37 2000/10/06 12:37:09 keisuke
* Add a vectorial Poisson equation Tenseur::poisson_vect_regu.
*
* Revision 2.36 2000/09/27 08:52:34 eric
* Modif commentaires.
*
* Revision 2.35 2000/09/13 12:21:36 eric
* Modif commentaires.
*
* Revision 2.34 2000/09/13 12:11:26 eric
* Ajout de la fonction allocate_all().
*
* Revision 2.33 2000/05/22 14:39:11 phil
* ajout de inc_dzpuis et dec_dzpuis
*
* Revision 2.32 2000/04/03 15:18:54 phil
* suppression de poisson_vect_dirichlet
*
* Revision 2.31 2000/03/31 13:29:43 phil
* poisson_vect_neumann devient poisson_vect_dirichlet
*
* Revision 2.30 2000/03/30 16:10:46 phil
* *** empty log message ***
*
* Revision 2.29 2000/02/18 10:45:05 eric
* Modif commentaires.
*
* Revision 2.28 2000/02/11 19:11:55 phil
* commentaires
*
* Revision 2.27 2000/02/10 16:26:48 eric
* Modif commentaires.
*
* Revision 2.26 2000/02/10 16:10:49 eric
* Ajout de la fonction change_triad.
*
* Revision 2.25 2000/02/09 19:29:25 eric
* MODIF IMPORTANTE: la triade de decomposition est desormais passee en
* argument des constructeurs.
*
* Revision 2.24 2000/02/09 09:53:56 phil
* ajout de poisson_vect_oohara
* ,
*
* Revision 2.23 2000/02/08 19:04:25 eric
* Les fonctions arithmetiques ne sont plus amies.
* Ajout de nouvelles fonctions arithmetiques.
*
* Revision 2.22 2000/02/01 15:40:19 eric
* Ajout de la fonction sqrt
*
* Revision 2.21 2000/02/01 14:13:56 eric
* Modif commentaires.
* Ajout de la fonction amie flat_scalar_prod.
*
* Revision 2.20 2000/01/21 12:48:55 phil
* changement prototypage de Tenseur::poisson_vect
*
* Revision 2.19 2000/01/20 16:02:20 eric
* Ajout des operator=(double ) et operator=(int ).
*
* Revision 2.18 2000/01/20 14:52:08 phil
* *** empty log message ***
*
* Revision 2.17 2000/01/20 14:50:56 phil
* *** empty log message ***
*
* Revision 2.16 2000/01/20 14:39:31 phil
* *** empty log message ***
*
* Revision 2.15 2000/01/20 13:10:56 phil
* Ajout de Tenseur::poisson_vect (double)
*
* Revision 2.14 2000/01/20 11:20:17 phil
* changement prototypage
*
* Revision 2.13 2000/01/20 10:31:30 phil
* ajout de xksk
*
* Revision 2.12 2000/01/14 14:17:47 eric
* Modif commentaires.
*
* Revision 2.11 2000/01/14 14:04:07 eric
* Ajout de la fonction annule.
* classe Tenseur: constructeurs pour les classes derivees et fonctions
* de gestion memoire (del_t(), ...) declarees protected et non plus
* private.
*
* Revision 2.10 2000/01/13 14:15:30 eric
* Modif commentaires.
*
* Revision 2.9 2000/01/13 14:10:18 eric
* Ajout du constructeur par copie d'un Cmp (pour un scalaire)
* ainsi que l'affectation a un Cmp.
*
* Revision 2.8 2000/01/13 13:45:56 eric
* Ajout du membre p_gradient_spher et des fonctions fait_gradient_spher(),
* gradient_spher() pour le calcul du gradient d'un scalaire en
* coordonnees spheriques sur la triade spherique associee.
*
* Revision 2.7 2000/01/12 13:17:49 eric
* Les operator::(...) renvoie desormais une reference const sur le c[...]
* correspondant et non plus un Cmp copie de c[...].
* (ceci grace a Map::cmp_zero()).
*
* Revision 2.6 2000/01/11 11:13:16 eric
* Changement de nom pour la base vectorielle : base --> triad
*
* Revision 2.5 2000/01/10 17:22:26 eric
* Modif des #include
*
* Revision 2.4 2000/01/10 15:14:58 eric
* Ajout du membre base (base vectorielle sur laquelle sont definies
* les composantes).
*
* Revision 2.3 1999/12/09 12:39:39 phil
* changement prototypage des derivees
*
* Revision 2.2 1999/12/07 15:24:17 phil
* ajout include
*
* Revision 2.1 1999/12/03 09:37:11 phil
* *** empty log message ***
*
* Revision 2.0 1999/12/02 17:15:32 phil
* *** empty log message ***
*
* Revision 1.1 1999/12/02 17:13:29 phil
* Initial revision
*
*
* $Header: /cvsroot/Lorene/C++/Include/tenseur.h,v 1.19 2016/09/19 15:26:22 j_novak Exp $
*
*/
#define COV -1
#define CON +1
#define N_MET_MAX 5
// Headers Lorene
#include "cmp.h"
#include "itbl.h"
#include "base_vect.h"
namespace Lorene {
class Metrique ;
class Tenseur_sym ;
//---------------------------------//
// class Tenseur //
//---------------------------------//
class Tenseur {
// Data :
// -----
protected:
const Map* const mp ;
int valence ;
const Base_vect* triad ;
Itbl type_indice ;
int n_comp ;
int etat ;
Cmp** c ;
double poids ;
const Metrique* metric ;
// Derived data :
// ------------
protected:
const Metrique** met_depend ;
mutable Tenseur* p_gradient ;
mutable Tenseur* p_gradient_spher ;
Tenseur** p_derive_cov ;
Tenseur** p_derive_con ;
Tenseur** p_carre_scal ;
// Constructors - Destructor :
// -------------------------
protected:
bool verif() const ;
void new_der_met() ;
public:
explicit Tenseur (const Map& map, const Metrique* met = 0x0,
double weight = 0) ;
explicit Tenseur (const Cmp& cmp, const Metrique* met = 0x0,
double weight = 0) ;
Tenseur (const Map& map, int val, const Itbl& tipe,
const Base_vect& triad_i, const Metrique* met = 0x0,
double weight = 0) ;
Tenseur (const Map& map, int val, const Itbl& tipe,
const Base_vect* triad_i, const Metrique* met = 0x0,
double weight = 0) ;
Tenseur (const Map& map, int val, int tipe, const
Base_vect& triad_i, const Metrique* met = 0x0,
double weight = 0) ;
Tenseur (const Tenseur&) ;
explicit Tenseur (const Tenseur_sym&) ;
Tenseur (const Map& map, const Base_vect& triad_i, FILE* fich,
const Metrique* met = 0x0) ;
Tenseur (const Map& map, FILE* fich, const Metrique* met = 0x0) ;
protected:
Tenseur (const Map& map, int val, const Itbl& tipe, int n_comp,
const Base_vect& triad_i, const Metrique* met = 0x0,
double weight = 0) ;
Tenseur (const Map&, int val, int tipe, int n_comp,
const Base_vect& triad_i, const Metrique* met = 0x0,
double weight = 0) ;
public:
virtual ~Tenseur() ;
// Memory management
// -----------------
protected:
void del_t() ;
void del_derive_met(int i) const ;
void del_derive() const ;
void set_der_met_0x0(int i) const ;
void set_der_0x0() const ;
// Mutators / assignment
// ---------------------
public:
void set_etat_nondef() ;
void set_etat_zero() ;
void set_etat_qcq() ;
void allocate_all() ;
void change_triad(const Base_vect& new_triad) ;
void set_triad(const Base_vect& new_triad) ;
void set_poids(double weight) ;
void set_metric(const Metrique& met) ;
virtual void operator=(const Tenseur& tens) ;
void operator=(const Cmp& field) ;
void operator=(double ) ;
void operator=(int ) ;
Cmp& set () ;
Cmp& set (int) ;
Cmp& set (int, int) ;
Cmp& set (int, int, int) ;
Cmp& set (const Itbl&) ;
void annule(int l) ;
void annule(int l_min, int l_max) ;
void set_std_base() ;
void dec_dzpuis() ;
void inc_dzpuis() ;
void dec2_dzpuis() ;
void inc2_dzpuis() ;
void mult_r_zec() ;
Tenseur inverse_poisson_vect (double lambda) const ;
// Accessors
// ---------
public:
virtual int donne_place (const Itbl& idx) const ;
virtual Itbl donne_indices (int place) const ;
const Map* get_mp() const {return mp ;} ;
const Base_vect* get_triad() const {return triad;} ;
int get_etat() const {return etat ;} ;
int get_valence() const {return valence ; } ;
int get_n_comp() const {return n_comp ;} ;
int get_type_indice (int i) const {return type_indice(i) ;};
Itbl get_type_indice () const {return type_indice ; } ;
double get_poids() const {return poids ; } ;
const Metrique* get_metric() const {return metric ; } ;
const Cmp& operator()() const ;
const Cmp& operator()(int) const ;
const Cmp& operator()(int, int) const ;
const Cmp& operator()(int, int, int) const ;
const Cmp& operator()(const Itbl&) const ;
// Outputs
// -------
public:
void sauve(FILE *) const ;
friend ostream& operator<<(ostream& , const Tenseur & ) ;
// Computation of derived members
// ------------------------------
protected:
virtual void fait_gradient () const ;
void fait_gradient_spher () const ;
virtual void fait_derive_cov (const Metrique& met, int i) const ;
virtual void fait_derive_con (const Metrique&, int i) const ;
void fait_carre_scal (const Metrique&, int i) const ;
void set_dependance (const Metrique& met) const ;
int get_place_met(const Metrique& metre) const ;
// Differential operators
// ----------------------
public:
const Tenseur& gradient() const ;
const Tenseur& gradient_spher() const ;
const Tenseur& derive_cov (const Metrique& met) const ;
const Tenseur& derive_con (const Metrique&) const ;
const Tenseur& carre_scal (const Metrique&) const ;
// Resolution d'EDP :
void poisson_vect(double lambda, Param& par, Tenseur& shift, Tenseur& vect
, Tenseur& scal) const ;
/*
* Same as poisson_vect with a Tau method
**/
void poisson_vect_tau(double lambda, Param& par, Tenseur& shift, Tenseur& vect
, Tenseur& scal) const ;
void poisson_vect_falloff(double lambda, Param& par, Tenseur& shift,
Tenseur& vect, Tenseur& scal, int* k_falloff) const ;
void poisson_vect_ylm(double lambda, Param& para, Tenseur& shift,
Tenseur& vecteur, Tenseur& scalaire, int nylm,
double* intvec) const ;
Tenseur poisson_vect(double lambda, Tenseur& vect , Tenseur& scal ) const ;
/*
* Same as poisson_vect with a Tau method
**/
Tenseur poisson_vect_tau(double lambda, Tenseur& vect , Tenseur& scal ) const ;
Tenseur poisson_vect_falloff(double lambda, Tenseur& vect ,
Tenseur& scal, int* k_falloff ) const ;
Tenseur poisson_vect_ylm(double lambda, Tenseur& vecteur,
Tenseur& scalaire, int nylm, double* intvec) const ;
void poisson_vect_oohara(double lambda, Param& par, Tenseur& shift,
Tenseur& scal) const ;
/*
* Same as poisson_vect_oohara with a Tau method
**/
void poisson_vect_oohara_tau(double lambda, Param& par, Tenseur& shift,
Tenseur& scal) const ;
Tenseur poisson_vect_oohara(double lambda, Tenseur& scal) const ;
/*
* Same as poisson_vect_oohara with a Tau method
**/
Tenseur poisson_vect_oohara_tau(double lambda, Tenseur& scal) const ;
void poisson_vect_regu(int k_div, int nzet, double unsgam1,
double lambda, Param& par, Tenseur& shift,
Tenseur& vect, Tenseur& scal) const ;
// Friend classes
// ---------------
friend class Tenseur_sym ;
friend class Metrique ;
// Mathematical operators
// ----------------------
friend Tenseur operator* (const Tenseur&, const Tenseur&) ;
friend Tenseur operator% (const Tenseur&, const Tenseur&) ;
friend Tenseur contract(const Tenseur&, int id1, int id2) ;
friend Tenseur contract(const Tenseur&, int id1, const Tenseur&,
int id2) ;
friend Tenseur contract_desal(const Tenseur&, int id1, const Tenseur&,
int id2) ;
friend Tenseur flat_scalar_prod(const Tenseur& t1, const Tenseur& t2) ;
friend Tenseur flat_scalar_prod_desal(const Tenseur& t1,
const Tenseur& t2) ;
friend Tenseur manipule(const Tenseur&, const Metrique&, int idx) ;
friend Tenseur manipule(const Tenseur&, const Metrique&) ;
friend Tenseur skxk (const Tenseur&) ;
friend Tenseur lie_derive(const Tenseur& , const Tenseur& ,
const Metrique* ) ;
};
Tenseur operator*(const Tenseur&, const Tenseur&) ;
Tenseur operator%(const Tenseur&, const Tenseur&) ;
Tenseur contract(const Tenseur&, int id1, int id2) ;
Tenseur contract(const Tenseur&, int id1, const Tenseur&, int id2) ;
Tenseur flat_scalar_prod(const Tenseur& t1, const Tenseur& t2) ;
Tenseur flat_scalar_prod_desal(const Tenseur& t1, const Tenseur& t2) ;
Tenseur manipule(const Tenseur&, const Metrique&, int idx) ;
Tenseur manipule(const Tenseur&, const Metrique&) ;
Tenseur skxk (const Tenseur&) ;
Tenseur lie_derive (const Tenseur& t, const Tenseur& x, const Metrique* = 0x0);
Tenseur sans_trace(const Tenseur& tens, const Metrique& metre) ;
Tenseur operator+(const Tenseur& ) ;
Tenseur operator-(const Tenseur& ) ;
Tenseur operator+(const Tenseur&, const Tenseur &) ;
Tenseur operator+(const Tenseur&, double ) ;
Tenseur operator+(double, const Tenseur& ) ;
Tenseur operator+(const Tenseur&, int ) ;
Tenseur operator+(int, const Tenseur& ) ;
Tenseur operator-(const Tenseur &, const Tenseur &) ;
Tenseur operator-(const Tenseur&, double ) ;
Tenseur operator-(double, const Tenseur& ) ;
Tenseur operator-(const Tenseur&, int ) ;
Tenseur operator-(int, const Tenseur& ) ;
Tenseur operator*(const Tenseur&, double ) ;
Tenseur operator*(double, const Tenseur& ) ;
Tenseur operator*(const Tenseur&, int ) ;
Tenseur operator*(int, const Tenseur& ) ;
Tenseur operator/(const Tenseur& a, const Tenseur& b) ;
Tenseur operator/(const Tenseur&, double ) ;
Tenseur operator/(double, const Tenseur &) ;
Tenseur operator/(const Tenseur&, int ) ;
Tenseur operator/(int, const Tenseur &) ;
Tenseur exp(const Tenseur& ) ;
Tenseur log(const Tenseur& ) ;
Tenseur sqrt(const Tenseur& ) ;
Tenseur abs(const Tenseur& ) ;
Tenseur pow(const Tenseur&, int ) ;
Tenseur pow(const Tenseur&, double ) ;
//---------------------------------//
// class Tenseur_sym //
//---------------------------------//
class Tenseur_sym : public Tenseur {
// Constructors - Destructor :
// -------------------------
public:
Tenseur_sym (const Map& map, int val, const Itbl& tipe,
const Base_vect& triad_i, const Metrique* met = 0x0,
double weight = 0) ;
Tenseur_sym (const Map& map, int val, int tipe,
const Base_vect& triad_i, const Metrique* met = 0x0,
double weight = 0) ;
Tenseur_sym (const Tenseur_sym&) ;
explicit Tenseur_sym (const Tenseur&) ;
Tenseur_sym (const Map& map, const Base_vect& triad_i, FILE* fich,
const Metrique* met = 0x0) ;
virtual ~Tenseur_sym() ;
// Mutators / assignment
// ---------------------
public:
virtual void operator= (const Tenseur&) ;
// Accessors
// ---------
public:
virtual int donne_place (const Itbl& idx) const ;
virtual Itbl donne_indices (int place) const ;
// Computation of derived members
// ------------------------------
protected:
virtual void fait_gradient () const ;
virtual void fait_derive_cov (const Metrique& met, int i) const ;
virtual void fait_derive_con (const Metrique&, int i) const ;
// Mathematical operators
// ----------------------
friend Tenseur_sym operator* (const Tenseur&, const Tenseur_sym&) ;
friend Tenseur_sym manipule(const Tenseur_sym&, const Metrique&) ;
friend Tenseur lie_derive (const Tenseur& , const Tenseur& ,
const Metrique* );
} ;
Tenseur_sym operator* (const Tenseur&, const Tenseur_sym&) ;
Tenseur_sym manipule(const Tenseur_sym&, const Metrique&) ;
Tenseur_sym lie_derive (const Tenseur_sym& t, const Tenseur& x,
const Metrique* = 0x0);
Tenseur_sym sans_trace(const Tenseur_sym& tens, const Metrique& metre) ;
Tenseur_sym operator+(const Tenseur_sym& ) ;
Tenseur_sym operator-(const Tenseur_sym& ) ;
Tenseur_sym operator+(const Tenseur_sym&, const Tenseur_sym &) ;
Tenseur_sym operator-(const Tenseur_sym &, const Tenseur_sym &) ;
Tenseur_sym operator*(const Tenseur_sym&, double ) ;
Tenseur_sym operator*(double, const Tenseur_sym& ) ;
Tenseur_sym operator*(const Tenseur_sym&, int ) ;
Tenseur_sym operator*(int, const Tenseur_sym& ) ;
Tenseur_sym operator/(const Tenseur_sym& a, const Tenseur& b) ;
Tenseur_sym operator/(const Tenseur_sym&, double ) ;
Tenseur_sym operator/(const Tenseur_sym&, int ) ;
}
#endif
Lorene::Base_vect
Vectorial bases (triads) with respect to which the tensorial components are defined.
Definition base_vect.h:105
Lorene::Cmp
Component of a tensorial field *** DEPRECATED : use class Scalar instead ***.
Definition cmp.h:446
Lorene::Itbl
Basic integer array class.
Definition itbl.h:122
Lorene::Param
Parameter storage.
Definition param.h:125
Lorene::Tenseur_sym
Class intended to describe tensors with a symmetry on the two last indices *** DEPRECATED : use class...
Definition tenseur.h:1253
Lorene::Tenseur_sym::fait_gradient
virtual void fait_gradient() const
Calculates, if needed, the gradient of *this .
Definition tenseur_sym.C:325
Lorene::Tenseur_sym::fait_derive_cov
virtual void fait_derive_cov(const Metrique &met, int i) const
Calculates, if needed, the covariant derivative of *this , with respect to met .
Definition tenseur_sym.C:372
Lorene::Tenseur_sym::operator*
friend Tenseur_sym operator*(const Tenseur &, const Tenseur_sym &)
Tensorial product.
Definition tenseur_sym_operateur.C:84
Lorene::Tenseur_sym::operator=
virtual void operator=(const Tenseur &)
Assignment from a Tenseur .
Definition tenseur_sym.C:283
Lorene::Tenseur_sym::donne_place
virtual int donne_place(const Itbl &idx) const
Returns the position in the Cmp 1-D array c of a component given by its indices.
Definition tenseur_sym.C:207
Lorene::Tenseur_sym::fait_derive_con
virtual void fait_derive_con(const Metrique &, int i) const
Calculates, if needed, the contravariant derivative of *this , with respect to met .
Definition tenseur_sym.C:441
Lorene::Tenseur_sym::manipule
friend Tenseur_sym manipule(const Tenseur_sym &, const Metrique &)
Raise or lower all the indices, depending on their type, using the given Metrique .
Definition tenseur_sym_operateur.C:137
Lorene::Tenseur_sym::Tenseur_sym
Tenseur_sym(const Map &map, int val, const Itbl &tipe, const Base_vect &triad_i, const Metrique *met=0x0, double weight=0)
Standard constructor.
Definition tenseur_sym.C:96
Lorene::Tenseur_sym::~Tenseur_sym
virtual ~Tenseur_sym()
Destructor.
Definition tenseur_sym.C:201
Lorene::Tenseur_sym::lie_derive
friend Tenseur lie_derive(const Tenseur &, const Tenseur &, const Metrique *)
Lie Derivative of t with respect to x .
Definition tenseur_operateur.C:816
Lorene::Tenseur_sym::donne_indices
virtual Itbl donne_indices(int place) const
Returns the indices of a component given by its position in the Cmp 1-D array c .
Definition tenseur_sym.C:252
Lorene::Tenseur
Tensor handling *** DEPRECATED : use class Tensor instead ***.
Definition tenseur.h:301
Lorene::Tenseur::get_triad
const Base_vect * get_triad() const
Returns the vectorial basis (triad) on which the components are defined.
Definition tenseur.h:704
Lorene::Tenseur::operator*
friend Tenseur operator*(const Tenseur &, const Tenseur &)
Tensorial product.
Definition tenseur_operateur.C:119
Lorene::Tenseur::donne_indices
virtual Itbl donne_indices(int place) const
Returns the indices of a component given by its position in the Cmp 1-D array c .
Definition tenseur.C:704
Lorene::Tenseur::operator%
friend Tenseur operator%(const Tenseur &, const Tenseur &)
Tensorial product with desaliasing.
Definition tenseur_operateur.C:199
Lorene::Tenseur::set
Cmp & set()
Read/write for a scalar (see also operator=(const Cmp&) ).
Definition tenseur.C:824
Lorene::Tenseur::p_derive_con
Tenseur ** p_derive_con
Array of pointers on the contravariant derivatives of *this with respect to the corresponding metri...
Definition tenseur.h:365
Lorene::Tenseur::fait_carre_scal
void fait_carre_scal(const Metrique &, int i) const
Calculates, if needed, the scalar square of *this , with respect to met .
Definition tenseur.C:1517
Lorene::Tenseur::allocate_all
void allocate_all()
Sets the logical state to ETATQCQ (ordinary state) and performs the memory allocation of all the elem...
Definition tenseur.C:657
Lorene::Tenseur::get_type_indice
int get_type_indice(int i) const
Returns the type of the index number i .
Definition tenseur.h:726
Lorene::Tenseur::metric
const Metrique * metric
For tensor densities: the metric defining the conformal factor.
Definition tenseur.h:325
Lorene::Tenseur::fait_gradient_spher
void fait_gradient_spher() const
Calculates, if needed, the gradient of *this in a spherical orthonormal basis (scalar field only).
Definition tenseur.C:1392
Lorene::Tenseur::sauve
void sauve(FILE *) const
Save in a file.
Definition tenseur.C:1325
Lorene::Tenseur::c
Cmp ** c
The components.
Definition tenseur.h:322
Lorene::Tenseur::operator=
virtual void operator=(const Tenseur &tens)
Assignment to another Tenseur.
Definition tenseur.C:718
Lorene::Tenseur::p_gradient_spher
Tenseur * p_gradient_spher
Pointer on the gradient of *this in a spherical orthonormal basis (scalar field only).
Definition tenseur.h:350
Lorene::Tenseur::mp
const Map *const mp
Reference mapping.
Definition tenseur.h:306
Lorene::Tenseur::inc_dzpuis
void inc_dzpuis()
dzpuis += 1 ;
Definition tenseur.C:1117
Lorene::Tenseur::fait_gradient
virtual void fait_gradient() const
Calculates, if needed, the gradient of *this .
Definition tenseur.C:1345
Lorene::Tenseur::met_depend
const Metrique ** met_depend
Array of pointers on the Metrique 's used to calculate derivatives members.
Definition tenseur.h:337
Lorene::Tenseur::donne_place
virtual int donne_place(const Itbl &idx) const
Returns the position in the Cmp 1-D array c of a component given by its indices.
Definition tenseur.C:690
Lorene::Tenseur::derive_cov
const Tenseur & derive_cov(const Metrique &met) const
Returns the covariant derivative of *this , with respect to met .
Definition tenseur.C:1554
Lorene::Tenseur::set_etat_qcq
void set_etat_qcq()
Sets the logical state to ETATQCQ (ordinary state).
Definition tenseur.C:636
Lorene::Tenseur::get_type_indice
Itbl get_type_indice() const
Returns the types of all the indices.
Definition tenseur.h:735
Lorene::Tenseur::poisson_vect_oohara
void poisson_vect_oohara(double lambda, Param &par, Tenseur &shift, Tenseur &scal) const
Solves the vectorial Poisson equation .
Definition tenseur_pde.C:218
Lorene::Tenseur::Tenseur
Tenseur(const Map &map, const Metrique *met=0x0, double weight=0)
Constructor for a scalar field.
Definition tenseur.C:209
Lorene::Tenseur::annule
void annule(int l)
Sets the Tenseur to zero in a given domain.
Definition tenseur.C:900
Lorene::Tenseur::flat_scalar_prod
friend Tenseur flat_scalar_prod(const Tenseur &t1, const Tenseur &t2)
Scalar product of two Tenseur when the metric is : performs the contraction of the last index of t1 w...
Definition tenseur_operateur.C:653
Lorene::Tenseur::get_mp
const Map * get_mp() const
Returns pointer on the mapping.
Definition tenseur.h:699
Lorene::Tenseur::dec2_dzpuis
void dec2_dzpuis()
dzpuis -= 2 ;
Definition tenseur.C:1130
Lorene::Tenseur::p_carre_scal
Tenseur ** p_carre_scal
Array of pointers on the scalar squares of *this with respect to the corresponding metric in *met_d...
Definition tenseur.h:372
Lorene::Tenseur::get_poids
double get_poids() const
Returns the weight.
Definition tenseur.h:738
Lorene::Tenseur::get_place_met
int get_place_met(const Metrique &metre) const
Returns the position of the pointer on metre in the array met_depend .
Definition tenseur.C:598
Lorene::Tenseur::dec_dzpuis
void dec_dzpuis()
dzpuis -= 1 ;
Definition tenseur.C:1104
Lorene::Tenseur::derive_con
const Tenseur & derive_con(const Metrique &) const
Returns the contravariant derivative of *this , with respect to met .
Definition tenseur.C:1568
Lorene::Tenseur::triad
const Base_vect * triad
Vectorial basis (triad) with respect to which the tensor components are defined.
Definition tenseur.h:312
Lorene::Tenseur::gradient
const Tenseur & gradient() const
Returns the gradient of *this (Cartesian coordinates).
Definition tenseur.C:1542
Lorene::Tenseur::~Tenseur
virtual ~Tenseur()
Destructor.
Definition tenseur.C:533
Lorene::Tenseur::contract
friend Tenseur contract(const Tenseur &, int id1, int id2)
Self contraction of two indices of a Tenseur .
Definition tenseur_operateur.C:279
Lorene::Tenseur::n_comp
int n_comp
Number of components, depending on the symmetry.
Definition tenseur.h:320
Lorene::Tenseur::new_der_met
void new_der_met()
Builds the arrays met_depend , p_derive_cov , p_derive_con and p_carre_scal and fills them with null ...
Definition tenseur.C:196
Lorene::Tenseur::set_std_base
void set_std_base()
Set the standard spectal basis of decomposition for each component.
Definition tenseur.C:1170
Lorene::Tenseur::set_etat_zero
void set_etat_zero()
Sets the logical state to ETATZERO (zero state).
Definition tenseur.C:645
Lorene::Tenseur::set_poids
void set_poids(double weight)
Sets the weight for a tensor density.
Definition tenseur.C:680
Lorene::Tenseur::mult_r_zec
void mult_r_zec()
Multiplication by r in the external zone.
Definition tenseur.C:1156
Lorene::Tenseur::type_indice
Itbl type_indice
Array of size valence contening the type of each index, COV for a covariant one and CON for a contrav...
Definition tenseur.h:318
Lorene::Tenseur::del_derive
void del_derive() const
Logical destructor of all the derivatives.
Definition tenseur.C:573
Lorene::Tenseur::verif
bool verif() const
Returns false for a tensor density without a defined metric.
Definition tenseur.C:192
Lorene::Tenseur::operator()
const Cmp & operator()() const
Read only for a scalar.
Definition tenseur.C:943
Lorene::Tenseur::set_metric
void set_metric(const Metrique &met)
Sets the pointer on the metric for a tensor density.
Definition tenseur.C:685
Lorene::Tenseur::flat_scalar_prod_desal
friend Tenseur flat_scalar_prod_desal(const Tenseur &t1, const Tenseur &t2)
Same as flat_scalar_prod but with desaliasing.
Definition tenseur_operateur.C:735
Lorene::Tenseur::valence
int valence
Valence.
Definition tenseur.h:307
Lorene::Tenseur::inc2_dzpuis
void inc2_dzpuis()
dzpuis += 2 ;
Definition tenseur.C:1143
Lorene::Tenseur::gradient_spher
const Tenseur & gradient_spher() const
Returns the gradient of *this (Spherical coordinates) (scalar field only).
Definition tenseur.C:1548
Lorene::Tenseur::del_derive_met
void del_derive_met(int i) const
Logical destructor of the derivatives depending on the i-th element of *met_depend .
Definition tenseur.C:554
Lorene::Tenseur::etat
int etat
Logical state ETATZERO , ETATQCQ or ETATNONDEF.
Definition tenseur.h:321
Lorene::Tenseur::poids
double poids
For tensor densities: the weight.
Definition tenseur.h:323
Lorene::Tenseur::skxk
friend Tenseur skxk(const Tenseur &)
Contraction of the last index of (*this) with or , depending on the type of S .
Definition tenseur_operateur.C:580
Lorene::Tenseur::set_der_met_0x0
void set_der_met_0x0(int i) const
Sets the pointers of the derivatives depending on the i-th element of *met_depend to zero (as well as...
Definition tenseur.C:583
Lorene::Tenseur::p_derive_cov
Tenseur ** p_derive_cov
Array of pointers on the covariant derivatives of *this with respect to the corresponding metric in...
Definition tenseur.h:358
Lorene::Tenseur::get_valence
int get_valence() const
Returns the valence.
Definition tenseur.h:710
Lorene::Tenseur::manipule
friend Tenseur manipule(const Tenseur &, const Metrique &, int idx)
Raise or lower the index idx depending on its type, using the given Metrique .
Definition tenseur_operateur.C:509
Lorene::Tenseur::get_metric
const Metrique * get_metric() const
Returns a pointer on the metric defining the conformal factor for tensor densities.
Definition tenseur.h:745
Lorene::Tenseur::set_der_0x0
void set_der_0x0() const
Sets the pointers of all the derivatives to zero.
Definition tenseur.C:591
Lorene::Tenseur::change_triad
void change_triad(const Base_vect &new_triad)
Sets a new vectorial basis (triad) of decomposition and modifies the components accordingly.
Definition tenseur.C:668
Lorene::Tenseur::lie_derive
friend Tenseur lie_derive(const Tenseur &, const Tenseur &, const Metrique *)
Lie Derivative of t with respect to x .
Definition tenseur_operateur.C:816
Lorene::Tenseur::set_triad
void set_triad(const Base_vect &new_triad)
Assigns a new vectorial basis (triad) of decomposition.
Definition tenseur.C:674
Lorene::Tenseur::inverse_poisson_vect
Tenseur inverse_poisson_vect(double lambda) const
Compute of *this , *this being of valence 1.
Definition tenseur_inv_pois_vect.C:59
Lorene::Tenseur::set_etat_nondef
void set_etat_nondef()
Sets the logical state to ETATNONDEF (undefined state).
Definition tenseur.C:650
Lorene::Tenseur::p_gradient
Tenseur * p_gradient
Pointer on the gradient of *this .
Definition tenseur.h:343
Lorene::Tenseur::get_n_comp
int get_n_comp() const
Returns the number of components.
Definition tenseur.h:713
Lorene::Tenseur::set_dependance
void set_dependance(const Metrique &met) const
To be used to describe the fact that the derivatives members have been calculated with met .
Definition tenseur.C:608
Lorene::Tenseur::del_t
void del_t()
Logical destructor.
Definition tenseur.C:545
Lorene::Tenseur::poisson_vect_regu
void poisson_vect_regu(int k_div, int nzet, double unsgam1, double lambda, Param &par, Tenseur &shift, Tenseur &vect, Tenseur &scal) const
Solves the vectorial Poisson equation : .
Definition tenseur_pde_regu.C:71
Lorene::Tenseur::poisson_vect
void poisson_vect(double lambda, Param &par, Tenseur &shift, Tenseur &vect, Tenseur &scal) const
Solves the vectorial Poisson equation : .
Definition tenseur_pde.C:118
Lorene::Tenseur::fait_derive_cov
virtual void fait_derive_cov(const Metrique &met, int i) const
Calculates, if needed, the covariant derivative of *this , with respect to met .
Definition tenseur.C:1429
Lorene::Tenseur::get_etat
int get_etat() const
Returns the logical state.
Definition tenseur.h:707
Lorene::Tenseur::carre_scal
const Tenseur & carre_scal(const Metrique &) const
Returns the scalar square of *this , with respect to met .
Definition tenseur.C:1577
Lorene::Tenseur::fait_derive_con
virtual void fait_derive_con(const Metrique &, int i) const
Calculates, if needed, the contravariant derivative of *this , with respect to met .
Definition tenseur.C:1499
Lorene::operator*
Base_val operator*(const Base_val &, const Base_val &)
This operator is used when calling multiplication or division of Valeur .
Definition base_val_mult.C:102
Lorene::operator-
Cmp operator-(const Cmp &)
- Cmp
Definition cmp_arithm.C:108
Lorene::sqrt
Cmp sqrt(const Cmp &)
Square root.
Definition cmp_math.C:220
Lorene::operator/
Cmp operator/(const Cmp &, const Cmp &)
Cmp / Cmp.
Definition cmp_arithm.C:457
Lorene::operator%
Cmp operator%(const Cmp &, const Cmp &)
Cmp * Cmp with desaliasing.
Definition cmp_arithm.C:364
Lorene::exp
Cmp exp(const Cmp &)
Exponential.
Definition cmp_math.C:270
Lorene::pow
Cmp pow(const Cmp &, int)
Power .
Definition cmp_math.C:348
Lorene::abs
Cmp abs(const Cmp &)
Absolute value.
Definition cmp_math.C:410
Lorene::operator+
Cmp operator+(const Cmp &)
Definition cmp_arithm.C:104
Lorene::log
Cmp log(const Cmp &)
Neperian logarithm.
Definition cmp_math.C:296
Lorene::skxk
Tenseur skxk(const Tenseur &)
Contraction of the last index of (*this) with or , depending on the type of S .
Definition tenseur_operateur.C:580
Lorene::flat_scalar_prod
Tenseur flat_scalar_prod(const Tenseur &t1, const Tenseur &t2)
Scalar product of two Tenseur when the metric is : performs the contraction of the last index of t1 w...
Definition tenseur_operateur.C:653
Lorene::sans_trace
Tenseur sans_trace(const Tenseur &tens, const Metrique &metre)
Computes the traceless part of a Tenseur of valence 2.
Definition tenseur_operateur.C:897
Lorene::contract
Tenseur contract(const Tenseur &, int id1, int id2)
Self contraction of two indices of a Tenseur .
Definition tenseur_operateur.C:279
Lorene::flat_scalar_prod_desal
Tenseur flat_scalar_prod_desal(const Tenseur &t1, const Tenseur &t2)
Same as flat_scalar_prod but with desaliasing.
Definition tenseur_operateur.C:735
Lorene::lie_derive
Tenseur lie_derive(const Tenseur &t, const Tenseur &x, const Metrique *=0x0)
Lie Derivative of t with respect to x .
Definition tenseur_operateur.C:816
Lorene::manipule
Tenseur manipule(const Tenseur &, const Metrique &, int idx)
Raise or lower the index idx depending on its type, using the given Metrique .
Definition tenseur_operateur.C:509
Lorene
Lorene prototypes.
Definition app_hor.h:64
Lorene::Map
Map(const Mg3d &)
Constructor from a multi-domain 3D grid.
Definition map.C:139
Lorene::x
Coord x
x coordinate centered on the grid
Definition map.h:726