UnaryOperator.hpp

#ifndef __FIBER_FIELD_UNARYY_OPERATOR_HPP
#define __FIBER_FIELD_UNARYY_OPERATOR_HPP
 
#include "ComputationalField.hpp"
 
namespace Fiber
{
 
 
template <class Operator, class PrimaryValueType,
          class ResultValueType = PrimaryValueType>
struct  UnaryOperator : public Operator
{
        typedef PrimaryValueType        primary_value_type;
 
        typedef ResultValueType         result_type;
 
        UnaryOperator()
        {}
 
        template <Dims_t Dims>
        bool computeFragment(MemArray<Dims, result_type>&Result,
                             const MemArray<Dims, primary_value_type>&PrimaryArray,
                             const RefPtr<FragmentID>&myFragmentID)
        {
        MultiIndex<Dims> Size = Result.Size();
        MultiIndex<Dims> Index;
 
                do
                {
                        this->Operator::binary(Result[ Index ], PrimaryArray[ Index ] );
                }
                while( Index.inc(Size) );
 
                return true;
        }
 
        template <Dims_t Dims>
        bool computeIndexedFragment(MemArray<Dims, result_type>&Result,
                             const MemArray<Dims, primary_value_type>&PrimaryArray,
                             const RefPtr<FragmentID>&myFragmentID)
        {
        MultiIndex<Dims> Size = Result.Size();
        MultiIndex<Dims> Index;
 
                do
                {
                        this->Operator::binary(Index, Result[ Index ], PrimaryArray[ Index ] );
                }
                while( Index.inc(Size) );
 
                return true;
        }
};
 
} // namespace Fiber
 
 
#endif // __FIBER_FIELD_UNARYY_OPERATOR_HPP