Mark matrix on host after clean_zeros

base/serial/impl/psb_c_csc_impl.f90

@@ -2500,6 +2500,7 @@ subroutine  psb_c_csc_clean_zeros(a, info)
   integer(psb_ipk_), allocatable :: ilcp(:) 
   
   info = 0
+  call a%sync()
   nc   = a%get_ncols()
   ilcp = a%icp(:) 
   a%icp(1) = 1
@@ -2515,6 +2516,7 @@ subroutine  psb_c_csc_clean_zeros(a, info)
     a%icp(i+1) = j
   end do
   call a%trim()
+  call a%set_host()
 end subroutine psb_c_csc_clean_zeros
 
 

base/serial/impl/psb_c_csr_impl.f90

@@ -1304,7 +1304,7 @@ function psb_c_csr_csnmi(a) result(res)
   if (a%is_dev())   call a%sync()
 
   do i = 1, a%get_nrows()
-    acc = dzero
+    acc = szero
     do j=a%irp(i),a%irp(i+1)-1  
       acc = acc + abs(a%val(j))
     end do
@@ -3332,6 +3332,7 @@ subroutine  psb_c_csr_clean_zeros(a, info)
   integer(psb_ipk_), allocatable :: ilrp(:) 
   
   info = 0
+  call a%sync()
   nr   = a%get_nrows()
   ilrp = a%irp(:) 
   a%irp(1) = 1
@@ -3347,6 +3348,7 @@ subroutine  psb_c_csr_clean_zeros(a, info)
     a%irp(i+1) = j
   end do
   call a%trim()
+  call a%set_host()
 end subroutine psb_c_csr_clean_zeros
 
 subroutine psb_ccsrspspmm(a,b,c,info)

base/serial/impl/psb_d_csc_impl.f90

@@ -2500,6 +2500,7 @@ subroutine  psb_d_csc_clean_zeros(a, info)
   integer(psb_ipk_), allocatable :: ilcp(:) 
   
   info = 0
+  call a%sync()
   nc   = a%get_ncols()
   ilcp = a%icp(:) 
   a%icp(1) = 1
@@ -2515,6 +2516,7 @@ subroutine  psb_d_csc_clean_zeros(a, info)
     a%icp(i+1) = j
   end do
   call a%trim()
+  call a%set_host()
 end subroutine psb_d_csc_clean_zeros
 
 

base/serial/impl/psb_d_csr_impl.f90

@@ -3332,6 +3332,7 @@ subroutine  psb_d_csr_clean_zeros(a, info)
   integer(psb_ipk_), allocatable :: ilrp(:) 
   
   info = 0
+  call a%sync()
   nr   = a%get_nrows()
   ilrp = a%irp(:) 
   a%irp(1) = 1
@@ -3347,6 +3348,7 @@ subroutine  psb_d_csr_clean_zeros(a, info)
     a%irp(i+1) = j
   end do
   call a%trim()
+  call a%set_host()
 end subroutine psb_d_csr_clean_zeros
 
 subroutine psb_dcsrspspmm(a,b,c,info)

base/serial/impl/psb_s_csc_impl.f90

@@ -2500,6 +2500,7 @@ subroutine  psb_s_csc_clean_zeros(a, info)
   integer(psb_ipk_), allocatable :: ilcp(:) 
   
   info = 0
+  call a%sync()
   nc   = a%get_ncols()
   ilcp = a%icp(:) 
   a%icp(1) = 1
@@ -2515,6 +2516,7 @@ subroutine  psb_s_csc_clean_zeros(a, info)
     a%icp(i+1) = j
   end do
   call a%trim()
+  call a%set_host()
 end subroutine psb_s_csc_clean_zeros
 
 

base/serial/impl/psb_s_csr_impl.f90

@@ -1304,7 +1304,7 @@ function psb_s_csr_csnmi(a) result(res)
   if (a%is_dev())   call a%sync()
 
   do i = 1, a%get_nrows()
-    acc = dzero
+    acc = szero
     do j=a%irp(i),a%irp(i+1)-1  
       acc = acc + abs(a%val(j))
     end do
@@ -3332,6 +3332,7 @@ subroutine  psb_s_csr_clean_zeros(a, info)
   integer(psb_ipk_), allocatable :: ilrp(:) 
   
   info = 0
+  call a%sync()
   nr   = a%get_nrows()
   ilrp = a%irp(:) 
   a%irp(1) = 1
@@ -3347,6 +3348,7 @@ subroutine  psb_s_csr_clean_zeros(a, info)
     a%irp(i+1) = j
   end do
   call a%trim()
+  call a%set_host()
 end subroutine psb_s_csr_clean_zeros
 
 subroutine psb_scsrspspmm(a,b,c,info)

base/serial/impl/psb_z_csc_impl.f90

@@ -2500,6 +2500,7 @@ subroutine  psb_z_csc_clean_zeros(a, info)
   integer(psb_ipk_), allocatable :: ilcp(:) 
   
   info = 0
+  call a%sync()
   nc   = a%get_ncols()
   ilcp = a%icp(:) 
   a%icp(1) = 1
@@ -2515,6 +2516,7 @@ subroutine  psb_z_csc_clean_zeros(a, info)
     a%icp(i+1) = j
   end do
   call a%trim()
+  call a%set_host()
 end subroutine psb_z_csc_clean_zeros
 
 

base/serial/impl/psb_z_csr_impl.f90

@@ -3332,6 +3332,7 @@ subroutine  psb_z_csr_clean_zeros(a, info)
   integer(psb_ipk_), allocatable :: ilrp(:) 
   
   info = 0
+  call a%sync()
   nr   = a%get_nrows()
   ilrp = a%irp(:) 
   a%irp(1) = 1
@@ -3347,6 +3348,7 @@ subroutine  psb_z_csr_clean_zeros(a, info)
     a%irp(i+1) = j
   end do
   call a%trim()
+  call a%set_host()
 end subroutine psb_z_csr_clean_zeros
 
 subroutine psb_zcsrspspmm(a,b,c,info)