Zeros.hpp

浏览该文件的文档.

 
 
#pragma once
 
#include "AstGlobal.h"
#include "zeros.h"
 
AST_NAMESPACE_BEGIN
 
namespace details{
    template<typename Func>
    A_ALWAYS_INLINE double callback_func(double x, void* func)
    {
        return (*(Func*)func)(x);
    }
}
 
template<typename Func>
double bisect(Func f, double xa, double xb, double xtol, double rtol, int iter, scipy_zeros_info &solver_stats)
{
    return ::bisect(details::callback_func<Func>, xa, xb, xtol, rtol, iter, &f, &solver_stats);
}
 
 
template<typename Func>
double ridder(Func f, double xa, double xb, double xtol, double rtol, int iter, scipy_zeros_info &solver_stats)
{
    return ::ridder(details::callback_func<Func>, xa, xb, xtol, rtol, iter, &f, &solver_stats);
}
 
 
template<typename Func>
double brenth(Func f, double xa, double xb, double xtol, double rtol, int iter, scipy_zeros_info &solver_stats)
{
    return ::brenth(details::callback_func<Func>, xa, xb, xtol, rtol, iter, &f, &solver_stats);
}
 
 
template<typename Func>
double brentq(Func f, double xa, double xb, double xtol, double rtol, int iter, scipy_zeros_info &solver_stats)
{
    return ::brentq(details::callback_func<Func>, xa, xb, xtol, rtol, iter, &f, &solver_stats);
}
 
 
template<typename Func>
double secant(Func f, double xa, double xb, double xtol, double rtol, int iter, scipy_zeros_info &solver_stats)
{
    return ::secant(details::callback_func<Func>, xa, xb, xtol, rtol, iter, &f, &solver_stats);
}
 
 
AST_NAMESPACE_END