util.h

#pragma once
#include <boost/numeric/ublas/vector.hpp>
#include <boost/numeric/ublas/vector_expression.hpp>

namespace algo { namespace qn {
    namespace ublas = boost::numeric::ublas;
    ublas::matrix<double> 
    initilizeQuasiNewtonHessian(const std::size_t size);
    ublas::matrix<double> 
    initilizeQuasiNewtonInverseHessian(const std::size_t size);
    template <typename E>
    double calculateSumOfSquares(const ublas::vector_expression<E>& x);
    double calculateDistanceNormL2(
        const ublas::vector<double>& x1,
        const ublas::vector<double>& x2);
    double calculateNormL2(const ublas::vector<double>& x);
} } // namespace algo { namespace qn {