28char connection_fcart_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Connection/connection_fcart.C,v 1.14 2014/10/13 08:52:50 j_novak Exp $" ;
97#include "connection.h"
138 cout <<
"Connection_fcart::operator= : not implemented yet !" << endl ;
158 int valence1 = valence0 + 1 ;
159 int valence1m1 = valence1 - 1 ;
178 Itbl tipe(valence1) ;
180 for (
int id = 0;
id<valence0;
id++) {
181 tipe.
set(
id) = tipeuu(
id) ;
183 tipe.
set(valence1m1) = COV ;
200 Itbl ind1(valence1) ;
201 Itbl ind0(valence0) ;
205 for (
int ic=0; ic<ncomp1; ic++) {
214 for (
int id = 0;
id < valence0;
id++) {
215 ind0.
set(
id) = ind1(
id) ;
219 int k = ind1(valence1m1) ;
223 cresu = (uu(ind0)).deriv(k) ;
244 int valence1 = valence0 - 1 ;
248 assert (valence0 >= 1) ;
264 Itbl tipe(valence1) ;
266 for (
int id = 0;
id<valence1;
id++) {
267 tipe.
set(
id) = tipeuu(
id) ;
305 Itbl ind0(valence0) ;
307 Itbl ind1(valence1) ;
310 for (
int ic=0; ic<ncomp1; ic++) {
316 for (
int k=1; k<=3; k++) {
319 for (
int id = 0;
id < valence1;
id++) {
320 ind0.
set(
id) = ind1(
id) ;
322 ind0.
set(valence0-1) = k ;
324 cresu += uu(ind0).deriv(k) ;
Cartesian vectorial bases (triads).
void operator=(const Connection_fcart &)
Assignment to another Connection_fcart.
Connection_fcart(const Map &, const Base_vect_cart &)
Contructor from a Cartesian flat-metric-orthonormal basis.
virtual ~Connection_fcart()
destructor
virtual Tensor * p_divergence(const Tensor &tens) const
Computes the divergence of a tensor (with respect to the current connection).
virtual Tensor * p_derive_cov(const Tensor &tens) const
Computes the covariant derivative of a tensor (with respect to the current connection).
const Base_vect *const triad
Triad with respect to which the connection coefficients are defined.
const Map *const mp
Reference mapping.
Basic integer array class.
int & set(int i)
Read/write of a particular element (index i ) (1D case)
Base class for coordinate mappings.
Tensor field of valence 0 (or component of a tensorial field).
virtual void set_etat_zero()
Sets the logical state to ETATZERO (zero).
Class intended to describe valence-2 symmetric tensors.
Symmetric tensors (with respect to two of their arguments).
Tensor field of valence 1.
int sym_index1() const
Number of the first symmetric index (0<= id_sym1 < valence )
int get_index_type(int i) const
Gives the type (covariant or contravariant) of the index number i .
int sym_index2() const
Number of the second symmetric index (id_sym1 < id_sym2 < valence )
int get_valence() const
Returns the valence.
virtual Itbl indices(int pos) const
Returns the indices of a component given by its position in the array cmp .
int get_n_comp() const
Returns the number of stored components.
const Base_vect * get_triad() const
Returns the vectorial basis (triad) on which the components are defined.
Scalar & set(const Itbl &ind)
Returns the value of a component (read/write version).