fdurastante/psblas3
1
0
Fork 0
mirror of https://github.com/sfilippone/psblas3.git synced 2026-10-20 04:42:08 +00:00
Code Issues Packages Projects Releases Wiki Activity

psblas:

Browse Source 
base/serial/sort/psb_c_qsort_impl.f90
 base/serial/sort/psb_d_qsort_impl.f90
 base/serial/sort/psb_i_qsort_impl.f90
 base/serial/sort/psb_s_qsort_impl.f90
 base/serial/sort/psb_z_qsort_impl.f90

New thresholds for qsort.
This commit is contained in:
Salvatore Filippone
2016-06-18 09:14:05 +00:00
parent 2489f25247
commit 3ca60328a2
5 changed files with 48 additions and 51 deletions
Download Patch File Download Diff File Expand all files Collapse all files
base/serial/sort/psb_c_qsort_impl.f90
+12 -12
View File
@@ -152,7 +152,7 @@ subroutine psi_clqsrx_up(n,x,idx)
integer(psb_ipk_) :: i, j, ilx, iux, istp, lpiv
integer(psb_ipk_) :: ixt, n1, n2
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=16
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=32
integer(psb_ipk_) :: istack(nparms,maxstack)
if (n > ithrs) then
@@ -308,7 +308,7 @@ subroutine psi_clqsrx_dw(n,x,idx)
integer(psb_ipk_) :: i, j, ilx, iux, istp, lpiv
integer(psb_ipk_) :: ixt, n1, n2
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=16
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=32
integer(psb_ipk_) :: istack(nparms,maxstack)
if (n > ithrs) then
@@ -463,7 +463,7 @@ subroutine psi_clqsr_up(n,x)
integer(psb_ipk_) :: i, j, ilx, iux, istp, lpiv
integer(psb_ipk_) :: n1, n2
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=16
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=32
integer(psb_ipk_) :: istack(nparms,maxstack)
@@ -604,7 +604,7 @@ subroutine psi_clqsr_dw(n,x)
integer(psb_ipk_) :: i, j, ilx, iux, istp, lpiv
integer(psb_ipk_) :: n1, n2
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=16
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=32
integer(psb_ipk_) :: istack(nparms,maxstack)
@@ -746,7 +746,7 @@ subroutine psi_calqsrx_up(n,x,idx)
integer(psb_ipk_) :: i, j, ilx, iux, istp, lpiv
integer(psb_ipk_) :: ixt, n1, n2
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=16
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=32
integer(psb_ipk_) :: istack(nparms,maxstack)
if (n > ithrs) then
@@ -901,7 +901,7 @@ subroutine psi_calqsrx_dw(n,x,idx)
integer(psb_ipk_) :: i, j, ilx, iux, istp, lpiv
integer(psb_ipk_) :: ixt, n1, n2
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=16
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=32
integer(psb_ipk_) :: istack(nparms,maxstack)
if (n > ithrs) then
@@ -1056,7 +1056,7 @@ subroutine psi_calqsr_up(n,x)
integer(psb_ipk_) :: i, j, ilx, iux, istp, lpiv
integer(psb_ipk_) :: n1, n2
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=16
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=32
integer(psb_ipk_) :: istack(nparms,maxstack)
@@ -1196,7 +1196,7 @@ subroutine psi_calqsr_dw(n,x)
integer(psb_ipk_) :: i, j, ilx, iux, istp, lpiv
integer(psb_ipk_) :: n1, n2
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=16
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=32
integer(psb_ipk_) :: istack(nparms,maxstack)
@@ -1337,7 +1337,7 @@ subroutine psi_caqsrx_up(n,x,idx)
integer(psb_ipk_) :: i, j, ilx, iux, istp, lpiv
integer(psb_ipk_) :: ixt, n1, n2
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=16
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=32
integer(psb_ipk_) :: istack(nparms,maxstack)
if (n > ithrs) then
@@ -1493,7 +1493,7 @@ subroutine psi_caqsrx_dw(n,x,idx)
integer(psb_ipk_) :: i, j, ilx, iux, istp, lpiv
integer(psb_ipk_) :: ixt, n1, n2
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=16
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=32
integer(psb_ipk_) :: istack(nparms,maxstack)
if (n > ithrs) then
!
@@ -1646,7 +1646,7 @@ subroutine psi_caqsr_up(n,x)
integer(psb_ipk_) :: i, j, ilx, iux, istp, lpiv
integer(psb_ipk_) :: ixt, n1, n2
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=16
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=32
integer(psb_ipk_) :: istack(nparms,maxstack)
if (n > ithrs) then
@@ -1786,7 +1786,7 @@ subroutine psi_caqsr_dw(n,x)
integer(psb_ipk_) :: i, j, ilx, iux, istp, lpiv
integer(psb_ipk_) :: ixt, n1, n2
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=16
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=32
integer(psb_ipk_) :: istack(nparms,maxstack)
if (n > ithrs) then
base/serial/sort/psb_d_qsort_impl.f90
+8 -9
View File
@@ -141,8 +141,7 @@ subroutine psi_dqsrx_up(n,x,idx)
real(psb_dpk_) :: piv, xk, xt
integer(psb_ipk_) :: i, j, ilx, iux, istp, lpiv
integer(psb_ipk_) :: ixt, n1, n2
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=16
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=60
integer(psb_ipk_) :: istack(nparms,maxstack)
if (n > ithrs) then
@@ -296,7 +295,7 @@ subroutine psi_dqsrx_dw(n,x,idx)
integer(psb_ipk_) :: i, j, ilx, iux, istp, lpiv
integer(psb_ipk_) :: ixt, n1, n2
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=16
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=60
integer(psb_ipk_) :: istack(nparms,maxstack)
if (n > ithrs) then
@@ -451,7 +450,7 @@ subroutine psi_dqsr_up(n,x)
integer(psb_ipk_) :: i, j, ilx, iux, istp, lpiv
integer(psb_ipk_) :: n1, n2
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=16
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=60
integer(psb_ipk_) :: istack(nparms,maxstack)
@@ -592,7 +591,7 @@ subroutine psi_dqsr_dw(n,x)
integer(psb_ipk_) :: i, j, ilx, iux, istp, lpiv
integer(psb_ipk_) :: n1, n2
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=16
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=60
integer(psb_ipk_) :: istack(nparms,maxstack)
@@ -734,7 +733,7 @@ subroutine psi_daqsrx_up(n,x,idx)
integer(psb_ipk_) :: i, j, ilx, iux, istp, lpiv
integer(psb_ipk_) :: ixt, n1, n2
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=16
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=60
integer(psb_ipk_) :: istack(nparms,maxstack)
if (n > ithrs) then
@@ -890,7 +889,7 @@ subroutine psi_daqsrx_dw(n,x,idx)
integer(psb_ipk_) :: i, j, ilx, iux, istp, lpiv
integer(psb_ipk_) :: ixt, n1, n2
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=16
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=60
integer(psb_ipk_) :: istack(nparms,maxstack)
if (n > ithrs) then
!
@@ -1043,7 +1042,7 @@ subroutine psi_daqsr_up(n,x)
integer(psb_ipk_) :: i, j, ilx, iux, istp, lpiv
integer(psb_ipk_) :: ixt, n1, n2
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=16
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=60
integer(psb_ipk_) :: istack(nparms,maxstack)
if (n > ithrs) then
@@ -1183,7 +1182,7 @@ subroutine psi_daqsr_dw(n,x)
integer(psb_ipk_) :: i, j, ilx, iux, istp, lpiv
integer(psb_ipk_) :: ixt, n1, n2
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=16
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=60
integer(psb_ipk_) :: istack(nparms,maxstack)
if (n > ithrs) then
base/serial/sort/psb_i_qsort_impl.f90
+8 -9
View File
@@ -141,8 +141,7 @@ subroutine psi_iqsrx_up(n,x,idx)
integer(psb_ipk_) :: piv, xk, xt
integer(psb_ipk_) :: i, j, ilx, iux, istp, lpiv
integer(psb_ipk_) :: ixt, n1, n2
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=16
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=72
integer(psb_ipk_) :: istack(nparms,maxstack)
if (n > ithrs) then
@@ -296,7 +295,7 @@ subroutine psi_iqsrx_dw(n,x,idx)
integer(psb_ipk_) :: i, j, ilx, iux, istp, lpiv
integer(psb_ipk_) :: ixt, n1, n2
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=16
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=72
integer(psb_ipk_) :: istack(nparms,maxstack)
if (n > ithrs) then
@@ -451,7 +450,7 @@ subroutine psi_iqsr_up(n,x)
integer(psb_ipk_) :: i, j, ilx, iux, istp, lpiv
integer(psb_ipk_) :: n1, n2
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=16
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=72
integer(psb_ipk_) :: istack(nparms,maxstack)
@@ -592,7 +591,7 @@ subroutine psi_iqsr_dw(n,x)
integer(psb_ipk_) :: i, j, ilx, iux, istp, lpiv
integer(psb_ipk_) :: n1, n2
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=16
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=72
integer(psb_ipk_) :: istack(nparms,maxstack)
@@ -734,7 +733,7 @@ subroutine psi_iaqsrx_up(n,x,idx)
integer(psb_ipk_) :: i, j, ilx, iux, istp, lpiv
integer(psb_ipk_) :: ixt, n1, n2
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=16
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=72
integer(psb_ipk_) :: istack(nparms,maxstack)
if (n > ithrs) then
@@ -890,7 +889,7 @@ subroutine psi_iaqsrx_dw(n,x,idx)
integer(psb_ipk_) :: i, j, ilx, iux, istp, lpiv
integer(psb_ipk_) :: ixt, n1, n2
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=16
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=72
integer(psb_ipk_) :: istack(nparms,maxstack)
if (n > ithrs) then
!
@@ -1043,7 +1042,7 @@ subroutine psi_iaqsr_up(n,x)
integer(psb_ipk_) :: i, j, ilx, iux, istp, lpiv
integer(psb_ipk_) :: ixt, n1, n2
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=16
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=72
integer(psb_ipk_) :: istack(nparms,maxstack)
if (n > ithrs) then
@@ -1183,7 +1182,7 @@ subroutine psi_iaqsr_dw(n,x)
integer(psb_ipk_) :: i, j, ilx, iux, istp, lpiv
integer(psb_ipk_) :: ixt, n1, n2
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=16
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=72
integer(psb_ipk_) :: istack(nparms,maxstack)
if (n > ithrs) then
base/serial/sort/psb_s_qsort_impl.f90
+8 -9
View File
@@ -141,8 +141,7 @@ subroutine psi_sqsrx_up(n,x,idx)
real(psb_spk_) :: piv, xk, xt
integer(psb_ipk_) :: i, j, ilx, iux, istp, lpiv
integer(psb_ipk_) :: ixt, n1, n2
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=16
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=72
integer(psb_ipk_) :: istack(nparms,maxstack)
if (n > ithrs) then
@@ -296,7 +295,7 @@ subroutine psi_sqsrx_dw(n,x,idx)
integer(psb_ipk_) :: i, j, ilx, iux, istp, lpiv
integer(psb_ipk_) :: ixt, n1, n2
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=16
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=72
integer(psb_ipk_) :: istack(nparms,maxstack)
if (n > ithrs) then
@@ -451,7 +450,7 @@ subroutine psi_sqsr_up(n,x)
integer(psb_ipk_) :: i, j, ilx, iux, istp, lpiv
integer(psb_ipk_) :: n1, n2
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=16
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=72
integer(psb_ipk_) :: istack(nparms,maxstack)
@@ -592,7 +591,7 @@ subroutine psi_sqsr_dw(n,x)
integer(psb_ipk_) :: i, j, ilx, iux, istp, lpiv
integer(psb_ipk_) :: n1, n2
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=16
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=72
integer(psb_ipk_) :: istack(nparms,maxstack)
@@ -734,7 +733,7 @@ subroutine psi_saqsrx_up(n,x,idx)
integer(psb_ipk_) :: i, j, ilx, iux, istp, lpiv
integer(psb_ipk_) :: ixt, n1, n2
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=16
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=72
integer(psb_ipk_) :: istack(nparms,maxstack)
if (n > ithrs) then
@@ -890,7 +889,7 @@ subroutine psi_saqsrx_dw(n,x,idx)
integer(psb_ipk_) :: i, j, ilx, iux, istp, lpiv
integer(psb_ipk_) :: ixt, n1, n2
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=16
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=72
integer(psb_ipk_) :: istack(nparms,maxstack)
if (n > ithrs) then
!
@@ -1043,7 +1042,7 @@ subroutine psi_saqsr_up(n,x)
integer(psb_ipk_) :: i, j, ilx, iux, istp, lpiv
integer(psb_ipk_) :: ixt, n1, n2
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=16
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=72
integer(psb_ipk_) :: istack(nparms,maxstack)
if (n > ithrs) then
@@ -1183,7 +1182,7 @@ subroutine psi_saqsr_dw(n,x)
integer(psb_ipk_) :: i, j, ilx, iux, istp, lpiv
integer(psb_ipk_) :: ixt, n1, n2
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=16
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=72
integer(psb_ipk_) :: istack(nparms,maxstack)
if (n > ithrs) then
base/serial/sort/psb_z_qsort_impl.f90
+12 -12
View File
@@ -152,7 +152,7 @@ subroutine psi_zlqsrx_up(n,x,idx)
integer(psb_ipk_) :: i, j, ilx, iux, istp, lpiv
integer(psb_ipk_) :: ixt, n1, n2
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=16
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=24
integer(psb_ipk_) :: istack(nparms,maxstack)
if (n > ithrs) then
@@ -308,7 +308,7 @@ subroutine psi_zlqsrx_dw(n,x,idx)
integer(psb_ipk_) :: i, j, ilx, iux, istp, lpiv
integer(psb_ipk_) :: ixt, n1, n2
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=16
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=24
integer(psb_ipk_) :: istack(nparms,maxstack)
if (n > ithrs) then
@@ -463,7 +463,7 @@ subroutine psi_zlqsr_up(n,x)
integer(psb_ipk_) :: i, j, ilx, iux, istp, lpiv
integer(psb_ipk_) :: n1, n2
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=16
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=24
integer(psb_ipk_) :: istack(nparms,maxstack)
@@ -604,7 +604,7 @@ subroutine psi_zlqsr_dw(n,x)
integer(psb_ipk_) :: i, j, ilx, iux, istp, lpiv
integer(psb_ipk_) :: n1, n2
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=16
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=24
integer(psb_ipk_) :: istack(nparms,maxstack)
@@ -746,7 +746,7 @@ subroutine psi_zalqsrx_up(n,x,idx)
integer(psb_ipk_) :: i, j, ilx, iux, istp, lpiv
integer(psb_ipk_) :: ixt, n1, n2
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=16
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=24
integer(psb_ipk_) :: istack(nparms,maxstack)
if (n > ithrs) then
@@ -901,7 +901,7 @@ subroutine psi_zalqsrx_dw(n,x,idx)
integer(psb_ipk_) :: i, j, ilx, iux, istp, lpiv
integer(psb_ipk_) :: ixt, n1, n2
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=16
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=24
integer(psb_ipk_) :: istack(nparms,maxstack)
if (n > ithrs) then
@@ -1056,7 +1056,7 @@ subroutine psi_zalqsr_up(n,x)
integer(psb_ipk_) :: i, j, ilx, iux, istp, lpiv
integer(psb_ipk_) :: n1, n2
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=16
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=24
integer(psb_ipk_) :: istack(nparms,maxstack)
@@ -1196,7 +1196,7 @@ subroutine psi_zalqsr_dw(n,x)
integer(psb_ipk_) :: i, j, ilx, iux, istp, lpiv
integer(psb_ipk_) :: n1, n2
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=16
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=24
integer(psb_ipk_) :: istack(nparms,maxstack)
@@ -1337,7 +1337,7 @@ subroutine psi_zaqsrx_up(n,x,idx)
integer(psb_ipk_) :: i, j, ilx, iux, istp, lpiv
integer(psb_ipk_) :: ixt, n1, n2
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=16
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=24
integer(psb_ipk_) :: istack(nparms,maxstack)
if (n > ithrs) then
@@ -1493,7 +1493,7 @@ subroutine psi_zaqsrx_dw(n,x,idx)
integer(psb_ipk_) :: i, j, ilx, iux, istp, lpiv
integer(psb_ipk_) :: ixt, n1, n2
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=16
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=24
integer(psb_ipk_) :: istack(nparms,maxstack)
if (n > ithrs) then
!
@@ -1646,7 +1646,7 @@ subroutine psi_zaqsr_up(n,x)
integer(psb_ipk_) :: i, j, ilx, iux, istp, lpiv
integer(psb_ipk_) :: ixt, n1, n2
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=16
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=24
integer(psb_ipk_) :: istack(nparms,maxstack)
if (n > ithrs) then
@@ -1786,7 +1786,7 @@ subroutine psi_zaqsr_dw(n,x)
integer(psb_ipk_) :: i, j, ilx, iux, istp, lpiv
integer(psb_ipk_) :: ixt, n1, n2
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=16
integer(psb_ipk_), parameter :: maxstack=64,nparms=3,ithrs=24
integer(psb_ipk_) :: istack(nparms,maxstack)
if (n > ithrs) then