|| || || || || || || ||

---

Apply.h

Go to the documentation of this file.

#ifndef Impala_Core_Vector_Apply_h
#define Impala_Core_Vector_Apply_h

#include "Core/Vector/VectorSet.h"
#include "Core/Matrix/Mat.h"
#include "Util/ProgressPrinter.h"
#include "Link/Mpi/MpiFuncs.h"

namespace Impala
{
namespace Core
{
namespace Vector
{

// is this a pattern?
Matrix::Mat* Apply(double (*fn)(const VectorTem<double>&, const VectorTem<double>&),
                   const VectorSet<Array::Array2dScalarReal64>* horizontal,
                   const VectorSet<Array::Array2dScalarReal64>* vertical)
{
    ILOG_VAR(Core.Vector.Apply);

    int size1 = horizontal->Size(), size2 = vertical->Size();
    ILOG_DEBUG("s1 = " << size1 << ", s2 = " << size2);
    Matrix::Mat* m = new Matrix::Mat(size1, size2, 0, 0);

    double *ptr = m->CPB();
    for(int row=0 ; row<size2 ; ++row)
    {
        for(int col=0 ; col<size1 ; ++col)
        {
            VectorTem<double> v1(horizontal->GetVector(col, true));
            VectorTem<double> v2(vertical->GetVector(row, true));
            *ptr = (*fn)(v1, v2);
            ++ptr;
        }
    }

    return m;
}

void Kernelise(VectorTem<double>* dst, VectorTem<double> src, 
               const VectorSet<Array::Array2dScalarReal64>* set,
               double (*fn)(const VectorTem<double>&, const VectorTem<double>&))
{
    ILOG_VAR(Core.Vector.Kernelise);

    int size = set->Size();
    ILOG_DEBUG("size = " << size);
    if(size != dst->Size())
        ILOG_ERROR("sizes of dst and set don't match: "<<dst->Size()<<"!="<<size);
    if(set->GetVectorLength(0) != src.Size())
        ILOG_ERROR("sizes of src and set feature length don't match: "<<src.Size()<<"!="<<set->GetVectorLength(0));

    double *ptr = dst->GetData();
    for(int row=0 ; row<size ; ++row)
    {
        VectorTem<double> v1(set->GetVector(row, true));
        *ptr = (*fn)(v1, src);
        ++ptr;
    }
}

} // namespace Vector
} // namespace Core
} // namespace Impala


#endif

---