Commit fa7bd33f authored by Tomáš Oberhuber's avatar Tomáš Oberhuber Committed by Tomáš Oberhuber
Browse files

Deleting deprecated method performSORIteration from matrices.

parent 6bcf11d2

src/TNL/Matrices/DenseMatrix.h

+0 −6
Original line number Diff line number Diff line
@@ -906,12 +906,6 @@ class DenseMatrix : public Matrix< Real, Device, Index, RealAllocator > 
      void getTransposition( const Matrix& matrix,

                             const RealType& matrixMultiplicator = 1.0 );



      template< typename Vector1, typename Vector2 >

      void performSORIteration( const Vector1& b,

                                const IndexType row,

                                Vector2& x,

                                const RealType& omega = 1.0 ) const;



      /**

       * \brief Assignment operator with exactly the same type of the dense matrix.

       *

src/TNL/Matrices/DenseMatrix.hpp

+0 −22
Original line number Diff line number Diff line
@@ -1037,28 +1037,6 @@ void DenseMatrix< Real, Device, Index, Organization, RealAllocator >::getTranspo 
   }

}



template< typename Real,

          typename Device,

          typename Index,

          ElementsOrganization Organization,

          typename RealAllocator >

   template< typename Vector1, typename Vector2 >

void DenseMatrix< Real, Device, Index, Organization, RealAllocator >::performSORIteration( const Vector1& b,

                                                        const IndexType row,

                                                        Vector2& x,

                                                        const RealType& omega ) const

{

   RealType sum( 0.0 ), diagonalValue;

   for( IndexType i = 0; i < this->getColumns(); i++ )

   {

      if( i == row )

         diagonalValue = this->getElement( row, row );

      else

         sum += this->getElement( row, i ) * x[ i ];

   }

   x[ row ] = ( 1.0 - omega ) * x[ row ] + omega / diagonalValue * ( b[ row ] - sum );

}



template< typename Real,

          typename Device,

          typename Index,

src/TNL/Matrices/DenseMatrixView.h

+0 −6
Original line number Diff line number Diff line
@@ -856,12 +856,6 @@ class DenseMatrixView : public MatrixView< Real, Device, Index > 
      void getTransposition( const Matrix& matrix,

                             const RealType& matrixMultiplicator = 1.0 );



      template< typename Vector1, typename Vector2 >

      void performSORIteration( const Vector1& b,

                                const IndexType row,

                                Vector2& x,

                                const RealType& omega = 1.0 ) const;



      /**

       * \brief Assignment operator with DenseMatrix.

       *

src/TNL/Matrices/DenseMatrixView.hpp

+0 −22
Original line number Diff line number Diff line
@@ -1006,28 +1006,6 @@ void DenseMatrixView< Real, Device, Index, Organization >::getTransposition( con 
   }

}



template< typename Real,

          typename Device,

          typename Index,

          ElementsOrganization Organization >

   template< typename Vector1, typename Vector2 >

void DenseMatrixView< Real, Device, Index, Organization >::performSORIteration( const Vector1& b,

                                                        const IndexType row,

                                                        Vector2& x,

                                                        const RealType& omega ) const

{

   RealType sum( 0.0 ), diagonalValue;

   for( IndexType i = 0; i < this->getColumns(); i++ )

   {

      if( i == row )

         diagonalValue = this->getElement( row, row );

      else

         sum += this->getElement( row, i ) * x[ i ];

   }

   x[ row ] = ( 1.0 - omega ) * x[ row ] + omega / diagonalValue * ( b[ row ] - sum );

}





template< typename Real,

          typename Device,

          typename Index,

src/TNL/Matrices/LambdaMatrix.h

+0 −7
Original line number Diff line number Diff line
@@ -490,13 +490,6 @@ class LambdaMatrix 
                          const IndexType begin = 0,

                          IndexType end = 0 ) const;





      template< typename Vector1, typename Vector2 >

      void performSORIteration( const Vector1& b,

                                const IndexType row,

                                Vector2& x,

                                const RealType& omega = 1.0 ) const;



      /**

       * \brief Method for printing the matrix to output stream.

       *