program test_sparse
use fvn_sparse
implicit none
integer(kind=ip_kind), parameter :: nz=12
integer(kind=ip_kind), parameter :: n=5
real(kind=dp_kind),dimension(nz) :: A
real(kind=dp_kind),dimension(n,n) :: As
integer(kind=ip_kind),dimension(nz) :: Ti,Tj
real(kind=dp_kind),dimension(n) :: B,x
integer(kind=ip_kind) :: status,i
! Description of the matrix in triplet form
A = (/ 2.,3.,3.,-1.,4.,4.,-3.,1.,2.,2.,6.,1. /)
B = (/ 8., 45., -3., 3., 19./)
Ti = (/ 1,2,1,3,5,2,3,4,5,3,2,5 /)
Tj = (/ 1,1,2,2,2,3,3,3,3,4,5,5 /)

! Reconstruction of the matrix in standard form
! just needed for printing the matrix here
As=0.
do i=1,nz
    As(Ti(i),Tj(i))=A(i)
end do
write(*,*) "Matrix in standard representation :"
do i=1,5
    write(*,'(5f8.4)') As(i,:)
end do
write(*,*)
write(*,'("Right hand side :",5f8.4)') B

!specific routine that will be used here
!call fvn_di_sparse_solve(n,nz,A,Ti,Tj,B,x,status)
call fvn_di_sparse_solve(n,nz,A,Ti,Tj,B,x,status)
write(*,'("Solution :",5f8.4)') x
write(*,'("Product matrix Solution :",5f8.4)') matmul(As,x)
end program