1 /* Copyright Jukka Jyl�nki
2
3    Licensed under the Apache License, Version 2.0 (the "License");
4    you may not use this file except in compliance with the License.
5    You may obtain a copy of the License at
6
7        http://www.apache.org/licenses/LICENSE-2.0
8
9    Unless required by applicable law or agreed to in writing, software
10    distributed under the License is distributed on an "AS IS" BASIS,
11    WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
12    See the License for the specific language governing permissions and
13    limitations under the License. */
14
15 /** @file SSEMath.h
16         @author Jukka Jyl�nki
17         @brief SIMD-enabled math helper utilities. */
18 #pragma once
19
20 #include "../MathBuildConfig.h"
21 #include "MathNamespace.h"
22 #include "../MathGeoLibFwd.h"
23 #include <stdint.h>
24 #include <cstddef>
25 #include "Reinterpret.h"
26 #include "simd.h"
27 #include <memory>
28 #include <stdexcept>
29
30 MATH_BEGIN_NAMESPACE
31
32 /// Allocates the given amount of memory at the given alignment.
33 void *AlignedMalloc(size_t size, size_t alignment);
34
35 /// \todo This is not really placement-new.
36 template<typename T>
37 inline T *AlignedNew(size_t numElements, size_t alignment) { return reinterpret_cast<T*>(AlignedMalloc(numElements*sizeof(T), alignment)); }
38
39 /// \todo This is not really placement-new.
40 template<typename T>
41 inline T *AlignedNew(size_t numElements) { return AlignedNew<T>(numElements, 16); }
42
43 /// Frees memory allocated by AlignedMalloc.
44 void AlignedFree(void *ptr);
45
46 template<class T, size_t Alignment>
47 struct AlignedAllocator : std::allocator<T>
48 {
49         template<class U>
50         struct rebind { typedef AlignedAllocator<U, Alignment> other; };
51
52         typedef std::allocator<T> base;
53
54         typedef typename base::const_pointer const_pointer;
55         typedef typename base::pointer pointer;
56         typedef typename base::reference reference;
57         typedef typename base::const_reference const_reference;
58         typedef typename base::value_type value_type;
59         typedef typename base::size_type size_type;
60         typedef typename base::difference_type difference_type;
61
62         AlignedAllocator(){}
63         AlignedAllocator(const AlignedAllocator &):std::allocator<T>(){}
64         template<typename U>
65         AlignedAllocator(const AlignedAllocator<U, Alignment> &){}
66
67         pointer allocate(size_type n, const void * = 0)
68         {
69                 return (pointer)AlignedMalloc(n*sizeof(T), Alignment);
70         }
71
72         void deallocate(pointer ptr, size_type)
73         {
74                 AlignedFree(ptr);
75         }
76 };
77
78 MATH_END_NAMESPACE