Added some unit tests for the kernel_matrix() functions.

--HG--
extra : convert_revision : svn%3Afdd8eb12-d10e-0410-9acb-85c331704f74/trunk%403330

dlib/test/CMakeLists.txt

@@ -41,6 +41,7 @@ set (tests
    hash_table.cpp
    image.cpp
    kcentroid.cpp
+   kernel_matrix.cpp
    lz77_buffer.cpp
    map.cpp
    matrix2.cpp

dlib/test/kernel_matrix.cpp

+// Copyright (C) 2009  Davis E. King (davisking@users.sourceforge.net)
+// License: Boost Software License   See LICENSE.txt for the full license.
+
+#include "tester.h"
+#include <dlib/svm.h>
+#include <vector>
+#include <sstream>
+
+namespace  
+{
+    using namespace test;
+    using namespace dlib;
+    using namespace std;
+    dlib::logger dlog("test.kernel_matrix");
+
+
+    class kernel_matrix_tester : public tester
+    {
+        /*!
+            WHAT THIS OBJECT REPRESENTS
+                This object represents a unit test.  When it is constructed
+                it adds itself into the testing framework.
+        !*/
+    public:
+        kernel_matrix_tester (
+        ) :
+            tester (
+                "test_kernel_matrix",       // the command line argument name for this test
+                "Run tests on the kernel_matrix functions.", // the command line argument description
+                0                     // the number of command line arguments for this test
+            )
+        {
+        }
+
+
+        void perform_test (
+        )
+        {
+            print_spinner();
+
+            typedef matrix<double,0,1> sample_type;
+            typedef radial_basis_kernel<sample_type> kernel_type;
+            kernel_type kern(0.1);
+
+            std::vector<sample_type> vect1;
+            std::vector<sample_type> vect2;
+
+            const sample_type samp = randm(4,1);
+
+            vect1.push_back(randm(4,1));
+            vect1.push_back(randm(4,1));
+            vect1.push_back(randm(4,1));
+            vect1.push_back(randm(4,1));
+
+            vect2.push_back(randm(4,1));
+            vect2.push_back(randm(4,1));
+            vect2.push_back(randm(4,1));
+            vect2.push_back(randm(4,1));
+            vect2.push_back(randm(4,1));
+
+            matrix<double> K;
+
+            K.set_size(vect1.size(), vect2.size());
+            for (long r = 0; r < K.nr(); ++r)
+            {
+                for (long c = 0; c < K.nc(); ++c)
+                {
+                    K(r,c) = kern(vect1[r], vect2[c]);
+                }
+            }
+            DLIB_TEST(equal(K, kernel_matrix(kern, vect1, vect2)));
+            DLIB_TEST(equal(K, kernel_matrix(kern, vector_to_matrix(vect1), vector_to_matrix(vect2))));
+
+
+            K.set_size(vect2.size(), vect1.size());
+            for (long r = 0; r < K.nr(); ++r)
+            {
+                for (long c = 0; c < K.nc(); ++c)
+                {
+                    K(r,c) = kern(vect2[r], vect1[c]);
+                }
+            }
+            DLIB_TEST(equal(K, kernel_matrix(kern, vect2, vect1)));
+            DLIB_TEST(equal(K, kernel_matrix(kern, vector_to_matrix(vect2), vector_to_matrix(vect1))));
+
+
+            K.set_size(vect1.size(), vect1.size());
+            for (long r = 0; r < K.nr(); ++r)
+            {
+                for (long c = 0; c < K.nc(); ++c)
+                {
+                    K(r,c) = kern(vect1[r], vect1[c]);
+                }
+            }
+            DLIB_TEST(equal(K, kernel_matrix(kern, vect1, vect1)));
+            DLIB_TEST(equal(K, kernel_matrix(kern, vect1)));
+            DLIB_TEST(equal(K, kernel_matrix(kern, vector_to_matrix(vect1), vector_to_matrix(vect1))));
+            DLIB_TEST(equal(K, kernel_matrix(kern, vector_to_matrix(vect1))));
+
+
+            K.set_size(vect1.size(),1);
+            for (long r = 0; r < K.nr(); ++r)
+            {
+                for (long c = 0; c < K.nc(); ++c)
+                {
+                    K(r,c) = kern(vect1[r], samp);
+                }
+            }
+            DLIB_TEST(equal(K, kernel_matrix(kern, vect1, samp)));
+            DLIB_TEST(equal(K, kernel_matrix(kern, vector_to_matrix(vect1), samp)));
+
+
+            K.set_size(1, vect1.size());
+            for (long r = 0; r < K.nr(); ++r)
+            {
+                for (long c = 0; c < K.nc(); ++c)
+                {
+                    K(r,c) = kern(samp, vect1[c]);
+                }
+            }
+            DLIB_TEST(equal(K, kernel_matrix(kern, samp, vect1)));
+            DLIB_TEST(equal(K, kernel_matrix(kern, samp, vector_to_matrix(vect1))));
+
+
+        }
+    };
+
+    // Create an instance of this object.  Doing this causes this test
+    // to be automatically inserted into the testing framework whenever this cpp file
+    // is linked into the project.  Note that since we are inside an unnamed-namespace 
+    // we won't get any linker errors about the symbol a being defined multple times. 
+    kernel_matrix_tester a;
+
+}
+
+

dlib/test/makefile

@@ -51,6 +51,7 @@ SRC += hash_set.cpp
 SRC += hash_table.cpp
 SRC += image.cpp
 SRC += kcentroid.cpp
+SRC += kernel_matrix.cpp
 SRC += lz77_buffer.cpp
 SRC += map.cpp
 SRC += matrix2.cpp