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，否則您正在做的事情會產生無意義的結果。