Difference between revisions of "cpp/numeric/math/remainder"
(Applied P1467R9.) |
(Added the missing additional overload notes.) |
||
Line 59: | Line 59: | ||
* {{lc|FE_INEXACT}} is never raised, the result is always exact. | * {{lc|FE_INEXACT}} is never raised, the result is always exact. | ||
* If {{c|x}} is ±∞ and {{c|y}} is not NaN, NaN is returned and {{lc|FE_INVALID}} is raised | * If {{c|x}} is ±∞ and {{c|y}} is not NaN, NaN is returned and {{lc|FE_INVALID}} is raised | ||
− | * If {{c|y}} is ±0 and {{c|x}} is not NaN, NaN is returned | + | * If {{c|y}} is ±0 and {{c|x}} is not NaN, NaN is returned and {{lc|FE_INVALID}} is raised |
* If either argument is NaN, NaN is returned | * If either argument is NaN, NaN is returned | ||
Line 66: | Line 66: | ||
{{lc|std::fmod}}, but not {{tt|std::remainder}} is useful for doing silent wrapping of floating-point types to unsigned integer types: {{c|1=(0.0 <= (y = std::fmod(std::rint(x), 65536.0)) ? y : 65536.0 + y)}} is in the range {{tt|[-0.0 .. 65535.0]}}, which corresponds to {{c/core|unsigned short}}, but {{c|std::remainder(std::rint(x), 65536.0)}} is in the range {{tt|[-32767.0, +32768.0]}}, which is outside of the range of {{c/core|signed short}}. | {{lc|std::fmod}}, but not {{tt|std::remainder}} is useful for doing silent wrapping of floating-point types to unsigned integer types: {{c|1=(0.0 <= (y = std::fmod(std::rint(x), 65536.0)) ? y : 65536.0 + y)}} is in the range {{tt|[-0.0 .. 65535.0]}}, which corresponds to {{c/core|unsigned short}}, but {{c|std::remainder(std::rint(x), 65536.0)}} is in the range {{tt|[-32767.0, +32768.0]}}, which is outside of the range of {{c/core|signed short}}. | ||
+ | |||
+ | {{cpp/numeric/math/additional overload note|remainder}} | ||
===Example=== | ===Example=== |
Revision as of 19:34, 14 March 2023
Defined in header <cmath>
|
||
(1) | ||
float remainder ( float x, float y ); double remainder ( double x, double y ); |
(since C++11) (until C++23) |
|
constexpr /* floating-point-type */ remainder ( /* floating-point-type */ x, |
(since C++23) | |
float remainderf( float x, float y ); |
(2) | (since C++11) (constexpr since C++23) |
long double remainderl( long double x, long double y ); |
(3) | (since C++11) (constexpr since C++23) |
Defined in header <cmath>
|
||
template< class Arithmetic1, class Arithmetic2 > /* common-floating-point-type */ |
(A) | (since C++11) (constexpr since C++23) |
std::remainder
for all cv-unqualified floating-point types as the type of the parameters x and y.(since C++23)The IEEE floating-point remainder of the division operation x / y calculated by this function is exactly the value x - quo * y, where the value quo is the integral value nearest the exact value x / y. When |quo-x/y| = ½, the value quo is chosen to be even.
In contrast to std::fmod, the returned value is not guaranteed to have the same sign as x.
If the returned value is zero, it will have the same sign as x.
Contents |
Parameters
x, y | - | floating-point or integer values |
Return value
If successful, returns the IEEE floating-point remainder of the division x / y as defined above.
If a domain error occurs, an implementation-defined value is returned (NaN where supported)
If a range error occurs due to underflow, the correct result is returned.
If y is zero, but the domain error does not occur, zero is returned.
Error handling
Errors are reported as specified in math_errhandling.
Domain error may occur if y is zero.
If the implementation supports IEEE floating-point arithmetic (IEC 60559),
- The current rounding mode has no effect.
- FE_INEXACT is never raised, the result is always exact.
- If x is ±∞ and y is not NaN, NaN is returned and FE_INVALID is raised
- If y is ±0 and x is not NaN, NaN is returned and FE_INVALID is raised
- If either argument is NaN, NaN is returned
Notes
POSIX requires that a domain error occurs if x is infinite or y is zero.
std::fmod, but not std::remainder
is useful for doing silent wrapping of floating-point types to unsigned integer types: (0.0 <= (y = std::fmod(std::rint(x), 65536.0)) ? y : 65536.0 + y) is in the range [-0.0 .. 65535.0]
, which corresponds to unsigned short, but std::remainder(std::rint(x), 65536.0) is in the range [-32767.0, +32768.0]
, which is outside of the range of signed short.
The additional overloads are not required to be provided exactly as (A). They only need to be sufficient to ensure that for their first argument num1 and second argument num2:
|
(until C++23) |
If num1 and num2 have arithmetic types, then std::remainder(num1, num2) has the same effect as std::remainder(static_cast</* common-floating-point-type */>(num1), If no such floating-point type with the greatest rank and subrank exists, then overload resolution does not result in a usable candidate from the overloads provided. |
(since C++23) |
Example
#include <cfenv> #include <cmath> #include <iostream> // #pragma STDC FENV_ACCESS ON int main() { std::cout << "remainder(+5.1, +3.0) = " << std::remainder(5.1, 3) << '\n' << "remainder(-5.1, +3.0) = " << std::remainder(-5.1, 3) << '\n' << "remainder(+5.1, -3.0) = " << std::remainder(5.1, -3) << '\n' << "remainder(-5.1, -3.0) = " << std::remainder(-5.1, -3) << '\n'; // special values std::cout << "remainder(-0.0, 1.0) = " << std::remainder(-0.0, 1) << '\n' << "remainder(5.1, Inf) = " << std::remainder(5.1, INFINITY) << '\n'; // error handling std::feclearexcept(FE_ALL_EXCEPT); std::cout << "remainder(+5.1, 0) = " << std::remainder(5.1, 0) << '\n'; if (fetestexcept(FE_INVALID)) std::cout << " FE_INVALID raised\n"; }
Possible output:
remainder(+5.1, +3.0) = -0.9 remainder(-5.1, +3.0) = 0.9 remainder(+5.1, -3.0) = -0.9 remainder(-5.1, -3.0) = 0.9 remainder(-0.0, 1.0) = -0 remainder(5.1, Inf) = 5.1 remainder(+5.1, 0) = -nan FE_INVALID raised
See also
(C++11) |
computes quotient and remainder of integer division (function) |
(C++11)(C++11) |
remainder of the floating point division operation (function) |
(C++11)(C++11)(C++11) |
signed remainder as well as the three last bits of the division operation (function) |
C documentation for remainder
|