Use generic DDOT kernel for WASM128_GENERIC

Author: teddygood

kernel/generic/dot.c

@@ -47,11 +47,46 @@ FLOAT CNAME(BLASLONG n, FLOAT *x, BLASLONG inc_x, FLOAT *y, BLASLONG inc_y)
 
 	if ( (inc_x == 1) && (inc_y == 1) )
 	{
-#if V_SIMD && !defined(DSDOT)
-        const int vstep = v_nlanes_f32;
-        const int unrollx4 = n & (-vstep * 4);
-        const int unrollx  = n &  -vstep;
-		v_f32 vsum0 = v_zero_f32();
+	#if defined(DOUBLE) && V_SIMD && V_SIMD_F64 && !defined(DSDOT)
+	        const int vstep = v_nlanes_f64;
+	        const int unrollx4 = n & (-vstep * 4);
+	        const int unrollx  = n &  -vstep;
+			v_f64 vsum0 = v_zero_f64();
+	        v_f64 vsum1 = v_zero_f64();
+	        v_f64 vsum2 = v_zero_f64();
+	        v_f64 vsum3 = v_zero_f64();
+			while(i < unrollx4)
+	        {
+	            vsum0 = v_muladd_f64(
+	                v_loadu_f64(x + i),           v_loadu_f64(y + i),           vsum0
+	            );
+	            vsum1 = v_muladd_f64(
+	                v_loadu_f64(x + i + vstep),   v_loadu_f64(y + i + vstep),   vsum1
+	            );
+	            vsum2 = v_muladd_f64(
+	                v_loadu_f64(x + i + vstep*2), v_loadu_f64(y + i + vstep*2), vsum2
+	            );
+	            vsum3 = v_muladd_f64(
+	                v_loadu_f64(x + i + vstep*3), v_loadu_f64(y + i + vstep*3), vsum3
+	            );
+	            i += vstep*4;
+	        }
+	        vsum0 = v_add_f64(
+	            v_add_f64(vsum0, vsum1), v_add_f64(vsum2 , vsum3)
+	        );
+			while(i < unrollx)
+	        {
+	            vsum0 = v_muladd_f64(
+	                v_loadu_f64(x + i), v_loadu_f64(y + i), vsum0
+	            );
+	            i += vstep;
+	        }
+	        dot = v_sum_f64(vsum0);
+	#elif V_SIMD && !defined(DSDOT)
+	        const int vstep = v_nlanes_f32;
+	        const int unrollx4 = n & (-vstep * 4);
+	        const int unrollx  = n &  -vstep;
+			v_f32 vsum0 = v_zero_f32();
         v_f32 vsum1 = v_zero_f32();
         v_f32 vsum2 = v_zero_f32();
         v_f32 vsum3 = v_zero_f32();
@@ -82,10 +117,10 @@ FLOAT CNAME(BLASLONG n, FLOAT *x, BLASLONG inc_x, FLOAT *y, BLASLONG inc_y)
             i += vstep;
         }
         dot = v_sum_f32(vsum0);
-#elif defined(DSDOT)
-        int n1 = n & -4;
-		for (; i < n1; i += 4)
-		{
+	#elif defined(DSDOT)
+	        int n1 = n & -4;
+			for (; i < n1; i += 4)
+			{
 			dot += (double) y[i] * (double) x[i]
 			    + (double) y[i+1] * (double) x[i+1]
 			    + (double) y[i+2] * (double) x[i+2]
@@ -133,5 +168,3 @@ FLOAT CNAME(BLASLONG n, FLOAT *x, BLASLONG inc_x, FLOAT *y, BLASLONG inc_y)
 	return(dot);
 
 }
-
-

kernel/wasm/KERNEL.WASM128_GENERIC

@@ -56,7 +56,7 @@ CCOPYKERNEL  = ../riscv64/zcopy.c
 ZCOPYKERNEL  = ../riscv64/zcopy.c
 
 SDOTKERNEL   = ../generic/dot.c
-DDOTKERNEL   = ../riscv64/dot.c
+DDOTKERNEL   = ../generic/dot.c
 CDOTKERNEL   = ../riscv64/zdot.c
 ZDOTKERNEL   = ../riscv64/zdot.c
 DSDOTKERNEL  = ../generic/dot.c