compiledDA_t.h

Go to the documentation of this file.

/******************************************************************************
*                                                                             *
* DIFFERENTIAL ALGEBRA CORE ENGINE                                            *
*                                                                             *
*******************************************************************************
*                                                                             *
* Copyright 2016 Politecnico di Milano (2014 Dinamica Srl)                    *
* Licensed under the Apache License, Version 2.0 (the "License");             *
* you may not use this file except in compliance with the License.            *
* You may obtain a copy of the License at                                     *
*                                                                             *
*    http://www.apache.org/licenses/LICENSE-2.0                               *
*                                                                             *
* Unless required by applicable law or agreed to in writing, software         *
* distributed under the License is distributed on an "AS IS" BASIS,           *
* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.    *
* See the License for the specific language governing permissions and         *
* limitations under the License.                                              *
*                                                                             *
*******************************************************************************/

/*
 * compiledDA_t.h
 *
 *  Created on: Apr 07, 2014
 *      Author: Dinamica Srl
 */

#ifndef DINAMICA_COMPILEDDA_T_H_
#define DINAMICA_COMPILEDDA_T_H_

// DACE classes
#include "dace/compiledDA.h"

namespace DACE{

/********************************************************************************
*     compiledDA evaluation routines
*********************************************************************************/
template<class V> V compiledDA::eval(const V &args) const{
    V res(dim);
    eval(args, res);

    return res;
}

template<class T> std::vector<T> compiledDA::eval(const std::initializer_list<T> l) const{
    std::vector<T> res(dim);
    eval(std::vector<T>(l), res);

    return res;
}

template<class T> std::vector<T> compiledDA::eval(const T args[], const unsigned int length) const{
    std::vector<T> arg(args,args+length);
    std::vector<T> res(dim);
    eval(arg, res);

    return res;
}

template<class T> std::vector<T> compiledDA::evalScalar(const T &arg) const{
    std::vector<T> args(1);
    std::vector<T> res(dim);
    args[0] = arg;
    eval(args, res);

    return res;
}

/* Actual evaluation with a std::vector<T> of algebraic type T.
   T must support at least:
     - default constructor    T::T()
     - assignment             T::operator=(const T& t)
     - assignment & addition  T::operator+=(const T& t)
     - double multiplication  T::operator*(const double d)
     - double addition        T::operator+(const double d)
 */
template<class T> void compiledDA::eval(const std::vector<T> &args, std::vector<T> &res) const{
    const unsigned int narg = args.size();
    unsigned int jlskip = ord+1;
    double *p = ac+2;
    T *xm = new T[ord+1];

    // prepare temporary powers
    xm[0] = args[0]*0.0 + 1.0;
    // constant part
    for(unsigned int i=0; i<dim; i++, p++)
        res[i] = args[0]*0.0 + (*p);
    // higher order terms
    for(unsigned int i=1; i<terms; i++){
        unsigned int jl = (unsigned int)(*p); p++;
        unsigned int jv = (unsigned int)(*p)-1; p++;
        if(jl > jlskip)
        {
            p += dim;
            continue;
        }
        if(jv >= narg)
        {
            jlskip = jl;
            p += dim;
            continue;
        }
        jlskip = ord+1;
        xm[jl] = xm[jl-1]*args[jv];
        for(unsigned int j=0; j<dim; j++, p++)
            if((*p)!=0.0) res[j] += xm[jl]*(*p);
    }

    delete[] xm;
}

}
#endif /* DINAMICA_COMPILEDDA_T_H_ */