1 /// @ref simd
2 /// @file glm/simd/integer.h
3 
4 #pragma once
5 
6 #if GLM_ARCH & GLM_ARCH_SSE2_BIT
7 
glm_i128_interleave(glm_uvec4 x)8 GLM_FUNC_QUALIFIER glm_uvec4 glm_i128_interleave(glm_uvec4 x)
9 {
10 	glm_uvec4 const Mask4 = _mm_set1_epi32(0x0000FFFF);
11 	glm_uvec4 const Mask3 = _mm_set1_epi32(0x00FF00FF);
12 	glm_uvec4 const Mask2 = _mm_set1_epi32(0x0F0F0F0F);
13 	glm_uvec4 const Mask1 = _mm_set1_epi32(0x33333333);
14 	glm_uvec4 const Mask0 = _mm_set1_epi32(0x55555555);
15 
16 	glm_uvec4 Reg1;
17 	glm_uvec4 Reg2;
18 
19 	// REG1 = x;
20 	// REG2 = y;
21 	//Reg1 = _mm_unpacklo_epi64(x, y);
22 	Reg1 = x;
23 
24 	//REG1 = ((REG1 << 16) | REG1) & glm::uint64(0x0000FFFF0000FFFF);
25 	//REG2 = ((REG2 << 16) | REG2) & glm::uint64(0x0000FFFF0000FFFF);
26 	Reg2 = _mm_slli_si128(Reg1, 2);
27 	Reg1 = _mm_or_si128(Reg2, Reg1);
28 	Reg1 = _mm_and_si128(Reg1, Mask4);
29 
30 	//REG1 = ((REG1 <<  8) | REG1) & glm::uint64(0x00FF00FF00FF00FF);
31 	//REG2 = ((REG2 <<  8) | REG2) & glm::uint64(0x00FF00FF00FF00FF);
32 	Reg2 = _mm_slli_si128(Reg1, 1);
33 	Reg1 = _mm_or_si128(Reg2, Reg1);
34 	Reg1 = _mm_and_si128(Reg1, Mask3);
35 
36 	//REG1 = ((REG1 <<  4) | REG1) & glm::uint64(0x0F0F0F0F0F0F0F0F);
37 	//REG2 = ((REG2 <<  4) | REG2) & glm::uint64(0x0F0F0F0F0F0F0F0F);
38 	Reg2 = _mm_slli_epi32(Reg1, 4);
39 	Reg1 = _mm_or_si128(Reg2, Reg1);
40 	Reg1 = _mm_and_si128(Reg1, Mask2);
41 
42 	//REG1 = ((REG1 <<  2) | REG1) & glm::uint64(0x3333333333333333);
43 	//REG2 = ((REG2 <<  2) | REG2) & glm::uint64(0x3333333333333333);
44 	Reg2 = _mm_slli_epi32(Reg1, 2);
45 	Reg1 = _mm_or_si128(Reg2, Reg1);
46 	Reg1 = _mm_and_si128(Reg1, Mask1);
47 
48 	//REG1 = ((REG1 <<  1) | REG1) & glm::uint64(0x5555555555555555);
49 	//REG2 = ((REG2 <<  1) | REG2) & glm::uint64(0x5555555555555555);
50 	Reg2 = _mm_slli_epi32(Reg1, 1);
51 	Reg1 = _mm_or_si128(Reg2, Reg1);
52 	Reg1 = _mm_and_si128(Reg1, Mask0);
53 
54 	//return REG1 | (REG2 << 1);
55 	Reg2 = _mm_slli_epi32(Reg1, 1);
56 	Reg2 = _mm_srli_si128(Reg2, 8);
57 	Reg1 = _mm_or_si128(Reg1, Reg2);
58 
59 	return Reg1;
60 }
61 
glm_i128_interleave2(glm_uvec4 x,glm_uvec4 y)62 GLM_FUNC_QUALIFIER glm_uvec4 glm_i128_interleave2(glm_uvec4 x, glm_uvec4 y)
63 {
64 	glm_uvec4 const Mask4 = _mm_set1_epi32(0x0000FFFF);
65 	glm_uvec4 const Mask3 = _mm_set1_epi32(0x00FF00FF);
66 	glm_uvec4 const Mask2 = _mm_set1_epi32(0x0F0F0F0F);
67 	glm_uvec4 const Mask1 = _mm_set1_epi32(0x33333333);
68 	glm_uvec4 const Mask0 = _mm_set1_epi32(0x55555555);
69 
70 	glm_uvec4 Reg1;
71 	glm_uvec4 Reg2;
72 
73 	// REG1 = x;
74 	// REG2 = y;
75 	Reg1 = _mm_unpacklo_epi64(x, y);
76 
77 	//REG1 = ((REG1 << 16) | REG1) & glm::uint64(0x0000FFFF0000FFFF);
78 	//REG2 = ((REG2 << 16) | REG2) & glm::uint64(0x0000FFFF0000FFFF);
79 	Reg2 = _mm_slli_si128(Reg1, 2);
80 	Reg1 = _mm_or_si128(Reg2, Reg1);
81 	Reg1 = _mm_and_si128(Reg1, Mask4);
82 
83 	//REG1 = ((REG1 <<  8) | REG1) & glm::uint64(0x00FF00FF00FF00FF);
84 	//REG2 = ((REG2 <<  8) | REG2) & glm::uint64(0x00FF00FF00FF00FF);
85 	Reg2 = _mm_slli_si128(Reg1, 1);
86 	Reg1 = _mm_or_si128(Reg2, Reg1);
87 	Reg1 = _mm_and_si128(Reg1, Mask3);
88 
89 	//REG1 = ((REG1 <<  4) | REG1) & glm::uint64(0x0F0F0F0F0F0F0F0F);
90 	//REG2 = ((REG2 <<  4) | REG2) & glm::uint64(0x0F0F0F0F0F0F0F0F);
91 	Reg2 = _mm_slli_epi32(Reg1, 4);
92 	Reg1 = _mm_or_si128(Reg2, Reg1);
93 	Reg1 = _mm_and_si128(Reg1, Mask2);
94 
95 	//REG1 = ((REG1 <<  2) | REG1) & glm::uint64(0x3333333333333333);
96 	//REG2 = ((REG2 <<  2) | REG2) & glm::uint64(0x3333333333333333);
97 	Reg2 = _mm_slli_epi32(Reg1, 2);
98 	Reg1 = _mm_or_si128(Reg2, Reg1);
99 	Reg1 = _mm_and_si128(Reg1, Mask1);
100 
101 	//REG1 = ((REG1 <<  1) | REG1) & glm::uint64(0x5555555555555555);
102 	//REG2 = ((REG2 <<  1) | REG2) & glm::uint64(0x5555555555555555);
103 	Reg2 = _mm_slli_epi32(Reg1, 1);
104 	Reg1 = _mm_or_si128(Reg2, Reg1);
105 	Reg1 = _mm_and_si128(Reg1, Mask0);
106 
107 	//return REG1 | (REG2 << 1);
108 	Reg2 = _mm_slli_epi32(Reg1, 1);
109 	Reg2 = _mm_srli_si128(Reg2, 8);
110 	Reg1 = _mm_or_si128(Reg1, Reg2);
111 
112 	return Reg1;
113 }
114 
115 #endif//GLM_ARCH & GLM_ARCH_SSE2_BIT
116