std::log1p, std::log1pf, std::log1pl
From cppreference.com
| ヘッダー <cmath> で定義 |
||
| (1) | ||
float log1p ( float num ); double log1p ( double num ); |
(C++23まで) | |
| /*浮動小数点数型*/ log1p ( /*浮動小数点型*/ num ); |
(C++23から) (C++26 以降 constexpr) |
|
float log1pf( float num ); |
(2) | (C++11以降) (C++26 以降 constexpr) |
long double log1pl( long double num ); |
(3) | (C++11以降) (C++26 以降 constexpr) |
| SIMDオーバーロード (C++26以降) |
||
| ヘッダー <simd> で定義 |
||
| template< /*math-floating-point*/ V > constexpr /*deduced-simd-t*/<V> |
(S) | (C++26以降) |
| 追加のオーバーロード (C++11以降) |
||
| ヘッダー <cmath> で定義 |
||
template< class Integer > double log1p ( Integer num ); |
(A) | (C++26 以降 constexpr) |
1-3) 1 + num の自然対数(底がe)を計算します。この関数は、num がゼロに近い場合、式 std::log(1 + num) よりも精度が高いです。ライブラリは、パラメータの型としてすべての cv 修飾されていない浮動小数点型に対する
std::log1p のオーバーロードを提供します。(since C++23)|
S) SIMD オーバーロードは、v_num に対して要素ごとの
std::log1p を実行します。
|
(C++26以降) |
|
A) すべての整数型に対する追加のオーバーロードが提供されます。これらは double として扱われます。
|
(C++11以降) |
目次 |
[edit] Parameters
| num | - | 浮動小数点数または整数値 |
[edit] Return value
エラーが発生しない場合、ln(1+num) が返されます。
領域エラーが発生した場合、実装定義の値が返される (サポートされている場合はNaN)。
If a pole error occurs, -HUGE_VAL, -HUGE_VALF, or -HUGE_VALL is returned.
アンダーフローによる範囲エラーが発生した場合、正確な結果 (丸め後) が返される。
[edit] Error handling
エラーは math_errhandling で指定された通りに報告される。
num が -1 未満の場合、ドメインエラーが発生します。
num が -1 の場合、極値エラーが発生する可能性があります。
実装がIEEE浮動小数点算術 (IEC 60559) をサポートしている場合、
- 引数が ±0 の場合、引数は変更されずに返されます。
- 引数が -1 の場合、-∞ が返され、FE_DIVBYZERO が発生します。
- 引数が -1 未満の場合、NaN が返され、FE_INVALID が発生します。
- 引数が +∞ の場合、+∞ が返されます。
- 引数が NaN の場合、NaN が返されます。
[edit] Notes
std::expm1 および std::log1p 関数は、例えば小さな日歩を計算する際の財務計算に役立ちます。 (1 + x)n
- 1 は std::expm1(n * std::log1p(x)) と表現できます。これらの関数は、逆双曲線関数の正確な記述を単純化します。
追加のオーバーロードは、(A) とまったく同じように提供される必要はありません。整数型の引数 num に対して、std::log1p(num) が std::log1p(static_cast<double>(num)) と同じ効果を持つことを保証するのに十分であればよいのです。
[edit] Example
このコードを実行
#include <cerrno> #include <cfenv> #include <cmath> #include <cstring> #include <iostream> // #pragma STDC FENV_ACCESS ON int main() { std::cout << "log1p(0) = " << log1p(0) << '\n' << "Interest earned in 2 days on $100, compounded daily at 1%\n" << " on a 30/360 calendar = " << 100 * expm1(2 * log1p(0.01 / 360)) << '\n' << "log(1+1e-16) = " << std::log(1 + 1e-16) << ", but log1p(1e-16) = " << std::log1p(1e-16) << '\n'; // special values std::cout << "log1p(-0) = " << std::log1p(-0.0) << '\n' << "log1p(+Inf) = " << std::log1p(INFINITY) << '\n'; // error handling errno = 0; std::feclearexcept(FE_ALL_EXCEPT); std::cout << "log1p(-1) = " << std::log1p(-1) << '\n'; if (errno == ERANGE) std::cout << " errno == ERANGE: " << std::strerror(errno) << '\n'; if (std::fetestexcept(FE_DIVBYZERO)) std::cout << " FE_DIVBYZERO raised\n"; }
実行結果の例
log1p(0) = 0
Interest earned in 2 days on $100, compounded daily at 1%
on a 30/360 calendar = 0.00555563
log(1+1e-16) = 0, but log1p(1e-16) = 1e-16
log1p(-0) = -0
log1p(+Inf) = inf
log1p(-1) = -inf
errno == ERANGE: Result too large
FE_DIVBYZERO raised[edit] See also
| (C++11)(C++11) |
自然対数 (e を底とする) を計算する (ln(x)) (関数) |
| (C++11)(C++11) |
常用対数 (10 を底とする) を計算する (log10(x)) (関数) |
| (C++11)(C++11)(C++11) |
与えられた数値の 2 を底とする対数 (log2(x)) (関数) |
| (C++11)(C++11)(C++11) |
与えられたべき乗に累乗した e から 1 を引いた値を返す (ex-1) (関数) |
| C ドキュメント for log1p
| |