Skip to content

Commit 535ec3a

Browse files
committed
refactor: improve float32 precision for addition assignment
--- 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: passed - task: lint_typescript_tests status: na - task: lint_license_headers status: passed ---
1 parent 5ad9f8e commit 535ec3a

1 file changed

Lines changed: 4 additions & 4 deletions

File tree

  • lib/node_modules/@stdlib/blas/base/ctrsv/lib

lib/node_modules/@stdlib/blas/base/ctrsv/lib/base.js

Lines changed: 4 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -214,8 +214,8 @@ function ctrsv( uplo, trans, diag, N, A, strideA1, strideA2, offsetA, x, strideX
214214
ima = sign * viewA[ ia + 1 ];
215215
rex = viewX[ ix0 ];
216216
imx = viewX[ ix0 + 1 ];
217-
retmp -= f32( ( rex * rea ) - ( imx * ima ) );
218-
imtmp -= f32( ( rex * ima ) + ( imx * rea ) );
217+
retmp = f32( retmp - f32( ( rex * rea ) - ( imx * ima ) ) );
218+
imtmp = f32( imtmp - f32( ( rex * ima ) + ( imx * rea ) ) );
219219
ix0 += sx;
220220
}
221221
if ( nonunit ) {
@@ -249,8 +249,8 @@ function ctrsv( uplo, trans, diag, N, A, strideA1, strideA2, offsetA, x, strideX
249249
ima = sign * viewA[ ia + 1 ];
250250
rex = viewX[ ix0 ];
251251
imx = viewX[ ix0 + 1 ];
252-
retmp -= f32( ( rex * rea ) - ( imx * ima ) );
253-
imtmp -= f32( ( rex * ima ) + ( imx * rea ) );
252+
retmp = f32( retmp - f32( ( rex * rea ) - ( imx * ima ) ) );
253+
imtmp = f32( imtmp - f32( ( rex * ima ) + ( imx * rea ) ) );
254254
ix0 -= sx;
255255
}
256256
if ( nonunit ) {

0 commit comments

Comments
 (0)