
  NOTE: this test should finish very quickly.
  Old versions of gfortran had a bug that made this loop forever,
  so if it seems to be stuck here, that might be the problem.

  test_square1

 -0.162154E+00   0.163913E+01  -0.384542E-01  -0.281585E+01

 rcond   2.0883392070603853E-002
 rcond   1.0000000000000000     
 test_square2

  0.48858D+00 -0.71219D-01  0.74908D+00 -0.10851D+01  0.10000D+01  0.17975D+02

 test_square_inv

        B1  0.10000D+01  0.20000D+01  0.30000D+01
 A1_init*X  0.10000D+01  0.20000D+01  0.30000D+01
  fac A1*X  0.10000D+01  0.20000D+01  0.30000D+01
     ipiv1            2            2            3

        B2  0.11000D+01  0.21000D+01  0.31000D+01
 A2_init*X  0.11000D+01  0.21000D+01  0.31000D+01
  fac A2*X  0.11000D+01  0.21000D+01  0.31000D+01
     ipiv2            2            3            3

  test_square_quad1

 -0.162154E+00   0.163913E+01  -0.384542E-01  -0.281585E+01

 test test_square_quad2

  0.48858D+00 -0.71219D-01  0.74908D+00 -0.10851D+01  0.10000D+01  0.17975D+02

 test_square_quad_inv

        B1  0.10000D+01  0.20000D+01  0.30000D+01
     ipiv1            2            2            3

        B2  0.11000D+01  0.21000D+01  0.31000D+01
 test_diffusion
               banded: init x, soln x(k)           1    3.3853553765573918D-01    4.9000759277093642D-01
               banded: init x, soln x(k)           2    3.3853553765572147D-01    4.9000759277093647D-01
               banded: init x, soln x(k)           3    4.5462481641995478D-01    4.9000759277093636D-01
               banded: init x, soln x(k)           4    4.9630359938855395D-01    4.9000759277093653D-01
               banded: init x, soln x(k)           5    4.9630359938855390D-01    4.9000759277093642D-01
               banded: init x, soln x(k)           6    4.9630359938855384D-01    4.9000759277093642D-01
               banded: init x, soln x(k)           7    4.9630359938855390D-01    4.9000759277093647D-01
               banded: init x, soln x(k)           8    5.8265967459668622D-01    4.9000759277093647D-01
               banded: init x, soln x(k)           9    6.1370782322384976D-01    4.9000759277093642D-01
               banded: init x, soln x(k)          10    6.1370782322384976D-01    4.9000759277093647D-01

                             init_xtotal    9.1839958264200943D+26
                           banded_xtotal    9.1839958264200943D+26

  test banded matrix2

 rcond_banded   1.0000000000000000     
 rcond_banded   1.4544324732842015E-009

 block_thomas_dble
 done block_thomas_dble

 test_block_dc_mt_dble
 done test_block_dc_mt_dble

 block_thomas_klu
 done block_thomas_klu

 test_block_dc_mt_klu
 done test_block_dc_mt_klu

 block_thomas_quad
 done block_thomas_quad

 test_block_dc_mt_quad
 done test_block_dc_mt_quad


 test PARDISO

 solve completed ... 
 the solution of the system is 
  b(           1 ) =  -0.52232142857142805     
  b(           2 ) =   -8.9285714285714159E-003
  b(           3 ) =    1.2209821428571399     
  b(           4 ) =  -0.50446428571428603     
  b(           5 ) =  -0.21428571428571400     
 the expected solution of the system is
 b(           1 ) =  -0.522321428571428
 b(           2 ) =  -8.928571428571416E-003
 b(           3 ) =    1.22098214285714
 b(           4 ) =  -0.504464285714286
 b(           5 ) =  -0.214285714285714

 test klu

                                    b(i)           1    1.8845682243899176D+02
                                    b(i)           2    1.3396799124371137D+02
                                    b(i)           3   -4.7023881375167809D+02
                                    b(i)           4   -2.7880340370786092D+02
                                    b(i)           5    1.9917307938954053D+01
                                    b(i)           6    2.7277269083325808D+02
                                    b(i)           7   -2.4780664233746421D+02
                                    b(i)           8   -3.1399766868523204D+02
                                    b(i)           9   -9.1277214458844540D+01
                                    b(i)          10    9.9759499750981945D+01

 
 test_format_conversion
 dense_to_row_sparse
 find_loc_in_row_sparse
 dense_to_column_sparse
 find_loc_in_column_sparse
 dense_to_band
 band_to_dense
 band_to_column_sparse
 column_sparse_to_band
 band_to_row_sparse
 row_sparse_to_band
 okay

 test_tridiag
           1  0.66666666666666663     
           2   2.0000000000000000     
           3   3.6666666666666670     
           4   5.0000000000000009     
           5   4.6666666666666670     

           2 lhs   1.0000000000000000     
           2 rhs   1.0000000000000000     

           3 lhs   1.9999999999999991     
           3 rhs   2.0000000000000000     

           4 lhs   3.0000000000000027     
           4 rhs   3.0000000000000000     


 test_quad_tridiag
           1  0.666666666666666666666666666666666635      
           2   2.00000000000000000000000000000000000      
           3   3.66666666666666666666666666666666692      
           4   5.00000000000000000000000000000000077      
           5   4.66666666666666666666666666666666692      


 test_tri_solve
           1   3.5238095238095242     
           2   5.2857142857142865     
           3   5.6666666666666679     
           4   4.8571428571428577     
           5   2.9523809523809526     

