// The template and inlines for the -*- C++ -*- gslice_array class.

// Copyright (C) 1997-1999, 2000 Free Software Foundation, Inc.
//
// This file is part of the GNU ISO C++ Library.  This library is free
// software; you can redistribute it and/or modify it under the
// terms of the GNU General Public License as published by the
// Free Software Foundation; either version 2, or (at your option)
// any later version.

// This library is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
// GNU General Public License for more details.

// You should have received a copy of the GNU General Public License along
// with this library; see the file COPYING.  If not, write to the Free
// Software Foundation, 59 Temple Place - Suite 330, Boston, MA 02111-1307,
// USA.

// As a special exception, you may use this file as part of a free software
// library without restriction.  Specifically, if other files instantiate
// templates or use macros or inline functions from this file, or you compile
// this file and link it with other files to produce an executable, this
// file does not by itself cause the resulting executable to be covered by
// the GNU General Public License.  This exception does not however
// invalidate any other reasons why the executable file might be covered by
// the GNU General Public License.

// Written by Gabriel Dos Reis <Gabriel.Dos-Reis@DPTMaths.ENS-Cachan.Fr>

#ifndef _CPP_BITS_GSLICE_ARRAY
#define _CPP_BITS_GSLICE_ARRAY 1

namespace std {

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

        void operator=  (const valarray<_Tp>&) const;
        void operator*= (const valarray<_Tp>&) const;
        void operator/= (const valarray<_Tp>&) const;
        void operator%= (const valarray<_Tp>&) const;
        void operator+= (const valarray<_Tp>&) const;
        void operator-= (const valarray<_Tp>&) const;
        void operator^= (const valarray<_Tp>&) const;
        void operator&= (const valarray<_Tp>&) const;
        void operator|= (const valarray<_Tp>&) const;
        void operator<<=(const valarray<_Tp>&) const;
        void operator>>=(const valarray<_Tp>&) const;
        void operator=(const _Tp&);

        template<class _Dom>
        void operator= (const _Expr<_Dom,_Tp>&) const;
        template<class _Dom>
        void operator*= (const _Expr<_Dom,_Tp>&) const;
        template<class _Dom>
        void operator/= (const _Expr<_Dom,_Tp>&) const;
        template<class _Dom>
        void operator%= (const _Expr<_Dom,_Tp>&) const;
        template<class _Dom>
        void operator+= (const _Expr<_Dom,_Tp>&) const;
        template<class _Dom>
        void operator-= (const _Expr<_Dom,_Tp>&) const;
        template<class _Dom>
        void operator^= (const _Expr<_Dom,_Tp>&) const;
        template<class _Dom>
        void operator&= (const _Expr<_Dom,_Tp>&) const;
        template<class _Dom>
        void operator|= (const _Expr<_Dom,_Tp>&) const;
        template<class _Dom>
        void operator<<= (const _Expr<_Dom,_Tp>&) const;
        template<class _Dom>
        void operator>>= (const _Expr<_Dom,_Tp>&) const;
        
    private:
        _Array<_Tp>    _M_array;
        const valarray<size_t>& _M_index;
        
        friend class valarray<_Tp>;
        
        gslice_array (_Array<_Tp>, const valarray<size_t>&);

        // this constructor needs to be implemented.
        gslice_array (const gslice_array&);

        // not implemented
        gslice_array();
        gslice_array& operator= (const gslice_array&);
    };

    template<typename _Tp>
    inline
    gslice_array<_Tp>::gslice_array (_Array<_Tp> __a,
                                     const valarray<size_t>& __i)
            : _M_array (__a), _M_index (__i) {}


    template<typename _Tp>
    inline
    gslice_array<_Tp>::gslice_array (const gslice_array<_Tp>& __a)
            : _M_array (__a._M_array), _M_index (__a._M_index) {}

    
    template<typename _Tp>
    inline void
    gslice_array<_Tp>::operator= (const _Tp& __t) 
    { 
        __valarray_fill (_M_array, _Array<size_t>(_M_index),
                         _M_index.size(), __t); 
    }

    template<typename _Tp>
    inline void
    gslice_array<_Tp>::operator= (const valarray<_Tp>& __v) const
    {
        __valarray_copy (_Array<_Tp> (__v), __v.size (),
                         _M_array, _Array<size_t>(_M_index));
    }

    template<typename _Tp>
    template<class E>
    inline void
    gslice_array<_Tp>::operator= (const _Expr<E, _Tp>& __e) const
    {
        __valarray_copy (__e, _M_index.size(), _M_array,
                         _Array<size_t>(_M_index));
    }

#undef _DEFINE_VALARRAY_OPERATOR
#define _DEFINE_VALARRAY_OPERATOR(op, name)				\
template<typename _Tp>							\
inline void								\
gslice_array<_Tp>::operator op##= (const valarray<_Tp>& __v) const	\
{									\
    _Array_augmented_##name (_M_array, _Array<size_t>(_M_index),	\
                              _Array<_Tp> (__v), __v.size ());		\
}									\
									\
template<typename _Tp> template<class E>        			\
inline void 								\
gslice_array<_Tp>::operator op##= (const _Expr<E, _Tp>& __e) const	\
{									\
    _Array_augmented_##name (_M_array, _Array<size_t>(_M_index), __e,	\
                              _M_index.size());				\
}

_DEFINE_VALARRAY_OPERATOR(*, multiplies)
_DEFINE_VALARRAY_OPERATOR(/, divides)    
_DEFINE_VALARRAY_OPERATOR(%, modulus)
_DEFINE_VALARRAY_OPERATOR(+, plus)    
_DEFINE_VALARRAY_OPERATOR(-, minus)
_DEFINE_VALARRAY_OPERATOR(^, xor)
_DEFINE_VALARRAY_OPERATOR(&, and)
_DEFINE_VALARRAY_OPERATOR(|, or)
_DEFINE_VALARRAY_OPERATOR(<<, shift_left)
_DEFINE_VALARRAY_OPERATOR(>>, shift_right)

#undef _DEFINE_VALARRAY_OPERATOR

} // std::

#endif /* _CPP_BITS_GSLICE_ARRAY */

// Local Variables:
// mode:c++
// End: