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