cubic.cpp File Reference

#include <iomanip>
#include "wcpplib/stream/prstream.h"
#include "wcpplib/math/cubic.h"
#include "wcpplib/math/parabol.h"
#include "wcpplib/math/DoubleAc.h"
#include "wcpplib/util/FunNameStack.h"
#include "wcpplib/safetl/AbsArr.h"
#include "wcpplib/matrix/multiply.h"
#include "wcpplib/matrix/inverse.h"

Go to the source code of this file.

Functions
std::ostream &	operator<< (std::ostream &file, const Cubic &f)

Function Documentation

operator<<()

std::ostream & operator<<	(	std::ostream &	file,
		const Cubic &	f
	)

Definition at line 150 of file cubic.cpp.

                                                         {
  Cubic::double_complex z1;
  Cubic::double_complex z2;
  Cubic::double_complex z3;
  Ifile << "Cubic: s_xzero=" << f.s_xzero() << '\n';
  indn.n += 2;
  f.find_zero(z1, z2, z3);
  Ifile << "Cubic: a=" << f.a() << " b=" << f.b() << " c=" << f.c()
        << " d=" << f.d() << '\n';
  file << " z1,2,3=" << z1 << ' ' << z2 << ' ' << z3 << '\n';
  double r[3];
  int q;
  q = f.find_real_zero(r);
  Ifile << "The number of real zeros =" << q << '\n';
  int n;
  Ifile << "Solutions=";
  for (n = 0; n < q; n++)
    file << ' ' << r[n];
  file << '\n';
  double xmm[2];
  double ymm[2];
  int s_mm[2];
  q = f.find_maxmin(xmm, ymm, s_mm);
  Ifile << "Max/Min, q=" << q << '\n';
  indn.n += 2;
  for (n = 0; n < q; n++) {
    Ifile << "n=" << n << " xmm[n]=" << std::setw(13) << xmm[n]
          << " ymm[n]=" << std::setw(13) << ymm[n]
          << " s_mm[n]=" << std::setw(13) << s_mm[n] << '\n';
  }
  indn.n -= 2;
  indn.n -= 2;
 
  return file;
}