Conversion.hpp

#ifndef __FIBEROPERATIONS_CONVERSION_HPP
#define __FIBEROPERATIONS_CONVERSION_HPP
 
#include "fish/fiber/field/Field.hpp"
#include "FiberHDF5/mtoolslib/MDefs.hpp"
 
namespace Fiber
{
 
template<int N, class T, class O>
RefPtr<MemArray<N, T> > convertData( RefPtr<Field>& field, RefPtr<FragmentID> fid )
{
RefPtr<MemArray<N, O>> original = field->getData( fid );
 
        if( !original)
                return NullPtr();
 
RefPtr<MemArray<N, T>> converted = new  MemArray<N, T>( original->Size() );
 
MultiArray<N, O>& o_data = *original;
MultiArray<N, T>& c_data = *converted;
 
MultiIndex<N> mi;
        for(unsigned i = 0; i < N ; ++i)
                mi[i] = 0;
 
        do
        {
                c_data[mi] = T( o_data[mi] );
//              std::cout << c_data[mi] << std::endl;
        }
        while( mi.inc( original->Size() ) );
 
 
        return converted;
}
 
template<int N, class T>
RefPtr<MemArray<N, T> > getConvertedData( RefPtr<Field>& field, RefPtr<FragmentID> fid )
{
        MTOOLS_RETURN_EV( !field, "Fiber::getConvertedData() no field found!", NullPtr() );
 
RefPtr<MemArray<N, T>> result = field->getData( fid );
 
 
        if(!result)
        {
        RefPtr<FiberTypeBase> ftb = field->getElementFiberType();
        RefPtr<FiberTypeBase> ftb_e = ftb->element_type();
        //              ftb_e.speak();
        //cout << "type:" << ftb_e.getType().name() << endl;
        //cout << "type:" << typeid(ftb_e.getType()).name() << endl;
 
                if( ftb_e->getType() == typeid( uint8_t ) )
                        result = convertData<N, T, uint8_t>( field, fid );
                else if( ftb_e->getType() == typeid( uint16_t ) )
                        result = convertData<N, T, uint16_t>( field, fid );
                else if( ftb_e->getType() == typeid( uint32_t ) )
                        result = convertData<N, T, uint32_t>( field, fid );
                else if( ftb_e->getType() == typeid( uint64_t ) )
                        result = convertData<N, T, uint64_t>( field, fid );
                else if( ftb_e->getType() == typeid( float ) )
                        result = convertData<N, T, float>( field, fid );
                else if( ftb_e->getType() == typeid( double ) )
                        result = convertData<N, T, double>( field, fid );
                else if( ftb_e->getType() == typeid( int8_t ) )
                        result = convertData<N, T, int8_t>( field, fid );
                else if( ftb_e->getType() == typeid( int16_t ) )
                        result = convertData<N, T, int16_t>( field, fid );
                else if( ftb_e->getType() == typeid( int32_t ) )
                        result = convertData<N, T, int32_t>( field, fid );
                else if( ftb_e->getType() == typeid( int64_t ) )
                        result = convertData<N, T, int64_t>( field, fid );
                else
                        MTOOLS_ERR_LOG( "Fiber::getConvertedData( no type matched for conversion" );
        }
 
        return result;
}
}
#endif