From de37e3602afe50cad4c7f9ca12e9f39e87dfe612 Mon Sep 17 00:00:00 2001
From: Salvatore Filippone <salvatore.filippone@uniroma2.it>
Date: Mon, 5 Dec 2022 10:52:26 +0100
Subject: [PATCH] Fix SV with CONJG

---
 base/serial/impl/psb_c_csr_impl.f90 | 4 ++--
 base/serial/impl/psb_d_csr_impl.f90 | 4 ++--
 base/serial/impl/psb_s_csr_impl.f90 | 4 ++--
 base/serial/impl/psb_z_csr_impl.f90 | 4 ++--
 4 files changed, 8 insertions(+), 8 deletions(-)

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)