wine中complex::pow 补丁的实现

来源：互联网发布：复杂网络算法与应用pdf 编辑：程序博客网时间：2024/06/05 18:32

我发了一个wine的补丁，实现complex::pow的部分重载：http://source.winehq.org/patches/data/104267

HongQian帮我在irc上问了Piotr对这个补丁的看法，他说：
“it generally looks ok but it's possible to implement it with better
precision, so I'm not sure what to think about it”

“我想了想发现，其实我们在实现complex::pow
在指数为实整数的情况下的重载的时候，没有必要先把实整数转换为复数的。一开始我是为了简化补丁的难度，就建议参考最像的函数进行实现，没有充分考虑到精度的问题。”

我去查了一些资料，发现，以前初学C\C++用pow函数的时候有点疑惑，为什么参数int要转化成double，以为是系统规定的，当时为了编译能过就用了转换了，也没有看源代码的想法。今天才知道，原来是考虑到那么几个问题（overflow、负数次幂的结果）等才用double的，而且还提高了复用，感觉设计得很好~（除了给人一点点confuse）

然后我做了测试，这是在win下的结果：

pow(complex_float= (3.14159f, 0), 3): (31.006203,0.000000)

complex_translation: (31.006201,0.000000)

*******************************************************

pow(complex_float= (32.1234f, 0), 6): (1098826752.000000,0.000000)

complex_translation: (1098826624.000000,0.000000)

*******************************************************

pow(complex_float= (7.12f, 0), 6): (130280.640625,0.000000)

complex_translation: (130280.648438,0.000000)

*******************Now is complex_double************************************

pow(complex_double= (3.14159/2, 0), 3): (3.875775, 0.000000)

double_complex_translation: (3.875775, 0.000000)

*******************************************************

pow(complex_double= (32.12345678, 0), 13): (38787300223012692000.000000, 0.000000)

double_complex_translation: (38787300223012684000.000000, 0.000000)

*******************************************************

pow(complex_double= (7.12, 0), 33): (13546895790785450000000000000.000000, 0.000000)

double_complex_translation: (13546895790785442000000000000.000000, 0.000000)

*******************************************************

二分求幂现在决定用用手写的底数为复数而指数为整数次幂的pow来解决，经zhenbo提醒，发现用二分求幂就可以高效率地实现循环相乘求值，于是使用实现之：

+complex_double* __cdecl complex_double_pow_ci(complex_double *ret, const complex_double *l, int r)+{+    complex_double c = *l;+    complex_double temp = { 1.0, 0 };+    complex_double_assign(ret, &temp);++    if(r < 0) {+        r = -r;+        complex_double_div(&c, &temp, l);+    }++    for(; r>0; r>>=1) {+        if(r & 1)+            complex_double_mult_assign(ret, &c);+        complex_double_assign(&temp, &c);+        complex_double_mult_assign(&c, &temp);+    }++    return ret;+}

0 0