锘??xml version="1.0" encoding="utf-8" standalone="yes"?>青青草原综合久久大伊人精品,国产精品免费福利久久,亚洲成色www久久网站夜月http://www.shnenglu.com/newplan/archive/2008/05/07/49127.html灞辨硥寮歡灞辨硥寮歡Wed, 07 May 2008 07:05:00 GMThttp://www.shnenglu.com/newplan/archive/2008/05/07/49127.htmlhttp://www.shnenglu.com/newplan/comments/49127.htmlhttp://www.shnenglu.com/newplan/archive/2008/05/07/49127.html#Feedback0http://www.shnenglu.com/newplan/comments/commentRss/49127.htmlhttp://www.shnenglu.com/newplan/services/trackbacks/49127.html
/*
Romberg Algorithm
寮鍙戣?newplan
寮鍙戞棩鏈?08.05.07
*/


/*=======================================*/
/*INCLUDES*/
#include <cstdlib>
#include <iostream>
#include <cmath>
/*=======================================*/
/*MACROS USED IN THIS FILE*/
#define MAX 20
#define PRECISION 0.000008
/*=======================================*/
/*DECLARE NAMES IN STL NAMESAPCE */
using std::cout;
using std::endl;

/*=======================================*/
/*CLASS FUNC (FUNCTION OBJECT): THE ORIGINAL FUNCTION WE WANT TO INTEGRAL*/
class func{
public:
func(double x=1.0):exp(x){}
double operator()(const double& dnum)const{return pow(dnum,exp);}
private:
double exp;
};
/*=======================================*/
/*CLASS ECHELONFUNC (FUNCTION OBJECT)姊艦娉曠殑閫掓帹鍏紡*/
class echelonFunc{
public:
echelonFunc(double begining,double ending,func & myfunc);
double operator()();
private:
double h;
int n;
double T;
func myfunc;
double begining;
double ending ;
};
/*=======================================*/
echelonFunc::echelonFunc(double begining,double ending,func & myfunc)
{
this->begining=begining;
this->ending=ending;
this->h=ending-begining;
this->n=0;
this->T=0;
this->myfunc=myfunc;//FUCNTION
}
/*------------------------------*/
/* INCREASE FUNCTION 閫掓帹鍑芥暟*/
double echelonFunc::operator()()
{ if(this->n==0)
{
this->T=h*0.5*(myfunc(this->begining)+myfunc(this->ending));
this->n=1;
return this->T;
}
double len=0.5*h;
double sum=0;
int k=0;
for(k=0;k<this->n;k++)
{
sum+=myfunc(len);
len=len+h;
}
this->T=0.5*this->T+0.5*h*sum;
this->h/=2;
this->n*=2;
return this->T;
}
/*=======================================*/
/*THE MAIN CLASS IN THIS PROGRAM*/
class Romberg{
public:
Romberg(double begining,double ending,double exp);
~Romberg();
private:
void RombergCPU();/*THE MOST IMPORTANT FUNCTION IN THIS PROGRAM*/
echelonFunc *echol;
double T[MAX][MAX];/*STO THE ROMBERG TABLE*/
};
/*------------------------------*/
Romberg::Romberg(double begining ,double ending ,double paraexp)
{
func myfunc(paraexp);
echol = new echelonFunc(begining,ending,myfunc);
RombergCPU();
}
/*------------------------------*/
Romberg::~Romberg()
{
delete echol;
}
/*------------------------------*/
void Romberg::RombergCPU()
{ clock_t Start; //TIME STARAT
clock_t End; //TIME END
double *p[MAX];//WE USE THIS POINTER ARRAY TO ACCELERATE ALGOTITHM
double **q;
int i,j;
Start = clock();//TIME START FROM HERE
for(i = 0,q = p; q < p+MAX; q++,i++)
*q= &T[i][0];
double a,b,pows;
cout<<"-----------------------Romberg Algorithm---------------------"<<endl;
*p[0]=(*echol)();
cout<<" "<<*p[0]<<endl;
p[0]++;
do{
*p[0]=(*echol)();
cout<<" "<<*p[0];
p[0]++;
for(i=1;;i++)
{
pows=pow(4.0,double(i));
a=pows/(pows-1);
b=1/(pows-1);
*p[i]=a*(*(p[i-1]-1))-b*(*(p[i-1]-2));//ROMBERG ALGORITHM
cout<<" "<<*p[i];
if(p[i]==&T[i][0])
{
p[i]++;
break;
}
p[i]++;
}
cout<<endl;//fabs(T[i][0]-T[i-1][0])
}while(fabs(T[i][0]-T[i-1][0])>PRECISION);
End = clock();//TIME END HERE
cout<<"-------------------------------------------------------------"<<
endl<<" TIME SPEND:"<<(double)(End-Start)/CLOCKS_PER_SEC<<endl;
}
/*=======================================*/
/*MAIN FUNCTION*/
int main(int argc, char *argv[])
{
Romberg romberg(0,1,1.5);//ROMBERG API :BEGIN(0) END(1) EXP(1.5)
system("PAUSE");
return EXIT_SUCCESS;
}

]]>
国产三级观看久久|
亚洲国产日韩欧美久久|
久久夜色精品国产欧美乱|
国产69精品久久久久久人妻精品|
久久精品国产免费观看三人同眠|
久久综合噜噜激激的五月天|
欧美精品一本久久男人的天堂|
久久久久国产一级毛片高清板
|
久久精品国产一区|
久久久国产精华液|
久久综合给合久久狠狠狠97色|
99久久精品免费看国产免费|
国产精品久久婷婷六月丁香|
精品久久久久久中文字幕|
97视频久久久|
久久er国产精品免费观看8|
亚洲av日韩精品久久久久久a|
国产精品无码久久四虎|
久久久久亚洲AV无码专区体验|
欧美日韩中文字幕久久久不卡|
av无码久久久久久不卡网站|
精品综合久久久久久98|
久久亚洲中文字幕精品一区四|
国产精品一区二区久久国产|
久久伊人五月丁香狠狠色|
久久久久亚洲AV成人网人人网站|
91精品国产色综合久久|
久久香蕉超碰97国产精品|
久久久久久久综合狠狠综合|
久久强奷乱码老熟女|
亚洲午夜精品久久久久久人妖|
久久久久久亚洲精品成人|
久久人人爽人人爽人人爽|
亚洲欧美日韩久久精品|
色综合久久久久综合99|
久久久久亚洲AV无码专区网站|
亚洲欧美日韩精品久久|
品成人欧美大片久久国产欧美...|
久久96国产精品久久久|
国产成人精品久久综合|
久久国产福利免费|