psblas3

base/modules/psb_glist_map_mod.f90
 test/kernel/d_file_spmv.f90
 test/kernel/runs/spmv.inp
 util/psb_mat_dist_impl.f90


Fixed idiotic bug in glist

base/modules/psb_glist_map_mod.f90

@@ -168,7 +168,7 @@ contains
 
     ngp = size(idxmap%vgp)
     do i=1, nv 
-      if ((1<=idx(i)).and.(idx(i)<ngp)) then
+      if ((1<=idx(i)).and.(idx(i)<=ngp)) then
         iprc(i) = idxmap%vgp(idx(i))
       else
         iprc(i) = -1

test/kernel/d_file_spmv.f90

@@ -195,7 +195,7 @@ program d_file_spmv
     call psb_matdist(aux_a, a, ictxt, &
          & desc_a,b_col_glob,b_col,info,fmt=afmt,v=ivg)
   else 
-    if (iam==psb_root_) write(psb_out_unit,'("Partition type: block")')
+    if (iam==psb_root_) write(psb_out_unit,'("Partition type: default block")')
     call psb_matdist(aux_a, a,  ictxt, &
          & desc_a,b_col_glob,b_col,info,fmt=afmt,parts=part_block)
   end if
@@ -270,6 +270,7 @@ program d_file_spmv
     write(psb_out_unit,'("MBYTES/S                         : ",F20.3)') bdwdth
     bdwdth = times*nbytes/(tt2*1.d6)
     write(psb_out_unit,'("MBYTES/S                  (trans): ",F20.3)') bdwdth
+    write(psb_out_unit,'("Storage type for DESC_A: ",a)') desc_a%indxmap%get_fmt()
     
   end if
 

test/kernel/runs/spmv.inp

@@ -1,3 +1,3 @@
-A_1M_gps.mtx   
+ASIC_100ks.mtx 
 MM
 0

util/psb_mat_dist_impl.f90

@@ -618,8 +618,10 @@ subroutine dmatdist(a_glob, a, ictxt, desc_a,&
   endif
   if (use_parts) then 
     call psb_cdall(ictxt,desc_a,info,mg=nrow,parts=parts)
-  else 
+  else if (use_v) then 
     call psb_cdall(ictxt,desc_a,info,vg=v)
+  else
+    info = -1
   end if
 
   if(info /= psb_success_) then
@@ -705,6 +707,7 @@ subroutine dmatdist(a_glob, a, ictxt, desc_a,&
           ll = ll + nz
         end do
 !!$        write(0,*) 'mat_dist: sending rows ',i_count,j_count-1,' to proc',iproc, ll
+!!$        write(0,*) 'mat_dist: sending ',irow(ll),icol(ll),val(ll )
         if (iproc == iam) then
           call psb_spins(ll,irow,icol,val,a,desc_a,info)
           if(info /= psb_success_) then