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;
        }
};