valarray_meta.h

俄罗斯高人Mamaich的Pocket gcc编译器（运行在PocketPC上）的全部源代码。

    //

  template<typename _Dom, typename _Arg>
    class _FunBase 
    {
    public:
      typedef typename _Dom::value_type value_type;
      
      _FunBase(const _Dom& __e, value_type __f(_Arg))
	: _M_expr(__e), _M_func(__f) {}

      value_type operator[](size_t __i) const
      { return _M_func (_M_expr[__i]); }

      size_t size() const { return _M_expr.size ();}

    private:
        const _Dom& _M_expr;
        value_type (*_M_func)(_Arg);
    };

  template<class _Dom>
    struct _ValFunClos<_Expr,_Dom> : _FunBase<_Dom, typename _Dom::value_type> 
    {
      typedef _FunBase<_Dom, typename _Dom::value_type> _Base;
      typedef typename _Base::value_type value_type;
      typedef value_type _Tp;
    
      _ValFunClos(const _Dom& __e, _Tp __f(_Tp)) : _Base(__e, __f) {}
    };

  template<typename _Tp>
    struct _ValFunClos<_ValArray,_Tp> : _FunBase<valarray<_Tp>, _Tp>
    {
      typedef _FunBase<valarray<_Tp>, _Tp> _Base;
      typedef _Tp value_type;
      
      _ValFunClos(const valarray<_Tp>& __v, _Tp __f(_Tp)) : _Base(__v, __f) {}
    };

  template<class _Dom>
    struct _RefFunClos<_Expr,_Dom> :
        _FunBase<_Dom, const typename _Dom::value_type&> 
    {
      typedef _FunBase<_Dom, const typename _Dom::value_type&> _Base;
      typedef typename _Base::value_type value_type;
      typedef value_type _Tp;
      
      _RefFunClos(const _Dom& __e, _Tp __f(const _Tp&))
	: _Base(__e, __f) {}
    };

  template<typename _Tp>
    struct _RefFunClos<_ValArray,_Tp> : _FunBase<valarray<_Tp>, const _Tp&> 
    {
      typedef _FunBase<valarray<_Tp>, const _Tp&> _Base;
      typedef _Tp value_type;
      
      _RefFunClos(const valarray<_Tp>& __v, _Tp __f(const _Tp&))
	: _Base(__v, __f) {}
    };
    
  //
  // Unary expression closure.
  //

  template<class _Oper, class _Arg>
    class _UnBase
    {
    public:
      typedef typename _Arg::value_type _Vt;
      typedef typename __fun<_Oper, _Vt>::result_type value_type;

      _UnBase(const _Arg& __e) : _M_expr(__e) {}

      value_type operator[](size_t __i) const
      { return _Oper()(_M_expr[__i]); }

      size_t size() const { return _M_expr.size(); }

    private:
      const _Arg& _M_expr;
    };

  template<class _Oper, class _Dom>
    struct _UnClos<_Oper, _Expr, _Dom> :  _UnBase<_Oper, _Dom>
    {
      typedef _Dom _Arg;
      typedef _UnBase<_Oper, _Dom> _Base;
      typedef typename _Base::value_type value_type;
      
      _UnClos(const _Arg& __e) : _Base(__e) {}
    };

  template<class _Oper, typename _Tp>
    struct _UnClos<_Oper, _ValArray, _Tp> : _UnBase<_Oper, valarray<_Tp> > 
    {
      typedef valarray<_Tp> _Arg;
      typedef _UnBase<_Oper, valarray<_Tp> > _Base;
      typedef typename _Base::value_type value_type;
      
      _UnClos(const _Arg& __e) : _Base(__e) {}
    };


  //
  // Binary expression closure.
  //

  template<class _Oper, class _FirstArg, class _SecondArg>
    class _BinBase 
    {
    public:
        typedef typename _FirstArg::value_type _Vt;
        typedef typename __fun<_Oper, _Vt>::result_type value_type;

      _BinBase(const _FirstArg& __e1, const _SecondArg& __e2)
	: _M_expr1(__e1), _M_expr2(__e2) {}
      
      value_type operator[](size_t __i) const
      { return _Oper()(_M_expr1[__i], _M_expr2[__i]); }

      size_t size() const { return _M_expr1.size(); }
        
    private:
      const _FirstArg& _M_expr1;
      const _SecondArg& _M_expr2;
    };


  template<class _Oper, class _Clos>
    class _BinBase2
    {
    public:
      typedef typename _Clos::value_type _Vt;
      typedef typename __fun<_Oper, _Vt>::result_type value_type;

      _BinBase2(const _Clos& __e, const _Vt& __t)
	: _M_expr1(__e), _M_expr2(__t) {}

      value_type operator[](size_t __i) const
      { return _Oper()(_M_expr1[__i], _M_expr2); }

      size_t size() const { return _M_expr1.size(); }

    private:
      const _Clos& _M_expr1;
      const _Vt& _M_expr2;
    };

  template<class _Oper, class _Clos>
    class _BinBase1
    {
    public:
      typedef typename _Clos::value_type _Vt;
      typedef typename __fun<_Oper, _Vt>::result_type value_type;

      _BinBase1(const _Vt& __t, const _Clos& __e)
	: _M_expr1(__t), _M_expr2(__e) {}

      value_type operator[](size_t __i) const
      { return _Oper()(_M_expr1, _M_expr2[__i]); }
      
      size_t size() const { return _M_expr2.size(); }

    private:
      const _Vt& _M_expr1;
      const _Clos& _M_expr2;
    };
    
  template<class _Oper, class _Dom1, class _Dom2>
    struct _BinClos<_Oper, _Expr, _Expr, _Dom1, _Dom2>
        : _BinBase<_Oper,_Dom1,_Dom2> 
    {
      typedef _BinBase<_Oper,_Dom1,_Dom2> _Base;
      typedef typename _Base::value_type value_type;
        
      _BinClos(const _Dom1& __e1, const _Dom2& __e2) : _Base(__e1, __e2) {}
    };

  template<class _Oper, typename _Tp>
    struct _BinClos<_Oper,_ValArray,_ValArray,_Tp,_Tp>
      : _BinBase<_Oper,valarray<_Tp>,valarray<_Tp> > 
    {
      typedef _BinBase<_Oper,valarray<_Tp>,valarray<_Tp> > _Base;
      typedef _Tp value_type;

      _BinClos(const valarray<_Tp>& __v, const valarray<_Tp>& __w)
	: _Base(__v, __w) {}
    };

  template<class _Oper, class _Dom>
    struct _BinClos<_Oper,_Expr,_ValArray,_Dom,typename _Dom::value_type>
      : _BinBase<_Oper,_Dom,valarray<typename _Dom::value_type> > 
    {
      typedef typename _Dom::value_type _Tp;
      typedef _BinBase<_Oper,_Dom,valarray<_Tp> > _Base;
      typedef typename _Base::value_type value_type;
      
      _BinClos(const _Dom& __e1, const valarray<_Tp>& __e2)
	: _Base(__e1, __e2) {}
    };

  template<class _Oper, class _Dom>
    struct  _BinClos<_Oper,_ValArray,_Expr,typename _Dom::value_type,_Dom>
      : _BinBase<_Oper,valarray<typename _Dom::value_type>,_Dom> 
    {
      typedef typename _Dom::value_type _Tp;
      typedef _BinBase<_Oper,valarray<_Tp>,_Dom> _Base;
      typedef typename _Base::value_type value_type;
      
      _BinClos(const valarray<_Tp>& __e1, const _Dom& __e2)
	: _Base(__e1, __e2) {}
    };

  template<class _Oper, class _Dom>
    struct _BinClos<_Oper,_Expr,_Constant,_Dom,typename _Dom::value_type>
      : _BinBase2<_Oper,_Dom> 
    {
      typedef typename _Dom::value_type _Tp;
      typedef _BinBase2<_Oper,_Dom> _Base;
      typedef typename _Base::value_type value_type;
      
      _BinClos(const _Dom& __e1, const _Tp& __e2) : _Base(__e1, __e2) {}
    };

  template<class _Oper, class _Dom>
    struct _BinClos<_Oper,_Constant,_Expr,typename _Dom::value_type,_Dom>
      : _BinBase1<_Oper,_Dom> 
    {
      typedef typename _Dom::value_type _Tp;
      typedef _BinBase1<_Oper,_Dom> _Base;
      typedef typename _Base::value_type value_type;
      
      _BinClos(const _Tp& __e1, const _Dom& __e2) : _Base(__e1, __e2) {}
    };
    
  template<class _Oper, typename _Tp>
    struct _BinClos<_Oper,_ValArray,_Constant,_Tp,_Tp>
      : _BinBase2<_Oper,valarray<_Tp> > 
    {
      typedef _BinBase2<_Oper,valarray<_Tp> > _Base;
      typedef typename _Base::value_type value_type;
      
      _BinClos(const valarray<_Tp>& __v, const _Tp& __t) : _Base(__v, __t) {}
    };

  template<class _Oper, typename _Tp>
    struct _BinClos<_Oper,_Constant,_ValArray,_Tp,_Tp>
      : _BinBase1<_Oper,valarray<_Tp> > 
    {
      typedef _BinBase1<_Oper,valarray<_Tp> > _Base;
      typedef typename _Base::value_type value_type;
      
      _BinClos(const _Tp& __t, const valarray<_Tp>& __v) : _Base(__t, __v) {}
    };
        

    //
    // slice_array closure.
    //
    template<typename _Dom>  class _SBase {
    public:
        typedef typename _Dom::value_type value_type;

        _SBase (const _Dom& __e, const slice& __s)
                : _M_expr (__e), _M_slice (__s) {}
        value_type operator[] (size_t __i) const
        { return _M_expr[_M_slice.start () + __i * _M_slice.stride ()]; }
        size_t size() const { return _M_slice.size (); }

    private:
        const _Dom& _M_expr;
        const slice& _M_slice;
    };

    template<typename _Tp> class _SBase<_Array<_Tp> > {
    public:
        typedef _Tp value_type;

        _SBase (_Array<_Tp> __a, const slice& __s)
                : _M_array (__a._M_data+__s.start()), _M_size (__s.size()),
                  _M_stride (__s.stride()) {}
        value_type operator[] (size_t __i) const
        { return _M_array._M_data[__i * _M_stride]; }
        size_t size() const { return _M_size; }

    private:
        const _Array<_Tp> _M_array;
        const size_t _M_size;
        const size_t _M_stride;
    };

    template<class _Dom> struct  _SClos<_Expr,_Dom> : _SBase<_Dom> {
        typedef _SBase<_Dom> _Base;
        typedef typename _Base::value_type value_type;
        
        _SClos (const _Dom& __e, const slice& __s) : _Base (__e, __s) {}
    };

    template<typename _Tp>
    struct _SClos<_ValArray,_Tp> : _SBase<_Array<_Tp> > {
        typedef  _SBase<_Array<_Tp> > _Base;
        typedef _Tp value_type;

        _SClos (_Array<_Tp> __a, const slice& __s) : _Base (__a, __s) {}
    };

    //
    // gslice_array closure.
    //
    template<class _Dom> class _GBase {
    public:
        typedef typename _Dom::value_type value_type;
        
        _GBase (const _Dom& __e, const valarray<size_t>& __i)
                : _M_expr (__e), _M_index(__i) {}
        value_type operator[] (size_t __i) const
        { return _M_expr[_M_index[__i]]; }
        size_t size () const { return _M_index.size(); }
        
    private:
        const _Dom&	 _M_expr;
        const valarray<size_t>& _M_index;
    };
    
    template<typename _Tp> class _GBase<_Array<_Tp> > {
    public:
        typedef _Tp value_type;
        
        _GBase (_Array<_Tp> __a, const valarray<size_t>& __i)
                : _M_array (__a), _M_index(__i) {}
        value_type operator[] (size_t __i) const
        { return _M_array._M_data[_M_index[__i]]; }
        size_t size () const { return _M_index.size(); }
        
    private:
        const _Array<_Tp>     _M_array;
        const valarray<size_t>& _M_index;
    };

    template<class _Dom> struct _GClos<_Expr,_Dom> : _GBase<_Dom> {
        typedef _GBase<_Dom> _Base;
        typedef typename _Base::value_type value_type;

        _GClos (const _Dom& __e, const valarray<size_t>& __i)
                : _Base (__e, __i) {}
    };

    template<typename _Tp>
    struct _GClos<_ValArray,_Tp> : _GBase<_Array<_Tp> > {
        typedef _GBase<_Array<_Tp> > _Base;
        typedef typename _Base::value_type value_type;

        _GClos (_Array<_Tp> __a, const valarray<size_t>& __i)
                : _Base (__a, __i) {}
    };

    //
    // indirect_array closure
    //

    template<class _Dom> class _IBase {
    public:
        typedef typename _Dom::value_type value_type;

        _IBase (const _Dom& __e, const valarray<size_t>& __i)
                : _M_expr (__e), _M_index (__i) {}
        value_type operator[] (size_t __i) const
        { return _M_expr[_M_index[__i]]; }
        size_t size() const { return _M_index.size(); }
        
    private:
        const _Dom& 	    _M_expr;
        const valarray<size_t>& _M_index;
    };

    template<class _Dom> struct _IClos<_Expr,_Dom> : _IBase<_Dom> {
        typedef _IBase<_Dom> _Base;
        typedef typename _Base::value_type value_type;

        _IClos (const _Dom& __e, const valarray<size_t>& __i)
                : _Base (__e, __i) {}