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