30char TBL_VAL_MATH_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Valencia/tbl_val_math.C,v 1.4 2014/10/13 08:53:48 j_novak Exp $" ;
70 assert(
ti.get_etat() != ETATNONDEF) ;
73 if (
ti.get_etat() == ETATZERO) {
81 int taille =
ti.get_taille() ;
82 for (
int i=0 ;
i<taille ;
i++) {
85 for (
int i=0;
i<
ti.get_taille_i(0);
i++)
87 if (
ti.txti != 0x0)
for (
int i=0;
i<
ti.get_taille_i(1);
i++)
89 if (
ti.typi != 0x0)
for (
int i=0;
i<
ti.get_taille_i(2);
i++)
101 assert(
ti.get_etat() != ETATNONDEF) ;
107 if (
ti.get_etat() == ETATZERO) {
114 int taille =
ti.get_taille() ;
115 for (
int i=0 ;
i<taille ;
i++) {
118 for (
int i=0;
i<
ti.get_taille_i(0);
i++)
120 if (
ti.txti != 0x0)
for (
int i=0;
i<
ti.get_taille_i(1);
i++)
122 if (
ti.typi != 0x0)
for (
int i=0;
i<
ti.get_taille_i(2);
i++)
134 assert(
ti.get_etat() != ETATNONDEF) ;
137 if (
ti.get_etat() == ETATZERO) {
145 int taille =
ti.get_taille() ;
146 for (
int i=0 ;
i<taille ;
i++) {
149 for (
int i=0;
i<
ti.get_taille_i(0);
i++)
151 if (
ti.txti != 0x0)
for (
int i=0;
i<
ti.get_taille_i(1);
i++)
153 if (
ti.typi != 0x0)
for (
int i=0;
i<
ti.get_taille_i(2);
i++)
165 assert(
ti.get_etat() != ETATNONDEF) ;
168 if (
ti.get_etat() == ETATZERO) {
176 int taille =
ti.get_taille() ;
177 for (
int i=0 ;
i<taille ;
i++) {
180 for (
int i=0;
i<
ti.get_taille_i(0);
i++)
182 if (
ti.txti != 0x0)
for (
int i=0;
i<
ti.get_taille_i(1);
i++)
184 if (
ti.typi != 0x0)
for (
int i=0;
i<
ti.get_taille_i(2);
i++)
197 assert(
ti.get_etat() != ETATNONDEF) ;
203 if (
ti.get_etat() == ETATZERO) {
210 int taille =
ti.get_taille() ;
211 for (
int i=0 ;
i<taille ;
i++) {
214 for (
int i=0;
i<
ti.get_taille_i(0);
i++)
216 if (
ti.txti != 0x0)
for (
int i=0;
i<
ti.get_taille_i(1);
i++)
218 if (
ti.typi != 0x0)
for (
int i=0;
i<
ti.get_taille_i(2);
i++)
230 assert(
ti.get_etat() != ETATNONDEF) ;
233 if (
ti.get_etat() == ETATZERO) {
241 int taille =
ti.get_taille() ;
242 for (
int i=0 ;
i<taille ;
i++) {
245 for (
int i=0;
i<
ti.get_taille_i(0);
i++)
247 if (
ti.txti != 0x0)
for (
int i=0;
i<
ti.get_taille_i(1);
i++)
249 if (
ti.typi != 0x0)
for (
int i=0;
i<
ti.get_taille_i(2);
i++)
261 assert(
ti.get_etat() != ETATNONDEF) ;
264 if (
ti.get_etat() == ETATZERO) {
272 int taille =
ti.get_taille() ;
273 for (
int i=0 ;
i<taille ;
i++) {
276 for (
int i=0;
i<
ti.get_taille_i(0);
i++)
278 if (
ti.txti != 0x0)
for (
int i=0;
i<
ti.get_taille_i(1);
i++)
280 if (
ti.typi != 0x0)
for (
int i=0;
i<
ti.get_taille_i(2);
i++)
292 assert(
ti.get_etat() != ETATNONDEF) ;
298 if (
ti.get_etat() == ETATZERO) {
305 int taille =
ti.get_taille() ;
306 for (
int i=0 ;
i<taille ;
i++) {
309 for (
int i=0;
i<
ti.get_taille_i(0);
i++)
311 if (
ti.txti != 0x0)
for (
int i=0;
i<
ti.get_taille_i(1);
i++)
313 if (
ti.typi != 0x0)
for (
int i=0;
i<
ti.get_taille_i(2);
i++)
325 assert(
ti.get_etat() != ETATNONDEF) ;
328 if (
ti.get_etat() == ETATZERO) {
329 cout <<
"Tbl_val log: log(ETATZERO) !" <<
endl ;
337 int taille =
ti.get_taille() ;
338 for (
int i=0 ;
i<taille ;
i++) {
341 for (
int i=0;
i<
ti.get_taille_i(0);
i++)
343 if (
ti.txti != 0x0)
for (
int i=0;
i<
ti.get_taille_i(1);
i++)
345 if (
ti.typi != 0x0)
for (
int i=0;
i<
ti.get_taille_i(2);
i++)
357 assert(
ti.get_etat() != ETATNONDEF) ;
360 if (
ti.get_etat() == ETATZERO) {
361 cout <<
"Tbl_val log10: log10(ETATZERO) !" <<
endl ;
369 int taille =
ti.get_taille() ;
370 for (
int i=0 ;
i<taille ;
i++) {
373 for (
int i=0;
i<
ti.get_taille_i(0);
i++)
375 if (
ti.txti != 0x0)
for (
int i=0;
i<
ti.get_taille_i(1);
i++)
377 if (
ti.typi != 0x0)
for (
int i=0;
i<
ti.get_taille_i(2);
i++)
389 assert(
ti.get_etat() != ETATNONDEF) ;
392 if (
ti.get_etat() == ETATZERO) {
397 cout <<
"Tbl_val pow: ETATZERO^n avec n<=0 ! "<<
endl ;
407 int taille =
ti.get_taille() ;
408 for (
int i=0 ;
i<taille ;
i++) {
411 for (
int i=0;
i<
ti.get_taille_i(0);
i++)
413 if (
ti.txti != 0x0)
for (
int i=0;
i<
ti.get_taille_i(1);
i++)
415 if (
ti.typi != 0x0)
for (
int i=0;
i<
ti.get_taille_i(2);
i++)
427 assert(
ti.get_etat() != ETATNONDEF) ;
430 if (
ti.get_etat() == ETATZERO) {
435 cout <<
"Tbl_val pow: ETATZERO^x avec x<=0 !" <<
endl ;
444 int taille =
ti.get_taille() ;
445 for (
int i=0 ;
i<taille ;
i++) {
448 for (
int i=0;
i<
ti.get_taille_i(0);
i++)
450 if (
ti.txti != 0x0)
for (
int i=0;
i<
ti.get_taille_i(1);
i++)
452 if (
ti.typi != 0x0)
for (
int i=0;
i<
ti.get_taille_i(2);
i++)
464 assert(
ti.get_etat() != ETATNONDEF) ;
467 if (
ti.get_etat() == ETATZERO) {
477 const double* xi =
ti.t ;
479 int taille =
ti.get_taille() ;
481 for (
int i=0 ;
i<taille ;
i++) {
484 for (
int i=0;
i<
ti.get_taille_i(0);
i++)
486 if (
ti.txti != 0x0)
for (
int i=0;
i<
ti.get_taille_i(1);
i++)
488 if (
ti.typi != 0x0)
for (
int i=0;
i<
ti.get_taille_i(2);
i++)
500 assert(
ti.get_etat() != ETATNONDEF) ;
503 if (
ti.get_etat() == ETATZERO) {
513 for (
int i=0 ;
i<
ti.get_taille() ;
i++)
516 for (
int i=0;
i<
ti.get_taille_i(0);
i++)
517 if (
ti.tzri[
i] < 0)
res.tzri[
i] *= -1 ;
518 if (
ti.txti != 0x0)
for (
int i=0;
i<
ti.get_taille_i(1);
i++)
519 if (
ti.txti[
i] < 0)
res.txti[
i] *= -1 ;
520 if (
ti.typi != 0x0)
for (
int i=0;
i<
ti.get_taille_i(2);
i++)
521 if (
ti.typi[
i] < 0)
res.typi[
i] *= -1 ;
533 assert(
ti.get_etat() != ETATNONDEF) ;
536 if (
ti.get_etat() == ETATZERO) {
543 const double* x =
ti.t ;
545 for (
int i=1;
i<
ti.get_taille();
i++) {
559 assert(
ti.get_etat() != ETATNONDEF) ;
562 if (
ti.get_etat() == ETATZERO) {
569 const double* x =
ti.t ;
571 for (
int i=1;
i<
ti.get_taille();
i++) {
585 assert(
ti.get_etat() != ETATNONDEF) ;
589 if (
ti.get_etat() != ETATZERO) {
592 const double* x =
ti.t ;
593 for (
int i=0;
i<
ti.get_taille();
i++) {
609 assert(
t1.get_etat() != ETATNONDEF) ;
610 assert(
t2.get_etat() != ETATNONDEF) ;
615 if (
norm2 ==
double(0) ) {
633 assert(
t1.get_etat() != ETATNONDEF) ;
634 assert(
t2.get_etat() != ETATNONDEF) ;
639 if (
max2 ==
double(0) ) {
Time evolution with partial storage (*** under development ***).
Finite-difference array intended to store field values.
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.