data.hh

//Copyright (C) Edward Rosten 2009, 2012

//All rights reserved.
//
//Redistribution and use in source and binary forms, with or without
//modification, are permitted provided that the following conditions
//are met:
//1. Redistributions of source code must retain the above copyright
//    notice, this list of conditions and the following disclaimer.
//2. Redistributions in binary form must reproduce the above copyright
//   notice, this list of conditions and the following disclaimer in the
//   documentation and/or other materials provided with the distribution.
//
//THIS SOFTWARE IS PROVIDED BY THE AUTHOR AND OTHER CONTRIBUTORS ``AS IS''
//AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
//IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
//ARE DISCLAIMED.  IN NO EVENT SHALL THE AUTHOR OR OTHER CONTRIBUTORS BE
//LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
//CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
//SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
//INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
//CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
//ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
//POSSIBILITY OF SUCH DAMAGE.

namespace TooN{
    namespace Internal{
        template<int N, class P> class Data;
    }


    template<int N, class P> struct Operator<Internal::Data<N, P> >
    {
        P vals[N];

        template<int R, int C, class T, class B>
        void eval(Matrix<R, C, T, B>& m) const
        {
            SizeMismatch<(R==-1?-1:(C==-1?-1:(R*C))), N>:: test(m.num_rows()*m.num_cols(), N);
            for(int r=0, n=0; r < R; r++)
                for(int c=0; c < C; c++, n++)
                    m[r][c] = vals[n];
        }

        template<int S, typename P2, typename B>
        void eval(Vector<S, P2, B>& v) const
        {
            SizeMismatch<S, N>::test(v.size(), N);
            for(int i=0; i <N; i++)
                v[i] = vals[i];
        }
    };

    #ifdef DOXYGEN_INCLUDE_ONLY_FOR_DOCS
    
        inline Operator<Internal::Data<N, double> > Data(double a, ...);

        template<typename Precision> inline Operator<Internal::Data<N, Precision> > Data(const Precision& a, ...);

    #endif

}