#pragma once

#include "./core/omp_config.h"
#include "detail/abs_serial.h"


namespace numerics{

	// ---------------- Elementwise ----------------
	template <typename T>
	inline void inplace_abs(T& c) {
		detail::inplace_abs_scalar_serial(c);
	}
	template <typename T>
	inline T abs(const T c) {
		T out = c;
		inplace_abs(out);
		return out;
	}

	template <typename T>
	inline void inplace_abs(utils::Matrix<T>& A) {
		detail::inplace_abs_elementwise_serial(A);
	}

	template <typename T>
	inline utils::Matrix<T> abs(const utils::Matrix<T>& A) {
	    utils::Matrix<T> out = A;
	    inplace_abs(out);
	    return out;
	}

	template <typename T>
	inline void inplace_abs(utils::Vector<T>& v) {
		detail::inplace_abs_elementwise_serial(v);
	}

	template <typename T>
	inline utils::Vector<T> abs(const utils::Vector<T>& v) {
	    utils::Vector<T> out = v;
	    inplace_abs(out);
	    return out;
	}


}