# HG changeset patch
# User alexc
# Date 1264201678 0
# Node ID abb3b23bda358166a03f9d116d101f76a96f04ba
# Parent  81033a08013654fc956b9680b4f1f72576f6a84a
Implement an sse version of scalarproduct_float().

diff -r 81033a080136 -r abb3b23bda35 x86/dsputil_mmx.c
--- a/x86/dsputil_mmx.c	Fri Jan 22 19:07:44 2010 +0000
+++ b/x86/dsputil_mmx.c	Fri Jan 22 23:07:58 2010 +0000
@@ -2510,6 +2510,8 @@
                                   int src_x, int src_y, int src_stride, slice_buffer * sb, int add, uint8_t * dst8);
 
 
+float ff_scalarproduct_float_sse(const float *v1, const float *v2, int order);
+
 void dsputil_init_mmx(DSPContext* c, AVCodecContext *avctx)
 {
     mm_flags = mm_support();
@@ -2965,6 +2967,9 @@
             c->vector_clipf = vector_clipf_sse;
             c->float_to_int16 = float_to_int16_sse;
             c->float_to_int16_interleave = float_to_int16_interleave_sse;
+#if HAVE_YASM
+            c->scalarproduct_float = ff_scalarproduct_float_sse;
+#endif
         }
         if(mm_flags & FF_MM_3DNOW)
             c->vector_fmul_add = vector_fmul_add_3dnow; // faster than sse
diff -r 81033a080136 -r abb3b23bda35 x86/dsputil_yasm.asm
--- a/x86/dsputil_yasm.asm	Fri Jan 22 19:07:44 2010 +0000
+++ b/x86/dsputil_yasm.asm	Fri Jan 22 23:07:58 2010 +0000
@@ -397,3 +397,27 @@
 .unaligned:
     ADD_HFYU_LEFT_LOOP 0
 
+
+; float ff_scalarproduct_float_sse(const float *v1, const float *v2, int len)
+cglobal scalarproduct_float_sse, 3,3,2, v1, v2, offset
+    neg offsetq
+    shl offsetq, 2
+    sub v1q, offsetq
+    sub v2q, offsetq
+    xorps xmm0, xmm0
+    .loop:
+        movaps   xmm1, [v1q+offsetq]
+        mulps    xmm1, [v2q+offsetq]
+        addps    xmm0, xmm1
+        add      offsetq, 16
+        js       .loop
+    movhlps xmm1, xmm0
+    addps   xmm0, xmm1
+    movss   xmm1, xmm0
+    shufps  xmm0, xmm0, 1
+    addss   xmm0, xmm1
+%ifndef ARCH_X86_64
+    movd    r0m,  xmm0
+    fld     dword r0m
+%endif
+    RET