30char tbl_math_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Tbl/tbl_math.C,v 1.4 2014/10/13 08:53:41 j_novak Exp $" ;
81 assert(
ti.get_etat() != ETATNONDEF) ;
84 if (
ti.get_etat() == ETATZERO) {
92 int taille =
ti.get_taille() ;
93 for (
int i=0 ;
i<taille ;
i++) {
106 assert(
ti.get_etat() != ETATNONDEF) ;
112 if (
ti.get_etat() == ETATZERO) {
113 int taille =
ti.get_taille() ;
114 for (
int i=0 ;
i<taille ;
i++) {
122 int taille =
ti.get_taille() ;
123 for (
int i=0 ;
i<taille ;
i++) {
136 assert(
ti.get_etat() != ETATNONDEF) ;
139 if (
ti.get_etat() == ETATZERO) {
147 int taille =
ti.get_taille() ;
148 for (
int i=0 ;
i<taille ;
i++) {
161 assert(
ti.get_etat() != ETATNONDEF) ;
164 if (
ti.get_etat() == ETATZERO) {
172 int taille =
ti.get_taille() ;
173 for (
int i=0 ;
i<taille ;
i++) {
186 assert(
ti.get_etat() != ETATNONDEF) ;
192 if (
ti.get_etat() == ETATZERO) {
193 int taille =
ti.get_taille() ;
194 for (
int i=0 ;
i<taille ;
i++) {
202 int taille =
ti.get_taille() ;
203 for (
int i=0 ;
i<taille ;
i++) {
216 assert(
ti.get_etat() != ETATNONDEF) ;
219 if (
ti.get_etat() == ETATZERO) {
227 int taille =
ti.get_taille() ;
228 for (
int i=0 ;
i<taille ;
i++) {
241 assert(
ti.get_etat() != ETATNONDEF) ;
244 if (
ti.get_etat() == ETATZERO) {
252 int taille =
ti.get_taille() ;
253 for (
int i=0 ;
i<taille ;
i++) {
266 assert(
ti.get_etat() != ETATNONDEF) ;
272 if (
ti.get_etat() == ETATZERO) {
273 int taille =
ti.get_taille() ;
274 for (
int i=0 ;
i<taille ;
i++) {
282 int taille =
ti.get_taille() ;
283 for (
int i=0 ;
i<taille ;
i++) {
296 assert(
ti.get_etat() != ETATNONDEF) ;
302 if (
ti.get_etat() == ETATZERO) {
303 int taille =
ti.get_taille() ;
304 for (
int i=0 ;
i<taille ;
i++) {
312 int taille =
ti.get_taille() ;
313 for (
int i=0 ;
i<taille ;
i++) {
329 assert(
ti.get_etat() != ETATNONDEF) ;
332 if (
ti.get_etat() == ETATZERO) {
333 cout <<
"Tbl log: log(ETATZERO) !" <<
endl ;
341 int taille =
ti.get_taille() ;
342 for (
int i=0 ;
i<taille ;
i++) {
355 assert(
ti.get_etat() != ETATNONDEF) ;
358 if (
ti.get_etat() == ETATZERO) {
359 cout <<
"Tbl log10: log10(ETATZERO) !" <<
endl ;
367 int taille =
ti.get_taille() ;
368 for (
int i=0 ;
i<taille ;
i++) {
381 assert(
ti.get_etat() != ETATNONDEF) ;
384 if (
ti.get_etat() == ETATZERO) {
389 cout <<
"Tbl pow: ETATZERO^n avec n<=0 ! "<<
endl ;
399 int taille =
ti.get_taille() ;
400 for (
int i=0 ;
i<taille ;
i++) {
413 assert(
ti.get_etat() != ETATNONDEF) ;
416 if (
ti.get_etat() == ETATZERO) {
421 cout <<
"Tbl pow: ETATZERO^x avec x<=0 !" <<
endl ;
430 int taille =
ti.get_taille() ;
431 for (
int i=0 ;
i<taille ;
i++) {
444 assert(
ti.get_etat() != ETATNONDEF) ;
447 if (
ti.get_etat() == ETATZERO) {
457 const double* xi =
ti.t ;
459 int taille =
ti.get_taille() ;
461 for (
int i=0 ;
i<taille ;
i++) {
474 assert(
ti.get_etat() != ETATNONDEF) ;
477 if (
ti.get_etat() == ETATZERO) {
487 for (
int i=0 ;
i<
ti.get_taille() ;
i++)
501 assert(
ti.get_etat() != ETATNONDEF) ;
504 if (
ti.get_etat() == ETATZERO) {
511 const double* x =
ti.t ;
513 for (
int i=1;
i<
ti.get_taille();
i++) {
527 assert(
ti.get_etat() != ETATNONDEF) ;
530 if (
ti.get_etat() == ETATZERO) {
537 const double* x =
ti.t ;
539 for (
int i=1;
i<
ti.get_taille();
i++) {
553 assert(
ti.get_etat() != ETATNONDEF) ;
557 if (
ti.get_etat() != ETATZERO) {
560 const double* x =
ti.t ;
561 for (
int i=0;
i<
ti.get_taille();
i++) {
577 assert(
t1.get_etat() != ETATNONDEF) ;
578 assert(
t2.get_etat() != ETATNONDEF) ;
583 if (
norm2 ==
double(0) ) {
601 assert(
t1.get_etat() != ETATNONDEF) ;
602 assert(
t2.get_etat() != ETATNONDEF) ;
607 if (
max2 ==
double(0) ) {
Time evolution with partial storage (*** under development ***).
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.