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 29namespace clover { 30 /// 31 /// N-dimensional coordinate array. 32 /// 33 template<typename T, int N> 34 class point { 35 public: 36 point() : a() { 37 } 38 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 44 point(const T *v) { 45 std::copy(v, v + N, a.begin()); 46 } 47 48 T &operator[](int i) { 49 return a[i]; 50 } 51 52 const T &operator[](int i) const { 53 return a[i]; 54 } 55 56 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 63 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