 do_test_newton using lapack
                                               expected    2.9347118120566711D-02
                                                 actual    2.9347118120566711D-02

 do_test_newton using block_thomas_dble
                                               expected    2.9347118120566711D-02
                                                 actual    2.9347118120566711D-02

 do_test_newton using block_thomas_quad
                                               expected    2.9347118120566711D-02
                                                 actual    2.9347118120566711D-02

 do_test_newton using bcyclic_dble
                                               expected    2.9347118120566711D-02
                                                 actual    2.9347118120566711D-02


 test_find0_quadratic
                                                      x    3.0000000000000000D+00
                                                      y   -0.0000000000000000D+00
                                                      x    8.0000000000000000D+00
                                                      y    0.0000000000000000D+00



 qsort into increasing order
    1   71    -0.99960
    2   72    -0.99929
    3   70    -0.99556
    4   73    -0.99463
    5   69    -0.98719
    6   74    -0.98564
    7   68    -0.97453
    8   75    -0.97237
    9   67    -0.95762
   10   76    -0.95486
   11   66    -0.93655
   12   77    -0.93320
   13   65    -0.91140
   14   78    -0.90748
   15   64    -0.88229
   16   79    -0.87782
   17   63    -0.84934
   18   80    -0.84433
   19   62    -0.81269
   20   81    -0.80717
   21   61    -0.77251
   22   82    -0.76649
   23   60    -0.72897
   24   83    -0.72248
   25   59    -0.68225
   26   84    -0.67533
   27   58    -0.63257
   28   85    -0.62524
   29   57    -0.58013
   30   86    -0.57243
   31   56    -0.52517
   32   87    -0.51713
   33   55    -0.46793
   34   88    -0.45958
   35   54    -0.40865
   36   89    -0.40003
   37   53    -0.34759
   38   90    -0.33874
   39   52    -0.28502
   40   91    -0.27597
   41   51    -0.22121
   42   92    -0.21201
   43   50    -0.15643
   44   93    -0.14712
   45   49    -0.09098
   46   94    -0.08159
   47   48    -0.02513
   48   95    -0.01571
   49   47     0.04083
   50   96     0.05024
   51    1     0.06593
   52   46     0.10661
   53   97     0.11598
   54    2     0.13156
   55   45     0.17193
   56   98     0.18121
   57    3     0.19663
   58   44     0.23650
   59   99     0.24565
   60    4     0.26084
   61   43     0.30004
   62  100     0.30902
   63    5     0.32392
   64   42     0.36228
   65    6     0.38558
   66   41     0.42293
   67    7     0.44557
   68   40     0.48175
   69    8     0.50362
   70   39     0.53848
   71    9     0.55948
   72   38     0.59286
   73   10     0.61291
   74   37     0.64466
   75   11     0.66367
   76   36     0.69365
   77   12     0.71154
   78   35     0.73963
   79   13     0.75631
   80   34     0.78239
   81   14     0.79779
   82   33     0.82175
   83   15     0.83581
   84   32     0.85753
   85   16     0.87018
   86   31     0.88958
   87   17     0.90077
   88   30     0.91775
   89   18     0.92745
   90   29     0.94194
   91   19     0.95008
   92   28     0.96203
   93   20     0.96858
   94   27     0.97793
   95   21     0.98287
   96   26     0.98958
   97   22     0.99288
   98   25     0.99692
   99   23     0.99857
  100   24     0.99992


    test NEWUOA with N = 2 and NPT =  5
    test NEWUOA with N = 4 and NPT =  9
    test NEWUOA with N = 6 and NPT = 13

    test BOBYQA with N = 2 and NPT =  5
    test BOBYQA with N = 4 and NPT =  9
    test BOBYQA with N = 6 and NPT = 13

 testing NM_simplex with FR4
                                  found f_final < 1d-10

 testing NM_simplex with ER
                                  found f_final < 1d-10

 testing NM_simplex with WD
                                  found f_final < 1d-10

 testing NM_simplex with BLE
                                  found f_final < 1d-10

 testing NM_simplex with PS
                                  found f_final < 1d-10

 testing NM_simplex with TR
                                  found f_final < 1d-10


 test brent
 
TEST_GLOMIN_ALL
  Test the Brent GLOMIN routine, which seeks
  a global minimizer of a function F(X)
  in an interval [A,B],
  given some upper bound M for F".
 
h_01(x) = 2 - x
 
      A                 X             B
    F(A)              F(X)          F(B)
 
      7.00000000      9.00000000      9.00000000
    -5.00000        -7.00000        -7.00000    
 
h_01(x) = 2 - x
 
      A                 X             B
    F(A)              F(X)          F(B)
 
      7.00000000      9.00000000      9.00000000
    -5.00000        -7.00000        -7.00000    
 
h_02(x) = x * x
 
      A                 X             B
    F(A)              F(X)          F(B)
 
     -1.00000000      0.00000000      2.00000000
     1.00000         0.00000         4.00000    
 
h_02(x) = x * x
 
      A                 X             B
    F(A)              F(X)          F(B)
 
     -1.00000000      0.00000000      2.00000000
     1.00000         0.00000         4.00000    
 
h_04(x) = ( x + sin(x) ) * exp(-x*x)
 
      A                 X             B
    F(A)              F(X)          F(B)
 
    -10.00000000     -0.67957866     10.00000000
   -0.351770E-42   -0.824239        0.351770E-42
 
h_05(x) = ( x - sin(x) ) * exp(-x*x)
 
      A                 X             B
    F(A)              F(X)          F(B)
 
    -10.00000000     -1.19513664     10.00000000
   -0.392246E-42   -0.634905E-01    0.392246E-42
 
TEST_LOCAL_MIN_ALL
  Test the Brent LOCAL_MIN routine, which seeks
  a local minimizer of a function F(X)
  in an interval [A,B].
 
g_01(x) = ( x - 2 ) * ( x - 2 ) + 1
 
      A                 X             B
    F(A)              F(X)          F(B)
 
      0.00000000      2.00000000      3.14159265
     5.00000         1.00000         2.30323    
 
g_02(x) = x * x + exp ( - x )
 
      A                 X             B
    F(A)              F(X)          F(B)
 
      0.00000000      0.35173370      1.00000000
     1.00000        0.827184         1.36788    
 
g_03(x) = x^4 + 2x^2 + x + 3
 
      A                 X             B
    F(A)              F(X)          F(B)
 
     -2.00000000     -0.23673290      2.00000000
     25.0000         2.87849         29.0000    
 
g_04(x) = exp ( x ) + 1 / ( 100 x )
 
      A                 X             B
    F(A)              F(X)          F(B)
 
      0.00010000      0.09534461      1.00000000
     101.000         1.20492         2.72828    
 
g_05(x) = exp ( x ) - 2x + 1/(100x) - 1/(1000000x^2)
 
      A                 X             B
    F(A)              F(X)          F(B)
 
      0.00020000      0.70320487      2.00000000
     25.9998        0.628026         3.39406    


 binary_search, increasing values
                val is less than vec(1):        0.000000        1.000000
 okay
                             vec(result)     1089.000000
                                     val     1091.000000
                           vec(result+1)     1156.000000
 okay
             val is greater than vec(n):    10001.000000    10000.000000
 okay

 binary_search, decreasing values
             val is greater than vec(n):        0.000000       -1.000000
 okay
                           vec(result+1)    -1156.000000
                                     val    -1091.000000
                             vec(result)    -1089.000000
 okay
                val is less than vec(1):   -10001.000000   -10000.000000
 okay


 root  0.74800610880593488     
 test_root2
 x1  -1.2000000000000000     
 x3   2.0000000000000000     
 y1  -1.2957717642721649     
 y3  0.50191042281580711     
 safe_root  0.50000000000000000     

 test_root3
 first safe_root_without_brackets  0.49999999999999356     
 second safe_root_without_brackets  0.50000000000000377     


 vdpol
 cash_karp

    i      calculated           reference    
    1    0.1763462599E+01    0.1763234540E+01
    2   -0.8354500992E+00   -0.8356886819E+00


 vdpol
 dopri5

    i      calculated           reference    
    1    0.1763536960E+01    0.1763234540E+01
    2   -0.8353884656E+00   -0.8356886819E+00


 vdpol
 dop853

    i      calculated           reference    
    1    0.1763177087E+01    0.1763234540E+01
    2   -0.8357624953E+00   -0.8356886819E+00

 ijob 1
 numerical jacobians
 vdpol
 ros2
 lapack_decsol

    i      calculated           reference    
    1    0.1763479136E+01    0.1763234540E+01
    2   -0.8354253151E+00   -0.8356886819E+00

 vdpol
 rose2
 lapack_decsol

    i      calculated           reference    
    1    0.1763259925E+01    0.1763234540E+01
    2   -0.8356492331E+00   -0.8356886819E+00

 vdpol
 ros3p
 lapack_decsol

    i      calculated           reference    
    1    0.1763213375E+01    0.1763234540E+01
    2   -0.8357203644E+00   -0.8356886819E+00

 vdpol
 ros3pl
 lapack_decsol

    i      calculated           reference    
    1    0.1763227588E+01    0.1763234540E+01
    2   -0.8356793386E+00   -0.8356886819E+00

 vdpol
 rodas3
 lapack_decsol

    i      calculated           reference    
    1    0.1763126408E+01    0.1763234540E+01
    2   -0.8357846804E+00   -0.8356886819E+00

 vdpol
 rodas4
 lapack_decsol

    i      calculated           reference    
    1    0.1763199425E+01    0.1763234540E+01
    2   -0.8357174395E+00   -0.8356886819E+00

 vdpol
 rodasp
 lapack_decsol

    i      calculated           reference    
    1    0.1763165261E+01    0.1763234540E+01
    2   -0.8357443988E+00   -0.8356886819E+00

 analytical jacobians
 vdpol
 ros2
 lapack_decsol

    i      calculated           reference    
    1    0.1763479136E+01    0.1763234540E+01
    2   -0.8354253142E+00   -0.8356886819E+00

 vdpol
 rose2
 lapack_decsol

    i      calculated           reference    
    1    0.1763259925E+01    0.1763234540E+01
    2   -0.8356492331E+00   -0.8356886819E+00

 vdpol
 ros3p
 lapack_decsol

    i      calculated           reference    
    1    0.1763213341E+01    0.1763234540E+01
    2   -0.8357074211E+00   -0.8356886819E+00

 vdpol
 ros3pl
 lapack_decsol

    i      calculated           reference    
    1    0.1763227587E+01    0.1763234540E+01
    2   -0.8356793412E+00   -0.8356886819E+00

 vdpol
 rodas3
 lapack_decsol

    i      calculated           reference    
    1    0.1763101295E+01    0.1763234540E+01
    2   -0.8358138392E+00   -0.8356886819E+00

 vdpol
 rodas4
 lapack_decsol

    i      calculated           reference    
    1    0.1763199416E+01    0.1763234540E+01
    2   -0.8357174474E+00   -0.8356886819E+00

 vdpol
 rodasp
 lapack_decsol

    i      calculated           reference    
    1    0.1763169482E+01    0.1763234540E+01
    2   -0.8357397820E+00   -0.8356886819E+00

 ijob 2
 numerical jacobians
 medakzo
 ros3pl
 lapack_decsol

    i      calculated           reference    
    1    0.5096063876E-05    0.5113983841E-05
    2    0.1049197001E-03    0.1052886196E-03
    3    0.2148970866E-03    0.2156525325E-03
    4    0.3158064483E-03    0.3169160481E-03
    5    0.3603617132E-03    0.3616265691E-03
    6    0.2573791644E-03    0.2582807380E-03
    7    0.7023531292E-06    0.7048158430E-06
    8    0.8691676699E-15    0.8721783456E-15
    9    0.1733538702E-28    0.1739533020E-28
   10    0.1242770039E-50    0.1247057765E-50
   11    0.6322901542E-88    0.6344600530E-88

 medakzo
 rodas3
 lapack_decsol

    i      calculated           reference    
    1    0.5104405672E-05    0.5113983841E-05
    2    0.1050914400E-03    0.1052886196E-03
    3    0.2152488158E-03    0.2156525325E-03
    4    0.3163232558E-03    0.3169160481E-03
    5    0.3609513128E-03    0.3616265691E-03
    6    0.2578002227E-03    0.2582807380E-03
    7    0.7035312127E-06    0.7048158430E-06
    8    0.8706032202E-15    0.8721783456E-15
    9    0.1736405602E-28    0.1739533020E-28
   10    0.1244849151E-50    0.1247057765E-50
   11    0.6333563088E-88    0.6344600530E-88

 medakzo
 rodas4
 lapack_decsol

    i      calculated           reference    
    1    0.5116551831E-05    0.5113983841E-05
    2    0.1053414765E-03    0.1052886196E-03
    3    0.2157606907E-03    0.2156525325E-03
    4    0.3170746572E-03    0.3169160481E-03
    5    0.3618068193E-03    0.3616265691E-03
    6    0.2584071467E-03    0.2582807380E-03
    7    0.7050702803E-06    0.7048158430E-06
    8    0.8723222968E-15    0.8721783456E-15
    9    0.1739862193E-28    0.1739533020E-28
   10    0.1247763053E-50    0.1247057765E-50
   11    0.6347882530E-88    0.6344600530E-88

 medakzo
 rodasp
 lapack_decsol

    i      calculated           reference    
    1    0.5117637983E-05    0.5113983841E-05
    2    0.1053638327E-03    0.1052886196E-03
    3    0.2158064371E-03    0.2156525325E-03
    4    0.3171417447E-03    0.3169160481E-03
    5    0.3618830662E-03    0.3616265691E-03
    6    0.2584601967E-03    0.2582807380E-03
    7    0.7051594409E-06    0.7048158430E-06
    8    0.8723112216E-15    0.8721783456E-15
    9    0.1739835083E-28    0.1739533020E-28
   10    0.1248053272E-50    0.1247057765E-50
   11    0.6349545606E-88    0.6344600530E-88

 analytical jacobians
 medakzo
 ros3pl
 lapack_decsol

    i      calculated           reference    
    1    0.5096063724E-05    0.5113983841E-05
    2    0.1049196970E-03    0.1052886196E-03
    3    0.2148970802E-03    0.2156525325E-03
    4    0.3158064390E-03    0.3169160481E-03
    5    0.3603617026E-03    0.3616265691E-03
    6    0.2573791569E-03    0.2582807380E-03
    7    0.7023531055E-06    0.7048158430E-06
    8    0.8691676439E-15    0.8721783456E-15
    9    0.1733538656E-28    0.1739533020E-28
   10    0.1242770008E-50    0.1247057765E-50
   11    0.6322901387E-88    0.6344600530E-88

 medakzo
 rodas3
 lapack_decsol

    i      calculated           reference    
    1    0.5104405773E-05    0.5113983841E-05
    2    0.1050914421E-03    0.1052886196E-03
    3    0.2152488200E-03    0.2156525325E-03
    4    0.3163232621E-03    0.3169160481E-03
    5    0.3609513200E-03    0.3616265691E-03
    6    0.2578002277E-03    0.2582807380E-03
    7    0.7035312240E-06    0.7048158430E-06
    8    0.8706032306E-15    0.8721783456E-15
    9    0.1736405623E-28    0.1739533020E-28
   10    0.1244849169E-50    0.1247057765E-50
   11    0.6333563205E-88    0.6344600530E-88

 medakzo
 rodas4
 lapack_decsol

    i      calculated           reference    
    1    0.5116551772E-05    0.5113983841E-05
    2    0.1053414753E-03    0.1052886196E-03
    3    0.2157606882E-03    0.2156525325E-03
    4    0.3170746537E-03    0.3169160481E-03
    5    0.3618068154E-03    0.3616265691E-03
    6    0.2584071440E-03    0.2582807380E-03
    7    0.7050702755E-06    0.7048158430E-06
    8    0.8723222945E-15    0.8721783456E-15
    9    0.1739862194E-28    0.1739533020E-28
   10    0.1247763052E-50    0.1247057765E-50
   11    0.6347882486E-88    0.6344600530E-88

 medakzo
 rodasp
 lapack_decsol

    i      calculated           reference    
    1    0.5117638912E-05    0.5113983841E-05
    2    0.1053638518E-03    0.1052886196E-03
    3    0.2158064761E-03    0.2156525325E-03
    4    0.3171418019E-03    0.3169160481E-03
    5    0.3618831310E-03    0.3616265691E-03
    6    0.2584602410E-03    0.2582807380E-03
    7    0.7051595172E-06    0.7048158430E-06
    8    0.8723112502E-15    0.8721783456E-15
    9    0.1739835032E-28    0.1739533020E-28
   10    0.1248053260E-50    0.1247057765E-50
   11    0.6349546149E-88    0.6344600530E-88

 ijob 3
 chemakzo
 M banded
 ros3pl
 lapack_decsol

                                    dense output for x=    1.0000000000000000D+02
                                                    val           1    1.4223489000473238D-01
                                                    val           2    1.1809782974441193D-03
                                                    val           3    1.4765482578382205D-01
                                                    val           4    5.1825658116203433D-04
                                                    val           5    1.6880751139654090D-02
                                                    val           6    8.5452781977769596D-03


    i      calculated           reference    
    1    0.1150794916E+00    0.1150794921E+00
    2    0.1203831472E-02    0.1203831472E-02
    3    0.1611562890E+00    0.1611562887E+00
    4    0.3656156340E-03    0.3656156421E-03
    5    0.1708010886E-01    0.1708010885E-01
    6    0.4873531240E-02    0.4873531310E-02

 chemakzo
 M banded
 rodas3
 lapack_decsol

                                    dense output for x=    1.0000000000000000D+02
                                                    val           1    1.4223488401777556D-01
                                                    val           2    1.1809784472428832D-03
                                                    val           3    1.4765482843510216D-01
                                                    val           4    5.1825589864234285D-04
                                                    val           5    1.6880750584724696D-02
                                                    val           6    8.5519891664236451D-03


    i      calculated           reference    
    1    0.1150794967E+00    0.1150794921E+00
    2    0.1203831514E-02    0.1203831472E-02
    3    0.1611562863E+00    0.1611562887E+00
    4    0.3656154156E-03    0.3656156421E-03
    5    0.1708010773E-01    0.1708010885E-01
    6    0.4873528478E-02    0.4873531310E-02

 chemakzo
 M banded
 rodas4
 lapack_decsol

                                    dense output for x=    1.0000000000000000D+02
                                                    val           1    1.4223489007110152D-01
                                                    val           2    1.1809784413771640D-03
                                                    val           3    1.4765482581841655D-01
                                                    val           4    5.1825671680107502D-04
                                                    val           5    1.6880752361952064D-02
                                                    val           6    8.5382931643355475D-03


    i      calculated           reference    
    1    0.1150794926E+00    0.1150794921E+00
    2    0.1203831471E-02    0.1203831472E-02
    3    0.1611562885E+00    0.1611562887E+00
    4    0.3656156925E-03    0.3656156421E-03
    5    0.1708011017E-01    0.1708010885E-01
    6    0.4873532008E-02    0.4873531310E-02

 chemakzo
 M banded
 rodasp
 lapack_decsol

                                    dense output for x=    1.0000000000000000D+02
                                                    val           1    1.4223489235744763D-01
                                                    val           2    1.1809788408432978D-03
                                                    val           3    1.4765482430359278D-01
                                                    val           4    5.1825597118094948D-04
                                                    val           5    1.6880754270063521D-02
                                                    val           6    8.5382150913009894D-03


    i      calculated           reference    
    1    0.1150794923E+00    0.1150794921E+00
    2    0.1203831452E-02    0.1203831472E-02
    3    0.1611562887E+00    0.1611562887E+00
    4    0.3656157579E-03    0.3656156421E-03
    5    0.1708011136E-01    0.1708010885E-01
    6    0.4873532869E-02    0.4873531310E-02

 ijob 5
 chemakzo
 M full
 ros3pl
 lapack_decsol

                                    dense output for x=    1.0000000000000000D+02
                                                    val           1    1.4223489000473238D-01
                                                    val           2    1.1809782974441193D-03
                                                    val           3    1.4765482578382205D-01
                                                    val           4    5.1825658116203433D-04
                                                    val           5    1.6880751139654090D-02
                                                    val           6    8.5452781977769596D-03


    i      calculated           reference    
    1    0.1150794916E+00    0.1150794921E+00
    2    0.1203831472E-02    0.1203831472E-02
    3    0.1611562890E+00    0.1611562887E+00
    4    0.3656156340E-03    0.3656156421E-03
    5    0.1708010886E-01    0.1708010885E-01
    6    0.4873531240E-02    0.4873531310E-02

 chemakzo
 M full
 rodas3
 lapack_decsol

                                    dense output for x=    1.0000000000000000D+02
                                                    val           1    1.4223488401777556D-01
                                                    val           2    1.1809784472428832D-03
                                                    val           3    1.4765482843510216D-01
                                                    val           4    5.1825589864234285D-04
                                                    val           5    1.6880750584724696D-02
                                                    val           6    8.5519891664236451D-03


    i      calculated           reference    
    1    0.1150794967E+00    0.1150794921E+00
    2    0.1203831514E-02    0.1203831472E-02
    3    0.1611562863E+00    0.1611562887E+00
    4    0.3656154156E-03    0.3656156421E-03
    5    0.1708010773E-01    0.1708010885E-01
    6    0.4873528478E-02    0.4873531310E-02

 chemakzo
 M full
 rodas4
 lapack_decsol

                                    dense output for x=    1.0000000000000000D+02
                                                    val           1    1.4223489007110152D-01
                                                    val           2    1.1809784413771640D-03
                                                    val           3    1.4765482581841655D-01
                                                    val           4    5.1825671680107502D-04
                                                    val           5    1.6880752361952064D-02
                                                    val           6    8.5382931643355475D-03


    i      calculated           reference    
    1    0.1150794926E+00    0.1150794921E+00
    2    0.1203831471E-02    0.1203831472E-02
    3    0.1611562885E+00    0.1611562887E+00
    4    0.3656156925E-03    0.3656156421E-03
    5    0.1708011017E-01    0.1708010885E-01
    6    0.4873532008E-02    0.4873531310E-02

 chemakzo
 M full
 rodasp
 lapack_decsol

                                    dense output for x=    1.0000000000000000D+02
                                                    val           1    1.4223489235744763D-01
                                                    val           2    1.1809788408432978D-03
                                                    val           3    1.4765482430359278D-01
                                                    val           4    5.1825597118094948D-04
                                                    val           5    1.6880754270063521D-02
                                                    val           6    8.5382150913009894D-03


    i      calculated           reference    
    1    0.1150794923E+00    0.1150794921E+00
    2    0.1203831452E-02    0.1203831472E-02
    3    0.1611562887E+00    0.1611562887E+00
    4    0.3656157579E-03    0.3656156421E-03
    5    0.1708011136E-01    0.1708010885E-01
    6    0.4873532869E-02    0.4873531310E-02

 ijob 11
 ijob 13
 vdpol_x
 M band
 rodas3
 lapack_decsol

    i      calculated           reference    
    1    0.1763224488E+01    0.1763234540E+01
    2   -0.8356888247E+00   -0.8356886819E+00

 vdpol_x
 M band
 rodas4
 lapack_decsol

    i      calculated           reference    
    1    0.1763213794E+01    0.1763234540E+01
    2   -0.8357033727E+00   -0.8356886819E+00

 vdpol_x
 M band
 rodasp
 lapack_decsol

    i      calculated           reference    
    1    0.1763198694E+01    0.1763234540E+01
    2   -0.8357162915E+00   -0.8356886819E+00

 ijob 14
 vdpol3_x
 ros3pl
 lapack_decsol

    i      calculated           reference    
    1    0.1763234137E+01    0.1763234540E+01
    2   -0.8356890527E+00   -0.8356886819E+00
    3    0.1763234137E+01    0.1763234540E+01
    4   -0.8356890527E+00   -0.8356886819E+00
    5    0.1763234137E+01    0.1763234540E+01
    6   -0.8356890527E+00   -0.8356886819E+00

 vdpol3_x
 rodas3
 lapack_decsol

    i      calculated           reference    
    1    0.1763232911E+01    0.1763234540E+01
    2   -0.8356901190E+00   -0.8356886819E+00
    3    0.1763232911E+01    0.1763234540E+01
    4   -0.8356901190E+00   -0.8356886819E+00
    5    0.1763232911E+01    0.1763234540E+01
    6   -0.8356901190E+00   -0.8356886819E+00

 vdpol3_x
 rodas4
 lapack_decsol

    i      calculated           reference    
    1    0.1763234447E+01    0.1763234540E+01
    2   -0.8356887070E+00   -0.8356886819E+00
    3    0.1763234447E+01    0.1763234540E+01
    4   -0.8356887070E+00   -0.8356886819E+00
    5    0.1763234447E+01    0.1763234540E+01
    6   -0.8356887070E+00   -0.8356886819E+00

 vdpol3_x
 rodasp
 lapack_decsol

    i      calculated           reference    
    1    0.1763233698E+01    0.1763234540E+01
    2   -0.8356893877E+00   -0.8356886819E+00
    3    0.1763233698E+01    0.1763234540E+01
    4   -0.8356893877E+00   -0.8356886819E+00
    5    0.1763233698E+01    0.1763234540E+01
    6   -0.8356893877E+00   -0.8356886819E+00

 ijob 15
 vdpol_x
 M full
 ros3pl
 lapack_decsol

    i      calculated           reference    
    1    0.1763227587E+01    0.1763234540E+01
    2   -0.8356793412E+00   -0.8356886819E+00

 vdpol_x
 M full
 rodas3
 lapack_decsol

    i      calculated           reference    
    1    0.1763224488E+01    0.1763234540E+01
    2   -0.8356888247E+00   -0.8356886819E+00

 vdpol_x
 M full
 rodas4
 lapack_decsol

    i      calculated           reference    
    1    0.1763213794E+01    0.1763234540E+01
    2   -0.8357033727E+00   -0.8356886819E+00

 vdpol_x
 M full
 rodasp
 lapack_decsol

    i      calculated           reference    
    1    0.1763198694E+01    0.1763234540E+01
    2   -0.8357162915E+00   -0.8356886819E+00

 diffusion
 ros2
 lapack_decsol

    i      calculated           reference    
    1    0.2583074334E+00    0.2583083212E+00
    2    0.2784691078E+00    0.2784700096E+00
    3    0.3169115295E+00    0.3169124424E+00
    4    0.3700188865E+00    0.3700196743E+00
    5    0.4327350319E+00    0.4327355325E+00
    6    0.4990030027E+00    0.4990030825E+00
    7    0.5623285013E+00    0.5623281002E+00
    8    0.6164218079E+00    0.6164209247E+00
    9    0.6558501736E+00    0.6558489203E+00
   10    0.6766211920E+00    0.6766197085E+00

 diffusion
 rose2
 lapack_decsol

    i      calculated           reference    
    1    0.2583083729E+00    0.2583083212E+00
    2    0.2784700727E+00    0.2784700096E+00
    3    0.3169124843E+00    0.3169124424E+00
    4    0.3700197135E+00    0.3700196743E+00
    5    0.4327355625E+00    0.4327355325E+00
    6    0.4990030826E+00    0.4990030825E+00
    7    0.5623280516E+00    0.5623281002E+00
    8    0.6164208531E+00    0.6164209247E+00
    9    0.6558488314E+00    0.6558489203E+00
   10    0.6766196422E+00    0.6766197085E+00

 diffusion
 ros3p
 lapack_decsol

    i      calculated           reference    
    1    0.2601107624E+00    0.2583083212E+00
    2    0.2801170348E+00    0.2784700096E+00
    3    0.3165286902E+00    0.3169124424E+00
    4    0.3695291167E+00    0.3700196743E+00
    5    0.4333315497E+00    0.4327355325E+00
    6    0.4991698879E+00    0.4990030825E+00
    7    0.5622495182E+00    0.5623281002E+00
    8    0.6180461065E+00    0.6164209247E+00
    9    0.6541840052E+00    0.6558489203E+00
   10    0.6718386088E+00    0.6766197085E+00

 diffusion
 ros3pl
 lapack_decsol

    i      calculated           reference    
    1    0.2583082375E+00    0.2583083212E+00
    2    0.2784702466E+00    0.2784700096E+00
    3    0.3169130787E+00    0.3169124424E+00
    4    0.3700205718E+00    0.3700196743E+00
    5    0.4327363566E+00    0.4327355325E+00
    6    0.4990034992E+00    0.4990030825E+00
    7    0.5623279410E+00    0.5623281002E+00
    8    0.6164202581E+00    0.6164209247E+00
    9    0.6558479094E+00    0.6558489203E+00
   10    0.6766185679E+00    0.6766197085E+00

 diffusion
 rodas3
 lapack_decsol

    i      calculated           reference    
    1    0.2583082622E+00    0.2583083212E+00
    2    0.2784701200E+00    0.2784700096E+00
    3    0.3169127477E+00    0.3169124424E+00
    4    0.3700201165E+00    0.3700196743E+00
    5    0.4327359432E+00    0.4327355325E+00
    6    0.4990032878E+00    0.4990030825E+00
    7    0.5623280104E+00    0.5623281002E+00
    8    0.6164205893E+00    0.6164209247E+00
    9    0.6558484221E+00    0.6558489203E+00
   10    0.6766191675E+00    0.6766197085E+00

 diffusion
 rodas4
 lapack_decsol

    i      calculated           reference    
    1    0.2583084842E+00    0.2583083212E+00
    2    0.2784700767E+00    0.2784700096E+00
    3    0.3169123482E+00    0.3169124424E+00
    4    0.3700194995E+00    0.3700196743E+00
    5    0.4327353817E+00    0.4327355325E+00
    6    0.4990030194E+00    0.4990030825E+00
    7    0.5623281154E+00    0.5623281002E+00
    8    0.6164209899E+00    0.6164209247E+00
    9    0.6558489780E+00    0.6558489203E+00
   10    0.6766197736E+00    0.6766197085E+00

 diffusion
 rodasp
 lapack_decsol

    i      calculated           reference    
    1    0.2583085103E+00    0.2583083212E+00
    2    0.2784700857E+00    0.2784700096E+00
    3    0.3169123343E+00    0.3169124424E+00
    4    0.3700194713E+00    0.3700196743E+00
    5    0.4327353556E+00    0.4327355325E+00
    6    0.4990030080E+00    0.4990030825E+00
    7    0.5623281204E+00    0.5623281002E+00
    8    0.6164210039E+00    0.6164209247E+00
    9    0.6558489923E+00    0.6558489203E+00
   10    0.6766197850E+00    0.6766197085E+00


 test mebdfi pollu

     1  5.6462559625294388E-02
     2  1.3424840790872916E-01
     4  5.5231408900399265E-03
     7  7.7842488273513560E-02
     8  3.2450753898265677E-01
     9  7.4940126390311366E-03
    12  2.2305066926244974E-03
    13  2.0871631301121632E-04
    14  1.3969206917218856E-05
    15  8.9648851080745932E-03
    17  6.8992196888701257E-03
    18  1.0078031112989498E-04
    19  1.7721468630894848E-06
    20  5.6829449013372444E-05

 test mebdfi medakzo

     1  5.1139896910174527E-06
    19  1.0528873944004347E-04
    37  2.1565277382420181E-04
    55  3.1691638935377929E-04
    73  3.6162693081377836E-04
    91  2.5828097182293799E-04
   109  7.0481666414038269E-07


