diff --git a/test/test_coordinate_matrix_sort.cpp b/test/test_coordinate_matrix_sort.cpp
new file mode 100644
index 00000000..48d34cac
--- /dev/null
+++ b/test/test_coordinate_matrix_sort.cpp
@@ -0,0 +1,65 @@
+#define BOOST_UBLAS_NO_ELEMENT_PROXIES
+
+#include<boost/numeric/ublas/matrix_sparse.hpp>
+#include<boost/numeric/ublas/io.hpp>
+
+#include "libs/numeric/ublas/test/utils.hpp"
+
+using std::cout;
+using std::endl;
+
+BOOST_UBLAS_TEST_DEF( test_coordinate_matrix_sort )
+{
+
+    boost::numeric::ublas::coordinate_matrix<double> matrix_mask(3, 3, 2);
+    cout << "Setting matrix(1,1) = 2.1" << endl;
+    matrix_mask(1,1) = 2.1;
+
+    cout << "Displaying matrix(1,1)" << endl;
+    std::cout << matrix_mask(1,1) << std::endl;
+
+    BOOST_UBLAS_DEBUG_TRACE( "Displaying matrix(1,1)" << matrix_mask(1,1) );
+    BOOST_UBLAS_TEST_CHECK( matrix_mask(1,1) == 2.1 );
+
+    BOOST_UBLAS_TEST_CHECK( matrix_mask.index1_data()[0] == 1 );
+    BOOST_UBLAS_TEST_CHECK( matrix_mask.index2_data()[0] == 1 );
+    BOOST_UBLAS_TEST_CHECK( matrix_mask.value_data()[0] == 2.1 );
+
+    BOOST_UBLAS_DEBUG_TRACE( "Setting matrix(0,1) = 1.1" );
+    matrix_mask(0, 1) = 1.1;
+
+    BOOST_UBLAS_TEST_CHECK( matrix_mask.index1_data()[0] == 1 );
+    BOOST_UBLAS_TEST_CHECK( matrix_mask.index2_data()[0] == 1 );
+    BOOST_UBLAS_TEST_CHECK( matrix_mask.value_data()[0] == 2.1 );
+
+    BOOST_UBLAS_TEST_CHECK( matrix_mask.index1_data()[1] == 0 );
+    BOOST_UBLAS_TEST_CHECK( matrix_mask.index2_data()[1] == 1 );
+    BOOST_UBLAS_TEST_CHECK( matrix_mask.value_data()[1] == 1.1 );
+
+    BOOST_UBLAS_DEBUG_TRACE( "Sort the matrix - this would be triggered by any element lookup." );
+    matrix_mask.sort();
+
+    BOOST_UBLAS_TEST_CHECK( matrix_mask.index1_data()[1] == 1 );
+    BOOST_UBLAS_TEST_CHECK( matrix_mask.index2_data()[1] == 1 );
+    BOOST_UBLAS_TEST_CHECK( matrix_mask.value_data()[1] == 2.1 );
+
+    BOOST_UBLAS_TEST_CHECK( matrix_mask.index1_data()[0] == 0 );
+    BOOST_UBLAS_TEST_CHECK( matrix_mask.index2_data()[0] == 1 );
+    BOOST_UBLAS_TEST_CHECK( matrix_mask.value_data()[0] == 1.1 );
+
+    BOOST_UBLAS_DEBUG_TRACE( "Displaying matrix(1,1)" << matrix_mask(1,1) );
+    BOOST_UBLAS_TEST_CHECK( matrix_mask(1,1) == 2.1 );
+
+    BOOST_UBLAS_DEBUG_TRACE( "Displaying matrix(0,1)" << matrix_mask(0,1) );
+    BOOST_UBLAS_TEST_CHECK( matrix_mask(1,1) == 1.1 );
+
+}
+
+int main()
+{
+    BOOST_UBLAS_TEST_BEGIN();
+
+    BOOST_UBLAS_TEST_DO( test_coordinate_matrix_sort );
+
+    BOOST_UBLAS_TEST_END();
+}