point.hh

Go to the documentation of this file.

/*
    @@@@@@@@  @@           @@@@@@   @@@@@@@@ @@
   /@@/////  /@@          @@////@@ @@////// /@@
   /@@       /@@  @@@@@  @@    // /@@       /@@
   /@@@@@@@  /@@ @@///@@/@@       /@@@@@@@@@/@@
   /@@////   /@@/@@@@@@@/@@       ////////@@/@@
   /@@       /@@/@@//// //@@    @@       /@@/@@
   /@@       @@@//@@@@@@ //@@@@@@  @@@@@@@@ /@@
   //       ///  //////   //////  ////////  //

   Copyright (c) 2016 Los Alamos National Laboratory, LLC
   All rights reserved
                                                                              */

#pragma once

#include "flecsi/util/common.hh"
#include "flecsi/util/dimensioned_array.hh"

#include <array>
#include <cmath>

namespace flecsi {
namespace util {

//----------------------------------------------------------------------------//
//----------------------------------------------------------------------------//

template<typename TYPE, size_t DIMENSION>
using point = dimensioned_array<TYPE, DIMENSION, 1>;

//----------------------------------------------------------------------------//
//----------------------------------------------------------------------------//

template<typename TYPE, size_t DIMENSION>
constexpr point<TYPE, DIMENSION> operator*(TYPE const val,
  point<TYPE, DIMENSION> const & p) {
  point<TYPE, DIMENSION> tmp(p);
  for(size_t d(0); d < DIMENSION; ++d) {
    tmp[d] *= val;
  } // for

  return tmp;
} // operator *

//----------------------------------------------------------------------------//
//----------------------------------------------------------------------------//

template<typename TYPE, size_t DIMENSION>
TYPE
distance(point<TYPE, DIMENSION> const & a, point<TYPE, DIMENSION> const & b) {
  TYPE sum(0);
  for(size_t d(0); d < DIMENSION; ++d) {
    sum += (square)(a[d] - b[d]);
  } // for

  return std::sqrt(sum);
} // distance

//----------------------------------------------------------------------------//
//----------------------------------------------------------------------------//

template<typename TYPE, size_t DIMENSION>
constexpr point<TYPE, DIMENSION>
midpoint(point<TYPE, DIMENSION> const & a, point<TYPE, DIMENSION> const & b) {
  return point<TYPE, DIMENSION>((a + b) / 2.0);
} // midpoint

//----------------------------------------------------------------------------//
//----------------------------------------------------------------------------//

template<template<typename...> class CONTAINER, typename TYPE, size_t DIMENSION>
constexpr auto
centroid(CONTAINER<point<TYPE, DIMENSION>> const & points) {
  point<TYPE, DIMENSION> tmp(0.0);

  for(auto p : points) {
    tmp += p;
  } // for

  tmp /= points.size();

  return tmp;
} // centroid

//----------------------------------------------------------------------------//
//----------------------------------------------------------------------------//

template<typename TYPE, size_t DIMENSION>
constexpr auto
centroid(std::initializer_list<point<TYPE, DIMENSION>> points) {
  point<TYPE, DIMENSION> tmp(0.0);

  for(auto p : points) {
    tmp += p;
  } // for

  tmp /= points.size();

  return tmp;
} // centroid

} // namespace util
} // namespace flecsi