質問
このサイトによると、エラー関数erf(x) math.hからしかし、実際にはmath.hを見ると、そこにはありません。g++はできるのに、gccは次のテストプログラムをコンパイルできません。
#include <math.h>
#include <stdio.h>
int main(int argc, char* argv[]) {
double x;
double erfX;
x = 1.0;
erfX = erf(x);
printf("erf(%f) = %f", x, erfX);
}
$ gcc mathHTest.c
/tmp/ccWfNox5.o: In function `main':
mathHTest.c:(.text+0x28): undefined reference to `erf'
collect2: ld returned 1 exit status
$ g++ mathHTest.c
g ++がgccでプルしないものは何ですか? / usr / includeを見ると、erf(x)を見つけることができる唯一の場所はtgmath.hでしたが、これは含めません。 g ++はgccとは異なるヘッダーを取得する必要がありますが、どのヘッダーを取得しますか?
編集:libmでgccとリンクしていなかったため、リンクエラーが発生しました。ただし、なぜerf()がmath.hに含まれていないのか、まだわかりません。どこから来たのですか?
解決
'erf'は実際にはbits / mathcalls.hで宣言されており、math.hに#includeされています。実際の宣言は、CとC ++の両方で正しいことをするためにマクロマジックによってかなり隠されています
他のヒント
同様の問題があり、 erf
の正確な定義を見つける必要があったので、これについて詳しく説明します。クリス・ドッドが言ったように、関数は maths.h
に含まれる bits / mathcalls.h
で宣言されています。
bits / mathcalls.h
:
...
#if defined __USE_MISC || defined __USE_XOPEN || defined __USE_ISOC99
__BEGIN_NAMESPACE_C99
/* Error and gamma functions. */
__MATHCALL (erf,, (_Mdouble_));
__MATHCALL (erfc,, (_Mdouble_));
__MATHCALL (lgamma,, (_Mdouble_));
__END_NAMESPACE_C99
#endif
...
マクロマジックは __ MATHCALL(erf ,,(_Mdouble _));
を
extern double erf (double) throw (); extern double __erf (double) throw ();
実際のコードは libm.a
または libm.so
( gcc -lm
)にあります:
$ nm /usr/lib/libm.a
...
s_erf.o:
00000400 T __erf
00000000 T __erfc
U __ieee754_exp
00000400 W erf
00000000 W erfc
...
ソースは gnu libc から取得できます。ウェブページ。実際の実装に関する大まかなアイデアについては、ソースの数行をここに示します。
sysdeps / ieee754 / dbl-64 / s_erf.c
:
/* double erf(double x)
* double erfc(double x)
* x
* 2 |\
* erf(x) = --------- | exp(-t*t)dt
* sqrt(pi) \|
* 0
*
* erfc(x) = 1-erf(x)
* Note that
* erf(-x) = -erf(x)
* erfc(-x) = 2 - erfc(x)
*
* Method:
* 1. For |x| in [0, 0.84375]
* erf(x) = x + x*R(x^2)
* erfc(x) = 1 - erf(x) if x in [-.84375,0.25]
* = 0.5 + ((0.5-x)-x*R) if x in [0.25,0.84375]
* where R = P/Q where P is an odd poly of degree 8 and
* Q is an odd poly of degree 10.
* -57.90
* | R - (erf(x)-x)/x | <= 2
*
*
* Remark. The formula is derived by noting
* erf(x) = (2/sqrt(pi))*(x - x^3/3 + x^5/10 - x^7/42 + ....)
* and that
* 2/sqrt(pi) = 1.128379167095512573896158903121545171688
* is close to one. The interval is chosen because the fix
* point of erf(x) is near 0.6174 (i.e., erf(x)=x when x is
* near 0.6174), and by some experiment, 0.84375 is chosen to
* guarantee the error is less than one ulp for erf.
*
* 2. For |x| in [0.84375,1.25], let s = |x| - 1, and
...
数学ライブラリ( libm
)もリンクする必要があります:
$ gcc mathHTest.c -lm
通常の数学ライブラリ関数はすべて実際にあり、標準Cライブラリ( libc
)にはありません。
私のテストによると、 g ++
には libm
が自動的に含まれますが、 gcc
には含まれません。
x86プロセッサでcygwinのgccを使用しても同じ問題が発生しました。 「-lm」ライブラリインクルードパラメータ(ファイルリストの後!)は完全に機能しました。