Skip to content

Commit 410012e

Browse files
committed
refactor: to use zero-based indexing
--- type: pre_commit_static_analysis_report description: Results of running static analysis checks when committing changes. report: - task: lint_filenames status: passed - task: lint_editorconfig status: passed - task: lint_markdown status: na - task: lint_package_json status: na - task: lint_repl_help status: na - task: lint_javascript_src status: passed - task: lint_javascript_cli status: na - task: lint_javascript_examples status: na - task: lint_javascript_tests status: na - task: lint_javascript_benchmarks status: na - task: lint_python status: na - task: lint_r status: na - task: lint_c_src status: na - task: lint_c_examples status: na - task: lint_c_benchmarks status: na - task: lint_c_tests_fixtures status: na - task: lint_shell status: na - task: lint_typescript_declarations status: na - task: lint_typescript_tests status: na - task: lint_license_headers status: passed ---
1 parent 7e66211 commit 410012e

File tree

1 file changed

+18
-20
lines changed
  • lib/node_modules/@stdlib/fft/base/fftpack/lib

1 file changed

+18
-20
lines changed

lib/node_modules/@stdlib/fft/base/fftpack/lib/radb2.js

Lines changed: 18 additions & 20 deletions
Original file line numberDiff line numberDiff line change
@@ -154,7 +154,7 @@ function optr( i, k, j, L, M, stride, offset ) {
154154
* @returns {void}
155155
*/
156156
function radb2( M, L, cc, oc, out, oo, twiddles, ot ) { // FIXME: support strides
157-
var idp2;
157+
var MP1;
158158
var tr2;
159159
var ti2;
160160
var ip1;
@@ -196,37 +196,35 @@ function radb2( M, L, cc, oc, out, oo, twiddles, ot ) { // FIXME: support stride
196196
return;
197197
}
198198

199-
// FIXME: everything below needs updating as `iptr` and `optr` have been converted to zero-based indexing
200-
201199
if ( M !== 2 ) {
202-
idp2 = M + 2;
200+
MP1 = M + 1;
203201

204202
// FIXME: describe what is happening in the loops below
205203
for ( k = 0; k < L; k++ ) {
206-
for ( i = 3; i <= M; i += 2 ) {
207-
ic = idp2 - i;
204+
for ( i = 2; i < M; i += 2 ) {
205+
ic = MP1 - i;
208206

209-
it1 = i - 2 + ot;
210-
it2 = i - 1 + ot;
207+
it1 = i - 1 + ot;
208+
it2 = i + ot;
211209

212-
ip1 = iptr( i-1, 1, k, M, oc );
213-
ip2 = iptr( ic-1, 2, k, M, oc );
210+
ip1 = iptr( i, 0, k, M, sc, oc );
211+
ip2 = iptr( ic-1, 1, k, M, sc, oc );
214212
tr2 = cc[ ip1 ] - cc[ ip2 ];
215213

216-
io = optr( i-1, k, 1, L, M, oo );
214+
io = optr( i, k, 0, L, M, so, oo );
217215
out[ io ] = cc[ ip1 ] + cc[ ip2 ];
218216

219-
ip1 = iptr( i, 1, k, M, oc );
220-
ip2 = iptr( ic, 2, k, M, oc );
217+
ip1 = iptr( i+1, 0, k, M, sc, oc );
218+
ip2 = iptr( ic, 1, k, M, sc, oc );
221219
ti2 = cc[ ip1 ] + cc[ ip2 ];
222220

223-
io = optr( i, k, 1, L, M, oo );
221+
io = optr( i+1, k, 0, L, M, so, oo );
224222
out[ io ] = cc[ ip1 ] - cc[ ip2 ];
225223

226-
io = optr( i-1, k, 2, L, M, oo );
224+
io = optr( i, k, 1, L, M, so, oo );
227225
out[ io ] = ( twiddles[ it1 ] * tr2 ) - ( twiddles[ it2 ] * ti2 );
228226

229-
io = optr( i, k, 2, L, M, oo );
227+
io = optr( i+1, k, 1, L, M, so, oo );
230228
out[ io ] = ( twiddles[ it1 ] * ti2 ) + ( twiddles[ it2 ] * tr2 );
231229
}
232230
}
@@ -237,12 +235,12 @@ function radb2( M, L, cc, oc, out, oo, twiddles, ot ) { // FIXME: support stride
237235
}
238236
// FIXME: describe what is happening here
239237
for ( k = 0; k < L; k++ ) {
240-
ip1 = iptr( M, 1, k, M, oc );
241-
io = optr( M, k, 1, L, M, oo );
238+
ip1 = iptr( M-1, 0, k, M, sc, oc );
239+
io = optr( M-1, k, 0, L, M, so, oo );
242240
out[ io ] = cc[ ip1 ] + cc[ ip1 ];
243241

244-
ip1 = iptr( 1, 2, k, M, oc );
245-
io = optr( M, k, 2, L, M, oo );
242+
ip1 = iptr( 0, 1, k, M, sc, oc );
243+
io = optr( M-1, k, 1, L, M, so, oo );
246244
out[ io ] = -( cc[ ip1 ] + cc[ ip1 ] );
247245
}
248246
}

0 commit comments

Comments
 (0)