問題タブ [quadruple-precision]

つい最近、GCC 4.6.0 がlibquadmathとともにリリースされました。残念ながら、GNU は Fortran をサポートしていますが、C や C++ はサポートしていません (含まれているのは .so だけです)。私はこれらの新機能を C++ で使用する方法を見つけていませんが、GNU C は__float128保証された 4 倍精度浮動小数点数の型をサポートしています。GNU C は、 libquadmathなど の数学関数をサポートしていないようですfabsq(絶対値、qquad の接尾辞)。

これらの関数を C++ で動作させる方法はありますか、または数学関数に使用できる代替ライブラリはあり__float128ますか? GCC で動作する 4 倍精度浮動小数点数を取得するための最良の方法は何ですか? 現在、それらを加算、減算、および乗算できますが、それらを文字列に変換したり、truncqおよびなどの関数fabsqを使用して独自の文字列関数を作成したりする方法がないことを考えると、これは役に立ちません。

コードを実行しようとしました

このコマンドで

しかし、私はエラーが発生しました：

バージョンは次のとおりです。

どうすれば問題を解決できますか？

を使用し<quadmath.h>ます。入力を正しく読み取ることができる引数の型はどれですか? double を使用すると、次のようになります。

たとえば、「l」の代わりにさまざまな可能性を試しました。

またはなくても

ただし、エラーが発生します。

標準では、IEEE754厳密に正の (非正規) 値の最小値は、4 倍精度の浮動小数点形式を使用して2 ⁻¹⁶⁴⁹³ ≈ 10 ^{−4965です。GCC が 10 -4949}未満のものを拒否するのはなぜですか? ^{制限を 10 -4965}ではなく10 ^-4949に決定する、その下で起こっている可能性のあるさまざまなことの説明を探しています。

オンラインで GNU GCC バージョン 4.8.1を使用しています。どのアーキテクチャで実行されているかわかりません (これが原因である可能性があります)。さまざまなアーキテクチャからの調査結果を自由に投稿してください。

の平方根を求めるにはどうすればよいReal(Real128)ですか?

C++ で記述された高精度 ODE (常微分方程式) ソルバーがあります。私はすべての計算をユーザー定義型で行いますreal_type。ヘッダーにこの型を宣言する typedef があります。

typedef long double real_type;

__float128より精度を上げるために、ロングダブルタイプに変更することにしました。これに加えて、quadmath.hすべての標準数学関数を含め、libquadmath のものに置き換えました。

「long double」バージョンが最適化フラグなしでビルドされた場合、一部の参照 ODE は 77 秒で解決されます。このバージョンが -O3 フラグでビルドされている場合、同じ ODE が 25 秒で解決されます。したがって、-O3 フラグは計算を 3 倍高速化します。

しかし、フラグなしでビルドされた "__float 128" バージョンでは、同様の ODE が 190 秒で解決され、-O3 を使用すると 160 秒で解決されます (~ 15% の差)。-O3 最適化が 4 倍精度の計算に弱い影響を与えるのはなぜですか? 他のコンパイラフラグを使用するか、他のライブラリを含める必要がありますか?

4 倍精度の quadmath ライブラリを使用しています。fwrite 関数を使用してバイナリ ファイルに書き込むことはできますか?