NewtonIter.cpp File Reference

#include <Newton_bundle.h>

Go to the source code of this file.

Classes
class	CppNoddy::Example::VCube_root_problem
	Defines the problem: here z^3 -1 = 0, in vector form using real and imaginary parts. More...

Namespaces
namespace	CppNoddy
	A collection of OO numerical routines aimed at simple (typical) applied problems in continuum mechanics.
namespace	CppNoddy::Example

Functions
int	main ()

Function Documentation

main()

int main

(

)

Definition at line 38 of file NewtonIter.cpp.

{
  cout << "\n";
  cout << "=== Newton: vector residual root test  ==============\n";
  cout << "\n";
 
  // Instantiate the problem
  Example::VCube_root_problem residual_problem;
  // A Newton object
  Newton<double> newton( &residual_problem );
 
  // Set an initial guess
  DenseVector<double> guess( 2, 0.0 );
  guess[ 0 ] = -0.5;
  guess[ 1 ] = 1.5;
 
  try
  {
    newton.iterate( guess );
  }
  catch (const std::runtime_error &error )
  {
    cout << " \033[1;31;48m  * FAILED THROUGH EXCEPTION BEING RAISED \033[0m\n";
    return 1;
  }
 
  const double tol = 1.e-7;
  if ( ( abs( guess[ 0 ] + 0.5 ) > tol ) ||
       ( abs( guess[ 1 ] - sqrt( 3. ) / 2. ) > tol ) )
  {
    cout << "\033[1;31;48m  * FAILED \033[0m\n";
    guess.dump();
    cout << abs( guess[ 1 ] - sqrt( 3. ) / 2. ) << "\n";
    cout << abs( guess[ 0 ] + 0.5 ) << "\n";
    return 1;
  }
  cout << "\033[1;32;48m  * PASSED \033[0m\n";
  return 0;
 
}

References CppNoddy::DenseVector< _Type >::dump(), and CppNoddy::Newton< _Type >::iterate().