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