arithmetic.hpp

Go to the documentation of this file.

/*************************************************************************
 *
 *  Project
 *                         _____ _____  __  __ _____
 *                        / ____|  __ \|  \/  |  __ \
 *  ___  _ __   ___ _ __ | |  __| |__) | \  / | |__) |
 * / _ \| '_ \ / _ \ '_ \| | |_ |  ___/| |\/| |  ___/
 *| (_) | |_) |  __/ | | | |__| | |    | |  | | |
 * \___/| .__/ \___|_| |_|\_____|_|    |_|  |_|_|
 *      | |
 *      |_|
 *
 * Copyright (C) Akiel Aries, <akiel@akiel.org>, et al.
 *
 * This software is licensed as described in the file LICENSE, which
 * you should have received as part of this distribution. The terms
 * among other details are referenced in the official documentation
 * seen here : https://akielaries.github.io/openGPMP/ along with
 * important files seen in this project.
 *
 * You may opt to use, copy, modify, merge, publish, distribute
 * and/or sell copies of the Software, and permit persons to whom
 * the Software is furnished to do so, under the terms of the
 * LICENSE file. As this is an Open Source effort, all implementations
 * must be of the same methodology.
 *
 *
 *
 * This software is distributed on an AS IS basis, WITHOUT
 * WARRANTY OF ANY KIND, either express or implied.
 *
 ************************************************************************/
 
#ifndef ARITHMETIC_HPP
#define ARITHMETIC_HPP
#include <cstdint>
 
namespace gpmp {
 
template <typename T> class Arith {
  public:
    T arr_add(T arr[], int64_t n) {
        T sum = 0;
 
        for (int64_t i = 0; i < n; i++) {
            sum += arr[i];
        }
        return sum;
    }
 
    T arr_sub(T arr[], int64_t n) {
        T diff = arr[0];
 
        for (int64_t i = 1; i < n; i++) {
            diff -= arr[i];
        }
        return diff;
    }
 
    T arr_mlt(T arr[], int64_t n) {
        T prod = arr[0];
 
        for (int64_t i = 0; i < n; i++) {
            prod *= arr[i];
        }
        return prod;
    }
};
 
class Basics {
  public:
    int64_t greatest_power(int64_t n, int64_t p);
 
    int64_t op_gcd(int64_t x, int64_t y);
 
    // TODO floating point operations for square root, logarithms,
    // exponentation, trigonometric functions of floating points
 
    // complex arithmetic, arithmetic on complex numbers
 
    // precision arithmetic
 
    // TODO FIXME XXX BUG CHECKME TESTME PENDING DOCME
};
 
} // namespace gpmp
 
//  public:
//    template <typename T> inline T add(T t) {
//        return t;
//    }
 
//    template <typename T, typename... Ts> inline auto add(T t, Ts... ts)
//    {
//        return t + add(ts...);
//    }
 
//    template <typename X> inline X sub(X x) {
//        return x;
//    }
//    template <typename X, typename... Xy> inline auto sub(X x, Xy... xy)
//    {
//        return x - sub(xy...);
//    }
 
//    template <typename W> inline W mult(W w) {
//        return w;
//    }
//    template <typename W, typename... Wv> inline auto mult(W w, Wv... wv)
//    {
//        return w * mult(wv...);
//    }
 
//    template <typename Z> inline Z exp(Z z) {
//        return z;
//    }
//    template <typename Z, typename... Zy> inline auto exp(Z z, Zy... zy)
//    {
//        return z *= exp(zy...);
//    }
//};
 
#endif