RootFind1dSecant.h

// 1-D Secant Root finding class
#pragma once
 
#include <functional>
#include <tuple>
 
#include "detail/RootFind1dBase.h"
 
namespace isce3 { namespace math {

class RootFind1dSecant : public detail::RootFind1dBase {
    // aliases
protected:
    using func_t = std::function<double(double)>;
    using tuple4_t = std::tuple<double, double, bool, int>;
 
public:
    // inherit all constructors from base class
    using detail::RootFind1dBase::RootFind1dBase;

    std::tuple<double, double, bool, int> root(
            const std::function<double(double)>& f, double x0, double x1) const;
 
private:
    tuple4_t _secant(const func_t& f, double x0, double x1) const;

    tuple4_t _secant_strict(const func_t& f, double x0, double x1) const;
};
 
}} // namespace isce3::math