DianaUtil.hpp

#include <CapeOpen.hpp>

#ifndef DIANA_UTIL_HPP
#define DIANA_UTIL_HPP
namespace Diana{
    namespace Util{

        using Common::Types::CapeArrayDouble;
        using Common::Types::CapeLong;
        using Common::Types::CapeDouble;
        inline void copyVec(double * dst, const Common::Types::CapeArrayDouble& src, 
                            CapeLong len=-1, CapeLong ssrc=0, CapeLong sdst=0)
        {
            if(len<0)
                len=src.size()-ssrc;
            Common::Types::CapeArrayDouble::const_iterator srci=src.begin()+ssrc,
                end=srci+len;
            dst+=sdst;
            for(;srci<end;srci++){
                *dst=*srci;dst++;
            }
        }
        inline void copyVec(CapeArrayDouble& dst, const CapeArrayDouble& src, 
                            CapeLong len=-1, CapeLong ssrc=0, CapeLong sdst=0)
        {
            if(len<0)
                len=src.size()-ssrc;
            CapeArrayDouble::const_iterator srci=src.begin()+ssrc,
                end=srci+len;
            CapeArrayDouble::iterator dsti=dst.begin()+sdst;
            for(;srci<end;srci++){
                *dsti=*srci;dsti++;
            }
        }

        inline CapeDouble* vecPtr(CapeArrayDouble& v){
            return static_cast<CapeDouble*>(&v[0]);
        }

    }
}
#endif //DIANA_UTIL_HPP