@@ -525,141 +525,141 @@ <h1><a href="../../../../index.html">All files</a> / <a href="index.html">blas/b
525525* cher2( 'lower', x.length, alpha, x, 1, 0, y, 1, 0, A, 2, 1, 0 );
526526* // A => <Complex64Array>[ 5.0, 0.0, 0.0, 0.0, 6.0, 3.0, 8.0, 0.0 ]
527527*/
528- function cher2( uplo, N, alpha, x, strideX, offsetX, y, strideY, offsetY, A, strideA1, strideA2, offsetA ) {
529- var alpha_re;
530- var alpha_im;
531- var tmp1_re;
532- var tmp1_im;
533- var tmp2_re;
534- var tmp2_im;
535- var re0_x;
536- var im0_x;
537- var re0_y;
538- var im0_y;
539- var re1_x;
540- var im1_x;
541- var re1_y;
542- var im1_y;
528+ function cher2( uplo, N, alpha, x, strideX, offsetX, y, strideY, offsetY, A, strideA1, strideA2, offsetA ) { // eslint-disable-line max-params, max-len
529+ var realpha;
530+ var imalpha;
531+ var retmp1;
532+ var imtmp1;
533+ var retmp2;
534+ var imtmp2;
543535 var viewX;
544- var viewY;
545- var viewA;
546- var isrm;
547- var idx;
548- var ix1;
549- var iy1;
550- var ix0;
551- var iy0;
552- var re0;
553- var im0;
554- var re1;
555- var im1;
556- var sa0;
557- var sa1;
558- var re;
559- var im;
560- var i1;
561- var i0;
562- var ix;
563- var iy;
564- var ia;
565- var sx;
566- var sy;
536+ var viewY;
537+ var viewA;
538+ var xre0;
539+ var xim0;
540+ var yre0;
541+ var yim0;
542+ var xre1;
543+ var xim1;
544+ var yre1;
545+ var yim1;
546+ var isrm;
547+ var idx;
548+ var ix1;
549+ var iy1;
550+ var ix0;
551+ var iy0;
552+ var re0;
553+ var im0;
554+ var re1;
555+ var im1;
556+ var sa0;
557+ var sa1;
558+ var re;
559+ var im;
560+ var i1;
561+ var i0;
562+ var ix;
563+ var iy;
564+ var ia;
565+ var sx;
566+ var sy;
567567
568- alpha_re = f32( alpha.re );
569- alpha_im = f32( alpha.im );
568+ realpha = f32( alpha.re );
569+ imalpha = f32( alpha.im );
570570
571- viewX = reinterpret( x, 0 );
572- viewY = reinterpret( y, 0 );
573- viewA = reinterpret( A, 0 );
571+ viewX = reinterpret( x, 0 );
572+ viewY = reinterpret( y, 0 );
573+ viewA = reinterpret( A, 0 );
574574
575- isrm = isRowMajor( [ strideA1, strideA2 ] );
576- if ( isrm ) {
577- sa0 = strideA2 * 2;
578- sa1 = strideA1 * 2;
579- } else {
580- sa0 = strideA1 * 2;
581- sa1 = strideA2 * 2;
582- }
583- ix = offsetX * 2;
584- iy = offsetY * 2;
585- ia = offsetA * 2;
586- sx = strideX * 2;
587- sy = strideY * 2;
588- if ( ( isrm && uplo === 'upper' ) || ( !isrm && uplo === 'lower' ) ) {
589- for ( i1 = 0; i1 < N; i1++ ) {
590- ix1 = ix + ( i1 * sx );
591- iy1 = iy + ( i1 * sy );
592- re0_x = f32( viewX[ ix1 ] );
593- im0_x = f32( viewX[ ix1 + 1 ] );
594- re0_y = f32( viewY[ iy1 ] );
595- im0_y = f32( viewY[ iy1 + 1 ] );
596- for ( i0 = i1; i0 < N; i0++ ) {
597- ix0 = ix + ( i0 * sx );
598- iy0 = iy + ( i0 * sy );
599- re1_x = f32( viewX[ ix0 ] );
600- im1_x = f32( viewX[ ix0 + 1 ] );
601- re1_y = f32( viewY[ iy0 ] );
602- im1_y = f32( viewY[ iy0 + 1 ] );
575+ isrm = isRowMajor( [ strideA1, strideA2 ] );
576+ if ( isrm ) {
577+ sa0 = strideA2 * 2;
578+ sa1 = strideA1 * 2;
579+ } else {
580+ sa0 = strideA1 * 2;
581+ sa1 = strideA2 * 2;
582+ }
583+ ix = offsetX * 2;
584+ iy = offsetY * 2;
585+ ia = offsetA * 2;
586+ sx = strideX * 2;
587+ sy = strideY * 2;
588+ if ( ( isrm && uplo === 'upper' ) || ( !isrm && uplo === 'lower' ) ) {
589+ for ( i1 = 0; i1 < N; i1++ ) {
590+ ix1 = ix + ( i1 * sx );
591+ iy1 = iy + ( i1 * sy );
592+ xre0 = f32( viewX[ ix1 ] );
593+ xim0 = f32( viewX[ ix1 + 1 ] );
594+ yre0 = f32( viewY[ iy1 ] );
595+ yim0 = f32( viewY[ iy1 + 1 ] );
596+ for ( i0 = i1; i0 < N; i0++ ) {
597+ ix0 = ix + ( i0 * sx );
598+ iy0 = iy + ( i0 * sy );
599+ xre1 = f32( viewX[ ix0 ] );
600+ xim1 = f32( viewX[ ix0 + 1 ] );
601+ yre1 = f32( viewY[ iy0 ] );
602+ yim1 = f32( viewY[ iy0 + 1 ] );
603603
604- re0 = f32( ( re0_x * re1_y ) + ( im0_x * im1_y ) );
605- im0 = f32( ( im0_x * re1_y ) - ( re0_x * im1_y ) );
606- tmp1_re = f32( ( alpha_re * re0 ) - ( alpha_im * im0 ) );
607- tmp1_im = f32( ( alpha_re * im0 ) + ( alpha_im * re0 ) );
604+ re0 = f32( ( xre0 * yre1 ) + ( xim0 * yim1 ) );
605+ im0 = f32( ( xim0 * yre1 ) - ( xre0 * yim1 ) );
606+ retmp1 = f32( ( realpha * re0 ) - ( imalpha * im0 ) );
607+ imtmp1 = f32( ( realpha * im0 ) + ( imalpha * re0 ) );
608608
609- re1 = f32( ( re0_y * re1_x ) + ( im0_y * im1_x ) );
610- im1 = f32( ( im0_y * re1_x ) - ( re0_y * im1_x ) );
611- tmp2_re = f32( ( alpha_re * re1 ) + ( alpha_im * im1 ) );
612- tmp2_im = f32( ( alpha_re * im1 ) - ( alpha_im * re1 ) );
609+ re1 = f32( ( yre0 * xre1 ) + ( yim0 * xim1 ) );
610+ im1 = f32( ( yim0 * xre1 ) - ( yre0 * xim1 ) );
611+ retmp2 = f32( ( realpha * re1 ) + ( imalpha * im1 ) );
612+ imtmp2 = f32( ( realpha * im1 ) - ( imalpha * re1 ) );
613613
614- re = f32( tmp1_re + tmp2_re );
615- im = f32( tmp1_im + tmp2_im );
616- idx = ia + ( i0 * sa0 ) + ( i1 * sa1 );
617- viewA[ idx ] = f32( viewA[ idx ] + re );
618- viewA[ idx + 1 ] = f32( viewA[ idx + 1 ] + im );
619- if ( i0 === i1 ) {
620- viewA[ idx + 1 ] = 0.0;
621- }
622- }
623- }
624- return A;
625- }
626- // ( isrm && uplo === 'lower' ) || ( !isrm && uplo === 'upper' )
627- for ( i1 = 0; i1 < N; i1++ ) {
628- ix1 = ix + ( i1 * sx );
629- iy1 = iy + ( i1 * sy );
630- re0_x = f32( viewX[ ix1 ] );
631- im0_x = f32( viewX[ ix1 + 1 ] );
632- re0_y = f32( viewY[ iy1 ] );
633- im0_y = f32( viewY[ iy1 + 1 ] );
634- for ( i0 = 0; i0 <= i1; i0++ ) {
635- ix0 = ix + ( i0 * sx );
636- iy0 = iy + ( i0 * sy );
637- re1_x = f32( viewX[ ix0 ] );
638- im1_x = f32( viewX[ ix0 + 1 ] );
639- re1_y = f32( viewY[ iy0 ] );
640- im1_y = f32( viewY[ iy0 + 1 ] );
614+ re = f32( retmp1 + retmp2 );
615+ im = f32( imtmp1 + imtmp2 );
616+ idx = ia + ( i0 * sa0 ) + ( i1 * sa1 );
617+ viewA[ idx ] = f32( viewA[ idx ] + re );
618+ viewA[ idx + 1 ] = f32( viewA[ idx + 1 ] + im );
619+ if ( i0 === i1 ) {
620+ viewA[ idx + 1 ] = 0.0;
621+ }
622+ }
623+ }
624+ return A;
625+ }
626+ // ( isrm && uplo === 'lower' ) || ( !isrm && uplo === 'upper' )
627+ for ( i1 = 0; i1 < N; i1++ ) {
628+ ix1 = ix + ( i1 * sx );
629+ iy1 = iy + ( i1 * sy );
630+ xre0 = f32( viewX[ ix1 ] );
631+ xim0 = f32( viewX[ ix1 + 1 ] );
632+ yre0 = f32( viewY[ iy1 ] );
633+ yim0 = f32( viewY[ iy1 + 1 ] );
634+ for ( i0 = 0; i0 <= i1; i0++ ) {
635+ ix0 = ix + ( i0 * sx );
636+ iy0 = iy + ( i0 * sy );
637+ xre1 = f32( viewX[ ix0 ] );
638+ xim1 = f32( viewX[ ix0 + 1 ] );
639+ yre1 = f32( viewY[ iy0 ] );
640+ yim1 = f32( viewY[ iy0 + 1 ] );
641641
642- re0 = f32( ( re0_x * re1_y ) + ( im0_x * im1_y ) );
643- im0 = f32( ( im0_x * re1_y ) - ( re0_x * im1_y ) );
644- tmp1_re = f32( ( alpha_re * re0 ) - ( alpha_im * im0 ) );
645- tmp1_im = f32( ( alpha_re * im0 ) + ( alpha_im * re0 ) );
642+ re0 = f32( ( xre0 * yre1 ) + ( xim0 * yim1 ) );
643+ im0 = f32( ( xim0 * yre1 ) - ( xre0 * yim1 ) );
644+ retmp1 = f32( ( realpha * re0 ) - ( imalpha * im0 ) );
645+ imtmp1 = f32( ( realpha * im0 ) + ( imalpha * re0 ) );
646646
647- re1 = f32( ( re0_y * re1_x ) + ( im0_y * im1_x ) );
648- im1 = f32( ( im0_y * re1_x ) - ( re0_y * im1_x ) );
649- tmp2_re = f32( ( alpha_re * re1 ) + ( alpha_im * im1 ) );
650- tmp2_im = f32( ( alpha_re * im1 ) - ( alpha_im * re1 ) );
647+ re1 = f32( ( yre0 * xre1 ) + ( yim0 * xim1 ) );
648+ im1 = f32( ( yim0 * xre1 ) - ( yre0 * xim1 ) );
649+ retmp2 = f32( ( realpha * re1 ) + ( imalpha * im1 ) );
650+ imtmp2 = f32( ( realpha * im1 ) - ( imalpha * re1 ) );
651651
652- re = f32( tmp1_re + tmp2_re );
653- im = f32( tmp1_im + tmp2_im );
654- idx = ia + ( i0 * sa0 ) + ( i1 * sa1 );
655- viewA[ idx ] = f32( viewA[ idx ] + re );
656- viewA[ idx + 1 ] = f32( viewA[ idx + 1 ] + im );
657- if ( i0 === i1 ) {
658- viewA[ idx + 1 ] = 0.0;
659- }
660- }
661- }
662- return A;
652+ re = f32( retmp1 + retmp2 );
653+ im = f32( imtmp1 + imtmp2 );
654+ idx = ia + ( i0 * sa0 ) + ( i1 * sa1 );
655+ viewA[ idx ] = f32( viewA[ idx ] + re );
656+ viewA[ idx + 1 ] = f32( viewA[ idx + 1 ] + im );
657+ if ( i0 === i1 ) {
658+ viewA[ idx + 1 ] = 0.0;
659+ }
660+ }
661+ }
662+ return A;
663663}
664664
665665
@@ -673,7 +673,7 @@ <h1><a href="../../../../index.html">All files</a> / <a href="index.html">blas/b
673673 < div class ='footer quiet pad2 space-top1 center small '>
674674 Code coverage generated by
675675 < a href ="https://istanbul.js.org/ " target ="_blank " rel ="noopener noreferrer "> istanbul</ a >
676- at 2025-06-06T11:45:25.659Z
676+ at 2025-06-06T11:57:50.795Z
677677 </ div >
678678 < script src ="../../../../prettify.js "> </ script >
679679 < script >
0 commit comments