c语言数学库math源码,math.h函数源代码

本文详细介绍了C语言math库中的几个关键函数的源代码实现,包括hypot函数用于计算直角三角形斜边长度,modf函数用于分解浮点数的整数和小数部分,fmod函数计算除法的余数,ldexp函数用于根据指数调整浮点数,frexp函数分解浮点数为尾数和指数,以及双曲函数tanh、sinh和cosh的实现,最后是acos函数的反余弦计算。通过源码分析,读者可以深入了解这些函数的工作原理和精度优化策略。
摘要由CSDN通过智能技术生成

hypot ( )/* -- C语言库函数源代码 - */

/*

hypot函数对于给定的直角三角形的两个直角边,

求其斜边的长度。

*/

//一般的常规算法:

double my_hypot01(double x, double y)

{

double hypotenuse;

x = fabs(x);

y = fabs(y);

IF (x < y)

{

double temp = x;

x = y;

y = temp;

}

if (x == 0.)

return 0.;

else

{

hypotenuse = y/x;

return x*sqrt(1.+hypotenuse*hypotenuse);

}

}

#define __SQRT_DBL_MAX 1.3e+154

#define __SQRT_DBL_MIN 2.3e-162

double my_hypot02(double x, double y)

{

double ratio;

double r, t, s, p, q;

x = fabs(x), y = fabs(y);

if (x < y)

{

double temp = x;

x = y;

y = temp;

}//保持x 是两个直角边中较长的边,y是较短的边。

if (y == 0.)

return x;

/*

主要考虑的是当x很大而y很小,那么它们的平方和将会造成

丢失小数的现象。首先要判断y是否是太小,x是不是太大。

如果出现这种情况则用,第一个公式来处理。其他的则用

这样可以让求出的斜边更加精确一些。

*/

if ((ratio = y / x) > __SQRT_DBL_MIN && x < __SQRT_DBL_MAX)

return x * sqrt(1. + ratio*ratio);

else

{//使用3次迭代是增加精确度。

r = ratio*ratio, p = x, q = y;

do

{

t = 4.+ r;

if (t == 4.)

break;

s = r / t;

p += 2. * s * p;

q *= s;

r = (q / p) * (q / p);

} while (1);

return p;

}

}

struct complex

{

double x;

double y;

}

double cabs(struct complex x)

{

return hypot(z.x,z.y);

}//hypot 函数的封装(这里不再举调用的例子了。)

# define DBL_MAX 1.79769313486231e+308

# define DBL_MIN 2.22507385850721e-308

int main(void)

{

printf("hypot(3, 4)                     =%25.17e\n",hypot(3., 4.));

printf("hypot(3*10^150,4*10^150)        =%25.17g\n",hypot(3.e+150, 4.e+150));

printf("hypot(3*10^306,4*10^306)        =%25.17g\n",hypot(3.e+306, 4.e+306));

printf("hypot(3*10^-320,4*10^-320)      =%25.17g\n",hypot(3.e-320, 4.e-320));

printf("hypot(0.7*DBL_MAX,0.7*DBL_MAX) =%25.17g\n",hypot(0.7*DBL_MAX,0.7*DBL_MAX));

printf("hypot(DBL_MAX, 1.0)             =%25.17g\n",hypot(DBL_MAX, 1.0));

printf("hypot(1.0, DBL_MAX)             =%25.17g\n",hypot(1.0, DBL_MAX));

printf("hypot(0.0, DBL_MAX)             =%25.17g\n",hypot(0.0, DBL_MAX));

printf("\n************************************************************\n");

printf("hypot(3, 4)                     =%25.17e\n",my_hypot01(3., 4.));

printf("hypot(3*10^150,4*10^150)        =%25.17g\n",my_hypot01(3.e+150, 4.e+150));

printf("hypot(3*10^306,4*10^306)      

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值