CLyapProblem.cpp

Go to the documentation of this file.

// Begin CVS Header
//   $Source: /Volumes/Home/Users/shoops/cvs/copasi_dev/copasi/lyap/CLyapProblem.cpp,v $
//   $Revision: 1.6 $
//   $Name:  $
//   $Author: ssahle $
//   $Date: 2012/04/22 14:54:54 $
// End CVS Header

// Copyright (C) 2012 by Pedro Mendes, Virginia Tech Intellectual
// Properties, Inc., University of Heidelberg, and The University
// of Manchester.
// All rights reserved.

// Copyright (C) 2001 - 2007 by Pedro Mendes, Virginia Tech Intellectual
// Properties, Inc. and EML Research, gGmbH.
// All rights reserved.

/**
 *  CLyapProblem class.
 */

#include <cmath>
#include <string>

#include "copasi.h"
#include "CLyapProblem.h"
#include "CLyapTask.h"
#include "model/CModel.h"
//#include "model/CState.h"
#include "CopasiDataModel/CCopasiDataModel.h"

/**
 *  Default constructor.
 */
CLyapProblem::CLyapProblem(const CCopasiContainer * pParent):
    CCopasiProblem(CCopasiTask::lyap, pParent),
    mpExponentNumber(NULL),
    //mpTimeSeriesRequested(NULL),
    mpDivergenceRequested(NULL),
    mpTransientTime(NULL)
{
  initializeParameter();
  CONSTRUCTOR_TRACE;
}

/**
 *  Copy constructor.
 *  @param "const CLyapProblem &" src
 */
CLyapProblem::CLyapProblem(const CLyapProblem & src,
                           const CCopasiContainer * pParent):
    CCopasiProblem(src, pParent),
    mpExponentNumber(NULL),
    //mpTimeSeriesRequested(NULL),
    mpDivergenceRequested(NULL),
    mpTransientTime(NULL)
{
  initializeParameter();
  CONSTRUCTOR_TRACE;
}

/**
 *  Destructor.
 */
CLyapProblem::~CLyapProblem()
{DESTRUCTOR_TRACE;}

bool CLyapProblem::elevateChildren()
{
  return true;
}

void CLyapProblem::initializeParameter()
{
  mpExponentNumber =
    assertParameter("ExponentNumber", CCopasiParameter::UINT, (unsigned C_INT32) 3)->getValue().pUINT;
  mpDivergenceRequested =
    assertParameter("DivergenceRequested", CCopasiParameter::BOOL, (bool) true)->getValue().pBOOL;
  mpTransientTime =
    assertParameter("TransientTime", CCopasiParameter::DOUBLE, (C_FLOAT64) 0.0)->getValue().pDOUBLE;

  elevateChildren();
}

void CLyapProblem::setExponentNumber(const unsigned C_INT32 & number)
{
  *mpExponentNumber = number;
}

const unsigned C_INT32 & CLyapProblem::getExponentNumber() const
{return *mpExponentNumber;}

void CLyapProblem::setTransientTime(const C_FLOAT64 & startTime)
{
  *mpTransientTime = startTime;
}

const C_FLOAT64 & CLyapProblem::getTransientTime() const
{return *mpTransientTime;}

/*void CLyapProblem::setTimeSeriesRequested(bool flag)
{
  *mpTimeSeriesRequested = flag;
}

bool CLyapProblem::timeSeriesRequested() const
  {return *mpTimeSeriesRequested;}
 */

void CLyapProblem::setDivergenceRequested(bool flag)
{
  *mpDivergenceRequested = flag;
}

bool CLyapProblem::divergenceRequested() const
{return *mpDivergenceRequested;}

void CLyapProblem::printResult(std::ostream * ostream) const
{
  CLyapTask* parent = dynamic_cast<CLyapTask*>(getObjectParent());

  if (!parent) return;

  parent->printResult(ostream);
}