summaryrefslogtreecommitdiff
path: root/modules/codec/c++/math.hpp
diff options
context:
space:
mode:
Diffstat (limited to 'modules/codec/c++/math.hpp')
-rw-r--r--modules/codec/c++/math.hpp18
1 files changed, 18 insertions, 0 deletions
diff --git a/modules/codec/c++/math.hpp b/modules/codec/c++/math.hpp
index ddeea3f..67d8ed3 100644
--- a/modules/codec/c++/math.hpp
+++ b/modules/codec/c++/math.hpp
@@ -21,5 +21,23 @@ data<schema::Scalar<T>, Encoding> dot(const data<schema::Vector<T,D>, Encoding>&
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;
+}
}
}