28 October 2014 8:42:17.773 AM POLYNOMIAL_PRB FORTRAN90 version. Test the POLYNOMIAL library. POLYNOMIAL_ADD_TEST POLYNOMIAL_ADD adds two polynomials. P1(X) = + 7.00000 * x^( 0, 0, 0) - 5.00000 * x^( 0, 0, 1) + 9.00000 * x^( 1, 0, 0) + 11.0000 * x^( 0, 0, 2) + 0.00000 * x^( 0, 1, 2) - 13.0000 * x^( 3, 0, 1). P2(X) = + 2.00000 * x^( 0, 0, 0) + 3.00000 * x^( 0, 1, 0) - 8.00000 * x^( 1, 0, 0) + 4.00000 * x^( 2, 0, 2) + 9.00000 * x^( 3, 0, 1). P1(X) + P2(X) = + 9.00000 * x^( 0, 0, 0) - 5.00000 * x^( 0, 0, 1) + 3.00000 * x^( 0, 1, 0) + 1.00000 * x^( 1, 0, 0) + 11.0000 * x^( 0, 0, 2) + 4.00000 * x^( 2, 0, 2) - 4.00000 * x^( 3, 0, 1). POLYNOMIAL_AXPY_TEST POLYNOMIAL_AXPY adds a multiple of one polynomial to another. P1(X) = + 7.00000 * x^( 0, 0, 0) - 5.00000 * x^( 0, 0, 1) + 9.00000 * x^( 1, 0, 0) + 11.0000 * x^( 0, 0, 2) + 0.00000 * x^( 0, 1, 2) - 13.0000 * x^( 3, 0, 1). P2(X) = + 2.00000 * x^( 0, 0, 0) + 3.00000 * x^( 0, 1, 0) - 8.00000 * x^( 1, 0, 0) + 4.00000 * x^( 2, 0, 2) + 9.00000 * x^( 3, 0, 1). 10.0000 * P1(X) + P2(X) = + 72.0000 * x^( 0, 0, 0) - 50.0000 * x^( 0, 0, 1) + 3.00000 * x^( 0, 1, 0) + 82.0000 * x^( 1, 0, 0) + 110.000 * x^( 0, 0, 2) + 4.00000 * x^( 2, 0, 2) - 121.000 * x^( 3, 0, 1). POLYNOMIAL_COMPRESS_TEST POLYNOMIAL_COMPRESS compresses a polynomial. Uncompressed P(X) = + 7.00000 * x^( 0, 0, 0) - 5.00000 * x^( 0, 0, 1) + 5.00000 * x^( 0, 0, 1) + 9.00000 * x^( 1, 0, 0) + 11.0000 * x^( 0, 0, 2) + 3.00000 * x^( 0, 0, 2) + 6.00000 * x^( 0, 0, 2) + 0.00000 * x^( 0, 1, 2) - 13.0000 * x^( 3, 0, 1) + 0.100000E-19 * x^( 4, 0, 0). Compressed P(X) = + 7.00000 * x^( 0, 0, 0) + 9.00000 * x^( 1, 0, 0) + 20.0000 * x^( 0, 0, 2) - 13.0000 * x^( 3, 0, 1). POLYNOMIAL_DIF_TEST POLYNOMIAL_DIF computes derivatives of a polynomial. P(X) = + 2.00000 * x^( 0, 0) + 3.00000 * x^( 3, 0) + 4.00000 * x^( 1, 3) + 5.00000 * x^( 3, 4). d3 P(X) dx1 dx1 dx2 = + 120.000 * x^( 1, 3). POLYNOMIAL_MUL_TEST POLYNOMIAL_MUL multiplies two polynomials. P1(X) = + 2.00000 * x^( 0, 0, 0) + 3.00000 * x^( 0, 1, 0) + 4.00000 * x^( 1, 0, 0) + 5.00000 * x^( 0, 1, 1). P2(X) = + 6.00000 * x^( 0, 0, 1) + 7.00000 * x^( 0, 0, 2). P1(X) * P2(X) = + 12.0000 * x^( 0, 0, 1) + 14.0000 * x^( 0, 0, 2) + 18.0000 * x^( 0, 1, 1) + 24.0000 * x^( 1, 0, 1) + 51.0000 * x^( 0, 1, 2) + 28.0000 * x^( 1, 0, 2) + 35.0000 * x^( 0, 1, 3). POLYNOMIAL_PRINT_TEST POLYNOMIAL_PRINT prints a polynomial. P1(X) = + 7.00000 * x^( 0, 0, 0) - 5.00000 * x^( 0, 0, 1) + 9.00000 * x^( 1, 0, 0) + 11.0000 * x^( 0, 0, 2) + 0.00000 * x^( 0, 1, 2) - 13.0000 * x^( 3, 0, 1). POLYNOMIAL_SCALE_TEST POLYNOMIAL_SCALE scales a polynomial by a multiplier S. P(X) = + 7.00000 * x^( 0, 0, 0) - 5.00000 * x^( 0, 0, 1) + 9.00000 * x^( 1, 0, 0) + 11.0000 * x^( 0, 0, 2) + 0.00000 * x^( 0, 1, 2) - 13.0000 * x^( 3, 0, 1). Apply scale factor S = -0.500000 S * P(X) = - 3.50000 * x^( 0, 0, 0) + 2.50000 * x^( 0, 0, 1) - 4.50000 * x^( 1, 0, 0) - 5.50000 * x^( 0, 0, 2) + 0.00000 * x^( 0, 1, 2) + 6.50000 * x^( 3, 0, 1). POLYNOMIAL_SORT_TEST POLYNOMIAL_SORT sorts a polynomial by exponent index.. Unsorted polynomial: + 0.00000 * x^( 0, 1, 2) + 9.00000 * x^( 1, 0, 0) - 5.00000 * x^( 0, 0, 1) - 13.0000 * x^( 3, 0, 1) + 7.00000 * x^( 0, 0, 0) + 11.0000 * x^( 0, 0, 2). Sorted polynomial: + 7.00000 * x^( 0, 0, 0) - 5.00000 * x^( 0, 0, 1) + 9.00000 * x^( 1, 0, 0) + 11.0000 * x^( 0, 0, 2) + 0.00000 * x^( 0, 1, 2) - 13.0000 * x^( 3, 0, 1). POLYNOMIAL_VALUE_TEST POLYNOMIAL_VALUE evaluates a polynomial. P(X) = + 7.00000 * x^( 0, 0, 0) - 5.00000 * x^( 0, 0, 1) + 9.00000 * x^( 1, 0, 0) + 11.0000 * x^( 0, 0, 2) + 0.00000 * x^( 0, 1, 2) - 13.0000 * x^( 3, 0, 1). P( 1.0000, 2.0000, 3.0000) = 61.0000 P( -2.0000, 4.0000, 1.0000) = 99.0000 POLYNOMIAL_PRB Normal end of execution. 28 October 2014 8:42:17.774 AM