test_sparse1.f90 running
 expect error
 sparse get error i=           0 , j=           1
 expect error
 sparse get error i=           1 , j=           0
 expect error
 sparse get error i=           5 , j=           1
 expect error
 sparse get error i=           1 , j=           6
(   1,   1)=   1.000000E+00
(   1,   3)=   1.300000E+00
(   1,   5)=   1.500000E+00
(   2,   2)=   2.000000E+00
(   2,   4)=   4.000000E+00
(   2,   5)=   2.500000E+00
 change 1,3 to 3.3
(   1,   1)=   1.000000E+00
(   1,   3)=   3.300000E+00
(   1,   5)=   1.500000E+00
(   2,   2)=   2.000000E+00
(   2,   4)=   4.000000E+00
(   2,   5)=   2.500000E+00
 change 1,1 to 1.1
(   1,   1)=   1.100000E+00
(   1,   3)=   3.300000E+00
(   1,   5)=   1.500000E+00
(   2,   2)=   2.000000E+00
(   2,   4)=   4.000000E+00
(   2,   5)=   2.500000E+00
 insert 1,2,2.2  between
(   1,   1)=   1.100000E+00
(   1,   2)=   2.200000E+00
(   1,   3)=   3.300000E+00
(   1,   5)=   1.500000E+00
(   2,   2)=   2.000000E+00
(   2,   4)=   4.000000E+00
(   2,   5)=   2.500000E+00
 insert at head 2,1,2.1
(   1,   1)=   1.100000E+00
(   1,   2)=   2.200000E+00
(   1,   3)=   3.300000E+00
(   1,   5)=   1.500000E+00
(   2,   1)=   2.100000E+00
(   2,   2)=   2.000000E+00
(   2,   4)=   4.000000E+00
(   2,   5)=   2.500000E+00
 get A(1,1) =    1.10000000000000     
 get A(2,5) =    2.50000000000000     
 fill in RHS
(   1,   1)=   1.100000E+00
(   1,   2)=   2.200000E+00
(   1,   3)=   3.300000E+00
(   1,   5)=   1.500000E+00
(   2,   1)=   2.100000E+00
(   2,   2)=   2.000000E+00
(   2,   4)=   4.000000E+00
(   2,   5)=   2.500000E+00
(   3,   3)=   1.000000E+00
(   3,   5)=   3.500000E+00
(   4,   4)=   1.000000E+00
(   4,   5)=   4.500000E+00
 get RHS  using getrhs
RHS=  1.500  2.500  3.500  4.500
 export matrix to non sparse form
row   1=  1.100  2.200  3.300  0.000  1.500
row   2=  2.100  2.000  0.000  4.000  2.500
row   3=  0.000  0.000  1.000  0.000  3.500
row   4=  0.000  0.000  0.000  1.000  4.500
 check putrhs
(   1,   1)=   1.100000E+00
(   1,   2)=   2.200000E+00
(   1,   3)=   3.300000E+00
(   1,   5)=  -1.500000E+00
(   2,   1)=   2.100000E+00
(   2,   2)=   2.000000E+00
(   2,   4)=   4.000000E+00
(   2,   5)=  -2.500000E+00
(   3,   3)=   1.000000E+00
(   3,   5)=  -3.500000E+00
(   4,   4)=   1.000000E+00
(   4,   5)=  -4.500000E+00
 check import
(   1,   1)=  -1.100000E+00
(   1,   2)=  -2.200000E+00
(   1,   3)=  -3.300000E+00
(   1,   5)=  -1.500000E+00
(   2,   1)=  -2.100000E+00
(   2,   2)=  -2.000000E+00
(   2,   4)=  -4.000000E+00
(   2,   5)=  -2.500000E+00
(   3,   3)=  -1.000000E+00
(   3,   5)=  -3.500000E+00
(   4,   4)=  -1.000000E+00
(   4,   5)=  -4.500000E+00
 check multiply
Yout= 18.700 26.150  3.500  4.500
 check add(A,1,1,98.0)    96.9000000000000     
 solve simultaneous equations
 row(k)=           1 , pivot=   96.9000000000000     
 simeq divide k=           1 , row(k)=           1 , j=           2
 simeq divide k=           1 , row(k)=           1 , j=           3
 simeq divide k=           1 , row(k)=           1 , j=           5
 kval=  -2.10000000000000     
 simeq reduce k=           1 , j=           2  kjval= -2.270381836945305E-002
 simeq reduce k=           1 , j=           4  kjval=   0.00000000000000     
 simeq reduce k=           1 , j=           5  kjval= -1.547987616099071E-002
 kval=   0.00000000000000     
 kval=   0.00000000000000     
 row(k)=           2 , pivot=  -2.04767801857585     
 simeq divide k=           2 , row(k)=           2 , j=           4
 simeq divide k=           2 , row(k)=           2 , j=           5
 kval= -2.270381836945305E-002
 simeq reduce k=           2 , j=           3  kjval=   0.00000000000000     
 simeq reduce k=           2 , j=           5  kjval=   1.23677048684608     
 kval=   0.00000000000000     
 kval=   0.00000000000000     
 row(k)=           3 , pivot=  -1.00000000000000     
 simeq divide k=           3 , row(k)=           3 , j=           5
 kval= -3.405572755417956E-002
 simeq reduce k=           3 , j=           5  kjval=   3.50000000000000     
 kval=   0.00000000000000     
 kval=   0.00000000000000     
 row(k)=           4 , pivot=  -1.00000000000000     
 simeq divide k=           4 , row(k)=           4 , j=           5
 kval=   0.00000000000000     
 kval=   1.95343211369822     
 simeq reduce k=           4 , j=           5  kjval=   4.50000000000000     
 kval=   0.00000000000000     
simeq=  0.132 -7.554  3.500  4.500
 check simeq A * X = Yout against Y
 check Y,Yout,err  -1.50000000000000        17.8389779256123        19.3389779256123     
 check Y,Yout,err  -2.50000000000000       -3.16942057423928       0.669420574239277     
 check Y,Yout,err  -3.50000000000000       -3.50000000000000        0.00000000000000     
 check Y,Yout,err  -4.50000000000000       -4.50000000000000        0.00000000000000     
simeq = -1.500 -2.500 -3.500 -4.500
 test_sparse1.f90 ending