--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/genericopenlibs/liboil/src/multsum_sse.c	Fri Apr 16 16:46:38 2010 +0300
@@ -0,0 +1,96 @@
+/*
+* Copyright (c) 2009 Nokia Corporation and/or its subsidiary(-ies).
+* All rights reserved.
+* This component and the accompanying materials are made available
+* under the terms of "Eclipse Public License v1.0"
+* which accompanies this distribution, and is available
+* at the URL "http://www.eclipse.org/legal/epl-v10.html".
+*
+* Initial Contributors:
+* Nokia Corporation - initial contribution.
+*
+* Contributors:
+*
+* Description: 
+*
+*/
+
+#ifdef HAVE_CONFIG_H
+#include "config.h"
+#endif
+#include <liboil/liboilclasses.h>
+#include <liboil/liboilfunction.h>
+#include <emmintrin.h>
+
+#define SSE_FUNCTION __attribute__((force_align_arg_pointer))
+
+#define MULTSUM_SSE2_NSTRIDED(i) { \
+  t1 = _mm_load_pd(&OIL_GET(src1, i, double)); \
+  t2 = _mm_load_pd(&OIL_GET(src2, i, double)); \
+  t1 = _mm_mul_pd(t1,t2); \
+  sum.reg = _mm_add_pd(sum.reg,t1); \
+}
+#define MULTSUM_SSE2_NSTRIDEDP(i) { \
+  t1 = _mm_load_pd(&OIL_GET(src1, i*sstr1, double)); \
+  t2 = _mm_loadl_pd(t2, &OIL_GET(src2, i*sstr2, double)); \
+  t2 = _mm_loadh_pd(t2, &OIL_GET(src2, (i+1)*sstr2, double)); \
+  t1 = _mm_mul_pd(t1,t2); \
+  sum.reg = _mm_add_pd(sum.reg,t1); \
+}
+#define MULTSUM_SSE2_STRIDED(i) { \
+  t1 = _mm_loadl_pd(t1, &OIL_GET(src1, i*sstr1, double)); \
+  t1 = _mm_loadh_pd(t1, &OIL_GET(src1, (i+1)*sstr1, double)); \
+  t2 = _mm_loadl_pd(t2, &OIL_GET(src2, i*sstr2, double)); \
+  t2 = _mm_loadh_pd(t2, &OIL_GET(src2, (i+1)*sstr2, double)); \
+  t1 = _mm_mul_pd(t1,t2); \
+  sum.reg = _mm_add_pd(sum.reg,t1); \
+}
+
+
+#ifdef ENABLE_BROKEN_IMPLS
+SSE_FUNCTION static void
+multsum_f64_sse2_unroll4(double *dest,
+     const double *src1, int sstr1,
+     const double *src2, int sstr2,
+     int n)
+{
+  __m128d t1, t2;
+  union {
+    __m128d reg;
+    double vals[2];
+  } sum;
+  int i = 0;
+
+  sum.reg = _mm_setzero_pd();
+  while (i < n-3) {
+    MULTSUM_SSE2_STRIDED(0);
+    MULTSUM_SSE2_STRIDED(2);
+
+    OIL_INCREMENT(src1, 4*sstr1);
+    OIL_INCREMENT(src2, 4*sstr2);
+    i += 4;
+  }
+  while (i < n-1) {
+    MULTSUM_SSE2_STRIDED(0);
+
+    OIL_INCREMENT(src1, 2*sstr1);
+    OIL_INCREMENT(src2, 2*sstr2);
+    i+=2;
+  }
+  *dest = sum.vals[0] + sum.vals[1];
+  if (i < n) {
+    *dest += (OIL_GET(src1,0,double)*OIL_GET(src2,0,double));
+  }
+}
+OIL_DEFINE_IMPL_FULL (multsum_f64_sse2_unroll4, multsum_f64, OIL_IMPL_FLAG_SSE2);
+#endif
+
+
+
+#ifdef	__SYMBIAN32__
+ 
+OilFunctionImpl* __oil_function_impl_multsum_f64_sse2_unroll4, multsum_f64() {
+		return &_oil_function_impl_multsum_f64_sse2_unroll4, multsum_f64;
+}
+#endif
+