/* ***************************************************** THIS IS AN AUTOMATICALLY GENERATED FILE. DO NOT EDIT. ***************************************************** Generated by: scitbx.source_generators.array_family.generate_algebras */ #ifndef SCITBX_ARRAY_FAMILY_SHARED_ALGEBRA_H #define SCITBX_ARRAY_FAMILY_SHARED_ALGEBRA_H #ifndef DOXYGEN_SHOULD_SKIP_THIS #include #include #include #include #include #include namespace scitbx { namespace af { template inline shared operator-(shared const& a) { typedef shared 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 shared operator!(shared const& a) { typedef shared 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 inline shared< typename binary_operator_traits< ElementType1, ElementType2>::arithmetic> operator+( shared const& a1, shared const& a2) { typedef shared< typename binary_operator_traits< ElementType1, ElementType2>::arithmetic> 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 shared operator+( shared const& a1, ElementType const& a2) { typedef shared 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 shared operator+( ElementType const& a1, shared const& a2) { typedef shared 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 inline shared& operator+=( shared& a1, shared 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 shared& operator+=( shared& a1, ElementType const& a2) { array_operation_in_place_a_s(fn::functor_ip_plus< ElementType, ElementType>(), a1.begin(), a2, a1.size()); return a1; } template inline shared< typename binary_operator_traits< ElementType1, ElementType2>::arithmetic> operator-( shared const& a1, shared const& a2) { typedef shared< typename binary_operator_traits< ElementType1, ElementType2>::arithmetic> 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 shared operator-( shared const& a1, ElementType const& a2) { typedef shared 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 shared operator-( ElementType const& a1, shared const& a2) { typedef shared 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 inline shared& operator-=( shared& a1, shared 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 shared& operator-=( shared& a1, ElementType const& a2) { array_operation_in_place_a_s(fn::functor_ip_minus< ElementType, ElementType>(), a1.begin(), a2, a1.size()); return a1; } template inline shared< typename binary_operator_traits< ElementType1, ElementType2>::arithmetic> operator*( shared const& a1, shared const& a2) { typedef shared< typename binary_operator_traits< ElementType1, ElementType2>::arithmetic> 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 shared operator*( shared const& a1, ElementType const& a2) { typedef shared 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 shared operator*( ElementType const& a1, shared const& a2) { typedef shared 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 inline shared& operator*=( shared& a1, shared 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 shared& operator*=( shared& a1, ElementType const& a2) { array_operation_in_place_a_s(fn::functor_ip_multiplies< ElementType, ElementType>(), a1.begin(), a2, a1.size()); return a1; } template inline shared< typename binary_operator_traits< ElementType1, ElementType2>::arithmetic> operator/( shared const& a1, shared const& a2) { typedef shared< typename binary_operator_traits< ElementType1, ElementType2>::arithmetic> 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 shared operator/( shared const& a1, ElementType const& a2) { typedef shared 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 shared operator/( ElementType const& a1, shared const& a2) { typedef shared 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 inline shared& operator/=( shared& a1, shared 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 shared& operator/=( shared& a1, ElementType const& a2) { array_operation_in_place_a_s(fn::functor_ip_divides< ElementType, ElementType>(), a1.begin(), a2, a1.size()); return a1; } template inline shared< typename binary_operator_traits< ElementType1, ElementType2>::arithmetic> operator%( shared const& a1, shared const& a2) { typedef shared< typename binary_operator_traits< ElementType1, ElementType2>::arithmetic> 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 shared operator%( shared const& a1, ElementType const& a2) { typedef shared 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 shared operator%( ElementType const& a1, shared const& a2) { typedef shared 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 inline shared& operator%=( shared& a1, shared 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 shared& operator%=( shared& a1, ElementType const& a2) { array_operation_in_place_a_s(fn::functor_ip_modulus< ElementType, ElementType>(), a1.begin(), a2, a1.size()); return a1; } template inline shared operator&&( shared const& a1, shared const& a2) { typedef shared 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 shared operator&&( shared const& a1, ElementType const& a2) { typedef shared 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 shared operator&&( ElementType const& a1, shared const& a2) { typedef shared 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 inline shared operator||( shared const& a1, shared const& a2) { typedef shared 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 shared operator||( shared const& a1, ElementType const& a2) { typedef shared 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 shared operator||( ElementType const& a1, shared const& a2) { typedef shared 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 inline shared operator==( shared const& a1, shared const& a2) { typedef shared 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 shared operator==( shared const& a1, ElementType const& a2) { typedef shared 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 shared operator==( ElementType const& a1, shared const& a2) { typedef shared 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 inline shared operator!=( shared const& a1, shared const& a2) { typedef shared 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 shared operator!=( shared const& a1, ElementType const& a2) { typedef shared 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 shared operator!=( ElementType const& a1, shared const& a2) { typedef shared 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 inline shared operator>( shared const& a1, shared const& a2) { typedef shared 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 shared operator>( shared const& a1, ElementType const& a2) { typedef shared 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 shared operator>( ElementType const& a1, shared const& a2) { typedef shared 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 inline shared operator<( shared const& a1, shared const& a2) { typedef shared 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 shared operator<( shared const& a1, ElementType const& a2) { typedef shared 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 shared operator<( ElementType const& a1, shared const& a2) { typedef shared 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 inline shared operator>=( shared const& a1, shared const& a2) { typedef shared 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 shared operator>=( shared const& a1, ElementType const& a2) { typedef shared 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 shared operator>=( ElementType const& a1, shared const& a2) { typedef shared 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 inline shared operator<=( shared const& a1, shared const& a2) { typedef shared 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 shared operator<=( shared const& a1, ElementType const& a2) { typedef shared 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 shared operator<=( ElementType const& a1, shared const& a2) { typedef shared 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 shared absolute(shared const& a) { typedef shared 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 shared pow2(shared const& a) { typedef shared 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 shared acos(shared const& a) { typedef shared 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 shared cos(shared const& a) { typedef shared 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 shared tan(shared const& a) { typedef shared 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 shared asin(shared const& a) { typedef shared 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 shared cosh(shared const& a) { typedef shared 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 shared tanh(shared const& a) { typedef shared 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 shared atan(shared const& a) { typedef shared 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 shared exp(shared const& a) { typedef shared 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 shared sin(shared const& a) { typedef shared 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 shared fabs(shared const& a) { typedef shared 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 shared log(shared const& a) { typedef shared 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 shared sinh(shared const& a) { typedef shared 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 shared ceil(shared const& a) { typedef shared 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 shared floor(shared const& a) { typedef shared 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 shared log10(shared const& a) { typedef shared 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 shared sqrt(shared const& a) { typedef shared 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 shared abs(shared const& a) { typedef shared 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 shared conj(shared const& a) { typedef shared 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 inline shared fmod_positive( shared const& a1, shared const& a2) { typedef shared 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 shared fmod_positive( shared const& a1, ElementType const& a2) { typedef shared 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 shared fmod_positive( ElementType const& a1, shared const& a2) { typedef shared 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 inline shared fmod( shared const& a1, shared const& a2) { typedef shared 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 shared fmod( shared const& a1, ElementType const& a2) { typedef shared 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 shared fmod( ElementType const& a1, shared const& a2) { typedef shared 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 inline shared pow( shared const& a1, shared const& a2) { typedef shared 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 shared pow( shared const& a1, ElementType const& a2) { typedef shared 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 shared pow( ElementType const& a1, shared const& a2) { typedef shared 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 inline shared atan2( shared const& a1, shared const& a2) { typedef shared 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 shared atan2( shared const& a1, ElementType const& a2) { typedef shared 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 shared atan2( ElementType const& a1, shared const& a2) { typedef shared 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 inline shared each_min( shared const& a1, shared const& a2) { typedef shared 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 shared each_min( shared const& a1, ElementType const& a2) { typedef shared 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 shared each_min( ElementType const& a1, shared const& a2) { typedef shared 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 inline shared each_max( shared const& a1, shared const& a2) { typedef shared 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 shared each_max( shared const& a1, ElementType const& a2) { typedef shared 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 shared