mirror/Sprout
1
0
Fork 0
mirror of https://github.com/bolero-MURAKAMI/Sprout.git synced 2025-01-25 20:56:38 +00:00
Code Issues Projects Releases Packages Wiki Activity Actions

sprout::proj should return an argument as it is.

Browse source 
sprout::proj should return the argument as it is unless the argument is infinite.
This commit is contained in:
Mitsuru Kariya 2014-07-13 02:49:50 +09:00
parent 37a8e6c6e7
commit 32f72555e1
1 changed files with 6 additions and 14 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

20
sprout/complex/values.hpp
View file

@ -15,6 +15,9 @@
#include <sprout/math/cos.hpp>
#include <sprout/math/atan2.hpp>
#include <sprout/math/sqrt.hpp>
#include <sprout/math/isinf.hpp>
#include <sprout/math/copysign.hpp>
#include <sprout/limits.hpp>
namespace sprout {
template<typename T>
@ -62,23 +65,12 @@ namespace sprout {
conj(sprout::complex<T> const& x) {
return sprout::complex<T>(x.real(), -x.imag());
}
namespace detail {
template<typename T>
inline SPROUT_CONSTEXPR sprout::complex<T>
proj_impl(sprout::complex<T> const& x, T const& den) {
return sprout::complex<T>(
T(2) * x.real() / den,
T(2) * x.imag() / den
);
}
} // detail
template<typename T>
inline SPROUT_CONSTEXPR sprout::complex<T>
proj(sprout::complex<T> const& x) {
return sprout::detail::proj_impl(
x,
sprout::norm(x) + T(1)
);
return sprout::isinf(x.real()) || sprout::isinf(x.imag())
? sprout::complex<T>(sprout::numeric_limits<T>::infinity(), sprout::copysign(T(), x.imag()))
: x;
}
template<typename T>
inline SPROUT_CONSTEXPR sprout::complex<T>