diff --git a/base/serial/impl/psb_c_base_mat_impl.F90 b/base/serial/impl/psb_c_base_mat_impl.F90
index 17f2cdc8..f4bc43cc 100644
--- a/base/serial/impl/psb_c_base_mat_impl.F90
+++ b/base/serial/impl/psb_c_base_mat_impl.F90
@@ -2006,8 +2006,8 @@ subroutine psb_c_base_vect_mv(alpha,a,x,beta,y,info,trans)
 
   ! For the time being we just throw everything back
   ! onto the normal routines.
-  call x%sync()
-  call y%sync()
+  if (x%is_dev()) call x%sync()
+  if (y%is_dev()) call y%sync()
   call a%spmm(alpha,x%v,beta,y%v,info,trans)
   call y%set_host()
 end subroutine psb_c_base_vect_mv
@@ -2060,8 +2060,6 @@ subroutine psb_c_base_vect_cssv(alpha,a,x,beta,y,info,trans,scale,d)
     goto 9999
   end if
 
-  call x%sync()
-  call y%sync()
   if (present(d)) then
     call d%sync()
     if (present(scale)) then
@@ -2161,8 +2159,11 @@ subroutine psb_c_base_inner_vect_sv(alpha,a,x,beta,y,info,trans)
 
   info  = psb_success_
   call psb_erractionsave(err_act)
+  if (x%is_dev()) call x%sync()
+  if (y%is_dev()) call y%sync()
 
   call a%inner_spsm(alpha,x%v,beta,y%v,info,trans)
+  call y%set_host()
 
   if (info /= psb_success_) then
     info = psb_err_from_subroutine_
diff --git a/base/serial/impl/psb_d_base_mat_impl.F90 b/base/serial/impl/psb_d_base_mat_impl.F90
index 69112529..1a8dc084 100644
--- a/base/serial/impl/psb_d_base_mat_impl.F90
+++ b/base/serial/impl/psb_d_base_mat_impl.F90
@@ -2006,8 +2006,8 @@ subroutine psb_d_base_vect_mv(alpha,a,x,beta,y,info,trans)
 
   ! For the time being we just throw everything back
   ! onto the normal routines.
-  call x%sync()
-  call y%sync()
+  if (x%is_dev()) call x%sync()
+  if (y%is_dev()) call y%sync()
   call a%spmm(alpha,x%v,beta,y%v,info,trans)
   call y%set_host()
 end subroutine psb_d_base_vect_mv
@@ -2060,8 +2060,6 @@ subroutine psb_d_base_vect_cssv(alpha,a,x,beta,y,info,trans,scale,d)
     goto 9999
   end if
 
-  call x%sync()
-  call y%sync()
   if (present(d)) then
     call d%sync()
     if (present(scale)) then
@@ -2161,8 +2159,11 @@ subroutine psb_d_base_inner_vect_sv(alpha,a,x,beta,y,info,trans)
 
   info  = psb_success_
   call psb_erractionsave(err_act)
+  if (x%is_dev()) call x%sync()
+  if (y%is_dev()) call y%sync()
 
   call a%inner_spsm(alpha,x%v,beta,y%v,info,trans)
+  call y%set_host()
 
   if (info /= psb_success_) then
     info = psb_err_from_subroutine_
diff --git a/base/serial/impl/psb_s_base_mat_impl.F90 b/base/serial/impl/psb_s_base_mat_impl.F90
index 4a99a684..0b47d472 100644
--- a/base/serial/impl/psb_s_base_mat_impl.F90
+++ b/base/serial/impl/psb_s_base_mat_impl.F90
@@ -2006,8 +2006,8 @@ subroutine psb_s_base_vect_mv(alpha,a,x,beta,y,info,trans)
 
   ! For the time being we just throw everything back
   ! onto the normal routines.
-  call x%sync()
-  call y%sync()
+  if (x%is_dev()) call x%sync()
+  if (y%is_dev()) call y%sync()
   call a%spmm(alpha,x%v,beta,y%v,info,trans)
   call y%set_host()
 end subroutine psb_s_base_vect_mv
@@ -2060,8 +2060,6 @@ subroutine psb_s_base_vect_cssv(alpha,a,x,beta,y,info,trans,scale,d)
     goto 9999
   end if
 
-  call x%sync()
-  call y%sync()
   if (present(d)) then
     call d%sync()
     if (present(scale)) then
@@ -2161,8 +2159,11 @@ subroutine psb_s_base_inner_vect_sv(alpha,a,x,beta,y,info,trans)
 
   info  = psb_success_
   call psb_erractionsave(err_act)
+  if (x%is_dev()) call x%sync()
+  if (y%is_dev()) call y%sync()
 
   call a%inner_spsm(alpha,x%v,beta,y%v,info,trans)
+  call y%set_host()
 
   if (info /= psb_success_) then
     info = psb_err_from_subroutine_
diff --git a/base/serial/impl/psb_z_base_mat_impl.F90 b/base/serial/impl/psb_z_base_mat_impl.F90
index 404027c5..2d97f698 100644
--- a/base/serial/impl/psb_z_base_mat_impl.F90
+++ b/base/serial/impl/psb_z_base_mat_impl.F90
@@ -2006,8 +2006,8 @@ subroutine psb_z_base_vect_mv(alpha,a,x,beta,y,info,trans)
 
   ! For the time being we just throw everything back
   ! onto the normal routines.
-  call x%sync()
-  call y%sync()
+  if (x%is_dev()) call x%sync()
+  if (y%is_dev()) call y%sync()
   call a%spmm(alpha,x%v,beta,y%v,info,trans)
   call y%set_host()
 end subroutine psb_z_base_vect_mv
@@ -2060,8 +2060,6 @@ subroutine psb_z_base_vect_cssv(alpha,a,x,beta,y,info,trans,scale,d)
     goto 9999
   end if
 
-  call x%sync()
-  call y%sync()
   if (present(d)) then
     call d%sync()
     if (present(scale)) then
@@ -2161,8 +2159,11 @@ subroutine psb_z_base_inner_vect_sv(alpha,a,x,beta,y,info,trans)
 
   info  = psb_success_
   call psb_erractionsave(err_act)
+  if (x%is_dev()) call x%sync()
+  if (y%is_dev()) call y%sync()
 
   call a%inner_spsm(alpha,x%v,beta,y%v,info,trans)
+  call y%set_host()
 
   if (info /= psb_success_) then
     info = psb_err_from_subroutine_