@@ -44,7 +44,7 @@ contains
44
44
call select(x_tmp, c, val)
45
45
46
46
if (mod(n, 2_int64) == 0) then
47
- call select(x_tmp, c+1, val1)
47
+ call select(x_tmp, c+1, val1, left = c )
48
48
#:if t1[0] == 'r'
49
49
res = (val + val1) / 2._${o1}$
50
50
#:else
@@ -111,7 +111,7 @@ contains
111
111
call select(x_tmp, c, val)
112
112
113
113
if (mod(n, 2) == 0) then
114
- call select(x_tmp, c+1, val1)
114
+ call select(x_tmp, c+1, val1, left = c )
115
115
res${reduce_subvector('j', rank, fi)}$ = &
116
116
#:if t1[0] == 'r'
117
117
(val + val1) / 2._${o1}$
@@ -171,7 +171,7 @@ contains
171
171
call select(x_tmp, c, val)
172
172
173
173
if (mod(n, 2_int64) == 0) then
174
- call select(x_tmp, c+1, val1)
174
+ call select(x_tmp, c+1, val1, left = c )
175
175
#:if t1[0] == 'r'
176
176
res = (val + val1) / 2._${o1}$
177
177
#:else
@@ -242,7 +242,7 @@ contains
242
242
call select(x_tmp, c, val)
243
243
244
244
if (mod(n, 2_int64) == 0) then
245
- call select(x_tmp, c+1, val1)
245
+ call select(x_tmp, c+1, val1, left = c )
246
246
res${reduce_subvector('j', rank, fi)}$ = &
247
247
#:if t1[0] == 'r'
248
248
(val + val1) / 2._${o1}$
0 commit comments