diff --git a/base/serial/impl/psb_c_csr_impl.f90 b/base/serial/impl/psb_c_csr_impl.f90
index bf5e85aa..55a91648 100644
--- a/base/serial/impl/psb_c_csr_impl.f90
+++ b/base/serial/impl/psb_c_csr_impl.f90
@@ -993,7 +993,7 @@ contains
           end do
         else if (.not.unit) then
           do i=n, 1, -1
-            y(i) = y(i)/val(irp(i+1)-1)
+            y(i) = y(i)/conjg(val(irp(i+1)-1))
             acc  = y(i)
             do j=irp(i), irp(i+1)-2
               jc    = ja(j)
@@ -1013,7 +1013,7 @@ contains
           end do
         else if (.not.unit) then
           do i=1, n
-            y(i) = y(i)/val(irp(i))
+            y(i) = y(i)/conjg(val(irp(i)))
             acc  = y(i)
             do j=irp(i)+1, irp(i+1)-1
               jc    = ja(j)
diff --git a/base/serial/impl/psb_d_csr_impl.f90 b/base/serial/impl/psb_d_csr_impl.f90
index 0bbd99a0..2c59c1a5 100644
--- a/base/serial/impl/psb_d_csr_impl.f90
+++ b/base/serial/impl/psb_d_csr_impl.f90
@@ -993,7 +993,7 @@ contains
           end do
         else if (.not.unit) then
           do i=n, 1, -1
-            y(i) = y(i)/val(irp(i+1)-1)
+            y(i) = y(i)/(val(irp(i+1)-1))
             acc  = y(i)
             do j=irp(i), irp(i+1)-2
               jc    = ja(j)
@@ -1013,7 +1013,7 @@ contains
           end do
         else if (.not.unit) then
           do i=1, n
-            y(i) = y(i)/val(irp(i))
+            y(i) = y(i)/(val(irp(i)))
             acc  = y(i)
             do j=irp(i)+1, irp(i+1)-1
               jc    = ja(j)
diff --git a/base/serial/impl/psb_s_csr_impl.f90 b/base/serial/impl/psb_s_csr_impl.f90
index 34f937a4..75358dbc 100644
--- a/base/serial/impl/psb_s_csr_impl.f90
+++ b/base/serial/impl/psb_s_csr_impl.f90
@@ -993,7 +993,7 @@ contains
           end do
         else if (.not.unit) then
           do i=n, 1, -1
-            y(i) = y(i)/val(irp(i+1)-1)
+            y(i) = y(i)/(val(irp(i+1)-1))
             acc  = y(i)
             do j=irp(i), irp(i+1)-2
               jc    = ja(j)
@@ -1013,7 +1013,7 @@ contains
           end do
         else if (.not.unit) then
           do i=1, n
-            y(i) = y(i)/val(irp(i))
+            y(i) = y(i)/(val(irp(i)))
             acc  = y(i)
             do j=irp(i)+1, irp(i+1)-1
               jc    = ja(j)
diff --git a/base/serial/impl/psb_z_csr_impl.f90 b/base/serial/impl/psb_z_csr_impl.f90
index 3084425f..4f2693c0 100644
--- a/base/serial/impl/psb_z_csr_impl.f90
+++ b/base/serial/impl/psb_z_csr_impl.f90
@@ -993,7 +993,7 @@ contains
           end do
         else if (.not.unit) then
           do i=n, 1, -1
-            y(i) = y(i)/val(irp(i+1)-1)
+            y(i) = y(i)/conjg(val(irp(i+1)-1))
             acc  = y(i)
             do j=irp(i), irp(i+1)-2
               jc    = ja(j)
@@ -1013,7 +1013,7 @@ contains
           end do
         else if (.not.unit) then
           do i=1, n
-            y(i) = y(i)/val(irp(i))
+            y(i) = y(i)/conjg(val(irp(i)))
             acc  = y(i)
             do j=irp(i)+1, irp(i+1)-1
               jc    = ja(j)