DenseVector.cpp File Reference

Some simple sanity checks for the NVector class with both double and complex types. More...

#include <algorithm>
#include <Types.h>
#include <Timer.h>
#include <Utility.h>
#include <Functors.h>
#include "../Utils_Fill.h"

Go to the source code of this file.

Functions
int	main ()

Detailed Description

Some simple sanity checks for the NVector class with both double and complex types.

Definition in file DenseVector.cpp.

Function Documentation

main()

int main

(

)

Definition at line 18 of file DenseVector.cpp.

{
 
  cout << "\n";
  cout << "=== Vector: (dense) double/complex Example  =========\n";
  cout << "\n";
 
  const unsigned N = 1000000; // size of vectors
 
  cout << " Using vectors of size  " << N << " \n";
 
  const double tol = 1.e-13;
  bool failed = false;
 
  Utils_Fill::time_seed();
  // a real vector of random entries
  DenseVector<double> V( N, 0.0 );
  Utils_Fill::fill_random( V );
  // a complex vector initialised with the same real vector
  DenseVector<D_complex> CExample( V );
 
  // another real vector of random entries
  DenseVector<double> DExample( N, 0.0 );
  Utils_Fill::fill_random( DExample );
 
  const unsigned M = 100;    // number of repeats
  for ( unsigned i = 0; i < M; ++i )
  {
    Utility::dot( V, DExample );
  }
 
  cout << " \nComplex vectors \n";
  cout << " Testing norms, nearest_index, max/minabs_index\n";
  if ( abs( CExample.inf_norm() - *min_element( CExample.begin(), CExample.end(), absDiff_predicate<D_complex>( 1.0 ) ) ) >
       tol )
  {
    std::cout << "  - Failed inf_norm/nearest_index Example" << "\n";
    failed = true;
  }
  if ( abs( CExample.inf_norm() - *max_element( CExample.begin(), CExample.end(), abs_predicate<D_complex>() ) )
       > tol )
  {
    std::cout << "  - Failed inf_norm/maxabs_index Example" << "\n";
    failed = true;
  }
 
  CExample.one_norm();
  CExample.two_norm();
  CExample.inf_norm();
 
  cout << "\n";
  cout << " Double vectors \n";
 
  cout << " Testing norms, nearest_index, max/minabs_index\n";
  if ( abs( DExample.inf_norm() - *min_element( DExample.begin(), DExample.end(), absDiff_predicate<double>( 1.0 ) ) )
       > tol )
  {
    std::cout << "  - Failed inf_norm/nearest_index Example" << "\n";
    failed = true;
  }
  if ( abs( DExample.inf_norm() - *max_element( DExample.begin(), DExample.end(), abs_predicate<double>() ) )
       > tol )
  {
    std::cout << "  - Failed inf_norm/maxabs_index Example" << "\n";
    failed = true;
  }
 
  DExample.one_norm();
  DExample.two_norm();
  DExample.inf_norm();
 
  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 CppNoddy::DenseVector< _Type >::begin(), CppNoddy::Utility::dot(), CppNoddy::DenseVector< _Type >::end(), Utils_Fill::fill_random(), CppNoddy::DenseVector< _Type >::inf_norm(), CppNoddy::DenseVector< _Type >::one_norm(), Utils_Fill::time_seed(), CppNoddy::DenseVector< _Type >::two_norm(), and V.