32char scalar_math_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Tensor/Scalar/scalar_math.C,v 1.6 2014/10/13 08:53:46 j_novak Exp $" ;
75 assert(
ci.get_etat() != ETATNONDEF) ;
78 if (
ci.get_etat() == ETATZERO) {
83 if (
ci.get_etat() == ETATUN) {
107 assert(
ci.get_etat() != ETATNONDEF) ;
112 if (
ci.get_etat() == ETATZERO) {
118 if (
ci.get_etat() == ETATUN) {
119 co =
cos(
double(1)) ;
139 assert(
ci.get_etat() != ETATNONDEF) ;
142 if (
ci.get_etat() == ETATZERO) {
147 if (
ci.get_etat() == ETATUN) {
170 assert(
ci.get_etat() != ETATNONDEF) ;
173 if (
ci.get_etat() == ETATZERO) {
178 if (
ci.get_etat() == ETATUN) {
202 assert(
ci.get_etat() != ETATNONDEF) ;
207 if (
ci.get_etat() == ETATZERO) {
213 if (
ci.get_etat() == ETATUN) {
234 assert(
ci.get_etat() != ETATNONDEF) ;
237 if (
ci.get_etat() == ETATZERO) {
242 if (
ci.get_etat() == ETATUN) {
266 assert(
ci.get_etat() != ETATNONDEF) ;
269 if (
ci.get_etat() == ETATZERO) {
274 if (
ci.get_etat() == ETATUN) {
296 assert(
ci.get_etat() != ETATNONDEF) ;
299 if (
ci.get_etat() == ETATZERO) {
304 if (
ci.get_etat() == ETATUN) {
326 assert(
ci.get_etat() != ETATNONDEF) ;
331 if (
ci.get_etat() == ETATZERO) {
336 if (
ci.get_etat() == ETATUN) {
358 assert(
ci.get_etat() != ETATNONDEF) ;
363 if (
ci.get_etat() == ETATZERO) {
368 if (
ci.get_etat() == ETATUN) {
388 assert(
ci.get_etat() != ETATNONDEF) ;
391 if (
ci.get_etat() == ETATZERO) {
392 cout <<
"Argument of log is ZERO in log(Scalar) !" <<
endl ;
397 if (
ci.get_etat() == ETATUN) {
399 resu.set_etat_zero() ;
421 assert(
ci.get_etat() != ETATNONDEF) ;
424 if (
ci.get_etat() == ETATZERO) {
425 cout <<
"Argument of log10 is ZERO in log10(Scalar) !" <<
endl ;
430 if (
ci.get_etat() == ETATUN) {
432 resu.set_etat_zero() ;
454 assert(
ci.get_etat() != ETATNONDEF) ;
457 if (
ci.get_etat() == ETATZERO) {
462 cout <<
"pow(Scalar, int) : ETATZERO^n with n <= 0 !" <<
endl ;
468 if (
ci.get_etat() == ETATUN) {
490 assert(
ci.get_etat() != ETATNONDEF) ;
493 if (
ci.get_etat() == ETATZERO) {
498 cout <<
"pow(Scalar, double) : ETATZERO^x with x <= 0 !" <<
endl ;
504 if (
ci.get_etat() == ETATUN) {
526 assert(
ci.get_etat() != ETATNONDEF) ;
529 if (
ci.get_etat() == ETATZERO) {
534 if (
ci.get_etat() == ETATUN) {
556 assert(
ci.get_etat() != ETATNONDEF) ;
561 if (
ci.get_etat() == ETATZERO) {
565 if (
ci.get_etat() == ETATUN) {
585 assert(
ci.get_etat() != ETATNONDEF) ;
590 if (
ci.get_etat() == ETATZERO) {
594 if (
ci.get_etat() == ETATUN) {
614 assert(
ci.get_etat() != ETATNONDEF) ;
616 Tbl resu(
ci.get_mp().get_mg()->get_nzone() ) ;
618 if (
ci.get_etat() == ETATZERO) {
622 if (
ci.get_etat() == ETATUN) {
642 assert(
ci.get_etat() != ETATNONDEF) ;
644 Tbl resu(
ci.get_mp().get_mg()->get_nzone() ) ;
646 if (
ci.get_etat() == ETATZERO) {
650 if (
ci.get_etat() == ETATUN) {
670 assert(
ci.get_etat() != ETATNONDEF) ;
672 Tbl resu(
ci.get_mp().get_mg()->get_nzone() ) ;
674 if (
ci.get_etat() == ETATZERO) {
678 if (
ci.get_etat() == ETATUN) {
698 assert(
c1.get_etat() != ETATNONDEF) ;
699 assert(
c2.get_etat() != ETATNONDEF) ;
701 int nz =
c1.get_mp().get_mg()->get_nzone() ;
712 resu.set_etat_qcq() ;
713 for (
int l=0;
l<nz;
l++) {
714 if (
norme2(
l) ==
double(0) ) {
733 assert(
c1.get_etat() != ETATNONDEF) ;
734 assert(
c2.get_etat() != ETATNONDEF) ;
736 int nz =
c1.get_mp().get_mg()->get_nzone() ;
748 resu.set_etat_qcq() ;
749 for (
int l=0;
l<nz;
l++) {
750 if (
max2(
l) ==
double(0) ) {
Time evolution with partial storage (*** under development ***).
Tensor field of valence 0 (or component of a tensorial field).
Cmp atan(const Cmp &)
Arctangent.
Cmp sqrt(const Cmp &)
Square root.
Cmp log10(const Cmp &)
Basis 10 logarithm.
Cmp exp(const Cmp &)
Exponential.
Cmp sin(const Cmp &)
Sine.
Tbl diffrel(const Cmp &a, const Cmp &b)
Relative difference between two Cmp (norme version).
Tbl norme(const Cmp &)
Sums of the absolute values of all the values of the Cmp in each domain.
Cmp acos(const Cmp &)
Arccosine.
Cmp asin(const Cmp &)
Arcsine.
Cmp racine_cubique(const Cmp &)
Cube root.
Tbl min(const Cmp &)
Minimum values of a Cmp in each domain.
Tbl max(const Cmp &)
Maximum values of a Cmp in each domain.
Cmp pow(const Cmp &, int)
Power .
Cmp cos(const Cmp &)
Cosine.
Cmp abs(const Cmp &)
Absolute value.
Tbl diffrelmax(const Cmp &a, const Cmp &b)
Relative difference between two Cmp (max version).
Cmp tan(const Cmp &)
Tangent.
Cmp log(const Cmp &)
Neperian logarithm.
Mtbl Heaviside(const Mtbl &)
Heaviside function.
double totalmin(const Mtbl &)
Minimum value of the Mtbl elements in all domain.
double totalmax(const Mtbl &)
Maximum value of the Mtbl elements in all domains.