Arithmetic.hpp

template <class ArrayType>
struct  Arithmetic
{
        enum { Dims = ArrayType::Dims };

        typedef typename ArrayType::value_type  value_type;

static  RefPtr<MemArrayBase<Dims> > create(const RefPtr<MemArrayBase<Dims> > &MB, double f)
        {
        typedef value_type result_t; 
        typedef MemArray<Dims, result_t> ResultArray_t; 

        RefPtr<TypedCreativeIterator<value_type> > TI = MB->getIterator();

        RefPtr<ResultArray_t> RA = new ResultArray_t( MB -> Size() );

        RefPtr<TypedCreativeIterator<result_t> > RI = RA->getIterator(); 

        const value_type*src = TI->ptr();
                result_t*dst  = RI->ptr();

                for(index_t i=0; i<TI->count(); i++)
                {
                        dst[i] = src[i] * f;
                }

                return RA;
        }
};