partialderivativetest.cpp File Reference

#include <print.h>
#include <partialderivativetest.h>
#include <exploreh.h>
#include <probsysequ01.h>
#include <partialderivative.h>

Include dependency graph for partialderivativetest.cpp:

Go to the source code of this file.

Functions
void	partialderivativetest01 ()

---

Function Documentation

void partialderivativetest01

(

)

Definition at line 12 of file partialderivativetest.cpp.

References partialderivative< EXP >::absvalue(), probsysequ01::counter, exploreh< FN, XI, T >::fn, exploreh< FN, XI, T >::print(), exploreh< FN, XI, T >::reset(), and SHOW.

Referenced by main().

{
  cout << "This test is to solve a non-linear system of equations" << endl;
  cout << " numerically and terminate the approximation when the" << endl;
  cout << " parital derivatives become small." << endl;
  exploreh<probsysequ01,double*,double> g(2,0.2,50000);

  double x0[2] = { 0.5, 0.5 };

  g.reset(x0);

  g.print(cout);

  double zero = 1e-5;
  cout << "Iterate until the partial derivative are less than" << endl;
  cout << " " << zero << " in magnitude." << endl;
  bool valid=true;

  partialderivative< exploreh<probsysequ01,double*,double> > partiald(g);

  uintc kmax=30;
  for (uint i=0; valid; ++i)
  {
    for (uint k=0; k<kmax; ++k)
      ++g;

    partiald();
    if (zero > partiald.absvalue() )
      valid=false;

cout << SHOW(i) << endl;
    //++g;

  }

  g.print(cout);
  cout << SHOW(g.fn.counter) << endl;

}