Matrix_Operators_T.cpp

Go to the documentation of this file.

//==============================================================================
//
//  This file is part of GNSSTk, the ARL:UT GNSS Toolkit.
//
//  The GNSSTk is free software; you can redistribute it and/or modify
//  it under the terms of the GNU Lesser General Public License as published
//  by the Free Software Foundation; either version 3.0 of the License, or
//  any later version.
//
//  The GNSSTk is distributed in the hope that it will be useful,
//  but WITHOUT ANY WARRANTY; without even the implied warranty of
//  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
//  GNU Lesser General Public License for more details.
//
//  You should have received a copy of the GNU Lesser General Public
//  License along with GNSSTk; if not, write to the Free Software Foundation,
//  Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110, USA
//
//  This software was developed by Applied Research Laboratories at the
//  University of Texas at Austin.
//  Copyright 2004-2022, The Board of Regents of The University of Texas System
//
//==============================================================================
 
//==============================================================================
//
//  This software was developed by Applied Research Laboratories at the
//  University of Texas at Austin, under contract to an agency or agencies
//  within the U.S. Department of Defense. The U.S. Government retains all
//  rights to use, duplicate, distribute, disclose, or release this software.
//
//  Pursuant to DoD Directive 523024
//
//  DISTRIBUTION STATEMENT A: This software has been approved for public
//                            release, distribution is unlimited.
//
//==============================================================================
 
#include "Matrix_T.hpp"
 
/* Test the add/subtract and set operators */
int Matrix_T::operatorTest (void)
{
        gnsstk::TestUtil testFramework("Matrix","Operators",__FILE__,__LINE__);
        std::string failMesg;
        int badCount = 0;
 
        gnsstk::Vector<int> v1 = gnsstk::Vector<int>(16);
        for(int i = 0; i < 16; i++) v1[i] = i+1;
        gnsstk::Vector<int> v2 = gnsstk::Vector<int>(16);
        for(int i = 0; i < 16; i++) v2[i] = 16-i;
        gnsstk::Vector<int> v3 = gnsstk::Vector<int>(4);
        for(int i = 0; i < 4; i++) v3[i] = i+1;
        gnsstk::Vector<int> v4 = gnsstk::Vector<int>(8);
        for(int i = 0; i < 4; i++) {v4[i] = i+1+4; v4[i+4] = i+1+4;}
 
        gnsstk::Matrix<int> a(2,2,1);
        gnsstk::Matrix<int> b(8, 2, 3);
        gnsstk::Matrix<int> c(4, 2, 5);
        gnsstk::Matrix<int> d(4, 4, 7);
        gnsstk::Matrix<int> e(8, 2, v1);
        gnsstk::Matrix<int> f(4, 4, v2);
        gnsstk::Matrix<int> g(2, 2, v3);
        gnsstk::Matrix<int> h(4, 2, v4);
 
//--------------Matrix_operatorTest_1 - Does 2x2 additon function as expected?
        a += g; // 2x2 addition
        for(int i = 0; i < a.rows(); i++)
                for(int j = 0; j < a.cols(); j++)
                        if (v3(i*g.cols()+j)+1 != a(i,j)) {badCount++;}
        failDescriptionStream << "Check to see if addition between two 2x2 gnsstk::Matrix objects works. " << badCount << " of the elements are incorrect.";
        failDescriptionString = failDescriptionStream.str(); failDescriptionStream.str("");
        testFramework.assert(badCount==0, failDescriptionString, __LINE__);
        badCount = 0; // Reset error counter;
        a -= g; // 2x2 reset a
 
//--------------Matrix_operatorTest_2 - Does 8x2 subtraction work as expected?
        b -= e; // 8x2 subtraction
        for(int i = 0; i < b.rows(); i++)
                for(int j = 0; j < b.cols(); j++)
                        //cout << b(i,j) << " " << 3-(1+i+j*8) << endl;
                        if (3-v1(i*e.cols()+j) != b(i,j)) {badCount++;}
        failDescriptionStream << "Check to see if subtraction between two 8x2 gnsstk::Matrix objects works. " << badCount << " of the elements are incorrect.";
        failDescriptionString = failDescriptionStream.str(); failDescriptionStream.str("");
        testFramework.assert(badCount==0, failDescriptionString, __LINE__);
        badCount = 0; // Reset error counter;
        b += e; // 8x2 reset b
 
//--------------Matrix_operatorTest_3 - Does 4x2 subtraction work as expected?
        c -= h; // 4x2 subtraction
        for(int i = 0; i < c.rows(); i++)
                        for(int j = 0; j < c.cols(); j++)
                                //cout << c(i,j) << " " << 5-(5.+i) << endl;
                                if (5-v4(i*h.cols()+j) != c(i,j)) {badCount++;}
        failDescriptionStream << "Check to see if subtraction between two 4x2 gnsstk::Matrix objects works. " << badCount << " of the elements are incorrect.";
        failDescriptionString = failDescriptionStream.str(); failDescriptionStream.str("");
        testFramework.assert(badCount==0, failDescriptionString, __LINE__);
        badCount = 0; // Reset error counter;
        c += h; // 4x2 reset c
 
//--------------Matrix_operatorTest_4 - Does 4x4 addition work as expected?
        d += f; // 4x4 addition
        for(int i = 0; i < d.rows(); i++)
                        for(int j = 0; j < d.cols(); j++)
                        //cout << d(i,j) << " " << 16-i-4*j+7 << endl;
                                if (v2(i*f.cols()+j)+7 != d(i,j)) {badCount++;}
        failDescriptionStream << "Check to see if addition between two 4x4 gnsstk::Matrix objects works. " << badCount << " of the elements are incorrect.";
        failDescriptionString = failDescriptionStream.str(); failDescriptionStream.str("");
        testFramework.assert(badCount==0, failDescriptionString, __LINE__);
        badCount = 0; // Reset error counter
        d -= f; // 4x4 reset d
 
        return testFramework.countFails();
}
 
int main()
{
        int check, errorCounter = 0;
        Matrix_T testClass;
 
        check = testClass.operatorTest();
        errorCounter += check;
 
        std::cout << "Total Failures for " << __FILE__ << ": " << errorCounter << std::endl;
 
        return errorCounter; //Return the total number of errors
}