Fix CUDA implementation of %set_scal and %zero

cuda/psb_c_cuda_vect_mod.F90

@@ -668,8 +668,7 @@ contains
     implicit none 
     class(psb_c_vect_cuda), intent(inout) :: x
     
-    if (allocated(x%v)) x%v=czero
+    call x%set_scal(czero)
-    call x%set_host()
   end subroutine c_cuda_zero
 
   subroutine c_cuda_asb_m(n, x, info)
@@ -807,7 +806,6 @@ contains
     if (present(first)) first_ = max(1,first)
     if (present(last))  last_  = min(last,last_)
     
-    if (x%is_host()) call x%sync()
     info = setScalDevice(val,first_,last_,1,x%deviceVect)
     call x%set_dev()
     

cuda/psb_d_cuda_vect_mod.F90

@@ -668,8 +668,7 @@ contains
     implicit none 
     class(psb_d_vect_cuda), intent(inout) :: x
     
-    if (allocated(x%v)) x%v=dzero
+    call x%set_scal(dzero)
-    call x%set_host()
   end subroutine d_cuda_zero
 
   subroutine d_cuda_asb_m(n, x, info)
@@ -807,7 +806,6 @@ contains
     if (present(first)) first_ = max(1,first)
     if (present(last))  last_  = min(last,last_)
     
-    if (x%is_host()) call x%sync()
     info = setScalDevice(val,first_,last_,1,x%deviceVect)
     call x%set_dev()
     

cuda/psb_i_cuda_vect_mod.F90

@@ -651,8 +651,7 @@ contains
     implicit none 
     class(psb_i_vect_cuda), intent(inout) :: x
     
-    if (allocated(x%v)) x%v=izero
+    call x%set_scal(izero)
-    call x%set_host()
   end subroutine i_cuda_zero
 
   subroutine i_cuda_asb_m(n, x, info)
@@ -790,7 +789,6 @@ contains
     if (present(first)) first_ = max(1,first)
     if (present(last))  last_  = min(last,last_)
     
-    if (x%is_host()) call x%sync()
     info = setScalDevice(val,first_,last_,1,x%deviceVect)
     call x%set_dev()
     

cuda/psb_s_cuda_vect_mod.F90

@@ -668,8 +668,7 @@ contains
     implicit none 
     class(psb_s_vect_cuda), intent(inout) :: x
     
-    if (allocated(x%v)) x%v=szero
+    call x%set_scal(szero)
-    call x%set_host()
   end subroutine s_cuda_zero
 
   subroutine s_cuda_asb_m(n, x, info)
@@ -807,7 +806,6 @@ contains
     if (present(first)) first_ = max(1,first)
     if (present(last))  last_  = min(last,last_)
     
-    if (x%is_host()) call x%sync()
     info = setScalDevice(val,first_,last_,1,x%deviceVect)
     call x%set_dev()
     

cuda/psb_z_cuda_vect_mod.F90

@@ -668,8 +668,7 @@ contains
     implicit none 
     class(psb_z_vect_cuda), intent(inout) :: x
     
-    if (allocated(x%v)) x%v=zzero
+    call x%set_scal(zzero)
-    call x%set_host()
   end subroutine z_cuda_zero
 
   subroutine z_cuda_asb_m(n, x, info)
@@ -807,7 +806,6 @@ contains
     if (present(first)) first_ = max(1,first)
     if (present(last))  last_  = min(last,last_)
     
-    if (x%is_host()) call x%sync()
     info = setScalDevice(val,first_,last_,1,x%deviceVect)
     call x%set_dev()