/*****************************************************************
 * Copyright (c) 2020 by Douglas Wilhelm Harder.
 * All rights reserved.
 *
 * To use, please contact dwharder@uwaterloo.ca
 *****************************************************************/

#pragma once

#include <utility>

template <typename T>
T lin_y0(    T *y, std::size_t n, unsigned int m, bool forwards = false );

template <typename T>
T lin_y1(    T *y, std::size_t n, unsigned int m, bool forwards = false ); 

template <typename T>
T lin_dy0(   T *y, std::size_t n, unsigned int m, bool forwards = false );

template <typename T>
T lin_iy0(   T *y, std::size_t n, unsigned int m, bool forwards = false );

template <typename T>
T lin_jit(   T *y, std::size_t n, unsigned int m, T epsilon, bool forwards = false );

template <typename T>
T lin_root(  T *y, std::size_t n, unsigned int m, bool forwards = false );

template <typename T>
T quad_y0(   T *y, std::size_t n, unsigned int m, bool forwards = false );

template <typename T>
T quad_y1(   T *y, std::size_t n, unsigned int m, bool forwards = false );

template <typename T>
T quad_dy0(  T *y, std::size_t n, unsigned int m, bool forwards = false );

template <typename T>
T quad_ddy0( T *y, std::size_t n, unsigned int m, bool forwards = false );

template <typename T>
T quad_iy0(  T *y, std::size_t n, unsigned int m, bool forwards = false );

template <typename T>
T quad_jit(  T *y, std::size_t n, unsigned int m, T epsilon, bool forwards = false );

template <typename T>
T quad_ext(  T *y, std::size_t n, unsigned int m, bool forwards = false );

template <typename T>
std::pair<T, T> quad_roots( T *y, std::size_t n, unsigned int m, bool forwards = false );

#include "least_sqrs_est.tpp"