使用虚拟容器的方法实现二分法求解方程的根(2)

#include <iostream>
#include <iterator>
#include <algorithm>
#include <functional>
#include <cmath>
#include "dummy.h"//在这个文件中实现了虚拟容器模板
//下面是求一元方程根的二分区间法
template< class Equation,long number=15 > struct Method
{
        typedef Equation equation_type;
        typedef std::pair<double,double> element_type;
        struct Generator:std::unary_function<const size_t&,element_type&>
        {
                element_type&operator()(const size_t&i)
                {
                        Method::process(Method::value());
                        return Method::value();
                }
        };
        struct Predicate:std::unary_function<const element_type&,bool>
        {
                bool operator()(const element_type&e)
                {
                        const double&eps = pow(0.1,number);
                        return fabs( e.first - e.second ) < eps;
                }
        };
public:
        //这个静态函数用来设置所求方程的根所在的初始区间
        static void range(const double&a,const double&b)
        {
                value() = std::make_pair(a,b);
        }
private:
        //下面的这个静态函数仅仅是为了提高效率而提供的静态分配的一个值变量
        static element_type&value()
        {
                static element_type _value;
                return _value;
        }
        //这个静态函数就是二分法求方程根的过程代码
        static void process(element_type&e)
        {
                Equation eqn;
                const double&eps = pow(0.1,number);
                double&a=e.first;
                double&b=e.second;
                double x,fa,fb,fmid;
                fa=eqn(a); fb=eqn(b);
                if(fabs(fa)<=eps) b = a;
                else if(fabs(fb)<=eps) a = b;
                else{
                        fmid=eqn(x=(a+b)/2);
                        if(fmid==0){ a=x,b=x; }
                        else if(fmid*fa<0){ b=(a+b)/2; }
                        else if(fmid*fb<0){ a=(a+b)/2; }
                }
        }
};
//方程为“x*x-2 = 0”很明显该方程有两个根分别为：+sqrt(2)和-sqrt(2)
struct Equation:std::unary_function<const double&,double>
{
        double operator()(const double&x)
        {
                return x*x - 2;
        }
};
//从上面的两个结构体，可以看出方法和方程已经完全的隔离开了，这样的代码组织是为了能够
//更加方便和快速的书写出数值算法。很显然如果我们将数值方法里面的其它的很多求解方程的
//根的方法都以上面的Method模板那样以库的形式给出，那么用户就可以很容易的以很多不同的
//方法来解同一个方程。更重要的是，这里采用了和STL兼容的方式实现代码，因此如果你熟悉
//STL的话，就可以很容易的使用这些数值方法。这样就可以减少使用该数值库的学习时间。
//当然到目前为止，这还是一个设想，但是随着时间的推移，我会慢慢的实现所有常用的算法的
//，当然是以STL兼容的方式提供，尽量减少重复代码。呵呵，忘了这里有一个假设就是你一定
//要会使用STL并熟悉STL的相关概念哦:)
int main()
{
        std::cout.precision(20);
        typedef Method<Equation,15> M;
        typedef dummy<M::element_type,M::Generator> SOLVE;
        {
                M::range(1.0,2.0);//设置根所在的区间
                SOLVE v(1,10000);
                SOLVE::iterator root = std::find_if(v.begin(),v.end(),M::Predicate());
                std::cout << "迭代次数：[" << static_cast<const size_t&>(root) << "]" << std::endl;
                std::cout << "左根：[" << (*root).first << "]" << std::endl;
                std::cout << "右根：[" << (*root).second << "]" << std::endl;
        }
        {
                M::range(-1.0,-2.0);//设置根所在的区间
                SOLVE v(1,10000);
                SOLVE::iterator root = std::find_if(v.begin(),v.end(),M::Predicate());
                std::cout << "迭代次数：[" << static_cast<const size_t&>(root) << "]" << std::endl;
                std::cout << "左根：[" << (*root).first << "]" << std::endl;
                std::cout << "右根：[" << (*root).second << "]" << std::endl;
        }
        return 0;
}

//该程序的输出如下：
/*******************************************************************************
迭代次数：[50]
左根：[1.4142135623730949]
右根：[1.4142135623730949]
迭代次数：[50]
左根：[-1.4142135623730949]
右根：[-1.4142135623730949]
*******************************************************************************/