1 // 2 // Copyright 2012 Francisco Jerez 3 // 4 // Permission is hereby granted, free of charge, to any person obtaining a 5 // copy of this software and associated documentation files (the "Software"), 6 // to deal in the Software without restriction, including without limitation 7 // the rights to use, copy, modify, merge, publish, distribute, sublicense, 8 // and/or sell copies of the Software, and to permit persons to whom the 9 // Software is furnished to do so, subject to the following conditions: 10 // 11 // The above copyright notice and this permission notice shall be included in 12 // all copies or substantial portions of the Software. 13 // 14 // THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR 15 // IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, 16 // FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL 17 // THE AUTHORS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, 18 // WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF 19 // OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE 20 // SOFTWARE. 21 // 22 23 #ifndef __CORE_GEOMETRY_HPP__ 24 #define __CORE_GEOMETRY_HPP__ 25 26 #include <array> 27 #include <algorithm> 28 29 namespace clover { 30 /// 31 /// N-dimensional coordinate array. 32 /// 33 template<typename T, int N> 34 class point { 35 public: point()36 point() : a() { 37 } 38 point(std::initializer_list<T> v)39 point(std::initializer_list<T> v) { 40 auto it = std::copy(v.begin(), v.end(), a.begin()); 41 std::fill(it, a.end(), 0); 42 } 43 point(const T * v)44 point(const T *v) { 45 std::copy(v, v + N, a.begin()); 46 } 47 operator [](int i)48 T &operator[](int i) { 49 return a[i]; 50 } 51 operator [](int i) const52 const T &operator[](int i) const { 53 return a[i]; 54 } 55 operator +(const point & p) const56 point operator+(const point &p) const { 57 point q; 58 std::transform(a.begin(), a.end(), p.a.begin(), 59 q.a.begin(), std::plus<T>()); 60 return q; 61 } 62 operator ()(const point & p) const63 T operator()(const point &p) const { 64 return std::inner_product(p.a.begin(), p.a.end(), a.begin(), 0); 65 } 66 67 protected: 68 std::array<T, N> a; 69 }; 70 } 71 72 #endif 73