fn_dpolyval.hpp

Go to the documentation of this file.

/*
 * BELFEM -- The Berkeley Lab Finite Element Framework
 * Copyright (c) 2026, The Regents of the University of California,
 * through Lawrence Berkeley National Laboratory (subject to receipt of any required
 * approvals from the U.S. Dept. of Energy).  All rights reserved.
 *
 * Developers: Christian Messe, Gregory Giard
 *
 * See the top-level LICENSE file for the complete license and disclaimer.
 */

 
#ifndef BELFEM_FN_DPOLYVAL_HPP
#define BELFEM_FN_DPOLYVAL_HPP
 
#include "typedefs.hpp"
#include "cl_Vector.hpp"
 
namespace belfem
{
    template < typename T >
    T
    dpolyval( const Vector< T > & aCoeffs, const T aX )
    {
        const index_t tN = aCoeffs.length() - 1 ;
        T tPow = ( T ) tN ;
        T aResult = tPow * aCoeffs( 0 );
 
        for( index_t k=1; k<tN; ++k )
        {
            aResult *= aX;
            tPow -= 1.0 ;
            aResult += tPow * aCoeffs( k );
        }
 
        return aResult;
    }
}
#endif //BELFEM_FN_DPOLYVAL_HPP