ltsa.cpp File Reference

#include "ltsa.h"

Include dependency graph for ltsa.cpp:

Go to the source code of this file.

Functions
void	computeWeightsVector (const Matrix2D< double > &A, Matrix1D< int > &weightVector)
void	getLessWeightNColumns (const Matrix2D< double > &A, const Matrix1D< int > &weightVector, Matrix2D< double > &B)

Function Documentation

computeWeightsVector()

void computeWeightsVector	(	const Matrix2D< double > &	A,
		Matrix1D< int > &	weightVector
	)

Definition at line 31 of file ltsa.cpp.

{
    weightVector.initZeros(MAT_XSIZE(A));
    for (size_t i = 0; i < MAT_XSIZE(A); ++i)
    {
        double diagelem = MAT_ELEM(A, i, i);
        for (size_t j = 0; j < MAT_XSIZE(A); ++j)
            if (((i != j) && (MAT_ELEM(A, j, j) > diagelem)) ||  ((MAT_ELEM(A, j, j) == diagelem) && (j < i)))
                VEC_ELEM(weightVector,i)++;
    }
}

getLessWeightNColumns()

void getLessWeightNColumns	(	const Matrix2D< double > &	A,
		const Matrix1D< int > &	weightVector,
		Matrix2D< double > &	B
	)

Definition at line 43 of file ltsa.cpp.

{
    size_t outputDim = MAT_XSIZE(B) - 1;
    for (size_t index = 0; index < outputDim; ++index)
        FOR_ALL_ELEMENTS_IN_MATRIX1D(weightVector)
            if (VEC_ELEM(weightVector,i) == index)
            {
                for (size_t j = 0; j < MAT_YSIZE(A); ++j)
                    MAT_ELEM(B, j, index + 1) = MAT_ELEM(A, j, i);
                break;
            }
}