CEigen.cpp File Reference

#include <cmath>
#include <complex>
#include "copasi.h"
#include "CEigen.h"
#include "report/CCopasiObjectReference.h"
#include "utilities/CReadConfig.h"
#include "utilities/CCopasiMessage.h"
#include "utilities/CSort.h"
#include "lapack/lapackwrap.h"

Include dependency graph for CEigen.cpp:

Go to the source code of this file.

Functions
std::ostream &	operator<< (std::ostream &os, const CEigen &A)

Function Documentation

std::ostream& operator<<	(	std::ostream &	os,
		const CEigen &	A
	)

Definition at line 693 of file CEigen.cpp.

References CEigen::mHierarchy, CEigen::mImagOfMaxComplex, CEigen::mMaximagpart, CEigen::mMaxRealOfComplex, CEigen::mMaxrealpart, CEigen::mNcplxconj, CEigen::mNimag, CEigen::mNnegreal, CEigen::mNposreal, CEigen::mNreal, CEigen::mNzero, CEigen::mResolution, and CEigen::mStiffness.

{
  os << std::endl;
  os << "KINETIC STABILITY ANALYSIS";
  os << std::endl;
  os << "The linear stability analysis based on the eigenvalues" << std::endl;
  os << "of the Jacobian matrix is only valid for steady states." << std::endl;
  os << std::endl;
  os << "Summary:" << std::endl;
  os << "This state ";

  // Output statistics

  if (A.mMaxrealpart > A.mResolution)
    os << "is unstable";
  else if (A.mMaxrealpart < -A.mResolution)
    os << "is asymptotically stable";
  else
    os << "'s stability is undetermined";

  if (A.mMaximagpart > A.mResolution)
    {
      os << "," << std::endl;
      os << "transient states in its vicinity have oscillatory components";
    }

  os << "." << std::endl;
  os << std::endl;

  os << "Eigenvalue statistics:" << std::endl;
  // Output Max Real Part
  os << " Largest real part: ";
  os << std::setprecision(6) << A.mMaxrealpart << std::endl;
  // Output Max imaginary Part
  os << " Largest absolute imaginary part:  ";
  os << std::setprecision(6) << A.mMaximagpart << std::endl;

  if (A.mImagOfMaxComplex > A.mResolution)
    os << " The complex eigenvalues with the largest real part are:  "
       << A.mMaxRealOfComplex << " +|- " << A.mImagOfMaxComplex << "i" << std::endl;

  // Output Eigen-nreal
  os.unsetf(std::ios_base::scientific);
  os.unsetf(std::ios_base::showpoint);
  os << " " << A.mNreal;
  os << " are purely real" << std::endl;
  // Output Eigen-nimage
  os << " " << A.mNimag;
  os << " are purely imaginary" << std::endl;
  // Output Eigen-ncplxconj
  os << " " << A.mNcplxconj;
  os << " are complex" << std::endl;
  // Output Eigen-nzero
  os << " " << A.mNzero;
  os << " are equal to zero" << std::endl;
  // Output Eigen-nposreal
  os << " " << A.mNposreal;
  os << " have positive real part" << std::endl;
  // Output Eigen-nnegreal
  os << " " << A.mNnegreal;
  os << " have negative real part" << std::endl;

  // Set point manipulators
  os.setf(std::ios_base::showpoint);
  // Output Eigne-stiffness
  os << " stiffness = " << A.mStiffness << std::endl;
  os << " time hierarchy = " << A.mHierarchy << std::endl;

  return os;
}