-
-template<typename T, unsigned N>
-inline LinAl::Vector<T, N+1> augment_vector(const LinAl::Vector<T, N> &v, T s)
-{
- LinAl::Vector<T, N+1> r;
- for(unsigned i=0; i<N; ++i)
- r[i] = v[i];
- r[N] = s;
- return r;
-}
-
-template<typename T, unsigned N>
-inline LinAl::Vector<T, N-1> reduce_vector(const LinAl::Vector<T, N> &v)
-{
- LinAl::Vector<T, N-1> r;
- for(unsigned i=0; i<N-1; ++i)
- r[i] = v[i];
- return r;
-}
-
-template<typename T, unsigned N>
-inline LinAl::Vector<T, N-1> divide_vector(const LinAl::Vector<T, N> &v)
-{
- LinAl::Vector<T, N-1> r;
- for(unsigned i=0; i<N-1; ++i)
- r[i] = v[i]/v[N-1];
- return r;
-}
-
-