MatrixBandedSolves_lapack.cpp File Reference

#include <cassert>
#include <Timer.h>
#include <Types.h>
#include <Utility.h>
#include <DenseLinearSystem.h>
#include <BandedLinearSystem.h>
#include "../Utils_Fill.h"

Go to the source code of this file.

Functions
int	main ()

Function Documentation

main()

int main

(

)

Definition at line 25 of file MatrixBandedSolves_lapack.cpp.

{
 
  cout << "\n";
  cout << "=== Matrix: Example linear banded solver  ============\n";
  cout << "\n";
 
  bool failed = false;
  // tolerance for the test
  const double tol = 1.e-10;
 
  //
  // SOLVE A BANDED REAL SYSTEM using LAPACK or native
  //
  const unsigned offdiag = 2;
  // N = size of some of the larger tests
  const unsigned N = 511;
  const double D = 12 * ( 1. / ( N - 1 ) ) * ( 1. / ( N - 1 ) );
 
  BandedMatrix<double> AB( N, offdiag, 0.0 );
  DenseVector<double> BB( N, D );
  Utils_Fill::fill_band( AB, 0, -30.0 );
  Utils_Fill::fill_band( AB, -1, 16.0 );
  Utils_Fill::fill_band( AB, 1, 16.0 );
  Utils_Fill::fill_band( AB, -2, -1.0 );
  Utils_Fill::fill_band( AB, 2, -1.0 );
  DenseMatrix<double> AD( AB );
  DenseVector<double> BD( BB );
 
  cout << "     Using dense matrix solver : " << N << "x" << N << " system \n";
  cout << "     Using the LAPACK LU dense routine\n";
  DenseLinearSystem<double> dense_system( &AD, &BD, "lapack" );
 
  try
  {
    dense_system.solve();
  }
  catch ( const std::runtime_error &error )
  {
    cout << " \033[1;31;48m  * FAILED THROUGH EXCEPTION BEING RAISED \033[0m\n";
    return 1;
  }
 
  cout << "  * Not checking.\n";
  cout << "\n";
  cout << "     Comparing the banded matrix solver solution : ";
  cout << N << "x" << 2 * AB.noffdiag() + 1 << " system \n";
 
  cout << "     Using the LAPACK LU banded routine\n";
  BandedLinearSystem<double> banded_system( &AB, &BB, "lapack" );
 
  try
  {
    banded_system.solve();
  }
  catch ( const std::runtime_error &error )
  {
    cout << " \033[1;31;48m * FAILED THROUGH EXCEPTION BEING RAISED \033[0m\n";
    return 1;
  }
 
  BB.sub( BD );
  if ( std::abs( BB.two_norm() ) > tol )
  {
    failed = true;
    cout << " \033[1;31;48m * Banded solver does not give same result as dense solver \033[0m\n";
  }
  else
  {
    cout << "     Banded solver agrees with the dense solver.\n";
  }
 
 
  //
  // CONCLUDING PASS/FAIL
  //
  if ( failed )
  {
    cout << "\033[1;31;48m  * FAILED \033[0m\n";
    return 1;
  }
  else
  {
    cout << "\033[1;32;48m  * PASSED \033[0m\n";
    return 0;
  }
}

References Utils_Fill::fill_band(), CppNoddy::BandedMatrix< _Type >::noffdiag(), CppNoddy::BandedLinearSystem< _Type >::solve(), CppNoddy::DenseLinearSystem< _Type >::solve(), CppNoddy::DenseVector< _Type >::sub(), and CppNoddy::DenseVector< _Type >::two_norm().