diff --git a/base/serial/impl/psb_c_csr_impl.f90 b/base/serial/impl/psb_c_csr_impl.f90
index 38c67839..51d45825 100644
--- 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 
diff --git a/base/serial/impl/psb_d_csr_impl.f90 b/base/serial/impl/psb_d_csr_impl.f90
index 72890e44..094f141e 100644
--- 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 
diff --git a/base/serial/impl/psb_s_csr_impl.f90 b/base/serial/impl/psb_s_csr_impl.f90
index 182da555..7fe08ee0 100644
--- 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 
diff --git a/base/serial/impl/psb_z_csr_impl.f90 b/base/serial/impl/psb_z_csr_impl.f90
index ab4878ff..4713da34 100644
--- 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