psblas3:

base/serial/impl/psb_c_csr_impl.f90 base/serial/impl/psb_d_csr_impl.f90 base/serial/impl/psb_s_csr_impl.f90 base/serial/impl/psb_z_csr_impl.f90 Try uysing qsort in csrspspmm
9 years ago · 12bf6b82f5
parent 2bd6693741
commit 12bf6b82f5
4 changed files with 4 additions and 4 deletions
--- a/base/serial/impl/psb_c_csr_impl.f90
+++ b/base/serial/impl/psb_c_csr_impl.f90
@ -3156,7 +3156,7 @@ contains
          if (info /= 0) return
        end if
          
-        call psb_msort(idxs(1:nrc))
+        call psb_qsort(idxs(1:nrc))
        do i=1, nrc
          irw        = idxs(i) 
          c%ja(nzc)  = irw 
--- a/base/serial/impl/psb_d_csr_impl.f90
+++ b/base/serial/impl/psb_d_csr_impl.f90
@ -3156,7 +3156,7 @@ contains
          if (info /= 0) return
        end if
          
-        call psb_msort(idxs(1:nrc))
+        call psb_qsort(idxs(1:nrc))
        do i=1, nrc
          irw        = idxs(i) 
          c%ja(nzc)  = irw 
--- a/base/serial/impl/psb_s_csr_impl.f90
+++ b/base/serial/impl/psb_s_csr_impl.f90
@ -3156,7 +3156,7 @@ contains
          if (info /= 0) return
        end if
          
-        call psb_msort(idxs(1:nrc))
+        call psb_qsort(idxs(1:nrc))
        do i=1, nrc
          irw        = idxs(i) 
          c%ja(nzc)  = irw 
--- a/base/serial/impl/psb_z_csr_impl.f90
+++ b/base/serial/impl/psb_z_csr_impl.f90
@ -3156,7 +3156,7 @@ contains
          if (info /= 0) return
        end if
          
-        call psb_msort(idxs(1:nrc))
+        call psb_qsort(idxs(1:nrc))
        do i=1, nrc
          irw        = idxs(i) 
          c%ja(nzc)  = irw