master
ブランチ master
VulkanExample / dep / include / glm / gtc / packing.inl
25.8 KB、952行 生 Blame 履歴
1/// @ref gtc_packing
2
3#include "../ext/scalar_relational.hpp"
4#include "../ext/vector_relational.hpp"
5#include "../common.hpp"
6#include "../vec2.hpp"
7#include "../vec3.hpp"
8#include "../vec4.hpp"
9#include "../detail/type_half.hpp"
10#include "type_ptr.hpp"
11#include <cstring>
12#include <limits>
13
14namespace glm{
15namespace detail
16{
17 GLM_FUNC_QUALIFIER glm::uint16 float2half(glm::uint32 f)
18 {
19 // 10 bits => EE EEEFFFFF
20 // 11 bits => EEE EEFFFFFF
21 // Half bits => SEEEEEFF FFFFFFFF
22 // Float bits => SEEEEEEE EFFFFFFF FFFFFFFF FFFFFFFF
23
24 // 0x00007c00 => 00000000 00000000 01111100 00000000
25 // 0x000003ff => 00000000 00000000 00000011 11111111
26 // 0x38000000 => 00111000 00000000 00000000 00000000
27 // 0x7f800000 => 01111111 10000000 00000000 00000000
28 // 0x00008000 => 00000000 00000000 10000000 00000000
29 return
30 ((f >> 16) & 0x8000) | // sign
31 ((((f & 0x7f800000) - 0x38000000) >> 13) & 0x7c00) | // exponential
32 ((f >> 13) & 0x03ff); // Mantissa
33 }
34
35 GLM_FUNC_QUALIFIER glm::uint32 float2packed11(glm::uint32 f)
36 {
37 // 10 bits => EE EEEFFFFF
38 // 11 bits => EEE EEFFFFFF
39 // Half bits => SEEEEEFF FFFFFFFF
40 // Float bits => SEEEEEEE EFFFFFFF FFFFFFFF FFFFFFFF
41
42 // 0x000007c0 => 00000000 00000000 00000111 11000000
43 // 0x00007c00 => 00000000 00000000 01111100 00000000
44 // 0x000003ff => 00000000 00000000 00000011 11111111
45 // 0x38000000 => 00111000 00000000 00000000 00000000
46 // 0x7f800000 => 01111111 10000000 00000000 00000000
47 // 0x00008000 => 00000000 00000000 10000000 00000000
48 return
49 ((((f & 0x7f800000) - 0x38000000) >> 17) & 0x07c0) | // exponential
50 ((f >> 17) & 0x003f); // Mantissa
51 }
52
53 GLM_FUNC_QUALIFIER glm::uint32 packed11ToFloat(glm::uint32 p)
54 {
55 // 10 bits => EE EEEFFFFF
56 // 11 bits => EEE EEFFFFFF
57 // Half bits => SEEEEEFF FFFFFFFF
58 // Float bits => SEEEEEEE EFFFFFFF FFFFFFFF FFFFFFFF
59
60 // 0x000007c0 => 00000000 00000000 00000111 11000000
61 // 0x00007c00 => 00000000 00000000 01111100 00000000
62 // 0x000003ff => 00000000 00000000 00000011 11111111
63 // 0x38000000 => 00111000 00000000 00000000 00000000
64 // 0x7f800000 => 01111111 10000000 00000000 00000000
65 // 0x00008000 => 00000000 00000000 10000000 00000000
66 return
67 ((((p & 0x07c0) << 17) + 0x38000000) & 0x7f800000) | // exponential
68 ((p & 0x003f) << 17); // Mantissa
69 }
70
71 GLM_FUNC_QUALIFIER glm::uint32 float2packed10(glm::uint32 f)
72 {
73 // 10 bits => EE EEEFFFFF
74 // 11 bits => EEE EEFFFFFF
75 // Half bits => SEEEEEFF FFFFFFFF
76 // Float bits => SEEEEEEE EFFFFFFF FFFFFFFF FFFFFFFF
77
78 // 0x0000001F => 00000000 00000000 00000000 00011111
79 // 0x0000003F => 00000000 00000000 00000000 00111111
80 // 0x000003E0 => 00000000 00000000 00000011 11100000
81 // 0x000007C0 => 00000000 00000000 00000111 11000000
82 // 0x00007C00 => 00000000 00000000 01111100 00000000
83 // 0x000003FF => 00000000 00000000 00000011 11111111
84 // 0x38000000 => 00111000 00000000 00000000 00000000
85 // 0x7f800000 => 01111111 10000000 00000000 00000000
86 // 0x00008000 => 00000000 00000000 10000000 00000000
87 return
88 ((((f & 0x7f800000) - 0x38000000) >> 18) & 0x03E0) | // exponential
89 ((f >> 18) & 0x001f); // Mantissa
90 }
91
92 GLM_FUNC_QUALIFIER glm::uint32 packed10ToFloat(glm::uint32 p)
93 {
94 // 10 bits => EE EEEFFFFF
95 // 11 bits => EEE EEFFFFFF
96 // Half bits => SEEEEEFF FFFFFFFF
97 // Float bits => SEEEEEEE EFFFFFFF FFFFFFFF FFFFFFFF
98
99 // 0x0000001F => 00000000 00000000 00000000 00011111
100 // 0x0000003F => 00000000 00000000 00000000 00111111
101 // 0x000003E0 => 00000000 00000000 00000011 11100000
102 // 0x000007C0 => 00000000 00000000 00000111 11000000
103 // 0x00007C00 => 00000000 00000000 01111100 00000000
104 // 0x000003FF => 00000000 00000000 00000011 11111111
105 // 0x38000000 => 00111000 00000000 00000000 00000000
106 // 0x7f800000 => 01111111 10000000 00000000 00000000
107 // 0x00008000 => 00000000 00000000 10000000 00000000
108 return
109 ((((p & 0x03E0) << 18) + 0x38000000) & 0x7f800000) | // exponential
110 ((p & 0x001f) << 18); // Mantissa
111 }
112
113 GLM_FUNC_QUALIFIER glm::uint half2float(glm::uint h)
114 {
115 return ((h & 0x8000) << 16) | ((( h & 0x7c00) + 0x1C000) << 13) | ((h & 0x03FF) << 13);
116 }
117
118 GLM_FUNC_QUALIFIER glm::uint floatTo11bit(float x)
119 {
120 if(x == 0.0f)
121 return 0u;
122 else if(glm::isnan(x))
123 return ~0u;
124 else if(glm::isinf(x))
125 return 0x1Fu << 6u;
126
127 uint Pack = 0u;
128 memcpy(&Pack, &x, sizeof(Pack));
129 return float2packed11(Pack);
130 }
131
132 GLM_FUNC_QUALIFIER float packed11bitToFloat(glm::uint x)
133 {
134 if(x == 0)
135 return 0.0f;
136 else if(x == ((1 << 11) - 1))
137 return ~0;//NaN
138 else if(x == (0x1f << 6))
139 return ~0;//Inf
140
141 uint Result = packed11ToFloat(x);
142
143 float Temp = 0;
144 memcpy(&Temp, &Result, sizeof(Temp));
145 return Temp;
146 }
147
148 GLM_FUNC_QUALIFIER glm::uint floatTo10bit(float x)
149 {
150 if(x == 0.0f)
151 return 0u;
152 else if(glm::isnan(x))
153 return ~0u;
154 else if(glm::isinf(x))
155 return 0x1Fu << 5u;
156
157 uint Pack = 0;
158 memcpy(&Pack, &x, sizeof(Pack));
159 return float2packed10(Pack);
160 }
161
162 GLM_FUNC_QUALIFIER float packed10bitToFloat(glm::uint x)
163 {
164 if(x == 0)
165 return 0.0f;
166 else if(x == ((1 << 10) - 1))
167 return ~0;//NaN
168 else if(x == (0x1f << 5))
169 return ~0;//Inf
170
171 uint Result = packed10ToFloat(x);
172
173 float Temp = 0;
174 memcpy(&Temp, &Result, sizeof(Temp));
175 return Temp;
176 }
177
178// GLM_FUNC_QUALIFIER glm::uint f11_f11_f10(float x, float y, float z)
179// {
180// return ((floatTo11bit(x) & ((1 << 11) - 1)) << 0) | ((floatTo11bit(y) & ((1 << 11) - 1)) << 11) | ((floatTo10bit(z) & ((1 << 10) - 1)) << 22);
181// }
182
183#if GLM_SILENT_WARNINGS == GLM_ENABLE
184# if defined(__clang__)
185# pragma clang diagnostic push
186# pragma clang diagnostic ignored "-Wpadded"
187# endif
188#endif
189
190 union u3u3u2
191 {
192 struct Data
193 {
194 uint x : 3;
195 uint y : 3;
196 uint z : 2;
197 } data;
198 uint8 pack;
199 };
200
201 union u4u4
202 {
203 struct Data
204 {
205 uint x : 4;
206 uint y : 4;
207 } data;
208 uint8 pack;
209 };
210
211 union u4u4u4u4
212 {
213 struct Data
214 {
215 uint x : 4;
216 uint y : 4;
217 uint z : 4;
218 uint w : 4;
219 } data;
220 uint16 pack;
221 };
222
223 union u5u6u5
224 {
225 struct Data
226 {
227 uint x : 5;
228 uint y : 6;
229 uint z : 5;
230 } data;
231 uint16 pack;
232 };
233
234 union u5u5u5u1
235 {
236 struct Data
237 {
238 uint x : 5;
239 uint y : 5;
240 uint z : 5;
241 uint w : 1;
242 } data;
243 uint16 pack;
244 };
245
246#if GLM_SILENT_WARNINGS == GLM_ENABLE
247# if defined(__clang__)
248# pragma clang diagnostic pop
249# endif
250#endif
251
252 union u10u10u10u2
253 {
254 struct Data
255 {
256 uint x : 10;
257 uint y : 10;
258 uint z : 10;
259 uint w : 2;
260 } data;
261 uint32 pack;
262 };
263
264 union i10i10i10i2
265 {
266 struct Data
267 {
268 int x : 10;
269 int y : 10;
270 int z : 10;
271 int w : 2;
272 } data;
273 uint32 pack;
274 };
275
276 union u9u9u9e5
277 {
278 struct Data
279 {
280 uint x : 9;
281 uint y : 9;
282 uint z : 9;
283 uint w : 5;
284 } data;
285 uint32 pack;
286 };
287
288 template<length_t L, qualifier Q>
289 struct compute_half
290 {};
291
292 template<qualifier Q>
293 struct compute_half<1, Q>
294 {
295 GLM_FUNC_QUALIFIER static vec<1, uint16, Q> pack(vec<1, float, Q> const& v)
296 {
297 int16 const Unpack(detail::toFloat16(v.x));
298 u16vec1 Packed;
299 memcpy(value_ptr(Packed), &Unpack, sizeof(Packed));
300 return Packed;
301 }
302
303 GLM_FUNC_QUALIFIER static vec<1, float, Q> unpack(vec<1, uint16, Q> const& v)
304 {
305 i16vec1 Unpack;
306 memcpy(value_ptr(Unpack), value_ptr(v), sizeof(Unpack));
307 return vec<1, float, Q>(detail::toFloat32(v.x));
308 }
309 };
310
311 template<qualifier Q>
312 struct compute_half<2, Q>
313 {
314 GLM_FUNC_QUALIFIER static vec<2, uint16, Q> pack(vec<2, float, Q> const& v)
315 {
316 vec<2, int16, Q> const Unpack(detail::toFloat16(v.x), detail::toFloat16(v.y));
317 u16vec2 Packed;
318 memcpy(value_ptr(Packed), value_ptr(Unpack), sizeof(Packed));
319 return Packed;
320 }
321
322 GLM_FUNC_QUALIFIER static vec<2, float, Q> unpack(vec<2, uint16, Q> const& v)
323 {
324 i16vec2 Unpack;
325 memcpy(value_ptr(Unpack), value_ptr(v), sizeof(Unpack));
326 return vec<2, float, Q>(detail::toFloat32(v.x), detail::toFloat32(v.y));
327 }
328 };
329
330 template<qualifier Q>
331 struct compute_half<3, Q>
332 {
333 GLM_FUNC_QUALIFIER static vec<3, uint16, Q> pack(vec<3, float, Q> const& v)
334 {
335 vec<3, int16, Q> const Unpack(detail::toFloat16(v.x), detail::toFloat16(v.y), detail::toFloat16(v.z));
336 u16vec3 Packed;
337 memcpy(value_ptr(Packed), value_ptr(Unpack), sizeof(Packed));
338 return Packed;
339 }
340
341 GLM_FUNC_QUALIFIER static vec<3, float, Q> unpack(vec<3, uint16, Q> const& v)
342 {
343 i16vec3 Unpack;
344 memcpy(value_ptr(Unpack), &v, sizeof(Unpack));
345 return vec<3, float, Q>(detail::toFloat32(v.x), detail::toFloat32(v.y), detail::toFloat32(v.z));
346 }
347 };
348
349 template<qualifier Q>
350 struct compute_half<4, Q>
351 {
352 GLM_FUNC_QUALIFIER static vec<4, uint16, Q> pack(vec<4, float, Q> const& v)
353 {
354 vec<4, int16, Q> const Unpack(detail::toFloat16(v.x), detail::toFloat16(v.y), detail::toFloat16(v.z), detail::toFloat16(v.w));
355 u16vec4 Packed;
356 memcpy(value_ptr(Packed), value_ptr(Unpack), sizeof(Packed));
357 return Packed;
358 }
359
360 GLM_FUNC_QUALIFIER static vec<4, float, Q> unpack(vec<4, uint16, Q> const& v)
361 {
362 i16vec4 Unpack;
363 memcpy(value_ptr(Unpack), &v, sizeof(Unpack));
364 return vec<4, float, Q>(detail::toFloat32(Unpack.x), detail::toFloat32(Unpack.y), detail::toFloat32(Unpack.z), detail::toFloat32(Unpack.w));
365 }
366 };
367}//namespace detail
368
369 GLM_FUNC_QUALIFIER uint8 packUnorm1x8(float v)
370 {
371 return static_cast<uint8>(round(clamp(v, 0.0f, 1.0f) * 255.0f));
372 }
373
374 GLM_FUNC_QUALIFIER float unpackUnorm1x8(uint8 p)
375 {
376 float const Unpack(p);
377 return Unpack * static_cast<float>(0.0039215686274509803921568627451); // 1 / 255
378 }
379
380 GLM_FUNC_QUALIFIER uint16 packUnorm2x8(vec2 const& v)
381 {
382 u8vec2 const Topack(round(clamp(v, 0.0f, 1.0f) * 255.0f));
383
384 uint16 Unpack = 0;
385 memcpy(&Unpack, &Topack, sizeof(Unpack));
386 return Unpack;
387 }
388
389 GLM_FUNC_QUALIFIER vec2 unpackUnorm2x8(uint16 p)
390 {
391 u8vec2 Unpack;
392 memcpy(value_ptr(Unpack), &p, sizeof(Unpack));
393 return vec2(Unpack) * float(0.0039215686274509803921568627451); // 1 / 255
394 }
395
396 GLM_FUNC_QUALIFIER uint8 packSnorm1x8(float v)
397 {
398 int8 const Topack(static_cast<int8>(round(clamp(v ,-1.0f, 1.0f) * 127.0f)));
399 uint8 Packed = 0;
400 memcpy(&Packed, &Topack, sizeof(Packed));
401 return Packed;
402 }
403
404 GLM_FUNC_QUALIFIER float unpackSnorm1x8(uint8 p)
405 {
406 int8 Unpack = 0;
407 memcpy(&Unpack, &p, sizeof(Unpack));
408 return clamp(
409 static_cast<float>(Unpack) * 0.00787401574803149606299212598425f, // 1.0f / 127.0f
410 -1.0f, 1.0f);
411 }
412
413 GLM_FUNC_QUALIFIER uint16 packSnorm2x8(vec2 const& v)
414 {
415 i8vec2 const Topack(round(clamp(v, -1.0f, 1.0f) * 127.0f));
416 uint16 Packed = 0;
417 memcpy(&Packed, value_ptr(Topack), sizeof(Packed));
418 return Packed;
419 }
420
421 GLM_FUNC_QUALIFIER vec2 unpackSnorm2x8(uint16 p)
422 {
423 i8vec2 Unpack;
424 memcpy(value_ptr(Unpack), &p, sizeof(Unpack));
425 return clamp(
426 vec2(Unpack) * 0.00787401574803149606299212598425f, // 1.0f / 127.0f
427 -1.0f, 1.0f);
428 }
429
430 GLM_FUNC_QUALIFIER uint16 packUnorm1x16(float s)
431 {
432 return static_cast<uint16>(round(clamp(s, 0.0f, 1.0f) * 65535.0f));
433 }
434
435 GLM_FUNC_QUALIFIER float unpackUnorm1x16(uint16 p)
436 {
437 float const Unpack(p);
438 return Unpack * 1.5259021896696421759365224689097e-5f; // 1.0 / 65535.0
439 }
440
441 GLM_FUNC_QUALIFIER uint64 packUnorm4x16(vec4 const& v)
442 {
443 u16vec4 const Topack(round(clamp(v , 0.0f, 1.0f) * 65535.0f));
444 uint64 Packed = 0;
445 memcpy(&Packed, &Topack, sizeof(Packed));
446 return Packed;
447 }
448
449 GLM_FUNC_QUALIFIER vec4 unpackUnorm4x16(uint64 p)
450 {
451 u16vec4 Unpack;
452 memcpy(value_ptr(Unpack), &p, sizeof(Unpack));
453 return vec4(Unpack) * 1.5259021896696421759365224689097e-5f; // 1.0 / 65535.0
454 }
455
456 GLM_FUNC_QUALIFIER uint16 packSnorm1x16(float v)
457 {
458 int16 const Topack = static_cast<int16>(round(clamp(v ,-1.0f, 1.0f) * 32767.0f));
459 uint16 Packed = 0;
460 memcpy(&Packed, &Topack, sizeof(Packed));
461 return Packed;
462 }
463
464 GLM_FUNC_QUALIFIER float unpackSnorm1x16(uint16 p)
465 {
466 int16 Unpack = 0;
467 memcpy(&Unpack, &p, sizeof(Unpack));
468 return clamp(
469 static_cast<float>(Unpack) * 3.0518509475997192297128208258309e-5f, //1.0f / 32767.0f,
470 -1.0f, 1.0f);
471 }
472
473 GLM_FUNC_QUALIFIER uint64 packSnorm4x16(vec4 const& v)
474 {
475 i16vec4 const Topack(round(clamp(v ,-1.0f, 1.0f) * 32767.0f));
476 uint64 Packed = 0;
477 memcpy(&Packed, value_ptr(Topack), sizeof(Packed));
478 return Packed;
479 }
480
481 GLM_FUNC_QUALIFIER vec4 unpackSnorm4x16(uint64 p)
482 {
483 i16vec4 Unpack;
484 memcpy(value_ptr(Unpack), &p, sizeof(Unpack));
485 return clamp(
486 vec4(Unpack) * 3.0518509475997192297128208258309e-5f, //1.0f / 32767.0f,
487 -1.0f, 1.0f);
488 }
489
490 GLM_FUNC_QUALIFIER uint16 packHalf1x16(float v)
491 {
492 int16 const Topack(detail::toFloat16(v));
493 uint16 Packed = 0;
494 memcpy(&Packed, &Topack, sizeof(Packed));
495 return Packed;
496 }
497
498 GLM_FUNC_QUALIFIER float unpackHalf1x16(uint16 v)
499 {
500 int16 Unpack = 0;
501 memcpy(&Unpack, &v, sizeof(Unpack));
502 return detail::toFloat32(Unpack);
503 }
504
505 GLM_FUNC_QUALIFIER uint64 packHalf4x16(glm::vec4 const& v)
506 {
507 i16vec4 const Unpack(
508 detail::toFloat16(v.x),
509 detail::toFloat16(v.y),
510 detail::toFloat16(v.z),
511 detail::toFloat16(v.w));
512 uint64 Packed = 0;
513 memcpy(&Packed, value_ptr(Unpack), sizeof(Packed));
514 return Packed;
515 }
516
517 GLM_FUNC_QUALIFIER glm::vec4 unpackHalf4x16(uint64 v)
518 {
519 i16vec4 Unpack;
520 memcpy(value_ptr(Unpack), &v, sizeof(Unpack));
521 return vec4(
522 detail::toFloat32(Unpack.x),
523 detail::toFloat32(Unpack.y),
524 detail::toFloat32(Unpack.z),
525 detail::toFloat32(Unpack.w));
526 }
527
528 GLM_FUNC_QUALIFIER uint32 packI3x10_1x2(ivec4 const& v)
529 {
530 detail::i10i10i10i2 Result;
531 Result.data.x = v.x;
532 Result.data.y = v.y;
533 Result.data.z = v.z;
534 Result.data.w = v.w;
535 return Result.pack;
536 }
537
538 GLM_FUNC_QUALIFIER ivec4 unpackI3x10_1x2(uint32 v)
539 {
540 detail::i10i10i10i2 Unpack;
541 Unpack.pack = v;
542 return ivec4(
543 Unpack.data.x,
544 Unpack.data.y,
545 Unpack.data.z,
546 Unpack.data.w);
547 }
548
549 GLM_FUNC_QUALIFIER uint32 packU3x10_1x2(uvec4 const& v)
550 {
551 detail::u10u10u10u2 Result;
552 Result.data.x = v.x;
553 Result.data.y = v.y;
554 Result.data.z = v.z;
555 Result.data.w = v.w;
556 return Result.pack;
557 }
558
559 GLM_FUNC_QUALIFIER uvec4 unpackU3x10_1x2(uint32 v)
560 {
561 detail::u10u10u10u2 Unpack;
562 Unpack.pack = v;
563 return uvec4(
564 Unpack.data.x,
565 Unpack.data.y,
566 Unpack.data.z,
567 Unpack.data.w);
568 }
569
570 GLM_FUNC_QUALIFIER uint32 packSnorm3x10_1x2(vec4 const& v)
571 {
572 ivec4 const Pack(round(clamp(v,-1.0f, 1.0f) * vec4(511.f, 511.f, 511.f, 1.f)));
573
574 detail::i10i10i10i2 Result;
575 Result.data.x = Pack.x;
576 Result.data.y = Pack.y;
577 Result.data.z = Pack.z;
578 Result.data.w = Pack.w;
579 return Result.pack;
580 }
581
582 GLM_FUNC_QUALIFIER vec4 unpackSnorm3x10_1x2(uint32 v)
583 {
584 detail::i10i10i10i2 Unpack;
585 Unpack.pack = v;
586
587 vec4 const Result(Unpack.data.x, Unpack.data.y, Unpack.data.z, Unpack.data.w);
588
589 return clamp(Result * vec4(1.f / 511.f, 1.f / 511.f, 1.f / 511.f, 1.f), -1.0f, 1.0f);
590 }
591
592 GLM_FUNC_QUALIFIER uint32 packUnorm3x10_1x2(vec4 const& v)
593 {
594 uvec4 const Unpack(round(clamp(v, 0.0f, 1.0f) * vec4(1023.f, 1023.f, 1023.f, 3.f)));
595
596 detail::u10u10u10u2 Result;
597 Result.data.x = Unpack.x;
598 Result.data.y = Unpack.y;
599 Result.data.z = Unpack.z;
600 Result.data.w = Unpack.w;
601 return Result.pack;
602 }
603
604 GLM_FUNC_QUALIFIER vec4 unpackUnorm3x10_1x2(uint32 v)
605 {
606 vec4 const ScaleFactors(1.0f / 1023.f, 1.0f / 1023.f, 1.0f / 1023.f, 1.0f / 3.f);
607
608 detail::u10u10u10u2 Unpack;
609 Unpack.pack = v;
610 return vec4(Unpack.data.x, Unpack.data.y, Unpack.data.z, Unpack.data.w) * ScaleFactors;
611 }
612
613 GLM_FUNC_QUALIFIER uint32 packF2x11_1x10(vec3 const& v)
614 {
615 return
616 ((detail::floatTo11bit(v.x) & ((1 << 11) - 1)) << 0) |
617 ((detail::floatTo11bit(v.y) & ((1 << 11) - 1)) << 11) |
618 ((detail::floatTo10bit(v.z) & ((1 << 10) - 1)) << 22);
619 }
620
621 GLM_FUNC_QUALIFIER vec3 unpackF2x11_1x10(uint32 v)
622 {
623 return vec3(
624 detail::packed11bitToFloat(v >> 0),
625 detail::packed11bitToFloat(v >> 11),
626 detail::packed10bitToFloat(v >> 22));
627 }
628
629 GLM_FUNC_QUALIFIER uint32 packF3x9_E1x5(vec3 const& v)
630 {
631 float const SharedExpMax = (pow(2.0f, 9.0f - 1.0f) / pow(2.0f, 9.0f)) * pow(2.0f, 31.f - 15.f);
632 vec3 const Color = clamp(v, 0.0f, SharedExpMax);
633 float const MaxColor = max(Color.x, max(Color.y, Color.z));
634
635 float const ExpSharedP = max(-15.f - 1.f, floor(log2(MaxColor))) + 1.0f + 15.f;
636 float const MaxShared = floor(MaxColor / pow(2.0f, (ExpSharedP - 15.f - 9.f)) + 0.5f);
637 float const ExpShared = equal(MaxShared, pow(2.0f, 9.0f), epsilon<float>()) ? ExpSharedP + 1.0f : ExpSharedP;
638
639 uvec3 const ColorComp(floor(Color / pow(2.f, (ExpShared - 15.f - 9.f)) + 0.5f));
640
641 detail::u9u9u9e5 Unpack;
642 Unpack.data.x = ColorComp.x;
643 Unpack.data.y = ColorComp.y;
644 Unpack.data.z = ColorComp.z;
645 Unpack.data.w = uint(ExpShared);
646 return Unpack.pack;
647 }
648
649 GLM_FUNC_QUALIFIER vec3 unpackF3x9_E1x5(uint32 v)
650 {
651 detail::u9u9u9e5 Unpack;
652 Unpack.pack = v;
653
654 return vec3(Unpack.data.x, Unpack.data.y, Unpack.data.z) * pow(2.0f, static_cast<float>(Unpack.data.w) - 15.f - 9.f);
655 }
656
657 // Based on Brian Karis http://graphicrants.blogspot.fr/2009/04/rgbm-color-encoding.html
658 template<typename T, qualifier Q>
659 GLM_FUNC_QUALIFIER vec<4, T, Q> packRGBM(vec<3, T, Q> const& rgb)
660 {
661 vec<3, T, Q> const Color(rgb * static_cast<T>(1.0 / 6.0));
662 T Alpha = clamp(max(max(Color.x, Color.y), max(Color.z, static_cast<T>(1e-6))), static_cast<T>(0), static_cast<T>(1));
663 Alpha = ceil(Alpha * static_cast<T>(255.0)) / static_cast<T>(255.0);
664 return vec<4, T, Q>(Color / Alpha, Alpha);
665 }
666
667 template<typename T, qualifier Q>
668 GLM_FUNC_QUALIFIER vec<3, T, Q> unpackRGBM(vec<4, T, Q> const& rgbm)
669 {
670 return vec<3, T, Q>(rgbm.x, rgbm.y, rgbm.z) * rgbm.w * static_cast<T>(6);
671 }
672
673 template<length_t L, qualifier Q>
674 GLM_FUNC_QUALIFIER vec<L, uint16, Q> packHalf(vec<L, float, Q> const& v)
675 {
676 return detail::compute_half<L, Q>::pack(v);
677 }
678
679 template<length_t L, qualifier Q>
680 GLM_FUNC_QUALIFIER vec<L, float, Q> unpackHalf(vec<L, uint16, Q> const& v)
681 {
682 return detail::compute_half<L, Q>::unpack(v);
683 }
684
685 template<typename uintType, length_t L, typename floatType, qualifier Q>
686 GLM_FUNC_QUALIFIER vec<L, uintType, Q> packUnorm(vec<L, floatType, Q> const& v)
687 {
688 static_assert(std::numeric_limits<uintType>::is_integer, "uintType must be an integer type");
689 static_assert(std::numeric_limits<floatType>::is_iec559, "floatType must be a floating point type");
690
691 return vec<L, uintType, Q>(round(clamp(v, static_cast<floatType>(0), static_cast<floatType>(1)) * static_cast<floatType>(std::numeric_limits<uintType>::max())));
692 }
693
694 template<typename floatType, length_t L, typename uintType, qualifier Q>
695 GLM_FUNC_QUALIFIER vec<L, floatType, Q> unpackUnorm(vec<L, uintType, Q> const& v)
696 {
697 static_assert(std::numeric_limits<uintType>::is_integer, "uintType must be an integer type");
698 static_assert(std::numeric_limits<floatType>::is_iec559, "floatType must be a floating point type");
699
700 return vec<L, floatType, Q>(v) * (static_cast<floatType>(1) / static_cast<floatType>(std::numeric_limits<uintType>::max()));
701 }
702
703 template<typename intType, length_t L, typename floatType, qualifier Q>
704 GLM_FUNC_QUALIFIER vec<L, intType, Q> packSnorm(vec<L, floatType, Q> const& v)
705 {
706 static_assert(std::numeric_limits<intType>::is_integer, "uintType must be an integer type");
707 static_assert(std::numeric_limits<floatType>::is_iec559, "floatType must be a floating point type");
708
709 return vec<L, intType, Q>(round(clamp(v , static_cast<floatType>(-1), static_cast<floatType>(1)) * static_cast<floatType>(std::numeric_limits<intType>::max())));
710 }
711
712 template<typename floatType, length_t L, typename intType, qualifier Q>
713 GLM_FUNC_QUALIFIER vec<L, floatType, Q> unpackSnorm(vec<L, intType, Q> const& v)
714 {
715 static_assert(std::numeric_limits<intType>::is_integer, "uintType must be an integer type");
716 static_assert(std::numeric_limits<floatType>::is_iec559, "floatType must be a floating point type");
717
718 return clamp(vec<L, floatType, Q>(v) * (static_cast<floatType>(1) / static_cast<floatType>(std::numeric_limits<intType>::max())), static_cast<floatType>(-1), static_cast<floatType>(1));
719 }
720
721 GLM_FUNC_QUALIFIER uint8 packUnorm2x4(vec2 const& v)
722 {
723 u32vec2 const Unpack(round(clamp(v, 0.0f, 1.0f) * 15.0f));
724 detail::u4u4 Result;
725 Result.data.x = Unpack.x;
726 Result.data.y = Unpack.y;
727 return Result.pack;
728 }
729
730 GLM_FUNC_QUALIFIER vec2 unpackUnorm2x4(uint8 v)
731 {
732 float const ScaleFactor(1.f / 15.f);
733 detail::u4u4 Unpack;
734 Unpack.pack = v;
735 return vec2(Unpack.data.x, Unpack.data.y) * ScaleFactor;
736 }
737
738 GLM_FUNC_QUALIFIER uint16 packUnorm4x4(vec4 const& v)
739 {
740 u32vec4 const Unpack(round(clamp(v, 0.0f, 1.0f) * 15.0f));
741 detail::u4u4u4u4 Result;
742 Result.data.x = Unpack.x;
743 Result.data.y = Unpack.y;
744 Result.data.z = Unpack.z;
745 Result.data.w = Unpack.w;
746 return Result.pack;
747 }
748
749 GLM_FUNC_QUALIFIER vec4 unpackUnorm4x4(uint16 v)
750 {
751 float const ScaleFactor(1.f / 15.f);
752 detail::u4u4u4u4 Unpack;
753 Unpack.pack = v;
754 return vec4(Unpack.data.x, Unpack.data.y, Unpack.data.z, Unpack.data.w) * ScaleFactor;
755 }
756
757 GLM_FUNC_QUALIFIER uint16 packUnorm1x5_1x6_1x5(vec3 const& v)
758 {
759 u32vec3 const Unpack(round(clamp(v, 0.0f, 1.0f) * vec3(31.f, 63.f, 31.f)));
760 detail::u5u6u5 Result;
761 Result.data.x = Unpack.x;
762 Result.data.y = Unpack.y;
763 Result.data.z = Unpack.z;
764 return Result.pack;
765 }
766
767 GLM_FUNC_QUALIFIER vec3 unpackUnorm1x5_1x6_1x5(uint16 v)
768 {
769 vec3 const ScaleFactor(1.f / 31.f, 1.f / 63.f, 1.f / 31.f);
770 detail::u5u6u5 Unpack;
771 Unpack.pack = v;
772 return vec3(Unpack.data.x, Unpack.data.y, Unpack.data.z) * ScaleFactor;
773 }
774
775 GLM_FUNC_QUALIFIER uint16 packUnorm3x5_1x1(vec4 const& v)
776 {
777 u32vec4 const Unpack(round(clamp(v, 0.0f, 1.0f) * vec4(31.f, 31.f, 31.f, 1.f)));
778 detail::u5u5u5u1 Result;
779 Result.data.x = Unpack.x;
780 Result.data.y = Unpack.y;
781 Result.data.z = Unpack.z;
782 Result.data.w = Unpack.w;
783 return Result.pack;
784 }
785
786 GLM_FUNC_QUALIFIER vec4 unpackUnorm3x5_1x1(uint16 v)
787 {
788 vec4 const ScaleFactor(1.f / 31.f, 1.f / 31.f, 1.f / 31.f, 1.f);
789 detail::u5u5u5u1 Unpack;
790 Unpack.pack = v;
791 return vec4(Unpack.data.x, Unpack.data.y, Unpack.data.z, Unpack.data.w) * ScaleFactor;
792 }
793
794 GLM_FUNC_QUALIFIER uint8 packUnorm2x3_1x2(vec3 const& v)
795 {
796 u32vec3 const Unpack(round(clamp(v, 0.0f, 1.0f) * vec3(7.f, 7.f, 3.f)));
797 detail::u3u3u2 Result;
798 Result.data.x = Unpack.x;
799 Result.data.y = Unpack.y;
800 Result.data.z = Unpack.z;
801 return Result.pack;
802 }
803
804 GLM_FUNC_QUALIFIER vec3 unpackUnorm2x3_1x2(uint8 v)
805 {
806 vec3 const ScaleFactor(1.f / 7.f, 1.f / 7.f, 1.f / 3.f);
807 detail::u3u3u2 Unpack;
808 Unpack.pack = v;
809 return vec3(Unpack.data.x, Unpack.data.y, Unpack.data.z) * ScaleFactor;
810 }
811
812 GLM_FUNC_QUALIFIER int16 packInt2x8(i8vec2 const& v)
813 {
814 int16 Pack = 0;
815 memcpy(&Pack, &v, sizeof(Pack));
816 return Pack;
817 }
818
819 GLM_FUNC_QUALIFIER i8vec2 unpackInt2x8(int16 p)
820 {
821 i8vec2 Unpack;
822 memcpy(value_ptr(Unpack), &p, sizeof(Unpack));
823 return Unpack;
824 }
825
826 GLM_FUNC_QUALIFIER uint16 packUint2x8(u8vec2 const& v)
827 {
828 uint16 Pack = 0;
829 memcpy(&Pack, &v, sizeof(Pack));
830 return Pack;
831 }
832
833 GLM_FUNC_QUALIFIER u8vec2 unpackUint2x8(uint16 p)
834 {
835 u8vec2 Unpack;
836 memcpy(value_ptr(Unpack), &p, sizeof(Unpack));
837 return Unpack;
838 }
839
840 GLM_FUNC_QUALIFIER int32 packInt4x8(i8vec4 const& v)
841 {
842 int32 Pack = 0;
843 memcpy(&Pack, &v, sizeof(Pack));
844 return Pack;
845 }
846
847 GLM_FUNC_QUALIFIER i8vec4 unpackInt4x8(int32 p)
848 {
849 i8vec4 Unpack;
850 memcpy(value_ptr(Unpack), &p, sizeof(Unpack));
851 return Unpack;
852 }
853
854 GLM_FUNC_QUALIFIER uint32 packUint4x8(u8vec4 const& v)
855 {
856 uint32 Pack = 0;
857 memcpy(&Pack, &v, sizeof(Pack));
858 return Pack;
859 }
860
861 GLM_FUNC_QUALIFIER u8vec4 unpackUint4x8(uint32 p)
862 {
863 u8vec4 Unpack;
864 memcpy(value_ptr(Unpack), &p, sizeof(Unpack));
865 return Unpack;
866 }
867
868 GLM_FUNC_QUALIFIER int packInt2x16(i16vec2 const& v)
869 {
870 int Pack = 0;
871 memcpy(&Pack, &v, sizeof(Pack));
872 return Pack;
873 }
874
875 GLM_FUNC_QUALIFIER i16vec2 unpackInt2x16(int p)
876 {
877 i16vec2 Unpack;
878 memcpy(value_ptr(Unpack), &p, sizeof(Unpack));
879 return Unpack;
880 }
881
882 GLM_FUNC_QUALIFIER int64 packInt4x16(i16vec4 const& v)
883 {
884 int64 Pack = 0;
885 memcpy(&Pack, &v, sizeof(Pack));
886 return Pack;
887 }
888
889 GLM_FUNC_QUALIFIER i16vec4 unpackInt4x16(int64 p)
890 {
891 i16vec4 Unpack;
892 memcpy(value_ptr(Unpack), &p, sizeof(Unpack));
893 return Unpack;
894 }
895
896 GLM_FUNC_QUALIFIER uint packUint2x16(u16vec2 const& v)
897 {
898 uint Pack = 0;
899 memcpy(&Pack, &v, sizeof(Pack));
900 return Pack;
901 }
902
903 GLM_FUNC_QUALIFIER u16vec2 unpackUint2x16(uint p)
904 {
905 u16vec2 Unpack;
906 memcpy(value_ptr(Unpack), &p, sizeof(Unpack));
907 return Unpack;
908 }
909
910 GLM_FUNC_QUALIFIER uint64 packUint4x16(u16vec4 const& v)
911 {
912 uint64 Pack = 0;
913 memcpy(&Pack, &v, sizeof(Pack));
914 return Pack;
915 }
916
917 GLM_FUNC_QUALIFIER u16vec4 unpackUint4x16(uint64 p)
918 {
919 u16vec4 Unpack;
920 memcpy(value_ptr(Unpack), &p, sizeof(Unpack));
921 return Unpack;
922 }
923
924 GLM_FUNC_QUALIFIER int64 packInt2x32(i32vec2 const& v)
925 {
926 int64 Pack = 0;
927 memcpy(&Pack, &v, sizeof(Pack));
928 return Pack;
929 }
930
931 GLM_FUNC_QUALIFIER i32vec2 unpackInt2x32(int64 p)
932 {
933 i32vec2 Unpack;
934 memcpy(value_ptr(Unpack), &p, sizeof(Unpack));
935 return Unpack;
936 }
937
938 GLM_FUNC_QUALIFIER uint64 packUint2x32(u32vec2 const& v)
939 {
940 uint64 Pack = 0;
941 memcpy(&Pack, &v, sizeof(Pack));
942 return Pack;
943 }
944
945 GLM_FUNC_QUALIFIER u32vec2 unpackUint2x32(uint64 p)
946 {
947 u32vec2 Unpack;
948 memcpy(value_ptr(Unpack), &p, sizeof(Unpack));
949 return Unpack;
950 }
951}//namespace glm
952