/* ***************************************************** THIS IS AN AUTOMATICALLY GENERATED FILE. DO NOT EDIT. ***************************************************** Generated by: scitbx.source_generators.array_family.generate_algebras */ #ifndef SCITBX_ARRAY_FAMILY_SMALL_ALGEBRA_H #define SCITBX_ARRAY_FAMILY_SMALL_ALGEBRA_H #ifndef DOXYGEN_SHOULD_SKIP_THIS #include #if (defined(BOOST_MSVC) && BOOST_MSVC <= 1300) // VC++ 7.0 #define SCITBX_ARRAY_FAMILY_SMALL_ALGEBRA_MIN_N1_N2 N1 #else #define SCITBX_ARRAY_FAMILY_SMALL_ALGEBRA_MIN_N1_N2 (N1 #include #include #include #include namespace scitbx { namespace af { template inline small operator-(small const& a) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a.size(), make_init_functor(make_array_functor_a( fn::functor_negate< return_element_type, ElementType>(), a.begin()))); } template inline small operator!(small const& a) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a.size(), make_init_functor(make_array_functor_a( fn::functor_logical_not< return_element_type, ElementType>(), a.begin()))); } template< typename ElementType1, std::size_t N1, typename ElementType2, std::size_t N2> inline small< typename binary_operator_traits< ElementType1, ElementType2>::arithmetic, SCITBX_ARRAY_FAMILY_SMALL_ALGEBRA_MIN_N1_N2> operator+( small const& a1, small const& a2) { typedef small< typename binary_operator_traits< ElementType1, ElementType2>::arithmetic, SCITBX_ARRAY_FAMILY_SMALL_ALGEBRA_MIN_N1_N2> return_array_type; typedef typename return_array_type::value_type return_element_type; if (a1.size() != a2.size()) throw_range_error(); return return_array_type(a1.size(), make_init_functor(make_array_functor_a_a( fn::functor_plus< return_element_type, ElementType1, ElementType2>(), a1.begin(), a2.begin()))); } template inline small operator+( small const& a1, ElementType const& a2) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a1.size(), make_init_functor(make_array_functor_a_s( fn::functor_plus< return_element_type, ElementType, ElementType>(), a1.begin(), a2))); } template inline small operator+( ElementType const& a1, small const& a2) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a2.size(), make_init_functor(make_array_functor_s_a( fn::functor_plus< return_element_type, ElementType, ElementType>(), a1, a2.begin()))); } template< typename ElementType1, std::size_t N1, typename ElementType2, std::size_t N2> inline small& operator+=( small& a1, small const& a2) { if (a1.size() != a2.size()) throw_range_error(); array_operation_in_place_a_a(fn::functor_ip_plus< ElementType1, ElementType2>(), a1.begin(), a2.begin(), a1.size()); return a1; } template inline small& operator+=( small& a1, ElementType const& a2) { array_operation_in_place_a_s(fn::functor_ip_plus< ElementType, ElementType>(), a1.begin(), a2, a1.size()); return a1; } template< typename ElementType1, std::size_t N1, typename ElementType2, std::size_t N2> inline small< typename binary_operator_traits< ElementType1, ElementType2>::arithmetic, SCITBX_ARRAY_FAMILY_SMALL_ALGEBRA_MIN_N1_N2> operator-( small const& a1, small const& a2) { typedef small< typename binary_operator_traits< ElementType1, ElementType2>::arithmetic, SCITBX_ARRAY_FAMILY_SMALL_ALGEBRA_MIN_N1_N2> return_array_type; typedef typename return_array_type::value_type return_element_type; if (a1.size() != a2.size()) throw_range_error(); return return_array_type(a1.size(), make_init_functor(make_array_functor_a_a( fn::functor_minus< return_element_type, ElementType1, ElementType2>(), a1.begin(), a2.begin()))); } template inline small operator-( small const& a1, ElementType const& a2) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a1.size(), make_init_functor(make_array_functor_a_s( fn::functor_minus< return_element_type, ElementType, ElementType>(), a1.begin(), a2))); } template inline small operator-( ElementType const& a1, small const& a2) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a2.size(), make_init_functor(make_array_functor_s_a( fn::functor_minus< return_element_type, ElementType, ElementType>(), a1, a2.begin()))); } template< typename ElementType1, std::size_t N1, typename ElementType2, std::size_t N2> inline small& operator-=( small& a1, small const& a2) { if (a1.size() != a2.size()) throw_range_error(); array_operation_in_place_a_a(fn::functor_ip_minus< ElementType1, ElementType2>(), a1.begin(), a2.begin(), a1.size()); return a1; } template inline small& operator-=( small& a1, ElementType const& a2) { array_operation_in_place_a_s(fn::functor_ip_minus< ElementType, ElementType>(), a1.begin(), a2, a1.size()); return a1; } template< typename ElementType1, std::size_t N1, typename ElementType2, std::size_t N2> inline small< typename binary_operator_traits< ElementType1, ElementType2>::arithmetic, SCITBX_ARRAY_FAMILY_SMALL_ALGEBRA_MIN_N1_N2> operator*( small const& a1, small const& a2) { typedef small< typename binary_operator_traits< ElementType1, ElementType2>::arithmetic, SCITBX_ARRAY_FAMILY_SMALL_ALGEBRA_MIN_N1_N2> return_array_type; typedef typename return_array_type::value_type return_element_type; if (a1.size() != a2.size()) throw_range_error(); return return_array_type(a1.size(), make_init_functor(make_array_functor_a_a( fn::functor_multiplies< return_element_type, ElementType1, ElementType2>(), a1.begin(), a2.begin()))); } template inline small operator*( small const& a1, ElementType const& a2) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a1.size(), make_init_functor(make_array_functor_a_s( fn::functor_multiplies< return_element_type, ElementType, ElementType>(), a1.begin(), a2))); } template inline small operator*( ElementType const& a1, small const& a2) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a2.size(), make_init_functor(make_array_functor_s_a( fn::functor_multiplies< return_element_type, ElementType, ElementType>(), a1, a2.begin()))); } template< typename ElementType1, std::size_t N1, typename ElementType2, std::size_t N2> inline small& operator*=( small& a1, small const& a2) { if (a1.size() != a2.size()) throw_range_error(); array_operation_in_place_a_a(fn::functor_ip_multiplies< ElementType1, ElementType2>(), a1.begin(), a2.begin(), a1.size()); return a1; } template inline small& operator*=( small& a1, ElementType const& a2) { array_operation_in_place_a_s(fn::functor_ip_multiplies< ElementType, ElementType>(), a1.begin(), a2, a1.size()); return a1; } template< typename ElementType1, std::size_t N1, typename ElementType2, std::size_t N2> inline small< typename binary_operator_traits< ElementType1, ElementType2>::arithmetic, SCITBX_ARRAY_FAMILY_SMALL_ALGEBRA_MIN_N1_N2> operator/( small const& a1, small const& a2) { typedef small< typename binary_operator_traits< ElementType1, ElementType2>::arithmetic, SCITBX_ARRAY_FAMILY_SMALL_ALGEBRA_MIN_N1_N2> return_array_type; typedef typename return_array_type::value_type return_element_type; if (a1.size() != a2.size()) throw_range_error(); return return_array_type(a1.size(), make_init_functor(make_array_functor_a_a( fn::functor_divides< return_element_type, ElementType1, ElementType2>(), a1.begin(), a2.begin()))); } template inline small operator/( small const& a1, ElementType const& a2) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a1.size(), make_init_functor(make_array_functor_a_s( fn::functor_divides< return_element_type, ElementType, ElementType>(), a1.begin(), a2))); } template inline small operator/( ElementType const& a1, small const& a2) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a2.size(), make_init_functor(make_array_functor_s_a( fn::functor_divides< return_element_type, ElementType, ElementType>(), a1, a2.begin()))); } template< typename ElementType1, std::size_t N1, typename ElementType2, std::size_t N2> inline small& operator/=( small& a1, small const& a2) { if (a1.size() != a2.size()) throw_range_error(); array_operation_in_place_a_a(fn::functor_ip_divides< ElementType1, ElementType2>(), a1.begin(), a2.begin(), a1.size()); return a1; } template inline small& operator/=( small& a1, ElementType const& a2) { array_operation_in_place_a_s(fn::functor_ip_divides< ElementType, ElementType>(), a1.begin(), a2, a1.size()); return a1; } template< typename ElementType1, std::size_t N1, typename ElementType2, std::size_t N2> inline small< typename binary_operator_traits< ElementType1, ElementType2>::arithmetic, SCITBX_ARRAY_FAMILY_SMALL_ALGEBRA_MIN_N1_N2> operator%( small const& a1, small const& a2) { typedef small< typename binary_operator_traits< ElementType1, ElementType2>::arithmetic, SCITBX_ARRAY_FAMILY_SMALL_ALGEBRA_MIN_N1_N2> return_array_type; typedef typename return_array_type::value_type return_element_type; if (a1.size() != a2.size()) throw_range_error(); return return_array_type(a1.size(), make_init_functor(make_array_functor_a_a( fn::functor_modulus< return_element_type, ElementType1, ElementType2>(), a1.begin(), a2.begin()))); } template inline small operator%( small const& a1, ElementType const& a2) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a1.size(), make_init_functor(make_array_functor_a_s( fn::functor_modulus< return_element_type, ElementType, ElementType>(), a1.begin(), a2))); } template inline small operator%( ElementType const& a1, small const& a2) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a2.size(), make_init_functor(make_array_functor_s_a( fn::functor_modulus< return_element_type, ElementType, ElementType>(), a1, a2.begin()))); } template< typename ElementType1, std::size_t N1, typename ElementType2, std::size_t N2> inline small& operator%=( small& a1, small const& a2) { if (a1.size() != a2.size()) throw_range_error(); array_operation_in_place_a_a(fn::functor_ip_modulus< ElementType1, ElementType2>(), a1.begin(), a2.begin(), a1.size()); return a1; } template inline small& operator%=( small& a1, ElementType const& a2) { array_operation_in_place_a_s(fn::functor_ip_modulus< ElementType, ElementType>(), a1.begin(), a2, a1.size()); return a1; } template< typename ElementType1, std::size_t N1, typename ElementType2, std::size_t N2> inline small operator&&( small const& a1, small const& a2) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; if (a1.size() != a2.size()) throw_range_error(); return return_array_type(a1.size(), make_init_functor(make_array_functor_a_a( fn::functor_logical_and< return_element_type, ElementType1, ElementType2>(), a1.begin(), a2.begin()))); } template inline small operator&&( small const& a1, ElementType const& a2) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a1.size(), make_init_functor(make_array_functor_a_s( fn::functor_logical_and< return_element_type, ElementType, ElementType>(), a1.begin(), a2))); } template inline small operator&&( ElementType const& a1, small const& a2) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a2.size(), make_init_functor(make_array_functor_s_a( fn::functor_logical_and< return_element_type, ElementType, ElementType>(), a1, a2.begin()))); } template< typename ElementType1, std::size_t N1, typename ElementType2, std::size_t N2> inline small operator||( small const& a1, small const& a2) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; if (a1.size() != a2.size()) throw_range_error(); return return_array_type(a1.size(), make_init_functor(make_array_functor_a_a( fn::functor_logical_or< return_element_type, ElementType1, ElementType2>(), a1.begin(), a2.begin()))); } template inline small operator||( small const& a1, ElementType const& a2) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a1.size(), make_init_functor(make_array_functor_a_s( fn::functor_logical_or< return_element_type, ElementType, ElementType>(), a1.begin(), a2))); } template inline small operator||( ElementType const& a1, small const& a2) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a2.size(), make_init_functor(make_array_functor_s_a( fn::functor_logical_or< return_element_type, ElementType, ElementType>(), a1, a2.begin()))); } template< typename ElementType1, std::size_t N1, typename ElementType2, std::size_t N2> inline small operator==( small const& a1, small const& a2) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; if (a1.size() != a2.size()) throw_range_error(); return return_array_type(a1.size(), make_init_functor(make_array_functor_a_a( fn::functor_equal_to< return_element_type, ElementType1, ElementType2>(), a1.begin(), a2.begin()))); } template inline small operator==( small const& a1, ElementType const& a2) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a1.size(), make_init_functor(make_array_functor_a_s( fn::functor_equal_to< return_element_type, ElementType, ElementType>(), a1.begin(), a2))); } template inline small operator==( ElementType const& a1, small const& a2) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a2.size(), make_init_functor(make_array_functor_s_a( fn::functor_equal_to< return_element_type, ElementType, ElementType>(), a1, a2.begin()))); } template< typename ElementType1, std::size_t N1, typename ElementType2, std::size_t N2> inline small operator!=( small const& a1, small const& a2) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; if (a1.size() != a2.size()) throw_range_error(); return return_array_type(a1.size(), make_init_functor(make_array_functor_a_a( fn::functor_not_equal_to< return_element_type, ElementType1, ElementType2>(), a1.begin(), a2.begin()))); } template inline small operator!=( small const& a1, ElementType const& a2) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a1.size(), make_init_functor(make_array_functor_a_s( fn::functor_not_equal_to< return_element_type, ElementType, ElementType>(), a1.begin(), a2))); } template inline small operator!=( ElementType const& a1, small const& a2) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a2.size(), make_init_functor(make_array_functor_s_a( fn::functor_not_equal_to< return_element_type, ElementType, ElementType>(), a1, a2.begin()))); } template< typename ElementType1, std::size_t N1, typename ElementType2, std::size_t N2> inline small operator>( small const& a1, small const& a2) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; if (a1.size() != a2.size()) throw_range_error(); return return_array_type(a1.size(), make_init_functor(make_array_functor_a_a( fn::functor_greater< return_element_type, ElementType1, ElementType2>(), a1.begin(), a2.begin()))); } template inline small operator>( small const& a1, ElementType const& a2) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a1.size(), make_init_functor(make_array_functor_a_s( fn::functor_greater< return_element_type, ElementType, ElementType>(), a1.begin(), a2))); } template inline small operator>( ElementType const& a1, small const& a2) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a2.size(), make_init_functor(make_array_functor_s_a( fn::functor_greater< return_element_type, ElementType, ElementType>(), a1, a2.begin()))); } template< typename ElementType1, std::size_t N1, typename ElementType2, std::size_t N2> inline small operator<( small const& a1, small const& a2) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; if (a1.size() != a2.size()) throw_range_error(); return return_array_type(a1.size(), make_init_functor(make_array_functor_a_a( fn::functor_less< return_element_type, ElementType1, ElementType2>(), a1.begin(), a2.begin()))); } template inline small operator<( small const& a1, ElementType const& a2) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a1.size(), make_init_functor(make_array_functor_a_s( fn::functor_less< return_element_type, ElementType, ElementType>(), a1.begin(), a2))); } template inline small operator<( ElementType const& a1, small const& a2) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a2.size(), make_init_functor(make_array_functor_s_a( fn::functor_less< return_element_type, ElementType, ElementType>(), a1, a2.begin()))); } template< typename ElementType1, std::size_t N1, typename ElementType2, std::size_t N2> inline small operator>=( small const& a1, small const& a2) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; if (a1.size() != a2.size()) throw_range_error(); return return_array_type(a1.size(), make_init_functor(make_array_functor_a_a( fn::functor_greater_equal< return_element_type, ElementType1, ElementType2>(), a1.begin(), a2.begin()))); } template inline small operator>=( small const& a1, ElementType const& a2) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a1.size(), make_init_functor(make_array_functor_a_s( fn::functor_greater_equal< return_element_type, ElementType, ElementType>(), a1.begin(), a2))); } template inline small operator>=( ElementType const& a1, small const& a2) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a2.size(), make_init_functor(make_array_functor_s_a( fn::functor_greater_equal< return_element_type, ElementType, ElementType>(), a1, a2.begin()))); } template< typename ElementType1, std::size_t N1, typename ElementType2, std::size_t N2> inline small operator<=( small const& a1, small const& a2) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; if (a1.size() != a2.size()) throw_range_error(); return return_array_type(a1.size(), make_init_functor(make_array_functor_a_a( fn::functor_less_equal< return_element_type, ElementType1, ElementType2>(), a1.begin(), a2.begin()))); } template inline small operator<=( small const& a1, ElementType const& a2) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a1.size(), make_init_functor(make_array_functor_a_s( fn::functor_less_equal< return_element_type, ElementType, ElementType>(), a1.begin(), a2))); } template inline small operator<=( ElementType const& a1, small const& a2) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a2.size(), make_init_functor(make_array_functor_s_a( fn::functor_less_equal< return_element_type, ElementType, ElementType>(), a1, a2.begin()))); } template inline small absolute(small const& a) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a.size(), make_init_functor(make_array_functor_a( fn::functor_absolute(), a.begin()))); } template inline small pow2(small const& a) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a.size(), make_init_functor(make_array_functor_a( fn::functor_pow2(), a.begin()))); } template inline small acos(small const& a) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a.size(), make_init_functor(make_array_functor_a( fn::functor_acos(), a.begin()))); } template inline small cos(small const& a) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a.size(), make_init_functor(make_array_functor_a( fn::functor_cos(), a.begin()))); } template inline small tan(small const& a) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a.size(), make_init_functor(make_array_functor_a( fn::functor_tan(), a.begin()))); } template inline small asin(small const& a) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a.size(), make_init_functor(make_array_functor_a( fn::functor_asin(), a.begin()))); } template inline small cosh(small const& a) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a.size(), make_init_functor(make_array_functor_a( fn::functor_cosh(), a.begin()))); } template inline small tanh(small const& a) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a.size(), make_init_functor(make_array_functor_a( fn::functor_tanh(), a.begin()))); } template inline small atan(small const& a) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a.size(), make_init_functor(make_array_functor_a( fn::functor_atan(), a.begin()))); } template inline small exp(small const& a) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a.size(), make_init_functor(make_array_functor_a( fn::functor_exp(), a.begin()))); } template inline small sin(small const& a) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a.size(), make_init_functor(make_array_functor_a( fn::functor_sin(), a.begin()))); } template inline small fabs(small const& a) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a.size(), make_init_functor(make_array_functor_a( fn::functor_fabs(), a.begin()))); } template inline small log(small const& a) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a.size(), make_init_functor(make_array_functor_a( fn::functor_log(), a.begin()))); } template inline small sinh(small const& a) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a.size(), make_init_functor(make_array_functor_a( fn::functor_sinh(), a.begin()))); } template inline small ceil(small const& a) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a.size(), make_init_functor(make_array_functor_a( fn::functor_ceil(), a.begin()))); } template inline small floor(small const& a) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a.size(), make_init_functor(make_array_functor_a( fn::functor_floor(), a.begin()))); } template inline small log10(small const& a) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a.size(), make_init_functor(make_array_functor_a( fn::functor_log10(), a.begin()))); } template inline small sqrt(small const& a) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a.size(), make_init_functor(make_array_functor_a( fn::functor_sqrt(), a.begin()))); } template inline small abs(small const& a) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a.size(), make_init_functor(make_array_functor_a( fn::functor_abs(), a.begin()))); } template inline small conj(small const& a) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a.size(), make_init_functor(make_array_functor_a( fn::functor_conj(), a.begin()))); } template< typename ElementType1, std::size_t N1, typename ElementType2, std::size_t N2> inline small fmod_positive( small const& a1, small const& a2) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; if (a1.size() != a2.size()) throw_range_error(); return return_array_type(a1.size(), make_init_functor(make_array_functor_a_a( fn::functor_fmod_positive(), a1.begin(), a2.begin()))); } template inline small fmod_positive( small const& a1, ElementType const& a2) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a1.size(), make_init_functor(make_array_functor_a_s( fn::functor_fmod_positive(), a1.begin(), a2))); } template inline small fmod_positive( ElementType const& a1, small const& a2) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a2.size(), make_init_functor(make_array_functor_s_a( fn::functor_fmod_positive(), a1, a2.begin()))); } template< typename ElementType1, std::size_t N1, typename ElementType2, std::size_t N2> inline small fmod( small const& a1, small const& a2) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; if (a1.size() != a2.size()) throw_range_error(); return return_array_type(a1.size(), make_init_functor(make_array_functor_a_a( fn::functor_fmod(), a1.begin(), a2.begin()))); } template inline small fmod( small const& a1, ElementType const& a2) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a1.size(), make_init_functor(make_array_functor_a_s( fn::functor_fmod(), a1.begin(), a2))); } template inline small fmod( ElementType const& a1, small const& a2) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a2.size(), make_init_functor(make_array_functor_s_a( fn::functor_fmod(), a1, a2.begin()))); } template< typename ElementType1, std::size_t N1, typename ElementType2, std::size_t N2> inline small pow( small const& a1, small const& a2) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; if (a1.size() != a2.size()) throw_range_error(); return return_array_type(a1.size(), make_init_functor(make_array_functor_a_a( fn::functor_pow(), a1.begin(), a2.begin()))); } template inline small pow( small const& a1, ElementType const& a2) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a1.size(), make_init_functor(make_array_functor_a_s( fn::functor_pow(), a1.begin(), a2))); } template inline small pow( ElementType const& a1, small const& a2) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a2.size(), make_init_functor(make_array_functor_s_a( fn::functor_pow(), a1, a2.begin()))); } template< typename ElementType1, std::size_t N1, typename ElementType2, std::size_t N2> inline small atan2( small const& a1, small const& a2) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; if (a1.size() != a2.size()) throw_range_error(); return return_array_type(a1.size(), make_init_functor(make_array_functor_a_a( fn::functor_atan2(), a1.begin(), a2.begin()))); } template inline small atan2( small const& a1, ElementType const& a2) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a1.size(), make_init_functor(make_array_functor_a_s( fn::functor_atan2(), a1.begin(), a2))); } template inline small atan2( ElementType const& a1, small const& a2) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a2.size(), make_init_functor(make_array_functor_s_a( fn::functor_atan2(), a1, a2.begin()))); } template< typename ElementType1, std::size_t N1, typename ElementType2, std::size_t N2> inline small each_min( small const& a1, small const& a2) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; if (a1.size() != a2.size()) throw_range_error(); return return_array_type(a1.size(), make_init_functor(make_array_functor_a_a( fn::functor_each_min(), a1.begin(), a2.begin()))); } template inline small each_min( small const& a1, ElementType const& a2) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a1.size(), make_init_functor(make_array_functor_a_s( fn::functor_each_min(), a1.begin(), a2))); } template inline small each_min( ElementType const& a1, small const& a2) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a2.size(), make_init_functor(make_array_functor_s_a( fn::functor_each_min(), a1, a2.begin()))); } template< typename ElementType1, std::size_t N1, typename ElementType2, std::size_t N2> inline small each_max( small const& a1, small const& a2) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; if (a1.size() != a2.size()) throw_range_error(); return return_array_type(a1.size(), make_init_functor(make_array_functor_a_a( fn::functor_each_max(), a1.begin(), a2.begin()))); } template inline small each_max( small const& a1, ElementType const& a2) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a1.size(), make_init_functor(make_array_functor_a_s( fn::functor_each_max(), a1.begin(), a2))); } template inline small each_max( ElementType const& a1, small const& a2) { typedef small return_array_type; typedef typename return_array_type::value_type return_element_type; return return_array_type(a2.size(), make_init_functor(make_array_functor_s_a( fn::functor_each_max(), a1, a2.begin()))); } template inline small real(small, N> const& a) { typedef small return_array_type; return return_array_type(a.size(), make_init_functor(make_array_functor_a( fn::functor_real< ElementType, std::complex >(), a.begin()))); } template inline small imag(small, N> const& a) { typedef small return_array_type; return return_array_type(a.size(), make_init_functor(make_array_functor_a( fn::functor_imag< ElementType, std::complex >(), a.begin()))); } template inline small abs(small, N> const& a) { typedef small return_array_type; return return_array_type(a.size(), make_init_functor(make_array_functor_a( fn::functor_abs< ElementType, std::complex >(), a.begin()))); } template inline small arg(small, N> const& a) { typedef small return_array_type; return return_array_type(a.size(), make_init_functor(make_array_functor_a( fn::functor_arg< ElementType, std::complex >(), a.begin()))); } template