@@ -70,7 +70,6 @@ int CNAME(BLASLONG rows, BLASLONG cols, FLOAT alpha_r, FLOAT alpha_i, FLOAT *a,
7070 FLOAT_VX2_T va , vb ;
7171 unsigned int gvl = 0 ;
7272
73-
7473 if ( rows <= 0 ) return (0 );
7574 if ( cols <= 0 ) return (0 );
7675
@@ -85,17 +84,13 @@ int CNAME(BLASLONG rows, BLASLONG cols, FLOAT alpha_r, FLOAT alpha_i, FLOAT *a,
8584 for (j = 0 ; j < rows ; j += gvl )
8685 {
8786 gvl = VSETVL (rows - j );
88- // bptr[ia + 0] = alpha_r * aptr[ia + 0] - alpha_i * aptr[ia+1];
89- // bptr[ia + 2] = alpha_r * aptr[ia + 2] - alpha_i * aptr[ia+3];
9087 va = VLSEG2_FLOAT (aptr + ia , gvl );
9188 aptr_v0 = VGET_VX2 (va , 0 );
9289 aptr_v1 = VGET_VX2 (va , 1 );
9390 bptr_v1 = VFMUL_VF_FLOAT ( aptr_v1 , alpha_r ,gvl );
9491 bptr_v1 = VFMACCVF_FLOAT (bptr_v1 , alpha_i , aptr_v0 , gvl );
9592 bptr_v0 = VFMUL_VF_FLOAT ( aptr_v0 ,alpha_r , gvl );
9693 bptr_v0 = VFNMSACVF_FLOAT (bptr_v0 , alpha_i , aptr_v1 , gvl );
97- // bptr[ia + 1] = alpha_r * aptr[ia + 1] + alpha_i * aptr[ia+0];
98- // bptr[ia + 3] = alpha_r * aptr[ia + 3] + alpha_i * aptr[ia+2];
9994 vb = VSET_VX2 (vb , 0 , bptr_v0 );
10095 vb = VSET_VX2 (vb , 1 , bptr_v1 );
10196 VSSEG2_FLOAT (& bptr [ia ], vb , gvl );
0 commit comments