文章/答案/技术大牛

发布

社区首页 >问答首页 >C预处理程序计算常量的sin()

问C预处理程序计算常量的sin()
EN

Stack Overflow用户

提问于 2012-03-20 21:53:44

回答 4查看 1.6K关注 0票数 6

有没有办法说服C预处理器在编译时评估一个常数的先验函数？

例如，将(int)256*sin(PI/4)替换为181。这将帮助我将魔法数字从我的代码中删除。

如果有差异，我使用的是MSPGCC4.5.3，运行时没有可用的sin()或cos()。

c-preprocessor

trigonometry

回答 4

Stack Overflow用户

回答已采纳

发布于 2012-03-21 00:17:44

只要您的参数都在范围内-π/4，+π/4，您就可以使用libm使用的相同公式标准实现来计算sin。它直到最后一个地方都是正确的(最多是1ulp错误)，就像IEEE标准所要求的那样：

static const double 
half =  5.00000000000000000000e-01, /* 0x3FE00000, 0x00000000 */
S1  = -1.66666666666666324348e-01, /* 0xBFC55555, 0x55555549 */
S2  =  8.33333333332248946124e-03, /* 0x3F811111, 0x1110F8A6 */
S3  = -1.98412698298579493134e-04, /* 0xBF2A01A0, 0x19C161D5 */
S4  =  2.75573137070700676789e-06, /* 0x3EC71DE3, 0x57B1FE7D */
S5  = -2.50507602534068634195e-08, /* 0xBE5AE5E6, 0x8A2B9CEB */
S6  =  1.58969099521155010221e-10; /* 0x3DE5D93A, 0x5ACFD57C */

double __kernel_sin(double x, double y, int iy)
{
    double z,r,v;
    int ix;
    ix = __HI(x)&0x7fffffff;    /* high word of x */
    if(ix<0x3e400000)           /* |x| < 2**-27 */
       {if((int)x==0) return x;}        /* generate inexact */
    z   =  x*x;
    v   =  z*x;
    r   =  S2+z*(S3+z*(S4+z*(S5+z*S6)));
    if(iy==0) return x+v*(S1+z*r);
    else      return x-((z*(half*y-v*r)-y)-v*S1);
}

来源：http://www.netlib.org/fdlibm/k_sin.c

虽然不是我所说的“愉快”，但您肯定可以将整个函数转换为一个宏，该宏将计算为(编译时)浮点常量表达式。(忽略一开始与值无关的位黑客行为，据我所知，您应该假定iy为0。)

票数 3

Stack Overflow用户

发布于 2012-03-20 21:58:30

C预处理程序不能提供sin()或cos()。

对于我的应用程序，我使用perl脚本创建一个单独的.h文件，其中包含所需的预计算。也许有更好的方法来做它，但这集成到我的工作流程足够好。

票数 4

Stack Overflow用户

发布于 2012-03-20 22:03:06

预处理器只能解析宏，这与执行函数完全不同。我能想到的最接近减少魔法数字的解决方案是创建一个标题，其中包含最常用的sins或共罪值：

#define SIN_PI  (-1)
#define SIN_PI2 0
#define SIN_PI4 0.707106781186548
...

然后你可以写：

256*SIN_PI2

并让编译器优化将其简化为单个常量。

票数 4

页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持

原文链接：

https://stackoverflow.com/questions/9795586

复制

相似问题

问C预处理程序计算常量的sin()
EN

回答 4

Stack Overflow用户

Stack Overflow用户

Stack Overflow用户

社区

活动

圈层

关于

腾讯云开发者

热门产品

热门推荐

更多推荐

问C预处理程序计算常量的sin()EN

回答 4

Stack Overflow用户

Stack Overflow用户

Stack Overflow用户

社区

活动

圈层

关于

腾讯云开发者

热门产品

热门推荐

更多推荐

问C预处理程序计算常量的sin()
EN