old-autovect-branch/libstdc++-v3/include/bits/valarray_after.h

   1 // The template and inlines for the -*- C++ -*- internal _Meta class.
   2
   3 // Copyright (C) 1997, 1998, 1999, 2000, 2001, 2002, 2003, 2004
   4 // Free Software Foundation, Inc.
   5 //
   6 // This file is part of the GNU ISO C++ Library.  This library is free
   7 // software; you can redistribute it and/or modify it under the
   8 // terms of the GNU General Public License as published by the
   9 // Free Software Foundation; either version 2, or (at your option)
  10 // any later version.
  11
  12 // This library is distributed in the hope that it will be useful,
  13 // but WITHOUT ANY WARRANTY; without even the implied warranty of
  14 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  15 // GNU General Public License for more details.
  16
  17 // You should have received a copy of the GNU General Public License along
  18 // with this library; see the file COPYING.  If not, write to the Free
  19 // Software Foundation, 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301,
  20 // USA.
  21
  22 // As a special exception, you may use this file as part of a free software
  23 // library without restriction.  Specifically, if other files instantiate
  24 // templates or use macros or inline functions from this file, or you compile
  25 // this file and link it with other files to produce an executable, this
  26 // file does not by itself cause the resulting executable to be covered by
  27 // the GNU General Public License.  This exception does not however
  28 // invalidate any other reasons why the executable file might be covered by
  29 // the GNU General Public License.
  30
  31 // Written by Gabriel Dos Reis <Gabriel.Dos-Reis@cmla.ens-cachan.fr>
  32
  33 /** @file valarray_after.h
  34  *  This is an internal header file, included by other library headers.
  35  *  You should not attempt to use it directly.
  36  */
  37
  38 #ifndef _VALARRAY_AFTER_H
  39 #define _VALARRAY_AFTER_H 1
  40
  41 #pragma GCC system_header
  42
  43 namespace std
  44 {
  45   //
  46   // gslice_array closure.
  47   //
  48   template<class _Dom>
  49     class _GBase
  50     {
  51     public:
  52       typedef typename _Dom::value_type value_type;
  53
  54       _GBase (const _Dom& __e, const valarray<size_t>& __i)
  55       : _M_expr (__e), _M_index(__i) {}
  56
  57       value_type
  58       operator[] (size_t __i) const
  59       { return _M_expr[_M_index[__i]]; }
  60
  61       size_t
  62       size () const
  63       { return _M_index.size(); }
  64
  65     private:
  66       const _Dom&             _M_expr;
  67       const valarray<size_t>& _M_index;
  68     };
  69
  70   template<typename _Tp>
  71     class _GBase<_Array<_Tp> >
  72     {
  73     public:
  74       typedef _Tp value_type;
  75
  76       _GBase (_Array<_Tp> __a, const valarray<size_t>& __i)
  77       : _M_array (__a), _M_index(__i) {}
  78
  79       value_type
  80       operator[] (size_t __i) const
  81       { return _M_array._M_data[_M_index[__i]]; }
  82
  83       size_t
  84       size () const
  85       { return _M_index.size(); }
  86
  87     private:
  88       const _Array<_Tp>       _M_array;
  89       const valarray<size_t>& _M_index;
  90     };
  91
  92   template<class _Dom>
  93     struct _GClos<_Expr, _Dom>
  94     : _GBase<_Dom>
  95     {
  96       typedef _GBase<_Dom> _Base;
  97       typedef typename _Base::value_type value_type;
  98
  99       _GClos (const _Dom& __e, const valarray<size_t>& __i)
 100       : _Base (__e, __i) {}
 101     };
 102
 103   template<typename _Tp>
 104     struct _GClos<_ValArray, _Tp>
 105     : _GBase<_Array<_Tp> >
 106     {
 107       typedef _GBase<_Array<_Tp> > _Base;
 108       typedef typename _Base::value_type value_type;
 109
 110       _GClos (_Array<_Tp> __a, const valarray<size_t>& __i)
 111       : _Base (__a, __i) {}
 112     };
 113
 114   //
 115   // indirect_array closure
 116   //
 117   template<class _Dom>
 118     class _IBase
 119     {
 120     public:
 121       typedef typename _Dom::value_type value_type;
 122
 123       _IBase (const _Dom& __e, const valarray<size_t>& __i)
 124       : _M_expr (__e), _M_index (__i) {}
 125
 126       value_type
 127       operator[] (size_t __i) const
 128       { return _M_expr[_M_index[__i]]; }
 129
 130       size_t
 131       size() const
 132       { return _M_index.size(); }
 133
 134     private:
 135       const _Dom&             _M_expr;
 136       const valarray<size_t>& _M_index;
 137     };
 138
 139   template<class _Dom>
 140     struct _IClos<_Expr, _Dom>
 141     : _IBase<_Dom>
 142     {
 143       typedef _IBase<_Dom> _Base;
 144       typedef typename _Base::value_type value_type;
 145
 146       _IClos (const _Dom& __e, const valarray<size_t>& __i)
 147       : _Base (__e, __i) {}
 148     };
 149
 150   template<typename _Tp>
 151     struct _IClos<_ValArray, _Tp>
 152     : _IBase<valarray<_Tp> >
 153     {
 154       typedef _IBase<valarray<_Tp> > _Base;
 155       typedef _Tp value_type;
 156
 157       _IClos (const valarray<_Tp>& __a, const valarray<size_t>& __i)
 158       : _Base (__a, __i) {}
 159     };
 160
 161   //
 162   // class _Expr
 163   //
 164   template<class _Clos, typename _Tp>
 165     class _Expr
 166     {
 167     public:
 168       typedef _Tp value_type;
 169
 170       _Expr(const _Clos&);
 171
 172       const _Clos& operator()() const;
 173
 174       value_type operator[](size_t) const;
 175       valarray<value_type> operator[](slice) const;
 176       valarray<value_type> operator[](const gslice&) const;
 177       valarray<value_type> operator[](const valarray<bool>&) const;
 178       valarray<value_type> operator[](const valarray<size_t>&) const;
 179
 180       _Expr<_UnClos<__unary_plus, std::_Expr, _Clos>, value_type>
 181       operator+() const;
 182
 183       _Expr<_UnClos<__negate, std::_Expr, _Clos>, value_type>
 184       operator-() const;
 185
 186       _Expr<_UnClos<__bitwise_not, std::_Expr, _Clos>, value_type>
 187       operator~() const;
 188
 189       _Expr<_UnClos<__logical_not, std::_Expr, _Clos>, bool>
 190       operator!() const;
 191
 192       size_t size() const;
 193       value_type sum() const;
 194
 195       valarray<value_type> shift(int) const;
 196       valarray<value_type> cshift(int) const;
 197
 198       value_type min() const;
 199       value_type max() const;
 200
 201       valarray<value_type> apply(value_type (*)(const value_type&)) const;
 202       valarray<value_type> apply(value_type (*)(value_type)) const;
 203
 204     private:
 205       const _Clos _M_closure;
 206     };
 207
 208   template<class _Clos, typename _Tp>
 209     inline
 210     _Expr<_Clos, _Tp>::_Expr(const _Clos& __c) : _M_closure(__c) {}
 211
 212   template<class _Clos, typename _Tp>
 213     inline const _Clos&
 214     _Expr<_Clos, _Tp>::operator()() const
 215     { return _M_closure; }
 216
 217   template<class _Clos, typename _Tp>
 218     inline _Tp
 219     _Expr<_Clos, _Tp>::operator[](size_t __i) const
 220     { return _M_closure[__i]; }
 221
 222   template<class _Clos, typename _Tp>
 223     inline valarray<_Tp>
 224     _Expr<_Clos, _Tp>::operator[](slice __s) const
 225     {
 226       valarray<_Tp> __v = valarray<_Tp>(*this)[__s];
 227       return __v;
 228     }
 229
 230   template<class _Clos, typename _Tp>
 231     inline valarray<_Tp>
 232     _Expr<_Clos, _Tp>::operator[](const gslice& __gs) const
 233     {
 234       valarray<_Tp> __v = valarray<_Tp>(*this)[__gs];
 235       return __v;
 236     }
 237
 238   template<class _Clos, typename _Tp>
 239     inline valarray<_Tp>
 240     _Expr<_Clos, _Tp>::operator[](const valarray<bool>& __m) const
 241     {
 242       valarray<_Tp> __v = valarray<_Tp>(*this)[__m];
 243       return __v;
 244     }
 245
 246   template<class _Clos, typename _Tp>
 247     inline valarray<_Tp>
 248     _Expr<_Clos, _Tp>::operator[](const valarray<size_t>& __i) const
 249     {
 250       valarray<_Tp> __v = valarray<_Tp>(*this)[__i];
 251       return __v;
 252     }
 253
 254   template<class _Clos, typename _Tp>
 255     inline size_t
 256     _Expr<_Clos, _Tp>::size() const
 257     { return _M_closure.size(); }
 258
 259   template<class _Clos, typename _Tp>
 260     inline valarray<_Tp>
 261     _Expr<_Clos, _Tp>::shift(int __n) const
 262     {
 263       valarray<_Tp> __v = valarray<_Tp>(*this).shift(__n);
 264       return __v;
 265     }
 266
 267   template<class _Clos, typename _Tp>
 268     inline valarray<_Tp>
 269     _Expr<_Clos, _Tp>::cshift(int __n) const
 270     {
 271       valarray<_Tp> __v = valarray<_Tp>(*this).cshift(__n);
 272       return __v;
 273     }
 274
 275   template<class _Clos, typename _Tp>
 276     inline valarray<_Tp>
 277     _Expr<_Clos, _Tp>::apply(_Tp __f(const _Tp&)) const
 278     {
 279       valarray<_Tp> __v = valarray<_Tp>(*this).apply(__f);
 280       return __v;
 281     }
 282
 283   template<class _Clos, typename _Tp>
 284     inline valarray<_Tp>
 285     _Expr<_Clos, _Tp>::apply(_Tp __f(_Tp)) const
 286     {
 287       valarray<_Tp> __v = valarray<_Tp>(*this).apply(__f);
 288       return __v;
 289     }
 290
 291   // XXX: replace this with a more robust summation algorithm.
 292   template<class _Clos, typename _Tp>
 293     inline _Tp
 294     _Expr<_Clos, _Tp>::sum() const
 295     {
 296       size_t __n = _M_closure.size();
 297       if (__n == 0)
 298         return _Tp();
 299       else
 300         {
 301           _Tp __s = _M_closure[--__n];
 302           while (__n != 0)
 303             __s += _M_closure[--__n];
 304           return __s;
 305         }
 306     }
 307
 308   template<class _Clos, typename _Tp>
 309     inline _Tp
 310     _Expr<_Clos, _Tp>::min() const
 311     { return __valarray_min(_M_closure); }
 312
 313   template<class _Clos, typename _Tp>
 314     inline _Tp
 315     _Expr<_Clos, _Tp>::max() const
 316     { return __valarray_max(_M_closure); }
 317
 318   template<class _Dom, typename _Tp>
 319     inline _Expr<_UnClos<__logical_not, _Expr, _Dom>, bool>
 320     _Expr<_Dom, _Tp>::operator!() const
 321     {
 322       typedef _UnClos<__logical_not, std::_Expr, _Dom> _Closure;
 323       return _Expr<_Closure, _Tp>(_Closure(this->_M_closure));
 324     }
 325
 326 #define _DEFINE_EXPR_UNARY_OPERATOR(_Op, _Name)                           \
 327   template<class _Dom, typename _Tp>                                      \
 328     inline _Expr<_UnClos<_Name, std::_Expr, _Dom>, _Tp>                   \
 329     _Expr<_Dom, _Tp>::operator _Op() const                                \
 330     {                                                                     \
 331       typedef _UnClos<_Name, std::_Expr, _Dom> _Closure;                  \
 332       return _Expr<_Closure, _Tp>(_Closure(this->_M_closure));            \
 333     }
 334
 335     _DEFINE_EXPR_UNARY_OPERATOR(+, __unary_plus)
 336     _DEFINE_EXPR_UNARY_OPERATOR(-, __negate)
 337     _DEFINE_EXPR_UNARY_OPERATOR(~, __bitwise_not)
 338
 339 #undef _DEFINE_EXPR_UNARY_OPERATOR
 340
 341 #define _DEFINE_EXPR_BINARY_OPERATOR(_Op, _Name)                        \
 342   template<class _Dom1, class _Dom2>                                    \
 343     inline _Expr<_BinClos<_Name, _Expr, _Expr, _Dom1, _Dom2>,           \
 344            typename __fun<_Name, typename _Dom1::value_type>::result_type> \
 345     operator _Op(const _Expr<_Dom1, typename _Dom1::value_type>& __v,   \
 346                  const _Expr<_Dom2, typename _Dom2::value_type>& __w)   \
 347     {                                                                   \
 348       typedef typename _Dom1::value_type _Arg;                          \
 349       typedef typename __fun<_Name, _Arg>::result_type _Value;          \
 350       typedef _BinClos<_Name, _Expr, _Expr, _Dom1, _Dom2> _Closure;     \
 351       return _Expr<_Closure, _Value>(_Closure(__v(), __w()));           \
 352     }                                                                   \
 353                                                                         \
 354   template<class _Dom>                                                  \
 355     inline _Expr<_BinClos<_Name, _Expr, _Constant, _Dom,                \
 356                           typename _Dom::value_type>,                   \
 357              typename __fun<_Name, typename _Dom::value_type>::result_type> \
 358     operator _Op(const _Expr<_Dom, typename _Dom::value_type>& __v,     \
 359                  const typename _Dom::value_type& __t)                  \
 360     {                                                                   \
 361       typedef typename _Dom::value_type _Arg;                           \
 362       typedef typename __fun<_Name, _Arg>::result_type _Value;          \
 363       typedef _BinClos<_Name, _Expr, _Constant, _Dom, _Arg> _Closure;   \
 364       return _Expr<_Closure, _Value>(_Closure(__v(), __t));             \
 365     }                                                                   \
 366                                                                         \
 367   template<class _Dom>                                                  \
 368     inline _Expr<_BinClos<_Name, _Constant, _Expr,                      \
 369                           typename _Dom::value_type, _Dom>,             \
 370              typename __fun<_Name, typename _Dom::value_type>::result_type> \
 371     operator _Op(const typename _Dom::value_type& __t,                  \
 372                  const _Expr<_Dom, typename _Dom::value_type>& __v)     \
 373     {                                                                   \
 374       typedef typename _Dom::value_type _Arg;                           \
 375       typedef typename __fun<_Name, _Arg>::result_type _Value;          \
 376       typedef _BinClos<_Name, _Constant, _Expr, _Arg, _Dom> _Closure;   \
 377       return _Expr<_Closure, _Value>(_Closure(__t, __v()));             \
 378     }                                                                   \
 379                                                                         \
 380   template<class _Dom>                                                  \
 381     inline _Expr<_BinClos<_Name, _Expr, _ValArray,                      \
 382                           _Dom, typename _Dom::value_type>,             \
 383              typename __fun<_Name, typename _Dom::value_type>::result_type> \
 384     operator _Op(const _Expr<_Dom,typename _Dom::value_type>& __e,      \
 385                  const valarray<typename _Dom::value_type>& __v)        \
 386     {                                                                   \
 387       typedef typename _Dom::value_type _Arg;                           \
 388       typedef typename __fun<_Name, _Arg>::result_type _Value;          \
 389       typedef _BinClos<_Name, _Expr, _ValArray, _Dom, _Arg> _Closure;   \
 390       return _Expr<_Closure, _Value>(_Closure(__e(), __v));             \
 391     }                                                                   \
 392                                                                         \
 393   template<class _Dom>                                                  \
 394     inline _Expr<_BinClos<_Name, _ValArray, _Expr,                      \
 395                  typename _Dom::value_type, _Dom>,                      \
 396              typename __fun<_Name, typename _Dom::value_type>::result_type> \
 397     operator _Op(const valarray<typename _Dom::value_type>& __v,        \
 398                  const _Expr<_Dom, typename _Dom::value_type>& __e)     \
 399     {                                                                   \
 400       typedef typename _Dom::value_type _Tp;                            \
 401       typedef typename __fun<_Name, _Tp>::result_type _Value;           \
 402       typedef _BinClos<_Name, _ValArray, _Expr, _Tp, _Dom> _Closure;    \
 403       return _Expr<_Closure, _Value>(_Closure(__v, __e ()));            \
 404     }
 405
 406     _DEFINE_EXPR_BINARY_OPERATOR(+, __plus)
 407     _DEFINE_EXPR_BINARY_OPERATOR(-, __minus)
 408     _DEFINE_EXPR_BINARY_OPERATOR(*, __multiplies)
 409     _DEFINE_EXPR_BINARY_OPERATOR(/, __divides)
 410     _DEFINE_EXPR_BINARY_OPERATOR(%, __modulus)
 411     _DEFINE_EXPR_BINARY_OPERATOR(^, __bitwise_xor)
 412     _DEFINE_EXPR_BINARY_OPERATOR(&, __bitwise_and)
 413     _DEFINE_EXPR_BINARY_OPERATOR(|, __bitwise_or)
 414     _DEFINE_EXPR_BINARY_OPERATOR(<<, __shift_left)
 415     _DEFINE_EXPR_BINARY_OPERATOR(>>, __shift_right)
 416     _DEFINE_EXPR_BINARY_OPERATOR(&&, __logical_and)
 417     _DEFINE_EXPR_BINARY_OPERATOR(||, __logical_or)
 418     _DEFINE_EXPR_BINARY_OPERATOR(==, __equal_to)
 419     _DEFINE_EXPR_BINARY_OPERATOR(!=, __not_equal_to)
 420     _DEFINE_EXPR_BINARY_OPERATOR(<, __less)
 421     _DEFINE_EXPR_BINARY_OPERATOR(>, __greater)
 422     _DEFINE_EXPR_BINARY_OPERATOR(<=, __less_equal)
 423     _DEFINE_EXPR_BINARY_OPERATOR(>=, __greater_equal)
 424
 425 #undef _DEFINE_EXPR_BINARY_OPERATOR
 426
 427 #define _DEFINE_EXPR_UNARY_FUNCTION(_Name)                               \
 428   template<class _Dom>                                                   \
 429     inline _Expr<_UnClos<__##_Name, _Expr, _Dom>,                        \
 430                  typename _Dom::value_type>                              \
 431     _Name(const _Expr<_Dom, typename _Dom::value_type>& __e)             \
 432     {                                                                    \
 433       typedef typename _Dom::value_type _Tp;                             \
 434       typedef _UnClos<__##_Name, _Expr, _Dom> _Closure;                  \
 435       return _Expr<_Closure, _Tp>(_Closure(__e()));                      \
 436     }                                                                    \
 437                                                                          \
 438   template<typename _Tp>                                                 \
 439     inline _Expr<_UnClos<__##_Name, _ValArray, _Tp>, _Tp>                \
 440     _Name(const valarray<_Tp>& __v)                                      \
 441     {                                                                    \
 442       typedef _UnClos<__##_Name, _ValArray, _Tp> _Closure;               \
 443       return _Expr<_Closure, _Tp>(_Closure(__v));                        \
 444     }
 445
 446     _DEFINE_EXPR_UNARY_FUNCTION(abs)
 447     _DEFINE_EXPR_UNARY_FUNCTION(cos)
 448     _DEFINE_EXPR_UNARY_FUNCTION(acos)
 449     _DEFINE_EXPR_UNARY_FUNCTION(cosh)
 450     _DEFINE_EXPR_UNARY_FUNCTION(sin)
 451     _DEFINE_EXPR_UNARY_FUNCTION(asin)
 452     _DEFINE_EXPR_UNARY_FUNCTION(sinh)
 453     _DEFINE_EXPR_UNARY_FUNCTION(tan)
 454     _DEFINE_EXPR_UNARY_FUNCTION(tanh)
 455     _DEFINE_EXPR_UNARY_FUNCTION(atan)
 456     _DEFINE_EXPR_UNARY_FUNCTION(exp)
 457     _DEFINE_EXPR_UNARY_FUNCTION(log)
 458     _DEFINE_EXPR_UNARY_FUNCTION(log10)
 459     _DEFINE_EXPR_UNARY_FUNCTION(sqrt)
 460
 461 #undef _DEFINE_EXPR_UNARY_FUNCTION
 462
 463 #define _DEFINE_EXPR_BINARY_FUNCTION(_Fun)                             \
 464   template<class _Dom1, class _Dom2>                                   \
 465     inline _Expr<_BinClos<__##_Fun, _Expr, _Expr, _Dom1, _Dom2>,       \
 466                  typename _Dom1::value_type>                           \
 467     _Fun(const _Expr<_Dom1, typename _Dom1::value_type>& __e1,         \
 468           const _Expr<_Dom2, typename _Dom2::value_type>& __e2)        \
 469     {                                                                  \
 470       typedef typename _Dom1::value_type _Tp;                          \
 471       typedef _BinClos<__##_Fun, _Expr, _Expr, _Dom1, _Dom2> _Closure; \
 472       return _Expr<_Closure, _Tp>(_Closure(__e1(), __e2()));           \
 473     }                                                                  \
 474                                                                        \
 475   template<class _Dom>                                                 \
 476     inline _Expr<_BinClos<__##_Fun, _Expr, _ValArray, _Dom,            \
 477                           typename _Dom::value_type>,                  \
 478                  typename _Dom::value_type>                            \
 479     _Fun(const _Expr<_Dom, typename _Dom::value_type>& __e,            \
 480          const valarray<typename _Dom::value_type>& __v)               \
 481     {                                                                  \
 482       typedef typename _Dom::value_type _Tp;                           \
 483       typedef _BinClos<__##_Fun, _Expr, _ValArray, _Dom, _Tp> _Closure; \
 484       return _Expr<_Closure, _Tp>(_Closure(__e(), __v));               \
 485     }                                                                  \
 486                                                                        \
 487   template<class _Dom>                                                 \
 488     inline _Expr<_BinClos<__##_Fun, _ValArray, _Expr,                  \
 489                           typename _Dom::value_type, _Dom>,            \
 490                  typename _Dom::value_type>                            \
 491     _Fun(const valarray<typename _Dom::valarray>& __v,                 \
 492          const _Expr<_Dom, typename _Dom::value_type>& __e)            \
 493     {                                                                  \
 494       typedef typename _Dom::value_type _Tp;                           \
 495       typedef _BinClos<__##_Fun, _ValArray, _Expr, _Tp, _Dom> _Closure; \
 496       return _Expr<_Closure, _Tp>(_Closure(__v, __e()));               \
 497     }                                                                  \
 498                                                                        \
 499   template<class _Dom>                                                 \
 500     inline _Expr<_BinClos<__##_Fun, _Expr, _Constant, _Dom,            \
 501                           typename _Dom::value_type>,                  \
 502                  typename _Dom::value_type>                            \
 503     _Fun(const _Expr<_Dom, typename _Dom::value_type>& __e,            \
 504          const typename _Dom::value_type& __t)                         \
 505     {                                                                  \
 506       typedef typename _Dom::value_type _Tp;                           \
 507       typedef _BinClos<__##_Fun, _Expr, _Constant, _Dom, _Tp> _Closure;\
 508       return _Expr<_Closure, _Tp>(_Closure(__e(), __t));               \
 509     }                                                                  \
 510                                                                        \
 511   template<class _Dom>                                                 \
 512     inline _Expr<_BinClos<__##_Fun, _Constant, _Expr,                  \
 513                           typename _Dom::value_type, _Dom>,            \
 514                  typename _Dom::value_type>                            \
 515     _Fun(const typename _Dom::value_type& __t,                         \
 516          const _Expr<_Dom, typename _Dom::value_type>& __e)            \
 517     {                                                                  \
 518       typedef typename _Dom::value_type _Tp;                           \
 519       typedef _BinClos<__##_Fun, _Constant, _Expr, _Tp, _Dom> _Closure; \
 520       return _Expr<_Closure, _Tp>(_Closure(__t, __e()));               \
 521     }                                                                  \
 522                                                                        \
 523   template<typename _Tp>                                               \
 524     inline _Expr<_BinClos<__##_Fun, _ValArray, _ValArray, _Tp, _Tp>, _Tp> \
 525     _Fun(const valarray<_Tp>& __v, const valarray<_Tp>& __w)           \
 526     {                                                                  \
 527       typedef _BinClos<__##_Fun, _ValArray, _ValArray, _Tp, _Tp> _Closure; \
 528       return _Expr<_Closure, _Tp>(_Closure(__v, __w));                 \
 529     }                                                                  \
 530                                                                        \
 531   template<typename _Tp>                                               \
 532     inline _Expr<_BinClos<__##_Fun, _ValArray, _Constant, _Tp, _Tp>, _Tp> \
 533     _Fun(const valarray<_Tp>& __v, const _Tp& __t)                     \
 534     {                                                                  \
 535       typedef _BinClos<__##_Fun, _ValArray, _Constant, _Tp, _Tp> _Closure; \
 536       return _Expr<_Closure, _Tp>(_Closure(__v, __t));                 \
 537     }                                                                  \
 538                                                                        \
 539   template<typename _Tp>                                               \
 540     inline _Expr<_BinClos<__##_Fun, _Constant, _ValArray, _Tp, _Tp>, _Tp> \
 541     _Fun(const _Tp& __t, const valarray<_Tp>& __v)                     \
 542     {                                                                  \
 543       typedef _BinClos<__##_Fun, _Constant, _ValArray, _Tp, _Tp> _Closure; \
 544       return _Expr<_Closure, _Tp>(_Closure(__t, __v));                 \
 545     }
 546
 547 _DEFINE_EXPR_BINARY_FUNCTION(atan2)
 548 _DEFINE_EXPR_BINARY_FUNCTION(pow)
 549
 550 #undef _DEFINE_EXPR_BINARY_FUNCTION
 551
 552 } // std::
 553
 554 #endif /* _CPP_VALARRAY_AFTER_H */
 555
 556 // Local Variables:
 557 // mode:c++
 558 // End: