78 assert(ci.
get_etat() != ETATNONDEF) ;
88 resu =
sin(
double(1)) ;
98 co.va =
sin( ci.
va ) ;
110 assert(ci.
get_etat() != ETATNONDEF) ;
122 co =
cos(
double(1)) ;
128 co.va =
cos( ci.
va ) ;
142 assert(ci.
get_etat() != ETATNONDEF) ;
152 resu =
tan(
double(1)) ;
161 co.va =
tan( ci.
va ) ;
173 assert(ci.
get_etat() != ETATNONDEF) ;
205 assert(ci.
get_etat() != ETATNONDEF) ;
212 co.va = double(0.5) * M_PI ;
237 assert(ci.
get_etat() != ETATNONDEF) ;
269 assert(ci.
get_etat() != ETATNONDEF) ;
299 assert(ci.
get_etat() != ETATNONDEF) ;
329 assert(ci.
get_etat() != ETATNONDEF) ;
347 co.va =
exp( ci.
va ) ;
361 assert(ci.
get_etat() != ETATNONDEF) ;
391 assert(ci.
get_etat() != ETATNONDEF) ;
395 cout <<
"Argument of log is ZERO in log(Scalar) !" << endl ;
412 co.va =
log( ci.
va ) ;
424 assert(ci.
get_etat() != ETATNONDEF) ;
428 cout <<
"Argument of log10 is ZERO in log10(Scalar) !" << endl ;
457 assert(ci.
get_etat() != ETATNONDEF) ;
465 cout <<
"pow(Scalar, int) : ETATZERO^n with n <= 0 !" << endl ;
481 co.va =
pow(ci.
va, n) ;
493 assert(ci.
get_etat() != ETATNONDEF) ;
501 cout <<
"pow(Scalar, double) : ETATZERO^x with x <= 0 !" << endl ;
517 co.va =
pow(ci.
va, x) ;
529 assert(ci.
get_etat() != ETATNONDEF) ;
547 co.va =
abs( ci.
va ) ;
559 assert(ci.
get_etat() != ETATNONDEF) ;
588 assert(ci.
get_etat() != ETATNONDEF) ;
617 assert(ci.
get_etat() != ETATNONDEF) ;
631 resu =
max( ci.
va ) ;
645 assert(ci.
get_etat() != ETATNONDEF) ;
659 resu =
min( ci.
va ) ;
673 assert(ci.
get_etat() != ETATNONDEF) ;
701 assert(c1.
get_etat() != ETATNONDEF) ;
702 assert(c2.
get_etat() != ETATNONDEF) ;
712 assert(normdiff.
get_etat() == ETATQCQ) ;
713 assert(norme2.
get_etat() == ETATQCQ) ;
716 for (
int l=0; l<nz; l++) {
717 if ( norme2(l) ==
double(0) ) {
718 resu.
set(l) = normdiff(l) ;
721 resu.
set(l) = normdiff(l) / norme2(l) ;
736 assert(c1.
get_etat() != ETATNONDEF) ;
737 assert(c2.
get_etat() != ETATNONDEF) ;
748 assert(maxdiff.
get_etat() == ETATQCQ) ;
749 assert(max2.
get_etat() == ETATQCQ) ;
752 for (
int l=0; l<nz; l++) {
753 if ( max2(l) ==
double(0) ) {
754 resu.
set(l) = maxdiff(l) ;
757 resu.
set(l) = maxdiff(l) / max2(l) ;
Cmp log(const Cmp &)
Neperian logarithm.
Cmp asin(const Cmp &)
Arcsine.
Cmp exp(const Cmp &)
Exponential.
Cmp sqrt(const Cmp &)
Square root.
virtual void set_etat_zero()
Sets the logical state to ETATZERO (zero).
const Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
Cmp racine_cubique(const Cmp &)
Cube root.
double & set(int i)
Read/write of a particular element (index i) (1D case)
Tensor field of valence 0 (or component of a tensorial field).
int get_etat() const
Gives the logical state.
int get_etat() const
Returns the logical state ETATNONDEF (undefined), ETATZERO (null) or ETATQCQ (ordinary).
Tbl min(const Cmp &)
Minimum values of a Cmp in each domain.
virtual void set_etat_qcq()
Sets the logical state to ETATQCQ (ordinary state).
Cmp cos(const Cmp &)
Cosine.
Tbl diffrel(const Cmp &a, const Cmp &b)
Relative difference between two Cmp (norme version).
void set_etat_qcq()
Sets the logical state to ETATQCQ (ordinary state).
Cmp tan(const Cmp &)
Tangent.
Tbl norme(const Cmp &)
Sums of the absolute values of all the values of the Cmp in each domain.
Cmp atan(const Cmp &)
Arctangent.
void set_etat_one()
Sets the logical state to ETATUN (one).
int get_nzone() const
Returns the number of domains.
Valeur va
The numerical value of the Scalar.
Tbl max(const Cmp &)
Maximum values of a Cmp in each domain.
double totalmax(const Mtbl &)
Maximum value of the Mtbl elements in all domains.
Cmp pow(const Cmp &, int)
Power .
double totalmin(const Mtbl &)
Minimum value of the Mtbl elements in all domain.
Cmp log10(const Cmp &)
Basis 10 logarithm.
Cmp acos(const Cmp &)
Arccosine.
Cmp abs(const Cmp &)
Absolute value.
Cmp sin(const Cmp &)
Sine.
Mtbl Heaviside(const Mtbl &)
Heaviside function.
const Map & get_mp() const
Returns the mapping.
void annule_hard()
Sets the Tbl to zero in a hard way.
Tbl diffrelmax(const Cmp &a, const Cmp &b)
Relative difference between two Cmp (max version).