C近似函数

Question

我正在为java中的C，exp，log和sqrt做一些近似函数。我对指针如何工作有些生疏 - 这个语法是否正确？

#include <math.h> 
#include "QDMath.h" 

JNIEXPORT jdouble JNICALL Java_QDMath_exp 
    (JNIEnv *env, jclass class, jdouble val) 
{ 
    jint tmp = (jint) (1512775 * val + 1072632447); 
    jdouble p = 0.0; 
    *(1 + (jint *) &p) = tmp; 
    return p; 
} 

JNIEXPORT jdouble JNICALL Java_QDMath_log 
    (JNIEnv *env, jclass class, jdouble val) 
{ 
    jint tmp = (*(1 + (jint *) &val)); 
    jdouble p = ((jdouble) tmp - 1072632447)/1512775; 
    return p; 
} 

JNIEXPORT jdouble JNICALL Java_QDMath_sqrt 
    (JNIEnv *env, jclass class, jdouble val) 
{ 
    jlong tmp = ((*(jlong *) &val) - 1065353216)>>1; 
    return *(jdouble *) &tmp; 
} 

Answer 1

一眼看去，语法看起来正确（尽管编译器会告诉你比我们更快）。

但是，你在做什么看起来相当严峻。重新解释double s为int s是平台相关的（考虑endianness和sizeof问题），并且也会落入“严格别名”规则的犯规。

Answer 2

一般来说，这看起来很不对，并为我提出了一些红旗，但你可能没问题。

使用Java类型是有帮助的。如果您投射到int，则在某些机器上会出现问题，因为ints可能是32位或64位。使用jint，它应该始终是32位，所以在这方面你是安全的。

我还不完全清楚你正在尝试什么。例如，该行：

*(1 + (jint *) &p) = tmp; 

你所说的做的是把P的地址（双），把它当作一个Java的整数地址，然后添加一个到它。这将意味着采取加倍的位置，并看看4个字节......将您置于IEEE encoded double的小数部分。

然后，您将该位置设置为tmp中的值...这是一个整数，而不是IEEE编码的值。

除非您通过直接操作double的位来做某件事incredibly clever，否则您正在做的事情会产生无意义的结果。