make trunk,branches,tag svn structure

git-svn-id: svn://mattst88.com/svn/alpha_mmintrin/trunk@2 f7d56953-e76f-4e43-a77e-20d50f6c004e

1 files changed, 92 insertions, 0 deletions

diff --git a/alpha_mmintrin.h b/alpha_mmintrin.h
new file mode 100644
index 0000000..34a2b63
--- /dev/null
+++ b/alpha_mmintrin.h
@@ -0,0 +1,92 @@
+#include <stdint.h>
+
+
+#if (defined(__GNUC__)						  \
+     && (__GNUC__ >= 4 || (__GNUC__ == 3 && __GNUC_MINOR__ >= 4)) \
+     && __alpha_max__)
+
+#define __bic		__builtin_alpha_bic
+#define __eqv		__builtin_alpha_eqv
+#define __ornot		__builtin_alpha_ornot
+#define __cmpbge	__builtin_alpha_cmpbge
+#define __extql		__builtin_alpha_extql	
+#define __extqh		__builtin_alpha_extqh	
+#define __extbl		__builtin_alpha_extbl
+#define __extwl		__builtin_alpha_extwl
+#define __extwh		__builtin_alpha_extwh
+#define __inswl		__builtin_alpha_inswl
+#define __zap		__builtin_alpha_zap	
+#define __zapnot	__builtin_alpha_zapnot
+#define __amask		__builtin_alpha_amask	
+#define __implver	__builtin_alpha_implver
+#define __rpcc		__builtin_alpha_rpcc	
+#define __minub8	__builtin_alpha_minub8
+#define __minsb8	__builtin_alpha_minsb8
+#define __minuw4	__builtin_alpha_minuw4
+#define __minsw4	__builtin_alpha_minsw4
+#define __maxub8	__builtin_alpha_maxub8
+#define __maxsb8	__builtin_alpha_maxsb8
+#define __maxuw4	__builtin_alpha_maxuw4
+#define __maxsw4	__builtin_alpha_maxsw4
+#define __perr		__builtin_alpha_perr	
+#define __pklb		__builtin_alpha_pklb	
+#define __pkwb		__builtin_alpha_pkwb	
+#define __unpkbl	__builtin_alpha_unpkbl
+#define __unpkbw	__builtin_alpha_unpkbw
+
+static __inline uint64_t
+__ctlz(uint64_t val) {
+	uint64_t ret;
+	asm volatile(
+		"ctlz	%1,%0"
+		: "=r" (ret)
+		: "r" (val)
+	);
+	return ret;
+}
+
+/* Add the 16-bit values in M1 to the 16-bit values in M2 (paddw).  */
+static __inline uint64_t
+__adduw4(uint64_t m1, uint64_t m2) {
+	uint64_t signs = (m1 ^ m2) & 0x8000800080008000;
+	m1 &= ~signs;             // ??? gcc doesn't use bic here
+	m2 &= ~signs;
+	m1 += m2;
+	m1 ^= signs;
+	return m1;
+}
+
+/* Add the 8-bit values in M1 to the 8-bit values in M2 using unsigned
+ * saturated arithmetic (paddusb) */
+static __inline uint64_t
+__addusb8(uint64_t m1, uint64_t m2) {
+	return m1 + __minub8(m2, ~m1);
+}
+
+/* Add the 16-bit values in M1 to the 16-bit values in M2 using unsigned
+ * saturating arithmetic (paddusw) */
+static __inline uint64_t
+__addusw4(uint64_t m1, uint64_t m2) {
+	return m1 + __minuw4(m2, ~m1);
+}
+
+/* Multiply four 16-bit values in M1 by four 16-bit values in M2 and produce
+   the low 16 bits of the results (pmullw).  */
+static __inline uint64_t
+__mulsw4sl(uint64_t m1, uint64_t m2) {
+	uint64_t t0, t2, t4, t6;
+
+	t0 = (int)         m1     * (int)         m2    ;
+	t2 = (int) __extwl(m1, 2) * (int) __extwl(m2, 2);
+	t4 = (int) __extwl(m1, 4) * (int) __extwl(m2, 4);
+	t6 = (int) __extwl(m1, 6) * (int) __extwl(m2, 6);
+
+	t0 = __inswl(t0, 0);
+	t2 = __inswl(t2, 2);
+	t4 = __inswl(t4, 4);
+	t6 = __inswl(t6, 6);
+
+	return t0 | t2 | t4 | t6;
+}
+
+#endif