test_sparse.f90 running
(   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
simeq= -0.071 -7.676  3.500  4.500
 check simeq A * X = Yout against Y
 check Y,Yout,err  -1.50000000000000       -1.50000000000000       4.440892098500626E-016
 check Y,Yout,err  -2.50000000000000       -2.50000000000000       1.776356839400250E-015
 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_sparse.f90 ending