Derivative.hpp

#ifndef __FIBER_DERIVATIVE_HPP
#define __FIBER_DERIVATIVE_HPP
 
#include "MemArray.hpp"
#include <vector/FiniteDifferences.hpp>
 
namespace Fiber
{
 
/***@ingroup field
   Class for N-dimensional MultiArrays with MemCore memory management
   that allow computing its derivative.
 
   @param  T The type of this array.
   @param dT The type of the derived quantity, such as a co-vector.
 
   @author Werner Benger
 */
template <class dT>
struct  Derivative
{
 
template <Dims_t N, class Type>
static  MemCore::RefPtr<MemArray<N, dT> >
        compute(const MemArray<N, Type>&Values)
        {
        MemCore::Ref<MemArray<N, dT> > DerivativeData( Values.Size() );
                ComputeDerivative( *DerivativeData, Values);
                return DerivativeData;
        }
};
 
} // Fiber
 
#endif // __FIBER_DERIVATIVE_HPP