-- FreeFem++ v4.14 (mer. 06 mars 2024 16:59:04 CET - git v4.14-1-g2b2052ae) file : math_test.edp Load: lg_fem lg_mesh lg_mesh3 eigenvalue 1 : // math_test.edp 2 : // 3 : // Discussion: 4 : // 5 : // Examine the FreeFem math functions: 6 : // 7 : // abs 8 : // acos 9 : // acosh 10 : // asin 11 : // asinh 12 : // atan 13 : // atan2 14 : // atanh 15 : // ceil 16 : // cos 17 : // cosh 18 : // erf 19 : // erfc 20 : // exp 21 : // floor 22 : // imag 23 : // j0 24 : // j1 25 : // jn 26 : // lgamma 27 : // log 28 : // log10 29 : // max 30 : // min 31 : // pow 32 : // rint 33 : // sin 34 : // sinh 35 : // sqrt 36 : // tan 37 : // tanh 38 : // tgamma 39 : // y0 40 : // y1 41 : // yn 42 : // 43 : // Licensing: 44 : // 45 : // This code is distributed under the MIT license. 46 : // 47 : // Modified: 48 : // 49 : // 15 February 2015 50 : // 51 : // Author: 52 : // 53 : // John Burkardt 54 : // 55 : // Reference: 56 : // 57 : // Frederic Hecht, 58 : // Freefem++, 59 : // Third Edition, version 3.22 60 : // 61 : cout << "\n"; 62 : cout << "math_test():\n"; 63 : cout << " FreeFem++ version\n"; 64 : cout << " Try out the FreeFem++ Math Library Functions.\n ... : "; 65 : // 66 : // Initialize the random number generator. 67 : // 68 : int seed; 69 : // 70 : // Y = ABS ( X ); 71 : // 72 : seed = 123456789; 73 : randinit ( seed ); 74 : cout << "\n"; 75 : cout << " Y = ABS ( X ):\n"; 76 : cout << "\n"; 77 : for ( int i = 0; i < 10; i++ ) 78 : { 79 : x = 10.0 * randreal1 ( ) - 5.0; 80 : y = abs ( x ); 81 : cout << " " << x << " " << y << "\n"; 82 : } 83 : // 84 : // Y = ACOS ( X ); 85 : // X <= 1.0 required. 86 : // 87 : seed = 123456789; 88 : randinit ( seed ); 89 : cout << "\n"; 90 : cout << " Y = ACOS ( X ):\n"; 91 : cout << "\n"; 92 : for ( int i = 0; i < 10; i++ ) 93 : { 94 : x = 2.0 * randreal1 ( ) - 1.0; 95 : y = acos ( x ); 96 : cout << " " << x << " " << y << "\n"; 97 : } 98 : // 99 : // Y = ACOSH ( X ); 100 : // 1 <= |X| required. 101 : // 102 : seed = 123456789; 103 : randinit ( seed ); 104 : cout << "\n"; 105 : cout << " Y = ACOSH ( X ):\n"; 106 : cout << "\n"; 107 : for ( int i = 0; i < 10; i++ ) 108 : { 109 : x = randreal1 ( ); 110 : if ( x < 0.0 ) 111 : { 112 : x = 10.0 * x - 1.0; 113 : } 114 : else 115 : { 116 : x = 10.0 * x + 1.0; 117 : } 118 : y = acosh ( x ); 119 : cout << " " << x << " " << y << "\n"; 120 : } 121 : // 122 : // Y = ASIN ( X ); 123 : // |X| <= 1. 124 : // 125 : seed = 123456789; 126 : randinit ( seed ); 127 : cout << "\n"; 128 : cout << " Y = ASIN ( X ):\n"; 129 : cout << "\n"; 130 : for ( int i = 0; i < 10; i++ ) 131 : { 132 : x = 2.0 * randreal1 ( ) - 1.0; 133 : y = asin ( x ); 134 : cout << " " << x << " " << y << "\n"; 135 : } 136 : // 137 : // Y = ASINH ( X ); 138 : // 139 : seed = 123456789; 140 : randinit ( seed ); 141 : cout << "\n"; 142 : cout << " Y = ASINH ( X ):\n"; 143 : cout << "\n"; 144 : for ( int i = 0; i < 10; i++ ) 145 : { 146 : x = 10.0 * randreal1 ( ) - 5.0; 147 : y = asinh ( x ); 148 : cout << " " << x << " " << y << "\n"; 149 : } 150 : // 151 : // Z = ATAN ( Y / X ); 152 : // 153 : seed = 123456789; 154 : randinit ( seed ); 155 : cout << "\n"; 156 : cout << " Z = ATAN ( Y / X ):\n"; 157 : cout << "\n"; 158 : for ( int i = 0; i < 10; i++ ) 159 : { 160 : x = 10.0 * randreal1 ( ) - 5.0; 161 : y = 10.0 * randreal1 ( ) - 5.0; 162 : z = atan ( y / x ); 163 : cout << " " << x << " " << y << " " << z << "\n"; 164 : } 165 : // 166 : // Z = ATAN2 ( Y, X ); 167 : // 168 : seed = 123456789; 169 : randinit ( seed ); 170 : cout << "\n"; 171 : cout << " Z = ATAN2 ( Y, X ):\n"; 172 : cout << "\n"; 173 : for ( int i = 0; i < 10; i++ ) 174 : { 175 : x = 10.0 * randreal1 ( ) - 5.0; 176 : y = 10.0 * randreal1 ( ) - 5.0; 177 : z = atan2 ( y, x ); 178 : cout << " " << x << " " << y << " " << z << "\n"; 179 : } 180 : // 181 : // Y = ATANH ( X ); 182 : // |X| < 1 required. 183 : // 184 : seed = 123456789; 185 : randinit ( seed ); 186 : cout << "\n"; 187 : cout << " Y = ATANH ( X ):\n"; 188 : cout << "\n"; 189 : for ( int i = 0; i < 10; i++ ) 190 : { 191 : x = 2.0 * randreal1 ( ) - 1.0; 192 : y = atanh ( x ); 193 : cout << " " << x << " " << y << "\n"; 194 : } 195 : // 196 : // Y = CEIL ( X ); 197 : // 198 : seed = 123456789; 199 : randinit ( seed ); 200 : cout << "\n"; 201 : cout << " Y = CEIL ( X ):\n"; 202 : cout << "\n"; 203 : for ( int i = 0; i < 10; i++ ) 204 : { 205 : x = 10.0 * randreal1 ( ) - 5.0; 206 : y = ceil ( x ); 207 : cout << " " << x << " " << y << "\n"; 208 : } 209 : // 210 : // Y = COS ( X ); 211 : // 212 : seed = 123456789; 213 : randinit ( seed ); 214 : cout << "\n"; 215 : cout << " Y = COS ( X ):\n"; 216 : cout << "\n"; 217 : for ( int i = 0; i < 10; i++ ) 218 : { 219 : x = 10.0 * randreal1 ( ) - 5.0; 220 : y = cos ( x ); 221 : cout << " " << x << " " << y << "\n"; 222 : } 223 : // 224 : // Y = COSH ( X ); 225 : // 226 : seed = 123456789; 227 : randinit ( seed ); 228 : cout << "\n"; 229 : cout << " Y = COSH ( X ):\n"; 230 : cout << "\n"; 231 : for ( int i = 0; i < 10; i++ ) 232 : { 233 : x = 10.0 * randreal1 ( ) - 5.0; 234 : y = cosh ( x ); 235 : cout << " " << x << " " << y << "\n"; 236 : } 237 : // 238 : // Y = ERF ( X ); 239 : // 240 : seed = 123456789; 241 : randinit ( seed ); 242 : cout << "\n"; 243 : cout << " Y = ERF ( X ):\n"; 244 : cout << "\n"; 245 : for ( int i = 0; i < 10; i++ ) 246 : { 247 : x = 10.0 * randreal1 ( ) - 5.0; 248 : y = erf ( x ); 249 : cout << " " << x << " " << y << "\n"; 250 : } 251 : // 252 : // Y = ERFC ( X ); 253 : // 254 : seed = 123456789; 255 : randinit ( seed ); 256 : cout << "\n"; 257 : cout << " Y = ERFC ( X ):\n"; 258 : cout << "\n"; 259 : for ( int i = 0; i < 10; i++ ) 260 : { 261 : x = 10.0 * randreal1 ( ) - 5.0; 262 : y = erfc ( x ); 263 : cout << " " << x << " " << y << "\n"; 264 : } 265 : // 266 : // Y = EXP ( X ); 267 : // 268 : seed = 123456789; 269 : randinit ( seed ); 270 : cout << "\n"; 271 : cout << " Y = EXP ( X ):\n"; 272 : cout << "\n"; 273 : for ( int i = 0; i < 10; i++ ) 274 : { 275 : x = 10.0 * randreal1 ( ) - 5.0; 276 : y = exp ( x ); 277 : cout << " " << x << " " << y << "\n"; 278 : } 279 : // 280 : // Y = FLOOR ( X ); 281 : // 282 : seed = 123456789; 283 : randinit ( seed ); 284 : cout << "\n"; 285 : cout << " Y = FLOOR ( X ):\n"; 286 : cout << "\n"; 287 : for ( int i = 0; i < 10; i++ ) 288 : { 289 : x = 10.0 * randreal1 ( ) - 5.0; 290 : y = floor ( x ); 291 : cout << " " << x << " " << y << "\n"; 292 : } 293 : // 294 : // Y = J0 ( X ); 295 : // 296 : seed = 123456789; 297 : randinit ( seed ); 298 : cout << "\n"; 299 : cout << " Y = J0 ( X ):\n"; 300 : cout << "\n"; 301 : for ( int i = 0; i < 10; i++ ) 302 : { 303 : x = 10.0 * randreal1 ( ) - 5.0; 304 : y = j0 ( x ); 305 : cout << " " << x << " " << y << "\n"; 306 : } 307 : // 308 : // Y = J1 ( X ); 309 : // 310 : seed = 123456789; 311 : randinit ( seed ); 312 : cout << "\n"; 313 : cout << " Y = J1 ( X ):\n"; 314 : cout << "\n"; 315 : for ( int i = 0; i < 10; i++ ) 316 : { 317 : x = 10.0 * randreal1 ( ) - 5.0; 318 : y = j1 ( x ); 319 : cout << " " << x << " " << y << "\n"; 320 : } 321 : // 322 : // Y = JN ( 2, X ); 323 : // 324 : seed = 123456789; 325 : randinit ( seed ); 326 : cout << "\n"; 327 : cout << " Y = JN ( 2, X ):\n"; 328 : cout << "\n"; 329 : for ( int i = 0; i < 10; i++ ) 330 : { 331 : x = 10.0 * randreal1 ( ) - 5.0; 332 : y = jn ( 2, x ); 333 : cout << " " << x << " " << y << "\n"; 334 : } 335 : // 336 : // Y = LGAMMA ( X ); 337 : // X must not be 0 or a negative integer. 338 : // 339 : seed = 123456789; 340 : randinit ( seed ); 341 : cout << "\n"; 342 : cout << " Y = LGAMMA ( X ):\n"; 343 : cout << "\n"; 344 : for ( int i = 0; i < 10; i++ ) 345 : { 346 : x = 10.0 * randreal1 ( ); 347 : y = lgamma ( x ); 348 : cout << " " << x << " " << y << "\n"; 349 : } 350 : // 351 : // Y = LOG ( X ); 352 : // 0 < X required. 353 : // 354 : seed = 123456789; 355 : randinit ( seed ); 356 : cout << "\n"; 357 : cout << " Y = LOG ( X ):\n"; 358 : cout << "\n"; 359 : for ( int i = 0; i < 10; i++ ) 360 : { 361 : x = 10.0 * randreal1 ( ); 362 : y = log ( x ); 363 : cout << " " << x << " " << y << "\n"; 364 : } 365 : // 366 : // Y = LOG10 ( X ); 367 : // 0 < X required. 368 : // 369 : seed = 123456789; 370 : randinit ( seed ); 371 : cout << "\n"; 372 : cout << " Y = LOG10 ( X ):\n"; 373 : cout << "\n"; 374 : for ( int i = 0; i < 10; i++ ) 375 : { 376 : x = 10.0 * randreal1 ( ); 377 : y = log10 ( x ); 378 : cout << " " << x << " " << y << "\n"; 379 : } 380 : // 381 : // Z = MAX ( X, Y ); 382 : // 383 : seed = 123456789; 384 : randinit ( seed ); 385 : cout << "\n"; 386 : cout << " Z = MAX ( X, Y ):\n"; 387 : cout << "\n"; 388 : for ( int i = 0; i < 10; i++ ) 389 : { 390 : x = 10.0 * randreal1 ( ) - 5.0; 391 : y = 10.0 * randreal1 ( ) - 5.0; 392 : z = max ( x, y ); 393 : cout << " " << x << " " << y << " " << z << "\n"; 394 : } 395 : // 396 : // Z = MIN ( X, Y ); 397 : // 398 : seed = 123456789; 399 : randinit ( seed ); 400 : cout << "\n"; 401 : cout << " Z = MIN ( X, Y ):\n"; 402 : cout << "\n"; 403 : for ( int i = 0; i < 10; i++ ) 404 : { 405 : x = 10.0 * randreal1 ( ) - 5.0; 406 : y = 10.0 * randreal1 ( ) - 5.0; 407 : z = min ( x, y ); 408 : cout << " " << x << " " << y << " " << z << "\n"; 409 : } 410 : // 411 : // Z = POW ( X, Y ); 412 : // If X negative, Y must be an integer. 413 : // If X 0, Y must be positive. 414 : // 415 : seed = 123456789; 416 : randinit ( seed ); 417 : cout << "\n"; 418 : cout << " Z = POW ( X, Y ):\n"; 419 : cout << "\n"; 420 : for ( int i = 0; i < 10; i++ ) 421 : { 422 : x = 10.0 * randreal1 ( ); 423 : y = 10.0 * randreal1 ( ) - 5.0; 424 : z = pow ( x, y ); 425 : cout << " " << x << " " << y << " " << z << "\n"; 426 : } 427 : // 428 : // Y = RINT ( X ); 429 : // 430 : seed = 123456789; 431 : randinit ( seed ); 432 : cout << "\n"; 433 : cout << " Y = RINT ( X ):\n"; 434 : cout << "\n"; 435 : for ( int i = 0; i < 10; i++ ) 436 : { 437 : x = 10.0 * randreal1 ( ) - 5.0; 438 : y = rint ( x ); 439 : cout << " " << x << " " << y << "\n"; 440 : } 441 : // 442 : // Y = SIN ( X ); 443 : // 444 : seed = 123456789; 445 : randinit ( seed ); 446 : cout << "\n"; 447 : cout << " Y = SIN ( X ):\n"; 448 : cout << "\n"; 449 : for ( int i = 0; i < 10; i++ ) 450 : { 451 : x = 10.0 * randreal1 ( ) - 5.0; 452 : y = sin ( x ); 453 : cout << " " << x << " " << y << "\n"; 454 : } 455 : // 456 : // Y = SINH ( X ); 457 : // 458 : seed = 123456789; 459 : randinit ( seed ); 460 : cout << "\n"; 461 : cout << " Y = SINH ( X ):\n"; 462 : cout << "\n"; 463 : for ( int i = 0; i < 10; i++ ) 464 : { 465 : x = 10.0 * randreal1 ( ) - 5.0; 466 : y = sinh ( x ); 467 : cout << " " << x << " " << y << "\n"; 468 : } 469 : // 470 : // Y = SQRT ( X ); 471 : // 0 <= X required. 472 : // 473 : seed = 123456789; 474 : randinit ( seed ); 475 : cout << "\n"; 476 : cout << " Y = SQRT ( X ):\n"; 477 : cout << "\n"; 478 : for ( int i = 0; i < 10; i++ ) 479 : { 480 : x = 10.0 * randreal1 ( ); 481 : y = sqrt ( x ); 482 : cout << " " << x << " " << y << "\n"; 483 : } 484 : // 485 : // Y = TAN ( X ); 486 : // 487 : seed = 123456789; 488 : randinit ( seed ); 489 : cout << "\n"; 490 : cout << " Y = TAN ( X ):\n"; 491 : cout << "\n"; 492 : for ( int i = 0; i < 10; i++ ) 493 : { 494 : x = 10.0 * randreal1 ( ) - 5.0; 495 : y = tan ( x ); 496 : cout << " " << x << " " << y << "\n"; 497 : } 498 : // 499 : // Y = TANH ( X ); 500 : // 501 : seed = 123456789; 502 : randinit ( seed ); 503 : cout << "\n"; 504 : cout << " Y = TANH ( X ):\n"; 505 : cout << "\n"; 506 : for ( int i = 0; i < 10; i++ ) 507 : { 508 : x = 10.0 * randreal1 ( ) - 5.0; 509 : y = tanh ( x ); 510 : cout << " " << x << " " << y << "\n"; 511 : } 512 : // 513 : // Y = TGAMMA ( X ); 514 : // X must not be 0 or a negative integer. 515 : // 516 : seed = 123456789; 517 : randinit ( seed ); 518 : cout << "\n"; 519 : cout << " Y = TGAMMA ( X ):\n"; 520 : cout << "\n"; 521 : for ( int i = 0; i < 10; i++ ) 522 : { 523 : x = 10.0 * randreal1 ( ) - 5.0; 524 : y = tgamma ( x ); 525 : cout << " " << x << " " << y << "\n"; 526 : } 527 : // 528 : // Y = Y0 ( X ); 529 : // 530 : seed = 123456789; 531 : randinit ( seed ); 532 : cout << "\n"; 533 : cout << " Y = Y0 ( X ):\n"; 534 : cout << "\n"; 535 : for ( int i = 0; i < 10; i++ ) 536 : { 537 : x = 10.0 * randreal1 ( ) - 5.0; 538 : y = y0 ( x ); 539 : cout << " " << x << " " << y << "\n"; 540 : } 541 : // 542 : // Y = Y1 ( X ); 543 : // 544 : seed = 123456789; 545 : randinit ( seed ); 546 : cout << "\n"; 547 : cout << " Y = Y1 ( X ):\n"; 548 : cout << "\n"; 549 : for ( int i = 0; i < 10; i++ ) 550 : { 551 : x = 10.0 * randreal1 ( ) - 5.0; 552 : y = y1 ( x ); 553 : cout << " " << x << " " << y << "\n"; 554 : } 555 : // 556 : // Y = YN ( 2, X ); 557 : // 558 : seed = 123456789; 559 : randinit ( seed ); 560 : cout << "\n"; 561 : cout << " Y = YN ( 2, X ):\n"; 562 : cout << "\n"; 563 : for ( int i = 0; i < 10; i++ ) 564 : { 565 : x = 10.0 * randreal1 ( ) - 5.0; 566 : y = yn ( 2, x ); 567 : cout << " " << x << " " << y << "\n"; 568 : } 569 : // 570 : // Terminate. 571 : // 572 : cout << "\n"; 573 : cout << "math_test():\n"; 574 : cout << " Normal end of execution.\n"; 575 : 576 : exit ( 0 ); 577 : 578 : sizestack + 1024 =1352 ( 328 ) math_test(): FreeFem++ version 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. current line = 576 exit(0) err code 0 , mpirank 0 CodeAlloc : nb ptr 5468, size :554904 mpirank: 0 Ok: Normal End