#pragma once

#include "data_math.hpp"

namespace saw {
namespace math {
/*
template<typename T, typename Encoding = encode::Native>
data<typename T::InnerType,Encoding> norm_2(const data<T,Encoding>& d){
	return {};
}
*/

template<typename T, uint64_t D, typename Encoding = encode::Native>
data<schema::Scalar<T>, Encoding> dot(const data<schema::Vector<T,D>, Encoding>& left, const data<schema::Vector<T,D>, Encoding>& right){
	data<schema::Scalar<T>,Encoding> val;
	auto& inner = val({});
	for(uint64_t i = 0u; i < D; ++i){
		inner = inner + left({{i}}) * right({{i}});
	}

	return val;
}

template<typename T,typename Encoding = encode::Native>
data<schema::Scalar<T>,Encoding> sqrt(const data<schema::Scalar<T>,Encoding>& inp){
	data<schema::Scalar<T>,Encoding> out;
	out.at({}).set(std::sqrt(inp.at({}).get()));
	return out;
}

template<typename T, uint64_t D, typename Encoding = encode::Native>
data<schema::Vector<T,D>, Encoding> normalize(const data<schema::Vector<T,D>>& input ) {
	auto inp_dot = dot<T,D,Encoding>(input,input);
	auto sqrt_inp_dot = sqrt<T>(inp_dot);

	saw::data<schema::Vector<T,D>, Encoding> out;
	out.at({{0u}}).set(out.at({{0u}}).get() / sqrt_inp_dot.at({}).get());
	out.at({{1u}}).set(out.at({{1u}}).get() / sqrt_inp_dot.at({}).get());
	return out;
}
}
}