-- FreeFem++ v4.6 (Thu Apr 2 15:47:38 CEST 2020 - git v4.6) Load: lg_fem lg_mesh lg_mesh3 eigenvalue 1 : // Discussion: 2 : // 3 : // Examine the math functions: 4 : // 5 : // abs 6 : // acos 7 : // acosh 8 : // asin 9 : // asinh 10 : // atan 11 : // atan2 12 : // atanh 13 : // ceil 14 : // cos 15 : // cosh 16 : // erf 17 : // erfc 18 : // exp 19 : // floor 20 : // imag 21 : // j0 22 : // j1 23 : // jn 24 : // lgamma 25 : // log 26 : // log10 27 : // max 28 : // min 29 : // pow 30 : // rint 31 : // sin 32 : // sinh 33 : // sqrt 34 : // tan 35 : // tanh 36 : // tgamma 37 : // y0 38 : // y1 39 : // yn 40 : // 41 : // Location: 42 : // 43 : // http://people.sc.fsu.edu/~jburkardt/freefem_src/math_test/math_test.edp 44 : // 45 : // Licensing: 46 : // 47 : // This code is distributed under the GNU LGPL license. 48 : // 49 : // Modified: 50 : // 51 : // 15 February 2015 52 : // 53 : // Reference: 54 : // 55 : // Frederic Hecht, 56 : // Freefem++, 57 : // Third Edition, version 3.22 58 : // 59 : cout << "\n"; 60 : cout << "math_test:\n"; 61 : cout << " Try out the FreeFem++ Math Library Functions.\n ... : "; 62 : // 63 : // Initialize the random number generator. 64 : // 65 : int seed; 66 : // 67 : // Y = ABS ( X ); 68 : // 69 : seed = 123456789; 70 : randinit ( seed ); 71 : cout << "\n"; 72 : cout << " Y = ABS ( X ):\n"; 73 : cout << "\n"; 74 : for ( int i = 0; i < 10; i++ ) 75 : { 76 : x = 10.0 * randreal1 ( ) - 5.0; 77 : y = abs ( x ); 78 : cout << " " << x << " " << y << "\n"; 79 : } 80 : // 81 : // Y = ACOS ( X ); 82 : // X <= 1.0 required. 83 : // 84 : seed = 123456789; 85 : randinit ( seed ); 86 : cout << "\n"; 87 : cout << " Y = ACOS ( X ):\n"; 88 : cout << "\n"; 89 : for ( int i = 0; i < 10; i++ ) 90 : { 91 : x = 2.0 * randreal1 ( ) - 1.0; 92 : y = acos ( x ); 93 : cout << " " << x << " " << y << "\n"; 94 : } 95 : // 96 : // Y = ACOSH ( X ); 97 : // 1 <= |X| required. 98 : // 99 : seed = 123456789; 100 : randinit ( seed ); 101 : cout << "\n"; 102 : cout << " Y = ACOSH ( X ):\n"; 103 : cout << "\n"; 104 : for ( int i = 0; i < 10; i++ ) 105 : { 106 : x = randreal1 ( ); 107 : if ( x < 0.0 ) 108 : { 109 : x = 10.0 * x - 1.0; 110 : } 111 : else 112 : { 113 : x = 10.0 * x + 1.0; 114 : } 115 : y = acosh ( x ); 116 : cout << " " << x << " " << y << "\n"; 117 : } 118 : // 119 : // Y = ASIN ( X ); 120 : // |X| <= 1. 121 : // 122 : seed = 123456789; 123 : randinit ( seed ); 124 : cout << "\n"; 125 : cout << " Y = ASIN ( X ):\n"; 126 : cout << "\n"; 127 : for ( int i = 0; i < 10; i++ ) 128 : { 129 : x = 2.0 * randreal1 ( ) - 1.0; 130 : y = asin ( x ); 131 : cout << " " << x << " " << y << "\n"; 132 : } 133 : // 134 : // Y = ASINH ( X ); 135 : // 136 : seed = 123456789; 137 : randinit ( seed ); 138 : cout << "\n"; 139 : cout << " Y = ASINH ( X ):\n"; 140 : cout << "\n"; 141 : for ( int i = 0; i < 10; i++ ) 142 : { 143 : x = 10.0 * randreal1 ( ) - 5.0; 144 : y = asinh ( x ); 145 : cout << " " << x << " " << y << "\n"; 146 : } 147 : // 148 : // Z = ATAN ( Y / X ); 149 : // 150 : seed = 123456789; 151 : randinit ( seed ); 152 : cout << "\n"; 153 : cout << " Z = ATAN ( Y / X ):\n"; 154 : cout << "\n"; 155 : for ( int i = 0; i < 10; i++ ) 156 : { 157 : x = 10.0 * randreal1 ( ) - 5.0; 158 : y = 10.0 * randreal1 ( ) - 5.0; 159 : z = atan ( y / x ); 160 : cout << " " << x << " " << y << " " << z << "\n"; 161 : } 162 : // 163 : // Z = ATAN2 ( Y, X ); 164 : // 165 : seed = 123456789; 166 : randinit ( seed ); 167 : cout << "\n"; 168 : cout << " Z = ATAN2 ( Y, X ):\n"; 169 : cout << "\n"; 170 : for ( int i = 0; i < 10; i++ ) 171 : { 172 : x = 10.0 * randreal1 ( ) - 5.0; 173 : y = 10.0 * randreal1 ( ) - 5.0; 174 : z = atan2 ( y, x ); 175 : cout << " " << x << " " << y << " " << z << "\n"; 176 : } 177 : // 178 : // Y = ATANH ( X ); 179 : // |X| < 1 required. 180 : // 181 : seed = 123456789; 182 : randinit ( seed ); 183 : cout << "\n"; 184 : cout << " Y = ATANH ( X ):\n"; 185 : cout << "\n"; 186 : for ( int i = 0; i < 10; i++ ) 187 : { 188 : x = 2.0 * randreal1 ( ) - 1.0; 189 : y = atanh ( x ); 190 : cout << " " << x << " " << y << "\n"; 191 : } 192 : // 193 : // Y = CEIL ( X ); 194 : // 195 : seed = 123456789; 196 : randinit ( seed ); 197 : cout << "\n"; 198 : cout << " Y = CEIL ( X ):\n"; 199 : cout << "\n"; 200 : for ( int i = 0; i < 10; i++ ) 201 : { 202 : x = 10.0 * randreal1 ( ) - 5.0; 203 : y = ceil ( x ); 204 : cout << " " << x << " " << y << "\n"; 205 : } 206 : // 207 : // Y = COS ( X ); 208 : // 209 : seed = 123456789; 210 : randinit ( seed ); 211 : cout << "\n"; 212 : cout << " Y = COS ( X ):\n"; 213 : cout << "\n"; 214 : for ( int i = 0; i < 10; i++ ) 215 : { 216 : x = 10.0 * randreal1 ( ) - 5.0; 217 : y = cos ( x ); 218 : cout << " " << x << " " << y << "\n"; 219 : } 220 : // 221 : // Y = COSH ( X ); 222 : // 223 : seed = 123456789; 224 : randinit ( seed ); 225 : cout << "\n"; 226 : cout << " Y = COSH ( X ):\n"; 227 : cout << "\n"; 228 : for ( int i = 0; i < 10; i++ ) 229 : { 230 : x = 10.0 * randreal1 ( ) - 5.0; 231 : y = cosh ( x ); 232 : cout << " " << x << " " << y << "\n"; 233 : } 234 : // 235 : // Y = ERF ( X ); 236 : // 237 : seed = 123456789; 238 : randinit ( seed ); 239 : cout << "\n"; 240 : cout << " Y = ERF ( X ):\n"; 241 : cout << "\n"; 242 : for ( int i = 0; i < 10; i++ ) 243 : { 244 : x = 10.0 * randreal1 ( ) - 5.0; 245 : y = erf ( x ); 246 : cout << " " << x << " " << y << "\n"; 247 : } 248 : // 249 : // Y = ERFC ( X ); 250 : // 251 : seed = 123456789; 252 : randinit ( seed ); 253 : cout << "\n"; 254 : cout << " Y = ERFC ( X ):\n"; 255 : cout << "\n"; 256 : for ( int i = 0; i < 10; i++ ) 257 : { 258 : x = 10.0 * randreal1 ( ) - 5.0; 259 : y = erfc ( x ); 260 : cout << " " << x << " " << y << "\n"; 261 : } 262 : // 263 : // Y = EXP ( X ); 264 : // 265 : seed = 123456789; 266 : randinit ( seed ); 267 : cout << "\n"; 268 : cout << " Y = EXP ( X ):\n"; 269 : cout << "\n"; 270 : for ( int i = 0; i < 10; i++ ) 271 : { 272 : x = 10.0 * randreal1 ( ) - 5.0; 273 : y = exp ( x ); 274 : cout << " " << x << " " << y << "\n"; 275 : } 276 : // 277 : // Y = FLOOR ( X ); 278 : // 279 : seed = 123456789; 280 : randinit ( seed ); 281 : cout << "\n"; 282 : cout << " Y = FLOOR ( X ):\n"; 283 : cout << "\n"; 284 : for ( int i = 0; i < 10; i++ ) 285 : { 286 : x = 10.0 * randreal1 ( ) - 5.0; 287 : y = floor ( x ); 288 : cout << " " << x << " " << y << "\n"; 289 : } 290 : // 291 : // Y = J0 ( X ); 292 : // 293 : seed = 123456789; 294 : randinit ( seed ); 295 : cout << "\n"; 296 : cout << " Y = J0 ( X ):\n"; 297 : cout << "\n"; 298 : for ( int i = 0; i < 10; i++ ) 299 : { 300 : x = 10.0 * randreal1 ( ) - 5.0; 301 : y = j0 ( x ); 302 : cout << " " << x << " " << y << "\n"; 303 : } 304 : // 305 : // Y = J1 ( X ); 306 : // 307 : seed = 123456789; 308 : randinit ( seed ); 309 : cout << "\n"; 310 : cout << " Y = J1 ( X ):\n"; 311 : cout << "\n"; 312 : for ( int i = 0; i < 10; i++ ) 313 : { 314 : x = 10.0 * randreal1 ( ) - 5.0; 315 : y = j1 ( x ); 316 : cout << " " << x << " " << y << "\n"; 317 : } 318 : // 319 : // Y = JN ( 2, X ); 320 : // 321 : seed = 123456789; 322 : randinit ( seed ); 323 : cout << "\n"; 324 : cout << " Y = JN ( 2, X ):\n"; 325 : cout << "\n"; 326 : for ( int i = 0; i < 10; i++ ) 327 : { 328 : x = 10.0 * randreal1 ( ) - 5.0; 329 : y = jn ( 2, x ); 330 : cout << " " << x << " " << y << "\n"; 331 : } 332 : // 333 : // Y = LGAMMA ( X ); 334 : // X must not be 0 or a negative integer. 335 : // 336 : seed = 123456789; 337 : randinit ( seed ); 338 : cout << "\n"; 339 : cout << " Y = LGAMMA ( X ):\n"; 340 : cout << "\n"; 341 : for ( int i = 0; i < 10; i++ ) 342 : { 343 : x = 10.0 * randreal1 ( ); 344 : y = lgamma ( x ); 345 : cout << " " << x << " " << y << "\n"; 346 : } 347 : // 348 : // Y = LOG ( X ); 349 : // 0 < X required. 350 : // 351 : seed = 123456789; 352 : randinit ( seed ); 353 : cout << "\n"; 354 : cout << " Y = LOG ( X ):\n"; 355 : cout << "\n"; 356 : for ( int i = 0; i < 10; i++ ) 357 : { 358 : x = 10.0 * randreal1 ( ); 359 : y = log ( x ); 360 : cout << " " << x << " " << y << "\n"; 361 : } 362 : // 363 : // Y = LOG10 ( X ); 364 : // 0 < X required. 365 : // 366 : seed = 123456789; 367 : randinit ( seed ); 368 : cout << "\n"; 369 : cout << " Y = LOG10 ( X ):\n"; 370 : cout << "\n"; 371 : for ( int i = 0; i < 10; i++ ) 372 : { 373 : x = 10.0 * randreal1 ( ); 374 : y = log10 ( x ); 375 : cout << " " << x << " " << y << "\n"; 376 : } 377 : // 378 : // Z = MAX ( X, Y ); 379 : // 380 : seed = 123456789; 381 : randinit ( seed ); 382 : cout << "\n"; 383 : cout << " Z = MAX ( X, Y ):\n"; 384 : cout << "\n"; 385 : for ( int i = 0; i < 10; i++ ) 386 : { 387 : x = 10.0 * randreal1 ( ) - 5.0; 388 : y = 10.0 * randreal1 ( ) - 5.0; 389 : z = max ( x, y ); 390 : cout << " " << x << " " << y << " " << z << "\n"; 391 : } 392 : // 393 : // Z = MIN ( X, Y ); 394 : // 395 : seed = 123456789; 396 : randinit ( seed ); 397 : cout << "\n"; 398 : cout << " Z = MIN ( X, Y ):\n"; 399 : cout << "\n"; 400 : for ( int i = 0; i < 10; i++ ) 401 : { 402 : x = 10.0 * randreal1 ( ) - 5.0; 403 : y = 10.0 * randreal1 ( ) - 5.0; 404 : z = min ( x, y ); 405 : cout << " " << x << " " << y << " " << z << "\n"; 406 : } 407 : // 408 : // Z = POW ( X, Y ); 409 : // If X negative, Y must be an integer. 410 : // If X 0, Y must be positive. 411 : // 412 : seed = 123456789; 413 : randinit ( seed ); 414 : cout << "\n"; 415 : cout << " Z = POW ( X, Y ):\n"; 416 : cout << "\n"; 417 : for ( int i = 0; i < 10; i++ ) 418 : { 419 : x = 10.0 * randreal1 ( ); 420 : y = 10.0 * randreal1 ( ) - 5.0; 421 : z = pow ( x, y ); 422 : cout << " " << x << " " << y << " " << z << "\n"; 423 : } 424 : // 425 : // Y = RINT ( X ); 426 : // 427 : seed = 123456789; 428 : randinit ( seed ); 429 : cout << "\n"; 430 : cout << " Y = RINT ( X ):\n"; 431 : cout << "\n"; 432 : for ( int i = 0; i < 10; i++ ) 433 : { 434 : x = 10.0 * randreal1 ( ) - 5.0; 435 : y = rint ( x ); 436 : cout << " " << x << " " << y << "\n"; 437 : } 438 : // 439 : // Y = SIN ( X ); 440 : // 441 : seed = 123456789; 442 : randinit ( seed ); 443 : cout << "\n"; 444 : cout << " Y = SIN ( X ):\n"; 445 : cout << "\n"; 446 : for ( int i = 0; i < 10; i++ ) 447 : { 448 : x = 10.0 * randreal1 ( ) - 5.0; 449 : y = sin ( x ); 450 : cout << " " << x << " " << y << "\n"; 451 : } 452 : // 453 : // Y = SINH ( X ); 454 : // 455 : seed = 123456789; 456 : randinit ( seed ); 457 : cout << "\n"; 458 : cout << " Y = SINH ( X ):\n"; 459 : cout << "\n"; 460 : for ( int i = 0; i < 10; i++ ) 461 : { 462 : x = 10.0 * randreal1 ( ) - 5.0; 463 : y = sinh ( x ); 464 : cout << " " << x << " " << y << "\n"; 465 : } 466 : // 467 : // Y = SQRT ( X ); 468 : // 0 <= X required. 469 : // 470 : seed = 123456789; 471 : randinit ( seed ); 472 : cout << "\n"; 473 : cout << " Y = SQRT ( X ):\n"; 474 : cout << "\n"; 475 : for ( int i = 0; i < 10; i++ ) 476 : { 477 : x = 10.0 * randreal1 ( ); 478 : y = sqrt ( x ); 479 : cout << " " << x << " " << y << "\n"; 480 : } 481 : // 482 : // Y = TAN ( X ); 483 : // 484 : seed = 123456789; 485 : randinit ( seed ); 486 : cout << "\n"; 487 : cout << " Y = TAN ( X ):\n"; 488 : cout << "\n"; 489 : for ( int i = 0; i < 10; i++ ) 490 : { 491 : x = 10.0 * randreal1 ( ) - 5.0; 492 : y = tan ( x ); 493 : cout << " " << x << " " << y << "\n"; 494 : } 495 : // 496 : // Y = TANH ( X ); 497 : // 498 : seed = 123456789; 499 : randinit ( seed ); 500 : cout << "\n"; 501 : cout << " Y = TANH ( X ):\n"; 502 : cout << "\n"; 503 : for ( int i = 0; i < 10; i++ ) 504 : { 505 : x = 10.0 * randreal1 ( ) - 5.0; 506 : y = tanh ( x ); 507 : cout << " " << x << " " << y << "\n"; 508 : } 509 : // 510 : // Y = TGAMMA ( X ); 511 : // X must not be 0 or a negative integer. 512 : // 513 : seed = 123456789; 514 : randinit ( seed ); 515 : cout << "\n"; 516 : cout << " Y = TGAMMA ( X ):\n"; 517 : cout << "\n"; 518 : for ( int i = 0; i < 10; i++ ) 519 : { 520 : x = 10.0 * randreal1 ( ) - 5.0; 521 : y = tgamma ( x ); 522 : cout << " " << x << " " << y << "\n"; 523 : } 524 : // 525 : // Y = Y0 ( X ); 526 : // 527 : seed = 123456789; 528 : randinit ( seed ); 529 : cout << "\n"; 530 : cout << " Y = Y0 ( X ):\n"; 531 : cout << "\n"; 532 : for ( int i = 0; i < 10; i++ ) 533 : { 534 : x = 10.0 * randreal1 ( ) - 5.0; 535 : y = y0 ( x ); 536 : cout << " " << x << " " << y << "\n"; 537 : } 538 : // 539 : // Y = Y1 ( X ); 540 : // 541 : seed = 123456789; 542 : randinit ( seed ); 543 : cout << "\n"; 544 : cout << " Y = Y1 ( X ):\n"; 545 : cout << "\n"; 546 : for ( int i = 0; i < 10; i++ ) 547 : { 548 : x = 10.0 * randreal1 ( ) - 5.0; 549 : y = y1 ( x ); 550 : cout << " " << x << " " << y << "\n"; 551 : } 552 : // 553 : // Y = YN ( 2, X ); 554 : // 555 : seed = 123456789; 556 : randinit ( seed ); 557 : cout << "\n"; 558 : cout << " Y = YN ( 2, X ):\n"; 559 : cout << "\n"; 560 : for ( int i = 0; i < 10; i++ ) 561 : { 562 : x = 10.0 * randreal1 ( ) - 5.0; 563 : y = yn ( 2, x ); 564 : cout << " " << x << " " << y << "\n"; 565 : } 566 : // 567 : // Terminate. 568 : // 569 : cout << "\n"; 570 : cout << "math_test:\n"; 571 : cout << " Normal end of execution.\n"; 572 : 573 : sizestack + 1024 =1352 ( 328 ) math_test: Try out the FreeFem++ Math Library Functions. Y = ABS ( X ): 0.32833 0.32833 4.90649 4.90649 0.341366 0.341366 -4.86846 4.86846 0.0955304 0.0955304 2.23866 2.23866 2.13564 2.13564 -3.64887 3.64887 -2.43001 2.43001 3.05856 3.05856 Y = ACOS ( X ): 0.065666 1.50508 0.981298 0.193702 0.0682732 1.50247 -0.973691 2.9117 0.0191061 1.55169 0.447731 1.10657 0.427128 1.12948 -0.729774 2.38879 -0.486002 2.07831 0.611712 0.912573 Y = ACOSH ( X ): 6.32833 2.53188 10.9065 3.0804 6.34137 2.53397 1.13154 0.50746 6.09553 2.49391 8.23866 2.79828 8.13564 2.7856 2.35113 1.4994 3.56999 1.94549 9.05856 2.8938 Y = ASIN ( X ): 0.065666 0.0657133 0.981298 1.37709 0.0682732 0.0683264 -0.973691 -1.3409 0.0191061 0.0191072 0.447731 0.464227 0.427128 0.441314 -0.729774 -0.817991 -0.486002 -0.507509 0.611712 0.658223 Y = ASINH ( X ): 0.32833 0.3227 4.90649 2.29393 0.341366 0.335061 -4.86846 -2.28631 0.0955304 0.0953857 2.23866 1.54554 2.13564 1.5027 -3.64887 -2.00583 -2.43001 -1.62092 3.05856 1.8368 Z = ATAN ( Y / X ): 0.32833 4.90649 1.50398 0.341366 -4.86846 -1.50079 0.0955304 2.23866 1.52815 2.13564 -3.64887 -1.04126 -2.43001 3.05856 -0.899421 2.52694 -4.37065 -1.04659 3.83879 3.61305 0.755114 -3.45101 4.26903 -0.890965 1.70546 -0.26315 -0.153091 1.43445 -2.21468 -0.996038 Z = ATAN2 ( Y, X ): 0.32833 4.90649 1.50398 0.341366 -4.86846 -1.50079 0.0955304 2.23866 1.52815 2.13564 -3.64887 -1.04126 -2.43001 3.05856 2.24217 2.52694 -4.37065 -1.04659 3.83879 3.61305 0.755114 -3.45101 4.26903 2.25063 1.70546 -0.26315 -0.153091 1.43445 -2.21468 -0.996038 Y = ATANH ( X ): 0.065666 0.0657607 0.981298 2.33145 0.0682732 0.0683796 -0.973691 -2.15887 0.0191061 0.0191084 0.447731 0.481859 0.427128 0.456379 -0.729774 -0.928244 -0.486002 -0.530813 0.611712 0.711653 Y = CEIL ( X ): 0.32833 1 4.90649 5 0.341366 1 -4.86846 -4 0.0955304 1 2.23866 3 2.13564 3 -3.64887 -3 -2.43001 -2 3.05856 4 Y = COS ( X ): 0.32833 0.946582 4.90649 0.192886 0.341366 0.942298 -4.86846 0.155433 0.0955304 0.99544 2.23866 -0.619308 2.13564 -0.535284 -3.64887 -0.87407 -2.43001 -0.75733 3.05856 -0.996555 Y = COSH ( X ): 0.32833 1.05439 4.90649 67.5859 0.341366 1.05883 -4.86846 65.0637 0.0955304 1.00457 2.23866 4.74366 2.13564 4.29032 -3.64887 19.2286 -2.43001 5.72352 3.05856 10.6719 Y = ERF ( X ): 0.32833 0.357588 4.90649 1 0.341366 0.370737 -4.86846 -1 0.0955304 0.107467 2.23866 0.998454 2.13564 0.997474 -3.64887 -1 -2.43001 -0.999411 3.05856 0.999985 Y = ERFC ( X ): 0.32833 0.642412 4.90649 3.9537e-12 0.341366 0.629263 -4.86846 2 0.0955304 0.892533 2.23866 0.00154583 2.13564 0.00252565 -3.64887 2 -2.43001 1.99941 3.05856 1.52216e-05 Y = EXP ( X ): 0.32833 1.38865 4.90649 135.164 0.341366 1.40687 -4.86846 0.00768523 0.0955304 1.10024 2.23866 9.38072 2.13564 8.46246 -3.64887 0.0260205 -2.43001 0.0880359 3.05856 21.2969 Y = FLOOR ( X ): 0.32833 0 4.90649 4 0.341366 0 -4.86846 -5 0.0955304 0 2.23866 2 2.13564 2 -3.64887 -4 -2.43001 -3 3.05856 3 Y = J0 ( X ): 0.32833 0.973231 4.90649 -0.207692 0.341366 0.971079 -4.86846 -0.219591 0.0955304 0.99772 2.23866 0.0889804 2.13564 0.146422 -3.64887 -0.395936 -2.43001 -0.0130053 3.05856 -0.279262 Y = J1 ( X ): 0.32833 0.161963 4.90649 -0.315632 0.341366 0.168209 -4.86846 0.30994 0.0955304 0.0477107 2.23866 0.550181 2.13564 0.564339 -3.64887 -0.0750638 -2.43001 -0.513604 3.05856 0.316919 Y = JN ( 2, X ): 0.32833 0.0133544 4.90649 0.0790334 0.341366 0.0144254 -4.86846 0.0922654 0.0955304 0.00113989 2.23866 0.402547 2.13564 0.382074 -3.64887 0.437079 -2.43001 0.435722 3.05856 0.486496 Y = LGAMMA ( X ): 5.32833 3.68421 9.90649 12.5917 5.34137 3.70477 0.131545 1.96587 5.09553 3.32294 7.23866 7.03052 7.13564 6.83468 1.35113 -0.115369 2.56999 0.335084 8.05856 8.64343 Y = LOG ( X ): 5.32833 1.67304 9.90649 2.29319 5.34137 1.67548 0.131545 -2.02841 5.09553 1.62836 7.23866 1.97944 7.13564 1.9651 1.35113 0.300941 2.56999 0.943902 8.05856 2.08673 Y = LOG10 ( X ): 5.32833 0.726591 9.90649 0.99592 5.34137 0.727652 0.131545 -0.880926 5.09553 0.707189 7.23866 0.859658 7.13564 0.853433 1.35113 0.130697 2.56999 0.409931 8.05856 0.906257 Z = MAX ( X, Y ): 0.32833 4.90649 4.90649 0.341366 -4.86846 0.341366 0.0955304 2.23866 2.23866 2.13564 -3.64887 2.13564 -2.43001 3.05856 3.05856 2.52694 -4.37065 2.52694 3.83879 3.61305 3.83879 -3.45101 4.26903 4.26903 1.70546 -0.26315 1.70546 1.43445 -2.21468 1.43445 Z = MIN ( X, Y ): 0.32833 4.90649 0.32833 0.341366 -4.86846 -4.86846 0.0955304 2.23866 0.0955304 2.13564 -3.64887 -3.64887 -2.43001 3.05856 -2.43001 2.52694 -4.37065 -4.37065 3.83879 3.61305 3.61305 -3.45101 4.26903 -3.45101 1.70546 -0.26315 -0.26315 1.43445 -2.21468 -2.21468 Z = POW ( X, Y ): 5.32833 4.90649 3672.93 5.34137 -4.86846 0.00028672 5.09553 2.23866 38.2961 7.13564 -3.64887 0.000769013 2.56999 3.05856 17.9391 7.52694 -4.37065 0.000147439 8.83879 3.61305 2626.44 1.54899 4.26903 6.47627 6.70546 -0.26315 0.606073 6.43445 -2.21468 0.016196 Y = RINT ( X ): 0.32833 0 4.90649 5 0.341366 0 -4.86846 -5 0.0955304 0 2.23866 2 2.13564 2 -3.64887 -4 -2.43001 -2 3.05856 3 Y = SIN ( X ): 0.32833 0.322463 4.90649 -0.981221 0.341366 0.334775 -4.86846 0.987846 0.0955304 0.0953851 2.23866 0.785148 2.13564 0.844672 -3.64887 0.4858 -2.43001 -0.653033 3.05856 0.0829369 Y = SINH ( X ): 0.32833 0.334261 4.90649 67.5785 0.341366 0.348035 -4.86846 -65.056 0.0955304 0.0956757 2.23866 4.63706 2.13564 4.17215 -3.64887 -19.2026 -2.43001 -5.63548 3.05856 10.625 Y = SQRT ( X ): 5.32833 2.30832 9.90649 3.14746 5.34137 2.31114 0.131545 0.362691 5.09553 2.25733 7.23866 2.69048 7.13564 2.67126 1.35113 1.16238 2.56999 1.60312 8.05856 2.83876 Y = TAN ( X ): 0.32833 0.34066 4.90649 -5.08706 0.341366 0.355275 -4.86846 6.35544 0.0955304 0.095822 2.23866 -1.26778 2.13564 -1.57799 -3.64887 -0.55579 -2.43001 0.862284 3.05856 -0.0832236 Y = TANH ( X ): 0.32833 0.31702 4.90649 0.999891 0.341366 0.328696 -4.86846 -0.999882 0.0955304 0.0952408 2.23866 0.977528 2.13564 0.972457 -3.64887 -0.998647 -2.43001 -0.984619 3.05856 0.9956 Y = TGAMMA ( X ): 0.32833 2.72159 4.90649 20.8675 0.341366 2.61322 -4.86846 -0.0814629 0.0955304 9.97764 2.23866 1.12571 2.13564 1.06515 -3.64887 0.245395 -2.43001 -1.04553 3.05856 2.11247 Y = Y0 ( X ): 0.32833 -0.74489 4.90649 -0.293224 0.341366 -0.717775 -4.86846 nan 0.0955304 -1.56376 2.23866 0.52034 2.13564 0.519792 -3.64887 nan -2.43001 nan 3.05856 0.357389 Y = Y1 ( X ): 0.32833 -2.11626 4.90649 0.179109 0.341366 -2.04479 -4.86846 nan 0.0955304 -6.75407 2.23866 0.0214141 2.13564 -0.0324696 -3.64887 nan -2.43001 nan 3.05856 0.339753 Y = YN ( 2, X ): 0.32833 -12.1461 4.90649 0.366233 0.341366 -11.2623 -4.86846 nan 0.0955304 -139.838 2.23866 -0.501209 2.13564 -0.5502 -3.64887 nan -2.43001 nan 3.05856 -0.135223 math_test: Normal end of execution. times: compile 0.007661s, execution 0.002942s, mpirank:0 CodeAlloc : nb ptr 5093, size :506928 mpirank: 0 Ok: Normal End