码迷,mamicode.com
首页 > 其他好文 > 详细

软件快速计算平方根与平方根的倒数

时间:2014-11-15 18:55:14      阅读:143      评论:0      收藏:0      [点我收藏+]

标签:算法   迭代   

#include<iostream>
#include<cmath>
using namespace std;



float Sqrt(float x);

float InvSqrt(float x);

int main(void)
{       
    system("color F0"); 
    cout.setf(ios::fixed);                     
    cout.setf(ios::showpoint);
    cout.precision(10);  


    //快速算法
    double result1,result2;
    result1=Sqrt(635);
    //系统方法
    double a=635;
    result2=sqrt(a);


    cout<<result1<<endl;
    cout<<result2<<endl;


    system("pause");
    return 0;
}




//浮点数平方根快速算法
float Sqrt(float x)
{
float xhalf = 0.5f*x;
int i = *(int*)&x;      // get bits for floating value
i = 0x5f375a86- (i>>1); // gives initial guess y0                     迅速找出第一次比较接近的迭代结果
x = *(float*)&i;        // convert bits BACK to float
x = x*(1.5f-xhalf*x*x); // Newton step, repeating increases accuracy  牛顿迭代法,次数越多越准确
x = x*(1.5f-xhalf*x*x); // Newton step, repeating increases accuracy
x = x*(1.5f-xhalf*x*x); // Newton step, repeating increases accuracy


return 1/x;

}

//实现平方根的倒数
float InvSqrt(float x)
{
float xhalf = 0.5f*x;
int i = *(int*)&x;      // get bits for floating value
i = 0x5f375a86- (i>>1); // gives initial guess y0                     迅速找出第一次比较接近的迭代结果
x = *(float*)&i;        // convert bits BACK to float
x = x*(1.5f-xhalf*x*x); // Newton step, repeating increases accuracy  牛顿迭代法,次数越多越准确
x = x*(1.5f-xhalf*x*x); // Newton step, repeating increases accuracy
x = x*(1.5f-xhalf*x*x); // Newton step, repeating increases accuracy


return x;

}

软件快速计算平方根与平方根的倒数

标签:算法   迭代   

原文地址:http://blog.csdn.net/u013467442/article/details/41147815

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!