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

---

TestMatNorm2Dist.h

Go to the documentation of this file.

#ifndef Impala_Core_Test_TestMatNorm2Dist_h
#define Impala_Core_Test_TestMatNorm2Dist_h

#include "Core/Test/TestBase.h"
#include "Core/Matrix/MatNorm2Dist.h"
#include "Core/Matrix/Mat.h"
#include "Core/Array/Abs.h"
#include "Core/Array/PixSum.h"
#include "Core/Array/Sub.h"
#include "Core/Array/PrintData.h"
#include "Core/Array/ReadRaw.h"
#include "Util/Database.h"

namespace Impala
{
namespace Core
{
namespace Test
{

using namespace Impala::Core::Array;
using namespace Impala::Core::Matrix;

class TestMatNorm2Dist : public TestBase
{
public:

    TestMatNorm2Dist(int testIndex)
    {
        mTestNameBase = "MatNorm2Dist";
        mTestIndex = testIndex;
    }

    virtual void
    Run()
    {
        Mat* a = 0;
        Mat* b = 0;
        Mat* c = 0;
        if(mTestIndex == 0) {
            /*
>> distance2([1,2,3,4,5; 4 3 2 9 0], [2.5 5 3; 0.9 0.4 3])

ans =

    3.4438    5.3814    2.2361
    2.1587    3.9699    1.0000
    1.2083    2.5612    1.0000
    8.2377    8.6579    6.0828
    2.6571    0.4000    3.6056*/

            a = MatCreate<Mat>(2, 5);
            *MatE(a, 0, 0) = 1;
            *MatE(a, 0, 1) = 2;
            *MatE(a, 0, 2) = 3;
            *MatE(a, 0, 3) = 4;
            *MatE(a, 0, 4) = 5;
            *MatE(a, 1, 0) = 4;
            *MatE(a, 1, 1) = 3;
            *MatE(a, 1, 2) = 2;
            *MatE(a, 1, 3) = 9;
            *MatE(a, 1, 4) = 0;
            b = MatCreate<Mat>(2, 3);
            *MatE(b, 0, 0) = 2.5;
            *MatE(b, 0, 1) = 5;
            *MatE(b, 0, 2) = 3;
            *MatE(b, 1, 0) = 0.9;
            *MatE(b, 1, 1) = 0.4;
            *MatE(b, 1, 2) = 3;
            c = MatCreate<Mat>(5, 3);
            *MatE(c, 0, 0) = 3.4438;
            *MatE(c, 0, 1) = 5.3814;
            *MatE(c, 0, 2) = 2.2361;
            *MatE(c, 1, 0) = 2.1587;
            *MatE(c, 1, 1) = 3.9699;
            *MatE(c, 1, 2) = 1.0;
            *MatE(c, 2, 0) = 1.2083;
            *MatE(c, 2, 1) = 2.5612;
            *MatE(c, 2, 2) = 1.0;
            *MatE(c, 3, 0) = 8.2377;
            *MatE(c, 3, 1) = 8.6579;
            *MatE(c, 3, 2) = 6.0828;
            *MatE(c, 4, 0) = 2.6571;
            *MatE(c, 4, 1) = 0.4000;
            *MatE(c, 4, 2) = 3.6056;
        } else if((mTestIndex == 1) || (mTestIndex == 2))
        {
            ReadRaw(a, "norm2dist_a.raw", &Util::Database::GetInstance());
            ReadRaw(b, "norm2dist_b.raw", &Util::Database::GetInstance());
            ReadRaw(c, "norm2dist_c.raw", &Util::Database::GetInstance());
            if(mTestIndex == 2)
            {
                Mat* temp1 = MatTranspose(a);
                Mat* temp2 = MatTranspose(b);
                delete a;
                delete b;
                a = temp1;
                b = temp2;
            }
        }
        Mat* result = 0;
        Timer speedTimer;
        if(mTestIndex <= 1) result = MatNorm2Dist(a,b);
        if(mTestIndex > 1) result = MatNorm2DistTransposed(a,b);
        std::cout << "Speed of " << mTestIndex << " " << speedTimer.SplitTimeStr() << std::endl;
        //PrintData(result);
        Sub(result, result, c);
        Abs(result, result);
        //PrintData(result);
        //PrintData(c);
        Real64 diff = PixSum(result);
        AddReport("Difference:" + MakeString(diff));
        std::cout << "Difference:" + MakeString(diff) << std::endl;
        if(diff > 0.0005) {
            // require small difference
            SetFailure("MatNorm2Dist too much difference failure!");
        }
        delete a;
        delete b;
        delete c;
        delete result;
    }

    virtual String
    Name() const
    {
        return FileNameCtor(mTestNameBase, MakeString(mTestIndex), "");
    }

private:
    String mTestNameBase;
    int    mTestIndex;

};


static TestBase* gMatNorm2DistSuite[] =
    {
        new TestMatNorm2Dist(0),
        new TestMatNorm2Dist(1),
        new TestMatNorm2Dist(2),
        (TestBase*)0
    };

} // namespace Test
} // namespace Core
} // namespace Impala

#endif

---