[ADD] Start of a spmm testing file, only supporting file read from MM format and randomly generating x vector

test/spmm/Makefile

@@ -0,0 +1,35 @@
+INSTALLDIR=../..
+INCDIR=$(INSTALLDIR)/include/
+MODDIR=$(INSTALLDIR)/modules/
+include $(INCDIR)/Make.inc.psblas
+#
+# Libraries used
+#
+LIBDIR=$(INSTALLDIR)/lib/
+PSBLAS_LIB= -L$(LIBDIR) -lpsb_util -lpsb_prec -lpsb_base
+LDLIBS=$(PSBLDLIBS)
+
+FINCLUDES=$(FMFLAG)$(MODDIR) $(FMFLAG).
+
+OBJS=psb_spmm_test.o   
+
+EXEDIR=./runs
+
+all: runsd psb_spmm_test
+
+runsd:
+	(if test ! -d runs ; then mkdir runs; fi)
+
+psb_spmm_test: $(OBJS)
+	$(FLINK) $(LOPT) $(OBJS) -o psb_spmm_test $(PSBLAS_LIB) $(LDLIBS) 
+	/bin/mv psb_spmm_test $(EXEDIR) 
+
+clean: 
+	/bin/rm -f $(OBJS)\
+        *$(.mod) $(EXEDIR)/psb_spmm_test 
+
+lib:
+	(cd ../../; make library)
+verycleanlib:
+	(cd ../../; make veryclean)
+

test/spmm/README.md

@@ -0,0 +1,4 @@
+# Introduction
+This is a directory developed by Luca Pepè Sciarria and Simone Staccone froma Tor Vergata University to start to create some unit tests for PSBLAS 3.9, in particular for psb_spmm routine.
+
+## Test Suite

test/spmm/psb_spmm_test.f90

@@ -0,0 +1,109 @@
+! Test program for y = AX spmm routine
+! Author: Luca Pepé Sciarria, Staccone Simone (Tor Vergata University)
+program psb_spmm_test
+    use psb_base_mod
+    use psb_prec_mod
+    use psb_util_mod
+
+    implicit none 
+
+    ! input parameters
+    character(len=40)           :: mtx_file, file_name, name
+
+    ! sparse matrices
+    type(psb_sspmat_type)       :: a
+
+    ! vectors
+    type(psb_s_vect_type)       :: x, y
+
+    ! matrix descriptor data structure
+    type(psb_desc_type)         :: desc_a
+
+    ! communication context
+    type(psb_ctxt_type)         :: ctxt
+    integer(psb_ipk_)           :: my_rank, np, info, err_act
+
+    ! matrix parameters
+    integer(psb_ipk_)           :: m, n, nnz
+
+
+    real(psb_spk_), allocatable :: x_global(:)
+    integer(psb_ipk_)           :: i
+
+    name     = "psb_spmm_test"    ! Name of the program to output in case of error
+    info     = psb_success_
+    mtx_file = "matrix/1138_bus.mtx"
+
+
+    call psb_init(ctxt)
+    call psb_info(ctxt,my_rank,np)
+
+    if (my_rank < 0) then 
+        ! This should not happen, but just in case
+        call psb_error(ctxt) 
+    endif
+
+    if (my_rank == psb_root_) then 
+        write(psb_out_unit,*) 'Welcome to PSBLAS version: ',psb_version_string_
+        write(psb_out_unit,*) 'This is the ',trim(name),' sample program'
+    end if
+
+
+    call mm_mat_read(a,info,filename=mtx_file)
+    
+    m   = a%get_nrows()
+    n   = a%get_ncols()
+    nnz = a%get_nzeros()
+
+    call psb_bcast(ctxt,m)
+    call psb_bcast(ctxt,n)
+    call psb_bcast(ctxt,nnz)
+
+    if(my_rank == psb_root_) then 
+        write(psb_out_unit,*) "Matrix stats"
+        write(psb_out_unit,*) "ROWS:", m
+        write(psb_out_unit,*) "COLS:", n
+        write(psb_out_unit,*) "NNZ: ", nnz
+    end if
+
+    ! Generate random array for b using always the same seed
+    if(my_rank == psb_root_) then
+        call random_init(repeatable=.true.,image_distinct=.false.) 
+        allocate(x_global(n))
+        do i=0,n
+            call random_number(x_global(i))
+        end do
+        write(psb_out_unit, *) x_global
+    end if
+
+    ! call psb_matdist to initialize desc_a
+
+    ! call psb_scatter(x_global,x,desc_a,info,root=psb_root_)
+    ! call psb_geall(y,desc_a,info)
+
+
+    ! call psb_spfree(a, desc_a,info)
+    ! if(info /= psb_success_) then
+    !     write(psb_out_unit,*) "Error in matrix A free routine"
+    !     goto 9999
+    ! end if
+    
+    ! call psb_cdfree(desc_a,info)
+    ! if(info /= psb_success_) then
+    !     write(psb_out_unit,*) "Error in matrix descriptor free routine"
+    !     goto 9999
+    ! end if
+
+    call psb_exit(ctxt)
+    stop
+
+
+    ! Error handling
+    9999 call psb_error(ctxt)
+    
+    call psb_errpush(info,name)
+    call psb_error_handler(ctxt,err_act)
+    call psb_exit(ctxt)
+    stop
+
+end program