AnimaL

Tutorial

Documentation

---

Documentation

• Overview
• Modules
• Namespace List
• Class Hierarchy
• Alphabetical List
• Compound List
• File List
• Namespace Members
• Compound Members
• File Members
• Examples

---

HermiteFunction.cpp

Go to the documentation of this file.

#include "HermiteFunction.h"
#include "assertions.h"
#include <math.h>
//#include <iostream.h>

namespace animal
{
namespace octree
{

HermiteFunction::HermiteFunction( const FloatingPointType *coeffs, const unsigned int nCoeffs ) :
    _nCoeffs(nCoeffs)
{
    Require( nCoeffs >= 2 );
    _coeffs = new FloatingPointType[_nCoeffs];
    _derivativeCoeffs = new FloatingPointType[_nCoeffs-1];
    
    for( unsigned int i=0 ; i<_nCoeffs ; ++i )
        _coeffs[i] = coeffs[i];
    for( unsigned int i=0 ; i<_nCoeffs-1 ; ++i )
    {
        _derivativeCoeffs[i] = (i+1)*_coeffs[i+1];
    }
}
HermiteFunction::HermiteFunction( const HermiteFunction& hf )
{
    _nCoeffs = hf._nCoeffs;
    
    _coeffs = new FloatingPointType[_nCoeffs];
    _derivativeCoeffs = new FloatingPointType[_nCoeffs-1];  
    
    for( unsigned int i=0 ; i<_nCoeffs ; ++i )
        _coeffs[i] = hf._coeffs[i];
    for( unsigned int i=0 ; i<_nCoeffs-1 ; ++i )
        _derivativeCoeffs[i] = hf._derivativeCoeffs[i];
}
HermiteFunction HermiteFunction::operator=( const HermiteFunction& hf )
{
    _nCoeffs = hf._nCoeffs;

    _coeffs = new FloatingPointType[_nCoeffs];
    _derivativeCoeffs = new FloatingPointType[_nCoeffs-1];
            
    for( unsigned int i=0 ; i<_nCoeffs ; ++i )
        _coeffs[i] = hf._coeffs[i];
    for( unsigned int i=0 ; i<_nCoeffs-1 ; ++i )
        _derivativeCoeffs[i] = hf._derivativeCoeffs[i];
        
    return *this;
}

HermiteFunction::~HermiteFunction()
{
    delete []_coeffs;
    delete []_derivativeCoeffs;
}
    


FloatingPointType HermiteFunction::compute( FloatingPointType x ) const
{
/*
    if( x < 0.0 )
        x = -x;
        
    Require( x<=2.0 );
    if( x > 1.0 )
        x = 2.0 - x;
*/
        
    FloatingPointType res = 0.0;
    //std::cerr << "x is " << x << "\n";
    for( int i=_nCoeffs-1 ; i>=0 ; --i )
    {
        //std::cerr << "we're adding " << res*x + _coeffs[i] << "( coeffs is " << _coeffs[i] << "\n";
        res = res*x + _coeffs[i];
    }
    return res;
}

FloatingPointType HermiteFunction::computeDerivative( FloatingPointType x ) const
{
/*  
    if( x < 0.0 )
        x = -x;
        
    Require( x<=2.0 );
    if( x > 1.0 )
        x = 2.0 - x;
*/
    
    FloatingPointType res = 0.0;
    //std::cerr << "Derivative : x is " << x << "\n";
    for( int i=_nCoeffs-2 ; i>=0 ; --i )
    {
        //std::cerr << "For i = " << i << "( coeffs is " << _derivativeCoeffs[i] << "\n";
        res = res*x + _derivativeCoeffs[i];
    }
    return res;
}

}
}

---