1 //:
2 // \file
3 // \author Peter Vanroose
4 // \date 6 September 2002
5 //
6 //-----------------------------------------------------------------------------
7
8 #include <complex>
9 #include "vnl_rational_traits.h"
10
11 const vnl_rational vnl_numeric_traits<vnl_rational>::zero = vnl_rational(0L, 1L);
12 const vnl_rational vnl_numeric_traits<vnl_rational>::one = vnl_rational(1L, 1L);
13 const vnl_rational vnl_numeric_traits<vnl_rational>::maxval = vnl_rational(vnl_numeric_traits<long>::maxval, 1L);
14
15 const std::complex<vnl_rational> vnl_numeric_traits<std::complex<vnl_rational>>::zero =
16 std::complex<vnl_rational>(vnl_rational(0L, 1L), vnl_rational(0L, 1L));
17 const std::complex<vnl_rational> vnl_numeric_traits<std::complex<vnl_rational>>::one =
18 std::complex<vnl_rational>(vnl_rational(1L, 1L), vnl_rational(0L, 1L));
19
20 namespace vnl_math
21 {
22 vnl_rational
squared_magnitude(std::complex<vnl_rational> const & x)23 squared_magnitude(std::complex<vnl_rational> const & x)
24 {
25 return x.real() * x.real() + x.imag() * x.imag();
26 }
27
28 vnl_rational
abs(std::complex<vnl_rational> const & x)29 abs(std::complex<vnl_rational> const & x)
30 {
31 return vnl_rational(std::sqrt(double(x.real() * x.real() + x.imag() * x.imag())));
32 }
33 } // namespace vnl_math
34
35 std::ostream &
operator <<(std::ostream & os,const std::complex<vnl_rational> & x)36 operator<<(std::ostream & os, const std::complex<vnl_rational> & x)
37 {
38 return os << x.real() << '+' << x.imag() << 'j';
39 }
40