#include <string>
#include <map>
#include "copasi/compareExpressions/CNormalFraction.h"

Include dependency graph for ConvertToCEvaluationNode.h:

This graph shows which files directly or indirectly include this file:

Functions
CEvaluationNode *	convertToCEvaluationNode (const CNormalFraction &fraction)

CEvaluationNode *	convertToCEvaluationNode (const CNormalItem &item)

CEvaluationNode *	convertToCEvaluationNode (const CNormalItemPower &itemPower)

CEvaluationNode *	convertToCEvaluationNode (const CNormalLcm &lcm)

CEvaluationNode *	convertToCEvaluationNode (const CNormalProduct &product)

CEvaluationNode *	convertToCEvaluationNode (const CNormalSum &sum)

CEvaluationNode *	convertToCEvaluationNode (const CNormalChoice &choice)

CEvaluationNode *	convertToCEvaluationNode (const CNormalChoiceLogical &choice)

CEvaluationNode *	convertToCEvaluationNode (const CNormalLogicalItem &choice)

CEvaluationNode *	convertToCEvaluationNode (const CNormalLogical &logical)

CEvaluationNode *	convertToCEvaluationNode (const CNormalBase &base)

CEvaluationNode *	convertToCEvaluationNode (const CNormalGeneralPower &pow)

CEvaluationNode *	convertToCEvaluationNode (const CNormalFunction &fun)

CEvaluationNode *	convertToCEvaluationNode (const CNormalCall &call)

CNormalCall *	createCall (const CEvaluationNode *node)

CNormalChoice *	createChoice (const CEvaluationNode *pNode)

CNormalFraction *	createFraction (const CEvaluationNode *node)

CNormalFunction *	createFunction (const CEvaluationNode *node)

CNormalGeneralPower *	createGeneralPower (const CEvaluationNode *node)

CNormalItem *	createItem (const CEvaluationNode *node)

CNormalItemPower *	createItemPower (const CEvaluationNode *node)

CNormalBase *	createItemPowerItem (const CEvaluationNode *pNode)

CNormalLogical *	createLogical (const CEvaluationNode *pNode)

CNormalChoiceLogical *	createLogicalChoice (const CEvaluationNode *pNode)

CNormalLogicalItem *	createLogicalItem (const CEvaluationNode *pNode)

CNormalFraction *	createNormalRepresentation (const CEvaluationNode *node)

CNormalProduct *	createProduct (const CEvaluationNode *node)

CNormalSum *	createSum (const CEvaluationNode *node)

bool	isLogical (const CEvaluationNode *pNode)

Variables
std::map< std::string, const CEvaluationNode * >	str2eval

Function Documentation

CEvaluationNode* convertToCEvaluationNode ( const CNormalFraction & fraction )

Definition at line 66 of file ConvertToCEvaluationNode.cpp.

References CCopasiNode< _Data >::addChild(), CNormalFraction::checkDenominatorOne(), convertToCEvaluationNode(), CEvaluationNodeOperator::DIVIDE, CNormalFraction::getDenominator(), CNormalFraction::getNumerator(), and pResult.

Referenced by convertToCEvaluationNode(), and CNormalTranslation::normAndSimplifyReptdly().

 {
   CEvaluationNode* pResult = NULL;
 
   if (fraction.checkDenominatorOne())
     {
       pResult = convertToCEvaluationNode(fraction.getNumerator());
     }
   else
     {
       CEvaluationNodeOperator* pDivision = new CEvaluationNodeOperator(CEvaluationNodeOperator::DIVIDE, "/");
       CEvaluationNode* pChild = convertToCEvaluationNode(fraction.getNumerator());
       pDivision->addChild(pChild);
       pChild = convertToCEvaluationNode(fraction.getDenominator());
       pDivision->addChild(pChild);
       pResult = pDivision;
     }
 
   return pResult;
 }

CEvaluationNode* convertToCEvaluationNode ( const CNormalItem & item )

Definition at line 87 of file ConvertToCEvaluationNode.cpp.

References CEvaluationNodeConstant::_INFINITY, CEvaluationNodeConstant::_NaN, CEvaluationNodeVariable::ANY, CNormalItem::CONSTANT, CEvaluationNodeConstant::EXPONENTIALE, CEvaluationNodeConstant::FALSE, CNormalItem::getName(), CNormalItem::getType(), CEvaluationNode::INVALID, CEvaluationNodeConstant::PI, CEvaluationNodeConstant::TRUE, and CNormalItem::VARIABLE.

 {
   CEvaluationNode* pNode = NULL;
   CEvaluationNode::Type type = CEvaluationNode::INVALID;
 
   switch (item.getType())
     {
       case CNormalItem::CONSTANT:
 
         // create a CEvaluationNodeNumber
         if (item.getName() == "pi" || item.getName() == "PI")
           {
             type = (CEvaluationNode::Type)(CEvaluationNodeConstant::PI);
           }
         else if (item.getName() == "EXPONENTIALE" || item.getName() == "exponentiale")
           {
             type = (CEvaluationNode::Type)(CEvaluationNodeConstant::EXPONENTIALE);
           }
         else if (item.getName() == "TRUE" || item.getName() == "true")
           {
             type = (CEvaluationNode::Type)(CEvaluationNodeConstant::TRUE);
           }
         else if (item.getName() == "FALSE" || item.getName() == "false")
           {
             type = (CEvaluationNode::Type)(CEvaluationNodeConstant::FALSE);
           }
         else if (item.getName() == "INFINITY" || item.getName() == "infinity")
           {
             type = (CEvaluationNode::Type)(CEvaluationNodeConstant::_INFINITY);
           }
         else if (item.getName() == "NAN" || item.getName() == "nan" || item.getName() == "Nan")
           {
             type = (CEvaluationNode::Type)(CEvaluationNodeConstant::_NaN);
           }
 
         pNode = new CEvaluationNodeConstant((CEvaluationNodeConstant::SubType)type, item.getName());
         break;
       case CNormalItem::VARIABLE:
         // create a CEvaluationNodeVariable
         pNode = new CEvaluationNodeVariable(CEvaluationNodeVariable::ANY, item.getName());
         break;
     }
 
   // TODO what about object nodes ?????
   return pNode;
 }

CEvaluationNode* convertToCEvaluationNode ( const CNormalItemPower & itemPower )

Definition at line 134 of file ConvertToCEvaluationNode.cpp.

References CCopasiNode< _Data >::addChild(), convertToCEvaluationNode(), CEvaluationNodeNumber::DOUBLE, CNormalItemPower::getExp(), CNormalItemPower::getItem(), CEvaluationNodeOperator::POWER, and pResult.

 {
   CEvaluationNode* pResult = NULL;
 
   if (fabs(itemPower.getExp() - 1.0) < 1e-12)
     {
       pResult = convertToCEvaluationNode(itemPower.getItem());
     }
   else
     {
       std::ostringstream sstream;
       CEvaluationNodeOperator* pPowerNode = new CEvaluationNodeOperator(CEvaluationNodeOperator::POWER, "^");
       CEvaluationNode* pChild = convertToCEvaluationNode(itemPower.getItem());
       pPowerNode->addChild(pChild);
       sstream << itemPower.getExp();
       pPowerNode->addChild(new CEvaluationNodeNumber(CEvaluationNodeNumber::DOUBLE, sstream.str()));
       pResult = pPowerNode;
     }
 
   return pResult;
 }

CEvaluationNode* convertToCEvaluationNode ( const CNormalLcm & lcm )

Definition at line 156 of file ConvertToCEvaluationNode.cpp.

References CCopasiNode< _Data >::addChild(), convertToCEvaluationNode(), CCopasiNode< _Data >::getChild(), CNormalLcm::getItemPowers(), CCopasiNode< _Data >::getParent(), CNormalLcm::getSums(), CEvaluationNodeOperator::MULTIPLY, pResult, and CCopasiNode< _Data >::removeChild().

 {
   CEvaluationNode* pResult = NULL;
   const std::set<CNormalItemPower*, compareItemPowers >& itemPowers = lcm.getItemPowers();
   std::set<CNormalItemPower*, compareItemPowers >::const_iterator it = itemPowers.begin();
   std::set<CNormalItemPower*, compareItemPowers >::const_iterator itEnd = itemPowers.end();
   CEvaluationNodeOperator* pMult = new CEvaluationNodeOperator(CEvaluationNodeOperator::MULTIPLY, "*");
   pResult = pMult;
   CEvaluationNode* pChild = NULL;
 
   while (it != itEnd)
     {
       assert(pMult != NULL);
       pChild = new CEvaluationNodeOperator(CEvaluationNodeOperator::MULTIPLY, "*");
       pMult->addChild(pChild);
       pChild = convertToCEvaluationNode(**it);
       pMult->addChild(pChild);
       pMult = dynamic_cast<CEvaluationNodeOperator*>(pMult->getChild());
       ++it;
     }
 
   const std::vector<CNormalSum*>& sums = lcm.getSums();
 
   std::vector<CNormalSum*>::const_iterator it2 = sums.begin();
 
   std::vector<CNormalSum*>::const_iterator itEnd2 = sums.end();
 
   while (it2 != itEnd2)
     {
       assert(pMult != NULL);
       pChild = new CEvaluationNodeOperator(CEvaluationNodeOperator::MULTIPLY, "*");
       pMult->addChild(pChild);
       pChild = convertToCEvaluationNode(**it);
       pMult->addChild(pChild);
       pMult = dynamic_cast<CEvaluationNodeOperator*>(pMult->getChild());
       ++it;
     }
 
   if (pMult->getParent() == pResult)
     {
       pResult->removeChild(pMult);
       delete pMult;
     }
   else
     {
       CEvaluationNode* pParent = dynamic_cast<CEvaluationNode*>(pMult->getParent());
       assert(pParent != NULL);
       pParent->removeChild(pMult);
       delete pMult;
       pMult = dynamic_cast<CEvaluationNodeOperator*>(pParent->getParent());
       assert(pMult != NULL);
       pMult->removeChild(pParent);
       pMult->addChild(dynamic_cast<CEvaluationNode*>(pParent->getChild())->copyBranch());
       delete pParent;
     }
 
   return pResult;
 }

CEvaluationNode* convertToCEvaluationNode ( const CNormalProduct & product )

Definition at line 215 of file ConvertToCEvaluationNode.cpp.

References convertToCEvaluationNode(), CNormalTranslation::createChain(), CEvaluationNodeNumber::DOUBLE, CNormalProduct::getFactor(), CNormalProduct::getItemPowers(), CNormalTranslation::ONE_NODE, pResult, and CNormalTranslation::TIMES_NODE.

 {
   CEvaluationNode* pResult = NULL;
   std::ostringstream sstream;
 
   if (product.getItemPowers().size() == 0)
     {
       sstream.precision(18);
       sstream << product.getFactor();
       pResult = new CEvaluationNodeNumber(CEvaluationNodeNumber::DOUBLE, sstream.str());
     }
   else
     {
       const std::set<CNormalItemPower*, compareItemPowers >& itemPowers = product.getItemPowers();
       std::set<CNormalItemPower*, compareItemPowers >::const_iterator it = itemPowers.begin();
       std::set<CNormalItemPower*, compareItemPowers >::const_iterator itEnd = itemPowers.end();
       CEvaluationNode* pChild = NULL;
       std::vector<CEvaluationNode*> products;
 
       while (it != itEnd)
         {
           pChild = convertToCEvaluationNode(**it);
           products.push_back(pChild);
           ++it;
         }
 
       if (fabs(product.getFactor() - 1.0) >= 1e-12)
         {
           sstream.precision(18);
           sstream << product.getFactor();
           products.push_back(new CEvaluationNodeNumber(CEvaluationNodeNumber::DOUBLE, sstream.str()));
         }
 
       pResult = CNormalTranslation::createChain(&CNormalTranslation::TIMES_NODE, &CNormalTranslation::ONE_NODE, products);
     }
 
   return pResult;
 }

CEvaluationNode* convertToCEvaluationNode ( const CNormalSum & sum )

Definition at line 254 of file ConvertToCEvaluationNode.cpp.

References convertToCEvaluationNode(), CEvaluationNode::copyBranch(), CNormalTranslation::createChain(), CNormalSum::getFractions(), CNormalSum::getProducts(), CNormalTranslation::PLUS_NODE, pResult, and CNormalTranslation::ZERO_NODE.

 {
   const std::set<CNormalFraction*>& fractions = sum.getFractions();
   std::set<CNormalFraction*>::const_iterator it = fractions.begin();
   std::set<CNormalFraction*>::const_iterator itEnd = fractions.end();
   std::vector<const CEvaluationNode*> summands;
   CEvaluationNode* pChild = NULL;
 
   while (it != itEnd)
     {
       pChild = convertToCEvaluationNode(**it);
       summands.push_back(pChild);
       ++it;
     }
 
   const std::set<CNormalProduct*, compareProducts >& products = sum.getProducts();
 
   std::set<CNormalProduct*, compareProducts >::const_iterator it2 = products.begin();
 
   std::set<CNormalProduct*, compareProducts >::const_iterator itEnd2 = products.end();
 
   while (it2 != itEnd2)
     {
       pChild = convertToCEvaluationNode(**it2);
       summands.push_back(pChild);
       ++it2;
     }
 
   CEvaluationNode* pResult = NULL;
 
   if (!summands.empty())
     {
       pResult = CNormalTranslation::createChain(&CNormalTranslation::PLUS_NODE, &CNormalTranslation::ZERO_NODE, summands);
     }
   else
     {
       pResult = CNormalTranslation::ZERO_NODE.copyBranch();
     }
 
   std::vector<const CEvaluationNode*>::iterator vIt = summands.begin(), vEndit = summands.end();
 
   while (vIt != vEndit)
     {
       delete *vIt;
       ++vIt;
     }
 
   return pResult;
 }

CEvaluationNode* convertToCEvaluationNode ( const CNormalChoice & choice )

Definition at line 1720 of file ConvertToCEvaluationNode.cpp.

References CCopasiNode< _Data >::addChild(), convertToCEvaluationNode(), CNormalChoice::getCondition(), CNormalChoice::getFalseExpression(), CNormalChoice::getTrueExpression(), and CEvaluationNodeChoice::IF.

 {
   // this will become a CEvaluationNodeChoice with an mpIf, an mpTrue and and
   // mpFalse
   CEvaluationNodeChoice* pChoiceNode = NULL;
   CEvaluationNode* pChild1 = convertToCEvaluationNode(choice.getCondition());
 
   if (pChild1 != NULL)
     {
       CEvaluationNode* pChild2 = convertToCEvaluationNode(choice.getTrueExpression());
 
       if (pChild2 != NULL)
         {
           CEvaluationNode* pChild3 = convertToCEvaluationNode(choice.getFalseExpression());
 
           if (pChild3 != NULL)
             {
               pChoiceNode = new CEvaluationNodeChoice(CEvaluationNodeChoice::IF, "IF");
               pChoiceNode->addChild(pChild1);
               pChoiceNode->addChild(pChild2);
               pChoiceNode->addChild(pChild3);
             }
         }
     }
 
   return pChoiceNode;
 }

CEvaluationNode* convertToCEvaluationNode ( const CNormalChoiceLogical & choice )

Definition at line 1748 of file ConvertToCEvaluationNode.cpp.

References CCopasiNode< _Data >::addChild(), convertToCEvaluationNode(), CNormalChoiceLogical::getCondition(), CNormalChoiceLogical::getFalseExpression(), CNormalChoiceLogical::getTrueExpression(), and CEvaluationNodeChoice::IF.

 {
   // this will become a CEvaluationNodeChoice with an mpIf, an mpTrue and and
   // mpFalse
   CEvaluationNodeChoice* pChoiceNode = NULL;
   CEvaluationNode* pChild1 = convertToCEvaluationNode(choice.getCondition());
 
   if (pChild1 != NULL)
     {
       CEvaluationNode* pChild2 = convertToCEvaluationNode(choice.getTrueExpression());
 
       if (pChild2 != NULL)
         {
           CEvaluationNode* pChild3 = convertToCEvaluationNode(choice.getFalseExpression());
 
           if (pChild3 != NULL)
             {
               pChoiceNode = new CEvaluationNodeChoice(CEvaluationNodeChoice::IF, "IF");
               pChoiceNode->addChild(pChild1);
               pChoiceNode->addChild(pChild2);
               pChoiceNode->addChild(pChild3);
             }
         }
     }
 
   return pChoiceNode;
 }

CEvaluationNode* convertToCEvaluationNode ( const CNormalLogicalItem & choice )

Definition at line 1776 of file ConvertToCEvaluationNode.cpp.

 {
   CEvaluationNode* pLogicalNode = NULL;
 
   switch (item.getType())
     {
       case CNormalLogicalItem::TRUE:
         pLogicalNode = new CEvaluationNodeConstant(CEvaluationNodeConstant::TRUE, "TRUE");
         pLogicalNode->compile(NULL);
         break;
       case CNormalLogicalItem::FALSE:
         pLogicalNode = new CEvaluationNodeConstant(CEvaluationNodeConstant::FALSE, "FALSE");
         pLogicalNode->compile(NULL);
         break;
       case CNormalLogicalItem::EQ:
         pLogicalNode = new CEvaluationNodeLogical(CEvaluationNodeLogical::EQ, "==");
         break;
       case CNormalLogicalItem::NE:
         pLogicalNode = new CEvaluationNodeLogical(CEvaluationNodeLogical::NE, "!=");
         break;
       case CNormalLogicalItem::LT:
         pLogicalNode = new CEvaluationNodeLogical(CEvaluationNodeLogical::LT, "<");
         break;
       case CNormalLogicalItem::GT:
         pLogicalNode = new CEvaluationNodeLogical(CEvaluationNodeLogical::GT, ">");
         break;
       case CNormalLogicalItem::GE:
         pLogicalNode = new CEvaluationNodeLogical(CEvaluationNodeLogical::GE, ">=");
         break;
       case CNormalLogicalItem::LE:
         pLogicalNode = new CEvaluationNodeLogical(CEvaluationNodeLogical::LE, "<=");
         break;
       case CNormalLogicalItem::INVALID:
         break;
     }
 
   if (item.getType() != CNormalLogicalItem::TRUE && item.getType() != CNormalLogicalItem::FALSE)
     {
       CEvaluationNode* pChild1 = convertToCEvaluationNode(item.getLeft());
 
       if (pChild1 == NULL)
         {
           delete pLogicalNode;
           pLogicalNode = NULL;
         }
       else
         {
           CEvaluationNode* pChild2 = convertToCEvaluationNode(item.getRight());
 
           if (pChild2 == NULL)
             {
               delete pLogicalNode;
               pLogicalNode = NULL;
             }
           else
             {
               pLogicalNode->addChild(pChild1);
               pLogicalNode->addChild(pChild2);
               pLogicalNode->compile(NULL);
             }
         }
     }
 
   return pLogicalNode;
 }

CEvaluationNode* convertToCEvaluationNode ( const CNormalLogical & logical )

Definition at line 1842 of file ConvertToCEvaluationNode.cpp.

References CCopasiNode< _Data >::addChild(), CEvaluationNodeLogical::AND, convertToCEvaluationNode(), CEvaluationNode::copyBranch(), CNormalTranslation::createChain(), CNormalLogical::getAndSets(), CNormalLogical::getChoices(), CNormalLogical::isNegated(), CNormalTranslation::NEUTRAL_ELEMENT_AND, CNormalTranslation::NEUTRAL_ELEMENT_OR, CEvaluationNodeFunction::NOT, CEvaluationNodeLogical::OR, and pResult.

 {
   CEvaluationNode* pResult = NULL;
   // go through all mAndSets and create AND combined chains which are then
   // combined with OR
   // do the same with all choice sets
   // combine the two results with OR
   CEvaluationNode* pNotNode = NULL;
   CEvaluationNode* pNode = NULL;
   CEvaluationNode* pOrNode = new CEvaluationNodeLogical(CEvaluationNodeLogical::OR, "OR");
   CEvaluationNode* pAndNode = new CEvaluationNodeLogical(CEvaluationNodeLogical::AND, "AND");
   std::vector<CEvaluationNode*> andElements;
   std::vector<CEvaluationNode*> orElements;
   CNormalLogical::ChoiceSetOfSets::const_iterator cIt = logical.getChoices().begin(), cEndit = logical.getChoices().end();
 
   while (cIt != cEndit)
     {
       CNormalLogical::ChoiceSet::const_iterator cInnerIt = (*cIt).first.begin(), cInnerEndit = (*cIt).first.end();
 
       while (cInnerIt != cInnerEndit)
         {
           pNode = convertToCEvaluationNode(*(*cInnerIt).first);
           assert(pNode != NULL);
 
           if ((*cInnerIt).second == true)
             {
               // only create the not node if it wouldn't result in the neutral
               // element
               if (*pNode != CNormalTranslation::NEUTRAL_ELEMENT_OR)
                 {
                   pNotNode = new CEvaluationNodeFunction(CEvaluationNodeFunction::NOT, "NOT");
                   pNotNode->addChild(pNode);
                   pNode = pNotNode;
                 }
               else
                 {
                   pNode = NULL;
                 }
             }
 
           ++cInnerIt;
 
           // if it is not the neutral element or it is the last element to be
           // inserted into an otherwise empty vector, insert the element
           if (((pNode != NULL) && ((*pNode) != CNormalTranslation::NEUTRAL_ELEMENT_AND)) || (andElements.empty() && cInnerIt == cInnerEndit))
             {
               andElements.push_back(pNode);
             }
         }
 
       // create the and chain
       pNode = CNormalTranslation::createChain(pAndNode, &CNormalTranslation::NEUTRAL_ELEMENT_AND, andElements);
       assert(pNode != NULL);
       andElements.clear();
 
       // check *cIt.second if it is true
       if ((*cIt).second == true)
         {
           // create a not node if it would not result in the creation of the
           // neutral element
           if (*pNode != CNormalTranslation::NEUTRAL_ELEMENT_AND)
             {
               pNotNode = new CEvaluationNodeFunction(CEvaluationNodeFunction::NOT, "NOT");
               pNode->addChild(pNode);
               pNode = pNotNode;
             }
           else
             {
               pNode = NULL;
             }
         }
 
       ++cIt;
 
       // if it is not the neutral element or it is the last element to be
       // inserted into an otherwise empty vector, insert the element
       if (((pNode != NULL) && ((*pNode) != CNormalTranslation::NEUTRAL_ELEMENT_OR)) || (orElements.empty() && cIt == cEndit))
         {
           orElements.push_back(pNode);
         }
     }
 
   // create the OR chain
   if (!orElements.empty())
     {
       pNode = CNormalTranslation::createChain(pOrNode, &CNormalTranslation::NEUTRAL_ELEMENT_OR, orElements);
       assert(pNode != NULL);
       orElements.clear();
     }
 
   pResult = pNode;
   pNode = NULL;
 
   CNormalLogical::ItemSetOfSets::const_iterator iIt = logical.getAndSets().begin(), iEndit = logical.getAndSets().end();
 
   while (iIt != iEndit)
     {
       CNormalLogical::ItemSet::const_iterator iInnerIt = (*iIt).first.begin(), iInnerEndit = (*iIt).first.end();
 
       while (iInnerIt != iInnerEndit)
         {
           pNode = convertToCEvaluationNode(*(*iInnerIt).first);
           assert(pNode != NULL);
 
           if ((*iInnerIt).second == true)
             {
               // only create the not node if it wouldn't result in the neutral
               // element
               if (*pNode != CNormalTranslation::NEUTRAL_ELEMENT_OR)
                 {
                   pNotNode = new CEvaluationNodeFunction(CEvaluationNodeFunction::NOT, "NOT");
                   pNotNode->addChild(pNode);
                   pNode = pNotNode;
                 }
               else
                 {
                   delete pNode;
                   pNode = NULL;
                 }
             }
 
           ++iInnerIt;
 
           // if it is not the neutral element or it is the last element to be
           // inserted into an otherwise empty vector, insert the element
           if (((pNode != NULL) && ((*pNode) != CNormalTranslation::NEUTRAL_ELEMENT_AND)) || (andElements.empty() && iInnerIt == iInnerEndit))
             {
               andElements.push_back(pNode);
             }
         }
 
       // create the and chain
       pNode = CNormalTranslation::createChain(pAndNode, &CNormalTranslation::NEUTRAL_ELEMENT_AND, andElements);
       assert(pNode != NULL);
       andElements.clear();
 
       // check *iIt.second if it is true
       if ((*iIt).second == true)
         {
 
           // create a not node if it does not result in the neutral node
           if (*pNode != CNormalTranslation::NEUTRAL_ELEMENT_AND)
             {
               pNotNode = new CEvaluationNodeFunction(CEvaluationNodeFunction::NOT, "NOT");
               pNotNode->addChild(pNode);
               pNode = pNotNode;
             }
           else
             {
               delete pNode;
               pNode = NULL;
             }
         }
 
       ++iIt;
 
       // if it is not the neutral element or it is the last element to be
       // inserted into an otherwise empty vector, insert the element
       if (((pNode != NULL) && ((*pNode) != CNormalTranslation::NEUTRAL_ELEMENT_OR)) || (orElements.empty() && iIt == iEndit))
         {
           orElements.push_back(pNode);
         }
     }
 
   // create the OR chain
   if (!orElements.empty())
     {
       pNode = CNormalTranslation::createChain(pOrNode, &CNormalTranslation::NEUTRAL_ELEMENT_OR, orElements);
       assert(pNode != NULL);
       orElements.clear();
     }
 
   if (pResult == NULL)
     {
       pResult = pNode;
       assert(pResult != NULL);
     }
   else
     {
       if (*pResult == CNormalTranslation::NEUTRAL_ELEMENT_OR)
         {
           pResult = pNode;
         }
       else if (*pNode == CNormalTranslation::NEUTRAL_ELEMENT_OR)
         {
           CEvaluationNode* pTmpNode = pOrNode->copyBranch();
           pTmpNode->addChild(pResult);
           pTmpNode->addChild(pNode);
           pResult = pTmpNode;
         }
     }
 
   // check if mNot is set
   if (logical.isNegated() == true)
     {
       pNotNode = new CEvaluationNodeFunction(CEvaluationNodeFunction::NOT, "NOT");
       pNode->addChild(pResult);
       pResult = pNotNode;
     }
 
   // clean up
   delete pOrNode;
   delete pAndNode;
   assert(pResult != NULL);
   return pResult;
 }

CEvaluationNode* convertToCEvaluationNode ( const CNormalBase & base )

Definition at line 1092 of file ConvertToCEvaluationNode.cpp.

References convertToCEvaluationNode().

 {
   CEvaluationNode* pNode = NULL;
 
   if (dynamic_cast<const CNormalItem*>(&base) != NULL)
     {
       pNode = convertToCEvaluationNode(dynamic_cast<const CNormalItem&>(base));
     }
   else if (dynamic_cast<const CNormalItemPower*>(&base) != NULL)
     {
       pNode = convertToCEvaluationNode(dynamic_cast<const CNormalItemPower&>(base));
     }
   else if (dynamic_cast<const CNormalGeneralPower*>(&base) != NULL)
     {
       pNode = convertToCEvaluationNode(dynamic_cast<const CNormalGeneralPower&>(base));
     }
   else if (dynamic_cast<const CNormalFunction*>(&base) != NULL)
     {
       pNode = convertToCEvaluationNode(dynamic_cast<const CNormalFunction&>(base));
     }
   else if (dynamic_cast<const CNormalCall*>(&base) != NULL)
     {
       pNode = convertToCEvaluationNode(dynamic_cast<const CNormalCall&>(base));
     }
   else if (dynamic_cast<const CNormalFraction*>(&base) != NULL)
     {
       pNode = convertToCEvaluationNode(dynamic_cast<const CNormalFraction&>(base));
     }
   else if (dynamic_cast<const CNormalProduct*>(&base) != NULL)
     {
       pNode = convertToCEvaluationNode(dynamic_cast<const CNormalProduct&>(base));
     }
   else if (dynamic_cast<const CNormalSum*>(&base) != NULL)
     {
       pNode = convertToCEvaluationNode(dynamic_cast<const CNormalSum&>(base));
     }
   else if (dynamic_cast<const CNormalLogical*>(&base) != NULL)
     {
       pNode = convertToCEvaluationNode(dynamic_cast<const CNormalLogical&>(base));
     }
   else if (dynamic_cast<const CNormalChoice*>(&base) != NULL)
     {
       pNode = convertToCEvaluationNode(dynamic_cast<const CNormalChoice&>(base));
     }
 
   return pNode;
 }

CEvaluationNode* convertToCEvaluationNode ( const CNormalGeneralPower & pow )

Definition at line 1140 of file ConvertToCEvaluationNode.cpp.

References CCopasiNode< _Data >::addChild(), CNormalFraction::checkIsOne(), convertToCEvaluationNode(), CNormalGeneralPower::getLeft(), CNormalGeneralPower::getRight(), CNormalGeneralPower::getType(), CNormalGeneralPower::INVALID, CNormalGeneralPower::MODULO, CEvaluationNodeOperator::MODULUS, CNormalGeneralPower::POWER, CEvaluationNodeOperator::POWER, and pResult.

 {
   CEvaluationNode* pResult = NULL;
 
   switch (pow.getType())
     {
       case CNormalGeneralPower::POWER:
         pResult = new CEvaluationNodeOperator(CEvaluationNodeOperator::POWER, "^");
         break;
       case CNormalGeneralPower::MODULO:
         pResult = new CEvaluationNodeOperator(CEvaluationNodeOperator::MODULUS, "%");
         break;
       case CNormalGeneralPower::INVALID:
         break;
     }
 
   if (pResult != NULL)
     {
       if (pow.getRight().checkIsOne())
         {
           delete pResult;
           pResult = convertToCEvaluationNode(pow.getLeft());
         }
       else
         {
           CEvaluationNode* pChild = convertToCEvaluationNode(pow.getLeft());
           pResult->addChild(pChild);
           pChild = convertToCEvaluationNode(pow.getRight());
           pResult->addChild(pChild);
         }
     }
 
   return pResult;
 }

CEvaluationNode* convertToCEvaluationNode ( const CNormalFunction & fun )

Definition at line 1220 of file ConvertToCEvaluationNode.cpp.

 {
   CEvaluationNode* pResult = NULL;
   std::string data;
   CEvaluationNodeFunction::SubType subType = CEvaluationNodeFunction::INVALID;
 
   switch (fun.getType())
     {
       case CNormalFunction::LOG:
         subType = CEvaluationNodeFunction::LOG;
         data = CNormalFunction::NAMES[fun.getType()];
         break;
       case CNormalFunction::LOG10:
         subType = CEvaluationNodeFunction::LOG10;
         data = CNormalFunction::NAMES[fun.getType()];
         break;
       case CNormalFunction::EXP:
         subType = CEvaluationNodeFunction::EXP;
         data = CNormalFunction::NAMES[fun.getType()];
         break;
       case CNormalFunction::SIN:
         subType = CEvaluationNodeFunction::SIN;
         data = CNormalFunction::NAMES[fun.getType()];
         break;
       case CNormalFunction:: COS:
         subType = CEvaluationNodeFunction::COS;
         data = CNormalFunction::NAMES[fun.getType()];
         break;
       case CNormalFunction::TAN:
         subType = CEvaluationNodeFunction::TAN;
         data = CNormalFunction::NAMES[fun.getType()];
         break;
       case CNormalFunction::SEC:
         subType = CEvaluationNodeFunction::SEC;
         data = CNormalFunction::NAMES[fun.getType()];
         break;
       case CNormalFunction::CSC:
         subType = CEvaluationNodeFunction::CSC;
         data = CNormalFunction::NAMES[fun.getType()];
         break;
       case CNormalFunction::COT:
         subType = CEvaluationNodeFunction::COT;
         data = CNormalFunction::NAMES[fun.getType()];
         break;
       case CNormalFunction::SINH:
         subType = CEvaluationNodeFunction::SINH;
         data = CNormalFunction::NAMES[fun.getType()];
         break;
       case CNormalFunction::COSH:
         subType = CEvaluationNodeFunction::COSH;
         data = CNormalFunction::NAMES[fun.getType()];
         break;
       case CNormalFunction::TANH:
         subType = CEvaluationNodeFunction::TANH;
         data = CNormalFunction::NAMES[fun.getType()];
         break;
       case CNormalFunction::SECH:
         subType = CEvaluationNodeFunction::SECH;
         data = CNormalFunction::NAMES[fun.getType()];
         break;
       case CNormalFunction::CSCH:
         subType = CEvaluationNodeFunction::CSCH;
         data = CNormalFunction::NAMES[fun.getType()];
         break;
       case CNormalFunction::COTH:
         subType = CEvaluationNodeFunction::COTH;
         data = CNormalFunction::NAMES[fun.getType()];
         break;
       case CNormalFunction::ARCSIN:
         subType = CEvaluationNodeFunction::ARCSIN;
         data = CNormalFunction::NAMES[fun.getType()];
         break;
       case CNormalFunction::ARCCOS:
         subType = CEvaluationNodeFunction::ARCCOS;
         data = CNormalFunction::NAMES[fun.getType()];
         break;
       case CNormalFunction::ARCTAN:
         subType = CEvaluationNodeFunction::ARCTAN;
         data = CNormalFunction::NAMES[fun.getType()];
         break;
       case CNormalFunction::ARCSEC:
         subType = CEvaluationNodeFunction::ARCSEC;
         data = CNormalFunction::NAMES[fun.getType()];
         break;
       case CNormalFunction::ARCCSC:
         subType = CEvaluationNodeFunction::ARCCSC;
         data = CNormalFunction::NAMES[fun.getType()];
         break;
       case CNormalFunction::ARCCOT:
         subType = CEvaluationNodeFunction::ARCCOT;
         data = CNormalFunction::NAMES[fun.getType()];
         break;
       case CNormalFunction::ARCSINH:
         subType = CEvaluationNodeFunction::ARCSINH;
         data = CNormalFunction::NAMES[fun.getType()];
         break;
       case CNormalFunction::ARCCOSH:
         subType = CEvaluationNodeFunction::ARCCOSH;
         data = CNormalFunction::NAMES[fun.getType()];
         break;
       case CNormalFunction::ARCTANH:
         subType = CEvaluationNodeFunction::ARCTANH;
         data = CNormalFunction::NAMES[fun.getType()];
         break;
       case CNormalFunction::ARCSECH:
         subType = CEvaluationNodeFunction::ARCSECH;
         data = CNormalFunction::NAMES[fun.getType()];
         break;
       case CNormalFunction::ARCCSCH:
         subType = CEvaluationNodeFunction::ARCCSCH;
         data = CNormalFunction::NAMES[fun.getType()];
         break;
       case CNormalFunction::ARCCOTH:
         subType = CEvaluationNodeFunction::ARCCOTH;
         data = CNormalFunction::NAMES[fun.getType()];
         break;
       case CNormalFunction::SQRT:
         subType = CEvaluationNodeFunction::SQRT;
         data = CNormalFunction::NAMES[fun.getType()];
         break;
       case CNormalFunction::ABS:
         subType = CEvaluationNodeFunction::ABS;
         data = CNormalFunction::NAMES[fun.getType()];
         break;
       case CNormalFunction::FLOOR:
         subType = CEvaluationNodeFunction::FLOOR;
         data = CNormalFunction::NAMES[fun.getType()];
         break;
       case CNormalFunction::CEIL:
         subType = CEvaluationNodeFunction::CEIL;
         data = CNormalFunction::NAMES[fun.getType()];
         break;
       case CNormalFunction::FACTORIAL:
         subType = CEvaluationNodeFunction::FACTORIAL;
         data = CNormalFunction::NAMES[fun.getType()];
         break;
       case CNormalFunction::INVALID:
         data = "@";
         break;
     }
 
   pResult = new CEvaluationNodeFunction(subType, data);
 
   if (subType != CEvaluationNodeFunction::INVALID)
     {
       CEvaluationNode* pChild = convertToCEvaluationNode(fun.getFraction());
       pResult->addChild(pChild);
     }
 
   return pResult;
 }

CEvaluationNode* convertToCEvaluationNode ( const CNormalCall & call )

Definition at line 1175 of file ConvertToCEvaluationNode.cpp.

References CCopasiNode< _Data >::addChild(), convertToCEvaluationNode(), CEvaluationNodeDelay::DELAY, CNormalCall::DELAY, CNormalCall::EXPRESSION, CEvaluationNodeCall::EXPRESSION, CNormalCall::FUNCTION, CEvaluationNodeCall::FUNCTION, CNormalCall::getFractions(), CNormalCall::getName(), CNormalCall::getType(), CEvaluationNodeCall::INVALID, and CNormalCall::INVALID.

 {
   CEvaluationNode* pCall = NULL;
 
   // check if the name contains any non-whitespace characters at all
   if (call.getName().find_first_not_of("\t\r\n ") != std::string::npos)
     {
       CEvaluationNodeCall::SubType type = CEvaluationNodeCall::INVALID;
 
       switch (call.getType())
         {
           case CNormalCall::EXPRESSION:
             type = CEvaluationNodeCall::EXPRESSION;
             break;
           case CNormalCall::FUNCTION:
             type = CEvaluationNodeCall::FUNCTION;
             break;
           case CNormalCall::DELAY:
             pCall = new CEvaluationNodeDelay(CEvaluationNodeDelay::DELAY, "delay");
             break;
           case CNormalCall::INVALID:
             break;
         }
 
       if (!pCall)
         {
           pCall = new CEvaluationNodeCall(type, call.getName());
         }
 
       const std::vector<CNormalFraction*>& children = call.getFractions();
 
       std::vector<CNormalFraction*>::const_iterator it = children.begin(), endit = children.end();
 
       while (it != endit)
         {
           CEvaluationNode* pChild = convertToCEvaluationNode(**it);
           assert(pChild != NULL);
           pCall->addChild(pChild);
           ++it;
         }
     }
 
   return pCall;
 }

CNormalCall* createCall ( const CEvaluationNode * node )

Create an function call from an evaluation node.

Returns: CNormalCall*, pointer to newly created call

Definition at line 1451 of file ConvertToCEvaluationNode.cpp.

References CNormalCall::add(), CEvaluationNode::CALL, createNormalRepresentation(), CNormalCall::DELAY, CEvaluationNode::DELAY, CNormalCall::EXPRESSION, CEvaluationNodeCall::EXPRESSION, CNormalCall::FUNCTION, CEvaluationNodeCall::FUNCTION, CCopasiNode< _Data >::getChild(), CCopasiNode< _Data >::getData(), CCopasiNode< _Data >::getSibling(), CEvaluationNode::getType(), CEvaluationNodeCall::INVALID, CNormalCall::INVALID, CNormalCall::setName(), CNormalCall::setType(), CEvaluationNode::subType(), and CEvaluationNode::type().

Referenced by createItemPower(), and createItemPowerItem().

 {
   CNormalCall* pCall = NULL;
   CEvaluationNode::Type type = CEvaluationNode::type(node->getType());
 
   if (type == CEvaluationNode::CALL || type == CEvaluationNode::DELAY)
     {
       // create a call object and add all children
       pCall = new CNormalCall();
       // set the name
       pCall->setName(node->getData());
       const CEvaluationNode* pChild = dynamic_cast<const CEvaluationNode*>(node->getChild());
 
       while (pChild != NULL)
         {
           CNormalFraction* pFraction = createNormalRepresentation(pChild);
           assert(pFraction != NULL);
           pCall->add(*pFraction);
           delete pFraction;
           pChild = dynamic_cast<const CEvaluationNode*>(pChild->getSibling());
         }
 
       if (type == CEvaluationNode::DELAY)
         {
           pCall->setType(CNormalCall::DELAY);
         }
       else
         {
           CEvaluationNodeCall::SubType subType = (CEvaluationNodeCall::SubType)CEvaluationNode::subType(node->getType());
 
           switch (subType)
             {
               case CEvaluationNodeCall::EXPRESSION:
                 pCall->setType(CNormalCall::EXPRESSION);
                 break;
               case CEvaluationNodeCall::FUNCTION:
                 pCall->setType(CNormalCall::FUNCTION);
                 break;
               case CEvaluationNodeCall::INVALID:
                 pCall->setType(CNormalCall::INVALID);
                 break;
             }
         }
     }
 
   return pCall;
 }

CNormalChoice* createChoice ( const CEvaluationNode * pNode )

This version of the method is wrong. CEvaluationNode* convertToCEvaluationNode(const CNormalLogical& logical) { CEvaluationNode* pResult = NULL; bool error = false; CEvaluationNodeFunction* pNot = NULL; CNormalLogical::ItemSetOfSets::const_iterator it = logical.getAndSets().begin(), endit = logical.getAndSets().end(); std::vector<CEvaluationNode*> andItems; while (it != endit && error == false) { CEvaluationNode* pTmpResult = NULL; if ((it).first.size() == 1) { pTmpResult = convertToCEvaluationNode(*(*(*it).first.begin()).first); if ((*(*it).first.begin()).second == true) { pNot = new CEvaluationNodeFunction(CEvaluationNodeFunction::NOT, "NOT"); pNot->addChild(pTmpResult); pTmpResult = pNot; } } else if ((*it).first.size() > 1) { CNormalLogical::ItemSet::const_iterator innerit = (*it).first.begin(), innerendit = (*it).first.end(); create the not CEvaluationNode pOrNode = new CEvaluationNodeLogical(CEvaluationNodeLogical::OR, "OR"); pTmpResult = pOrNode; while (innerit != innerendit && error == false) { CEvaluationNode* pItem = convertToCEvaluationNode(*(innerit).first); if (pItem != NULL) { if ((*innerit).second == true) { pNot = new CEvaluationNodeFunction(CEvaluationNodeFunction::NOT, "NOT"); pNot->addChild(pItem); pItem = pNot; } pOrNode->addChild(pItem); pOrNode->addChild(new CEvaluationNodeLogical(CEvaluationNodeLogical::OR, "OR")); pOrNode = dynamic_cast<CEvaluationNodeLogical>(pOrNode->getChild()->getSibling()); } else { error = true; delete pTmpResult; pTmpResult = NULL; } ++innerit; } if (error == false) { replace the parent of the current pOrNode with its first child; CEvaluationNode * pGrandParent = dynamic_cast<CEvaluationNode*>(pOrNode->getParent()); pGrandParent->removeChild(pOrNode->getParent()); pOrNode = dynamic_cast<CEvaluationNode*>(pOrNode->getParent()); CEvaluationNode* pChild = dynamic_cast<CEvaluationNode*>(pOrNode->getChild()); pOrNode->removeChild(pChild); pGrandParent->addChild(pChild); delete pOrNode; } } if (pTmpResult != NULL) { if ((it).second == true) { pNot = new CEvaluationNodeFunction(CEvaluationNodeFunction::NOT, "NOT"); pNot->addChild(pTmpResult); } andItems.push_back(pTmpResult); } ++it; } CNormalLogical::ChoiceSetOfSets::const_iterator it2 = logical.getChoices().begin(), endit2 = logical.getChoices().end(); while (it2 != endit2 && error == false) { CEvaluationNode pTmpResult = NULL; if ((it2).first.size() == 1) { pTmpResult = convertToCEvaluationNode(*(*(*it2).first.begin()).first); if ((*(*it2).first.begin()).second == true) { pNot = new CEvaluationNodeFunction(CEvaluationNodeFunction::NOT, "NOT"); pNot->addChild(pTmpResult); pTmpResult = pNot; } } else if ((*it2).first.size() > 1) { CNormalLogical::ChoiceSet::const_iterator innerit = (*it2).first.begin(), innerendit = (*it2).first.end(); create the not CEvaluationNode pOrNode = new CEvaluationNodeLogical(CEvaluationNodeLogical::OR, "OR"); pTmpResult = pOrNode; while (innerit != innerendit && error != true) { CEvaluationNode* pItem = convertToCEvaluationNode(*(innerit).first); if (pItem != NULL) { if ((*innerit).second == true) { pNot = new CEvaluationNodeFunction(CEvaluationNodeFunction::NOT, "NOT"); pNot->addChild(pItem); pItem = pNot; } pOrNode->addChild(pItem); pOrNode->addChild(new CEvaluationNodeLogical(CEvaluationNodeLogical::OR, "OR")); pOrNode = dynamic_cast<CEvaluationNodeLogical>(pOrNode->getChild()->getSibling()); } else { delete pTmpResult; pTmpResult = NULL; error = true; } ++innerit; } replace the parent of the current pOrNode with its first child; if (error == false) { CEvaluationNode * pGrandParent = dynamic_cast<CEvaluationNode*>(pOrNode->getParent()); pGrandParent->removeChild(pOrNode->getParent()); pOrNode = dynamic_cast<CEvaluationNode*>(pOrNode->getParent()); CEvaluationNode* pChild = dynamic_cast<CEvaluationNode*>(pOrNode->getChild()); pOrNode->removeChild(pChild); pGrandParent->addChild(pChild); delete pOrNode; } } if (pTmpResult != NULL) { if ((it2).second == true) { pNot = new CEvaluationNodeFunction(CEvaluationNodeFunction::NOT, "NOT"); pNot->addChild(pTmpResult); } andItems.push_back(pTmpResult); } ++it2; } if (!andItems.empty()) { if (error == false) { if (andItems.size() == 1) { pResult = *andItems.begin(); } else { std::vector<CEvaluationNode>::iterator andit = andItems.begin(), andendit = andItems.end(); CEvaluationNode* pAndNode = new CEvaluationNodeLogical(CEvaluationNodeLogical::AND, "AND"); pResult = pAndNode; while (andit != andendit) { pAndNode->addChild(<em>andit); pAndNode->addChild(new CEvaluationNodeLogical(CEvaluationNodeLogical::AND, "AND")); pAndNode = dynamic_cast<CEvaluationNodeLogical>(pAndNode->getChild()->getSibling()); ++andit; } replace the parent of the current pAndNode with its first child; CEvaluationNode* pGrandParent = dynamic_cast<CEvaluationNode*>(pAndNode->getParent()); pGrandParent->removeChild(pAndNode->getParent()); pAndNode = dynamic_cast<CEvaluationNode*>(pAndNode->getParent()); CEvaluationNode* pChild = dynamic_cast<CEvaluationNode*>(pAndNode->getChild()); pAndNode->removeChild(pChild); pGrandParent->addChild(pChild); delete pAndNode; } } else { an error has occured, we need to free the memory taken by all andItems std::vector<CEvaluationNode*>::iterator andit = andItems.begin(), andendit = andItems.end(); while (andit != andendit) { delete (*andit); ++andit; } } } if (pResult != NULL) { if (logical.isNegated() == true) { pNot = new CEvaluationNodeFunction(CEvaluationNodeFunction::NOT, "NOT"); pNot->addChild(pResult); pResult = pNot; } pResult->compile(NULL); } return pResult; }

Definition at line 2245 of file ConvertToCEvaluationNode.cpp.

References CEvaluationNode::CHOICE, createLogical(), createNormalRepresentation(), CCopasiNode< _Data >::getChild(), CEvaluationNode::getType(), pResult, CNormalChoice::setCondition(), CNormalChoice::setFalseExpression(), CNormalChoice::setTrueExpression(), and CEvaluationNode::type().

Referenced by createItemPower(), and createItemPowerItem().

 {
   CNormalChoice* pResult = NULL;
 
   if (pNode != NULL && CEvaluationNode::type(pNode->getType()) == CEvaluationNode::CHOICE)
     {
       CNormalLogical* pLogical = createLogical(dynamic_cast<const CEvaluationNode*>(pNode->getChild()));
 
       if (pLogical != NULL)
         {
           CNormalFraction* pTrueFraction = createNormalRepresentation(dynamic_cast<const CEvaluationNode*>(pNode->getChild()->getSibling()));
 
           if (pTrueFraction != NULL)
             {
               // the false branch is optional
               const CEvaluationNode* pFalseBranch = dynamic_cast<const CEvaluationNode*>(pNode->getChild()->getSibling()->getSibling());
 
               if (pFalseBranch != NULL)
                 {
                   CNormalFraction* pFalseFraction = createNormalRepresentation(pFalseBranch);
 
                   if (pFalseFraction != NULL)
                     {
                       pResult = new CNormalChoice();
                       pResult->setCondition(*pLogical);
                       pResult->setTrueExpression(*pTrueFraction);
                       pResult->setFalseExpression(*pFalseFraction);
                       delete pFalseFraction;
                     }
                 }
               else
                 {
                   CNormalFraction* pFalseFraction = new CNormalFraction();
                   pResult = new CNormalChoice();
                   pResult->setCondition(*pLogical);
                   pResult->setTrueExpression(*pTrueFraction);
                   pResult->setFalseExpression(*pFalseFraction);
                   delete pFalseFraction;
                 }
 
               delete pTrueFraction;
             }
 
           delete pLogical;
         }
     }
 
   return pResult;
 }

CNormalFraction* createFraction ( const CEvaluationNode * node )

Create a fraction from a node. node is necessarily a DIVIDE operator if not the root node of the tree.

Returns: CNormalFraction*, pointer to newly created fraction.

Definition at line 308 of file ConvertToCEvaluationNode.cpp.

References CNormalSum::add(), CNormalFraction::cancel(), CNormalTranslation::createChain(), createSum(), CCopasiNode< _Data >::getData(), CNormalTranslation::ONE_NODE, CNormalFraction::setDenominator(), CNormalFraction::setNumerator(), CNormalTranslation::splitProduct(), and CNormalTranslation::TIMES_NODE.

Referenced by createItemPowerItem(), and createNormalRepresentation().

 {
   CNormalFraction* pFraction = new CNormalFraction();
 
   if (node->getData() == "/")
     {// always executed except on root node possibly not
       // find a product chain and create new temporary nodes for the
       // numerator and denominator which are then converted to sums
       std::vector<const CEvaluationNode*> multiplications, divisions;
       CNormalTranslation::splitProduct(node, multiplications, divisions, false);
       std::vector<CEvaluationNode*> tmp;
       std::vector<const CEvaluationNode*>::const_iterator it = multiplications.begin(), endit = multiplications.end();
 
       while (it != endit)
         {
           tmp.push_back((*it)->copyBranch());
           ++it;
         }
 
       CEvaluationNode* pTmpNode = CNormalTranslation::createChain(&CNormalTranslation::TIMES_NODE, &CNormalTranslation::ONE_NODE, tmp);
       assert(pTmpNode != NULL);
       CNormalSum* pNum = createSum(pTmpNode);
       assert(pNum != NULL);
       delete pTmpNode;
       tmp.clear();
       it = divisions.begin();
       endit = divisions.end();
 
       while (it != endit)
         {
           tmp.push_back((*it)->copyBranch());
           ++it;
         }
 
       pTmpNode = CNormalTranslation::createChain(&CNormalTranslation::TIMES_NODE, &CNormalTranslation::ONE_NODE, tmp);
       assert(pTmpNode != NULL);
       CNormalSum* pDenom = createSum(pTmpNode);
       assert(pDenom != NULL);
       delete pTmpNode;
       pFraction->setNumerator(*pNum);
       pFraction->setDenominator(*pDenom);
       //pFraction->cancel();
       delete pNum;
       delete pDenom;
     }
   else
     {// only possible for root node
       CNormalSum* pNum = createSum(node);
       CNormalSum* pDenom = new CNormalSum();
       CNormalProduct* pProduct = new CNormalProduct();
       pDenom->add(*pProduct);
       pFraction->setNumerator(*pNum);
       pFraction->setDenominator(*pDenom);
       pFraction->cancel();
       delete pProduct;
       delete pNum;
       delete pDenom;
     }
 
   return pFraction;
 };

CNormalFunction* createFunction ( const CEvaluationNode * node )

Create an function from an evaluation node.

Returns: CNormalFunction*, pointer to newly created function

Definition at line 1503 of file ConvertToCEvaluationNode.cpp.

Referenced by createItemPower(), and createItemPowerItem().

 {
   CNormalFunction* pFun = NULL;
   CEvaluationNode::Type type = CEvaluationNode::type(node->getType());
 
   if (type == CEvaluationNode::FUNCTION)
     {
       CNormalFunction::Type type = CNormalFunction::INVALID;
 
       switch ((CEvaluationNodeFunction::SubType)CEvaluationNode::subType(node->getType()))
         {
           case CEvaluationNodeFunction::LOG:
             type = CNormalFunction::LOG;
             break;
           case CEvaluationNodeFunction::LOG10:
             type = CNormalFunction::LOG10;
             break;
           case CEvaluationNodeFunction::EXP:
             type = CNormalFunction::EXP;
             break;
           case CEvaluationNodeFunction::SIN:
             type = CNormalFunction::SIN;
             break;
           case CEvaluationNodeFunction:: COS:
             type = CNormalFunction::COS;
             break;
           case CEvaluationNodeFunction::TAN:
             type = CNormalFunction::TAN;
             break;
           case CEvaluationNodeFunction::SEC:
             type = CNormalFunction::SEC;
             break;
           case CEvaluationNodeFunction::CSC:
             type = CNormalFunction::CSC;
             break;
           case CEvaluationNodeFunction::COT:
             type = CNormalFunction::COT;
             break;
           case CEvaluationNodeFunction::SINH:
             type = CNormalFunction::SINH;
             break;
           case CEvaluationNodeFunction::COSH:
             type = CNormalFunction::COSH;
             break;
           case CEvaluationNodeFunction::TANH:
             type = CNormalFunction::TANH;
             break;
           case CEvaluationNodeFunction::SECH:
             type = CNormalFunction::SECH;
             break;
           case CEvaluationNodeFunction::CSCH:
             type = CNormalFunction::CSCH;
             break;
           case CEvaluationNodeFunction::COTH:
             type = CNormalFunction::COTH;
             break;
           case CEvaluationNodeFunction::ARCSIN:
             type = CNormalFunction::ARCSIN;
             break;
           case CEvaluationNodeFunction::ARCCOS:
             type = CNormalFunction::ARCCOS;
             break;
           case CEvaluationNodeFunction::ARCTAN:
             type = CNormalFunction::ARCTAN;
             break;
           case CEvaluationNodeFunction::ARCSEC:
             type = CNormalFunction::ARCSEC;
             break;
           case CEvaluationNodeFunction::ARCCSC:
             type = CNormalFunction::ARCCSC;
             break;
           case CEvaluationNodeFunction::ARCCOT:
             type = CNormalFunction::ARCCOT;
             break;
           case CEvaluationNodeFunction::ARCSINH:
             type = CNormalFunction::ARCSINH;
             break;
           case CEvaluationNodeFunction::ARCCOSH:
             type = CNormalFunction::ARCCOSH;
             break;
           case CEvaluationNodeFunction::ARCTANH:
             type = CNormalFunction::ARCTANH;
             break;
           case CEvaluationNodeFunction::ARCSECH:
             type = CNormalFunction::ARCSECH;
             break;
           case CEvaluationNodeFunction::ARCCSCH:
             type = CNormalFunction::ARCCSCH;
             break;
           case CEvaluationNodeFunction::ARCCOTH:
             type = CNormalFunction::ARCCOTH;
             break;
           case CEvaluationNodeFunction::SQRT:
             type = CNormalFunction::SQRT;
             break;
           case CEvaluationNodeFunction::ABS:
             type = CNormalFunction::ABS;
             break;
           case CEvaluationNodeFunction::FLOOR:
             type = CNormalFunction::FLOOR;
             break;
           case CEvaluationNodeFunction::CEIL:
             type = CNormalFunction::CEIL;
             break;
           case CEvaluationNodeFunction::FACTORIAL:
             type = CNormalFunction::FACTORIAL;
             break;
           case CEvaluationNodeFunction::INVALID:
           default:
             assert(false);
             break;
         }
 
       if (type != CNormalFunction::INVALID)
         {
           CNormalFraction* pFrac = createNormalRepresentation(dynamic_cast<const CEvaluationNode*>(node->getChild()));
 
           if (pFrac != NULL)
             {
               pFun = new CNormalFunction();
               pFun->setType(type);
               pFun->setFraction(*pFrac);
               delete pFrac;
             }
         }
     }
 
   return pFun;
 }

CNormalGeneralPower* createGeneralPower ( const CEvaluationNode * node )

Create an general power from an evaluation node. Node MUST BE a POWER operator with second child (exponent) of type not NUMBER !!

Returns: CNormalItemPower*, pointer to newly created general power (or modulo).

Create an general power from an evaluation node. if it is a power node, create a for each child and set base and exponent, else only create the base from the node and set the exponent to 1

Returns: CNormalItemPower*, pointer to newly created general power (or modulo).

Definition at line 1378 of file ConvertToCEvaluationNode.cpp.

References createNormalRepresentation(), CEvaluationNodeNumber::DOUBLE, CCopasiNode< _Data >::getChild(), CEvaluationNode::getType(), CNormalGeneralPower::MODULO, CEvaluationNodeOperator::MODULUS, CEvaluationNode::OPERATOR, CNormalGeneralPower::POWER, CEvaluationNodeOperator::POWER, CNormalGeneralPower::setLeft(), CNormalGeneralPower::setRight(), CNormalGeneralPower::setType(), CEvaluationNode::subType(), and CEvaluationNode::type().

Referenced by createItemPower(), and createItemPowerItem().

 {
   CNormalGeneralPower* pPow = NULL;
 
   if (node != NULL)
     {
       if (CEvaluationNode::type(node->getType()) == CEvaluationNode::OPERATOR)
         {
           if (((CEvaluationNodeOperator::SubType)CEvaluationNode::subType(node->getType())) == CEvaluationNodeOperator::POWER)
             {
               pPow = new CNormalGeneralPower();
               pPow->setType(CNormalGeneralPower::POWER);
             }
           else if (((CEvaluationNodeOperator::SubType)CEvaluationNode::subType(node->getType())) == CEvaluationNodeOperator::MODULUS)
             {
               pPow = new CNormalGeneralPower();
               pPow->setType(CNormalGeneralPower::MODULO);
             }
 
           if (pPow != NULL)
             {
               // add the left and the right side
               CNormalFraction* pBase = createNormalRepresentation(dynamic_cast<const CEvaluationNode*>(node->getChild()));
               CNormalFraction* pExponent = createNormalRepresentation(dynamic_cast<const CEvaluationNode*>(node->getChild()->getSibling()));
               assert(pBase != NULL);
               assert(pExponent != NULL);
               pPow->setLeft(*pBase);
               pPow->setRight(*pExponent);
               delete pBase;
               delete pExponent;
             }
           else
             {
               // create a fraction for the base and a unit fraction for the exponent
               pPow = new CNormalGeneralPower();
               pPow->setType(CNormalGeneralPower::POWER);
               CNormalFraction* pBase = createNormalRepresentation(dynamic_cast<const CEvaluationNode*>(node));
               CEvaluationNode* pTmpNode = new CEvaluationNodeNumber(CEvaluationNodeNumber::DOUBLE, "1.0");
               CNormalFraction* pExponent = createNormalRepresentation(pTmpNode);
               delete pTmpNode;
               assert(pBase != NULL);
               assert(pExponent != NULL);
               pPow->setLeft(*pBase);
               pPow->setRight(*pExponent);
               delete pBase;
               delete pExponent;
             }
         }
       else
         {
           // create a fraction for the base and a unit fraction for the exponent
           pPow = new CNormalGeneralPower();
           pPow->setType(CNormalGeneralPower::POWER);
           CNormalFraction* pBase = createNormalRepresentation(dynamic_cast<const CEvaluationNode*>(node));
           CEvaluationNode* pTmpNode = new CEvaluationNodeNumber(CEvaluationNodeNumber::DOUBLE, "1.0");
           CNormalFraction* pExponent = createNormalRepresentation(pTmpNode);
           delete pTmpNode;
           assert(pBase != NULL);
           assert(pExponent != NULL);
           pPow->setLeft(*pBase);
           pPow->setRight(*pExponent);
           delete pBase;
           delete pExponent;
         }
     }
 
   return pPow;
 }

CNormalItem* createItem ( const CEvaluationNode * pNode )

Create an item from an evaluation node that need not be of specific type.

Returns: CNormalItem*, pointer to newly created item.

Definition at line 374 of file ConvertToCEvaluationNode.cpp.

References CEvaluationNode::buildInfix(), CNormalItem::CONSTANT, CEvaluationNode::CONSTANT, CEvaluationNode::getType(), CEvaluationNode::type(), CNormalItem::VARIABLE, and CEvaluationNode::VARIABLE.

Referenced by createItemPower(), and createItemPowerItem().

 {
   CNormalItem* pItem = NULL;
 
   switch (CEvaluationNode::type(pNode->getType()))
     {
       case CEvaluationNode::VARIABLE:
         pItem = new CNormalItem(pNode->buildInfix(), CNormalItem::VARIABLE);
         break;
       case CEvaluationNode::CONSTANT:
         pItem = new CNormalItem(pNode->buildInfix(), CNormalItem::CONSTANT);
         break;
         /*
         case CEvaluationNode::OPERATOR:
         if (((CEvaluationNodeOperator::SubType)CEvaluationNode::subType(pNode->getType()))==CEvaluationNodeOperator::POWER
             || ((CEvaluationNodeOperator::SubType)CEvaluationNode::subType(pNode->getType()))==CEvaluationNodeOperator::MODULUS)
         {
             // calling createItem will add the wrong pNode for the
             // string into the lookup table
             CNormalItem * child1 = createItem(dynamic_cast<const CEvaluationNode*>(pNode->getChild()));
             CNormalItem * child2 = createItem(dynamic_cast<const CEvaluationNode*>(pNode->getChild()->getSibling()));
 
             std::stringstream tmp;
             if ((child1->getType() == CNormalItem::VARIABLE) || (child1->getType() == CNormalItem::CONSTANT))
                 tmp << *child1;
             else
             {
                 tmp << "(" << *child1 << ")";
             }
             tmp << pNode->getData();
             if ((child2->getType() == CNormalItem::VARIABLE) || (child2->getType() == CNormalItem::CONSTANT))
                 tmp << *child2;
             else
             {
                 tmp << "(" << *child2 << ")";
             }
 
             pItem = new CNormalItem(tmp.str(), CNormalItem::FUNCTION);
             std::string s=tmp.str();
             std::map<std::string,const CEvaluationNode*>::const_iterator pos=str2eval.find(tmp.str());
             if(pos==str2eval.end())
             {
                 CEvaluationNodeOperator::SubType type=CEvaluationNodeOperator::POWER;
                 std::string s("^");
                 if(pNode->getData()=="%")
                 {
                     type=CEvaluationNodeOperator::MODULUS;
                     s="%";
                 }
                 CEvaluationNodeOperator* pNode=new CEvaluationNodeOperator(type,s);
                 // when the child pNode is converted, the wrongly added
                 // pNode from above is returned
 
                 pNode->addChild(convertToCEvaluationNode(*child1));
                 pNode->addChild(convertToCEvaluationNode(*child2));
                 //pNode->printRecursively(std::cerr);
                 str2eval[tmp.str()]=pNode;
             }
             delete child1;
             delete child2;
         }
         else // can be called only by createItem('OPERATOR pNode')
         {
             CNormalFraction * normedNode = createFraction(pNode);
             normedNode->simplify();
             CEvaluationNode* pTmpNode=convertToCEvaluationNode(*normedNode);
             std::stringstream tmp;
             tmp << *normedNode;
             pItem = new CNormalItem(tmp.str(), CNormalItem::FUNCTION);
             std::map<std::string,const CEvaluationNode*>::const_iterator pos=str2eval.find(normedNode->toString());
             if(pos==str2eval.end())
             {
                 str2eval[tmp.str()]=pTmpNode;
             }
 
             delete normedNode;
         }
         break;
         */
       default:   //cases CALL, CHOICE, LOGICAL, OBJECT, VECTOR.  NUMBER should not occur!
         throw std::exception();
         break;
     }
 
   return pItem;
 }

CNormalItemPower* createItemPower ( const CEvaluationNode * node )

Create an item power from an evaluation node. Node MUST BE a POWER operator with second child (exponent) of type NUMBER !!

Returns: CNormalItemPower*, pointer to newly created item power.

Create an item power from an evaluation node. Depending on the given node, we have to create different items. Either we have to create a CNormalItem, a CNormalGeneralPower, a function or a choice.

Definition at line 467 of file ConvertToCEvaluationNode.cpp.

References CCopasiNode< _Data >::addChild(), C_FLOAT64, CEvaluationNode::CALL, CEvaluationNode::CHOICE, CEvaluationNode::CONSTANT, createCall(), createChoice(), createFunction(), createGeneralPower(), createItem(), createItemPower(), createItemPowerItem(), CEvaluationNode::DELAY, CEvaluationNodeNumber::DOUBLE, CEvaluationNode::FUNCTION, CCopasiNode< _Data >::getChild(), CEvaluationNode::getType(), CEvaluationNode::LOGICAL, CEvaluationNodeFunction::MINUS, CEvaluationNodeOperator::MULTIPLY, CEvaluationNode::NUMBER, CEvaluationNode::OBJECT, CEvaluationNode::OPERATOR, CEvaluationNodeOperator::POWER, CNormalItemPower::setExp(), CNormalItemPower::setItem(), CEvaluationNode::subType(), CEvaluationNode::type(), and CEvaluationNode::VARIABLE.

Referenced by createItemPower(), createItemPowerItem(), and createProduct().

 {
   // make sure we create the correct item for the given node
   CNormalItemPower * pItemPower = new CNormalItemPower();
 
   if (CEvaluationNode::type(node->getType()) == CEvaluationNode::OPERATOR && (CEvaluationNodeOperator::SubType)CEvaluationNode::subType(node->getType()) == CEvaluationNodeOperator::POWER)
     {
       // check if the second child is a number
       if (CEvaluationNode::type(dynamic_cast<const CEvaluationNode*>(node->getChild()->getSibling())->getType()) == CEvaluationNode::NUMBER)
         {
           // we set the exponent to that number
           pItemPower->setExp((C_FLOAT64)dynamic_cast<const CEvaluationNodeNumber*>(node->getChild()->getSibling())->getValue());
           // check if we can create a CNormalItem object for the fist child, else we
           // create a general power with exponent 1
           CEvaluationNode::Type type = CEvaluationNode::type(dynamic_cast<const CEvaluationNode*>(node->getChild())->getType());
 
           if (type == CEvaluationNode::CONSTANT || type == CEvaluationNode::OBJECT || type == CEvaluationNode::VARIABLE || type == CEvaluationNode::FUNCTION || type == CEvaluationNode::CHOICE || type == CEvaluationNode::CALL || type == CEvaluationNode::LOGICAL)
             {
               CNormalBase* pItem = createItemPowerItem(dynamic_cast<const CEvaluationNode*>(node->getChild()));
               assert(pItem != NULL);
               pItemPower->setItem(*pItem);
               delete pItem;
             }
           else
             {
               // create a general power for the child node
               CNormalGeneralPower* pGeneralPower = createGeneralPower(dynamic_cast<const CEvaluationNode*>(node->getChild()));
               assert(pGeneralPower != NULL);
               pItemPower->setItem(*pGeneralPower);
               delete pGeneralPower;
             }
         }
       else
         {
           // create a general power for the node
           CNormalGeneralPower* pGeneralPower = createGeneralPower(node);
           assert(pGeneralPower != NULL);
           pItemPower->setItem(*pGeneralPower);
           pItemPower->setExp(1.0);
           delete pGeneralPower;
         }
     }
   else if (CEvaluationNode::type(node->getType()) == CEvaluationNode::FUNCTION)
     {
       if (((CEvaluationNodeFunction::SubType)CEvaluationNode::subType(node->getType())) == CEvaluationNodeFunction::MINUS)
         {
           // multiply the second node by -1 and call createItempower with the
           // result
           CEvaluationNodeNumber* pNumberNode = new CEvaluationNodeNumber(CEvaluationNodeNumber::DOUBLE, "-1.0");
           CEvaluationNodeOperator* pOperatorNode = new CEvaluationNodeOperator(CEvaluationNodeOperator::MULTIPLY, "*");
           pOperatorNode->addChild(pNumberNode);
           pOperatorNode->addChild(dynamic_cast<const CEvaluationNode*>(node->getChild())->copyBranch());
           // delete the old item power object
           delete pItemPower;
           // create a new item power object by calling this method again with
           // the newly created multiplication node
           pItemPower = createItemPower(pOperatorNode);
           // delete the multiplication node
           delete pOperatorNode;
         }
       else
         {
           CNormalFunction* pFunction = createFunction(node);
           assert(pFunction != NULL);
           pItemPower->setItem(*pFunction);
           delete pFunction;
           pItemPower->setExp(1.0);
         }
     }
   else if (CEvaluationNode::type(node->getType()) == CEvaluationNode::CALL || CEvaluationNode::type(node->getType()) == CEvaluationNode::DELAY)
     {
       CNormalCall* pCall = createCall(node);
       assert(pCall != NULL);
       pItemPower->setItem(*pCall);
       delete pCall;
       pItemPower->setExp(1.0);
     }
   else if (CEvaluationNode::type(node->getType()) == CEvaluationNode::CHOICE)
     {
       CNormalChoice* pChoice = createChoice(node);
       assert(pChoice != NULL);
       pItemPower->setItem(*pChoice);
       delete pChoice;
       pItemPower->setExp(1.0);
     }
   else if (CEvaluationNode::type(node->getType()) == CEvaluationNode::CONSTANT || CEvaluationNode::type(node->getType()) == CEvaluationNode::OBJECT || CEvaluationNode::type(node->getType()) == CEvaluationNode::VARIABLE)
     {
       CNormalItem* pItem = createItem(node);
       assert(pItem != NULL);
       pItemPower->setItem(*pItem);
       delete pItem;
       pItemPower->setExp(1.0);
     }
   else if (CEvaluationNode::type(node->getType()) == CEvaluationNode::LOGICAL)
     {
       CNormalBase* pItem = createItemPowerItem(node);
       assert(pItem != NULL);
       pItemPower->setItem(*pItem);
       delete pItem;
       pItemPower->setExp(1.0);
     }
   else
     {
       // create a general power for the node
       CNormalGeneralPower* pGeneralPower = createGeneralPower(node);
       assert(pGeneralPower != NULL);
       pItemPower->setItem(*pGeneralPower);
       pItemPower->setExp(1.0);
       delete pGeneralPower;
     }
 
   return pItemPower;
 }

CNormalBase* createItemPowerItem ( const CEvaluationNode * pNode )

Definition at line 1660 of file ConvertToCEvaluationNode.cpp.

References CEvaluationNode::CALL, CEvaluationNode::CHOICE, CEvaluationNode::CONSTANT, createCall(), createChoice(), createFraction(), createFunction(), createGeneralPower(), createItem(), createItemPower(), createLogical(), CEvaluationNode::FUNCTION, CCopasiNode< _Data >::getChild(), CEvaluationNode::getType(), CEvaluationNode::LOGICAL, CEvaluationNodeOperator::MODULUS, CEvaluationNode::NUMBER, CEvaluationNode::OPERATOR, CEvaluationNodeOperator::POWER, pResult, CEvaluationNode::subType(), CEvaluationNode::type(), and CEvaluationNode::VARIABLE.

Referenced by createItemPower().

 {
   CNormalBase* pResult = NULL;
   CEvaluationNode::Type type = CEvaluationNode::type(pNode->getType());
   CEvaluationNodeOperator::SubType subType;
 
   switch (type)
     {
       case CEvaluationNode::OPERATOR:
         subType = (CEvaluationNodeOperator::SubType)CEvaluationNode::subType(pNode->getType());
 
         if (subType == CEvaluationNodeOperator::POWER)
           {
             if ((CEvaluationNode::type(dynamic_cast<const CEvaluationNode*>(pNode->getChild()->getSibling())->getType())) == CEvaluationNode::NUMBER)
               {
                 pResult = createItemPower(pNode);
               }
             else
               {
                 pResult = createGeneralPower(pNode);
               }
           }
         else if (CEvaluationNodeOperator::MODULUS)
           {
             pResult = createGeneralPower(pNode);
           }
         else
           {
             pResult = createFraction(pNode);
           }
 
         break;
       case CEvaluationNode::NUMBER:
         pResult = createItem(pNode);
         break;
       case CEvaluationNode::FUNCTION:
         pResult = createFunction(pNode);
         break;
       case CEvaluationNode::CALL:
         pResult = createCall(pNode);
         break;
       case CEvaluationNode::CHOICE:
         pResult = createChoice(pNode);
         break;
       case CEvaluationNode::LOGICAL:
         pResult = createLogical(pNode);
         break;
       case CEvaluationNode::VARIABLE:
         pResult = createItem(pNode);
         break;
       case CEvaluationNode::CONSTANT:
         pResult = createItem(pNode);
         break;
       default:
         break;
     }
 
   return pResult;
 }

CNormalLogical* createLogical ( const CEvaluationNode * pNode )

Definition at line 2345 of file ConvertToCEvaluationNode.cpp.

Referenced by createChoice(), createItemPowerItem(), createLogical(), and createLogicalChoice().

 {
   CNormalLogical* pResult = NULL;
 
   if (pNode != NULL)
     {
       CEvaluationNode::Type type = pNode->getType();
 
       if (CEvaluationNode::type(type) == CEvaluationNode::CONSTANT)
         {
           if (((CEvaluationNodeConstant::SubType)CEvaluationNode::subType(type)) == CEvaluationNodeConstant::TRUE || ((CEvaluationNodeConstant::SubType)CEvaluationNode::subType(type)) == CEvaluationNodeConstant::FALSE)
             {
               CNormalLogicalItem* pLogicalItem = createLogicalItem(pNode);
 
               if (pLogicalItem != NULL)
                 {
                   pResult = new CNormalLogical();
                   CNormalLogical::ItemSet tmp;
                   tmp.insert(std::make_pair(pLogicalItem, false));
                   pResult->getAndSets().insert(std::make_pair(tmp, false));
                 }
             }
         }
       else if (CEvaluationNode::type(type) == CEvaluationNode::LOGICAL)
         {
           CNormalLogicalItem* pLogicalItem = NULL;
           CEvaluationNode* pOrNode = NULL;
           CEvaluationNode* pAndNode = NULL;
           const CEvaluationNode* pA = NULL;
           const CEvaluationNode* pB = NULL;
           CEvaluationNode* pNotNode = NULL;
           CNormalLogical* pLeftLogical = NULL;
           CNormalLogical::ChoiceSetOfSets::const_iterator it, endit;
           CNormalLogical::ItemSetOfSets::const_iterator it2, endit2;
 
           switch ((CEvaluationNodeLogical::SubType)CEvaluationNode::subType(type))
             {
               case CEvaluationNodeLogical::EQ:
               case CEvaluationNodeLogical::NE:
               case CEvaluationNodeLogical::GT:
               case CEvaluationNodeLogical::LT:
               case CEvaluationNodeLogical::GE:
               case CEvaluationNodeLogical::LE:
                 pLogicalItem = createLogicalItem(pNode);
 
                 if (pLogicalItem != NULL)
                   {
                     pResult = new CNormalLogical();
                     CNormalLogical::ItemSet tmp;
                     tmp.insert(std::make_pair(pLogicalItem, false));
                     pResult->getAndSets().insert(std::make_pair(tmp, false));
                   }
 
                 break;
               case CEvaluationNodeLogical::XOR:
                 // replace A xor B by A OR B AND NOT(A AND B)
                 pA = dynamic_cast<const CEvaluationNode*>(pNode->getChild());
                 pB = dynamic_cast<const CEvaluationNode*>(pA->getSibling());
                 pAndNode = new CEvaluationNodeLogical(CEvaluationNodeLogical::AND, "AND");
                 pAndNode->addChild(pA->copyBranch());
                 pAndNode->addChild(pB->copyBranch());
                 pNotNode = new CEvaluationNodeFunction(CEvaluationNodeFunction::NOT, "NOT");
                 pNotNode->addChild(pAndNode);
                 pOrNode = new CEvaluationNodeLogical(CEvaluationNodeLogical::OR, "OR");
                 pOrNode->addChild(pA->copyBranch());
                 pOrNode->addChild(pB->copyBranch());
                 pAndNode = new CEvaluationNodeLogical(CEvaluationNodeLogical::AND, "AND");
                 pAndNode->addChild(pOrNode);
                 pAndNode->addChild(pNotNode);
                 pResult = createLogical(pAndNode);
                 delete pAndNode;
                 break;
               case CEvaluationNodeLogical::AND:
                 pLeftLogical = createLogical(dynamic_cast<const CEvaluationNode*>(pNode->getChild()));
 
                 if (pLeftLogical != NULL)
                   {
                     pLeftLogical->simplify();
                     CNormalLogical* pRightLogical = createLogical(dynamic_cast<const CEvaluationNode*>(pNode->getChild()->getSibling()));
 
                     if (pRightLogical != NULL)
                       {
                         pRightLogical->simplify();
                         // integrate the items into a new logical item
                         // we have to create an AND set for every
                         // combination between two AND Sets from the
                         // two logicals
                         pResult = new CNormalLogical();
                         CNormalLogical::ChoiceSetOfSets::const_iterator it = pLeftLogical->getChoices().begin(), endit = pLeftLogical->getChoices().end();
 
                         while (it != endit)
                           {
                             CNormalLogical::ChoiceSetOfSets::const_iterator it2 = pRightLogical->getChoices().begin(), endit2 = pRightLogical->getChoices().end();
 
                             if ((*it).second == true) throw std::exception();
 
                             while (it2 != endit2)
                               {
                                 // create an AND set from the two, but only
                                 // if both don't have the NOT flag set
                                 // which should always be the case since we
                                 // called simplify on the items.
                                 if ((*it2).second == true) throw std::exception();
 
                                 CNormalLogical::ChoiceSet tmp;
                                 CNormalLogical::ChoiceSet::const_iterator innerit = (*it).first.begin(), innerendit = (*it).first.end();
 
                                 while (innerit != innerendit)
                                   {
                                     tmp.insert(std::make_pair(new CNormalChoiceLogical(*(*innerit).first), (*innerit).second));
                                     ++innerit;
                                   }
 
                                 innerit = (*it2).first.begin(), innerendit = (*it2).first.end();
 
                                 while (innerit != innerendit)
                                   {
                                     tmp.insert(std::make_pair(new CNormalChoiceLogical(*(*innerit).first), (*innerit).second));
                                     ++innerit;
                                   }
 
                                 pResult->getChoices().insert(std::make_pair(tmp, false));
                                 ++it2;
                               }
 
                             ++it;
                           }
 
                         CNormalLogical::ItemSetOfSets::const_iterator it2 = pLeftLogical->getAndSets().begin(), endit2 = pLeftLogical->getAndSets().end();
 
                         while (it2 != endit2)
                           {
                             CNormalLogical::ItemSetOfSets::const_iterator it3 = pRightLogical->getAndSets().begin(), endit3 = pRightLogical->getAndSets().end();
 
                             if ((*it2).second == true) throw std::exception();
 
                             while (it3 != endit3)
                               {
                                 // create an AND set from the two, but only
                                 // if both don't have the NOT flag set
                                 // which should always be the case since we
                                 // called simplify on the items.
                                 if ((*it3).second == true) throw std::exception();
 
                                 CNormalLogical::ItemSet tmp;
                                 CNormalLogical::ItemSet::const_iterator innerit = (*it2).first.begin(), innerendit = (*it2).first.end();
 
                                 while (innerit != innerendit)
                                   {
                                     tmp.insert(std::make_pair(new CNormalLogicalItem(*(*innerit).first), (*innerit).second));
                                     ++innerit;
                                   }
 
                                 innerit = (*it3).first.begin(), innerendit = (*it3).first.end();
 
                                 while (innerit != innerendit)
                                   {
                                     tmp.insert(std::make_pair(new CNormalLogicalItem(*(*innerit).first), (*innerit).second));
                                     ++innerit;
                                   }
 
                                 pResult->getAndSets().insert(std::make_pair(tmp, false));
                                 ++it3;
                               }
 
                             ++it2;
                           }
 
                         delete pLeftLogical;
                         delete pRightLogical;
                       }
                     else
                       {
                         delete pLeftLogical;
                       }
                   }
 
                 break;
               case CEvaluationNodeLogical::OR:
                 // 1. create two logicals, one for the left side and one for the
                 //    right side
                 pResult = createLogical(dynamic_cast<const CEvaluationNode*>(pNode->getChild()));
                 pLeftLogical = createLogical(dynamic_cast<const CEvaluationNode*>(pNode->getChild()->getSibling()));
                 // 2. add all or items from the second logical to the first and
                 //    delete the second
                 it = pLeftLogical->getChoices().begin(), endit = pLeftLogical->getChoices().end();
 
                 while (it != endit)
                   {
                     CNormalLogical::ChoiceSet tmpSet;
                     CNormalLogical::copySet((*it).first, tmpSet);
                     pResult->getChoices().insert(std::make_pair(tmpSet, (*it).second));
                     ++it;
                   }
 
                 it2 = pLeftLogical->getAndSets().begin(), endit2 = pLeftLogical->getAndSets().end();
 
                 while (it2 != endit2)
                   {
                     CNormalLogical::ItemSet tmpSet;
                     CNormalLogical::copySet((*it2).first, tmpSet);
                     pResult->getAndSets().insert(std::make_pair(tmpSet, (*it2).second));
                     ++it2;
                   }
 
                 delete pLeftLogical;
                 break;
               default:
                 break;
             }
         }
       else if (CEvaluationNode::type(type) == CEvaluationNode::FUNCTION && ((CEvaluationNodeFunction::SubType)CEvaluationNode::subType(type)) == CEvaluationNodeFunction::NOT)
         {
           pResult = createLogical(dynamic_cast<const CEvaluationNode*>(pNode->getChild()));
 
           if (pResult != NULL)
             {
               pResult->negate();
             }
         }
       else if (CEvaluationNode::type(type) == CEvaluationNode::CHOICE)
         {
           CNormalChoiceLogical* pLogicalChoice = createLogicalChoice(pNode);
 
           if (pLogicalChoice != NULL)
             {
               CNormalLogical::ChoiceSet tmp;
               tmp.insert(std::make_pair(pLogicalChoice, false));
               pResult = new CNormalLogical();
               pResult->getChoices().insert(std::make_pair(tmp, false));
             }
         }
     }
 
   return pResult;
 }

CNormalChoiceLogical* createLogicalChoice ( const CEvaluationNode * pNode )

Definition at line 2295 of file ConvertToCEvaluationNode.cpp.

References CEvaluationNode::CHOICE, createLogical(), CCopasiNode< _Data >::getChild(), CEvaluationNode::getType(), pResult, CNormalChoiceLogical::setCondition(), CNormalChoiceLogical::setFalseExpression(), CNormalChoiceLogical::setTrueExpression(), and CEvaluationNode::type().

Referenced by createLogical().

 {
   CNormalChoiceLogical* pResult = NULL;
 
   if (pNode != NULL && CEvaluationNode::type(pNode->getType()) == CEvaluationNode::CHOICE)
     {
       CNormalLogical* pLogical = createLogical(dynamic_cast<const CEvaluationNode*>(pNode->getChild()));
 
       if (pLogical != NULL)
         {
           CNormalLogical* pTrueLogical = createLogical(dynamic_cast<const CEvaluationNode*>(pNode->getChild()->getSibling()));
 
           if (pTrueLogical != NULL)
             {
               // the false branch is optional
               const CEvaluationNode* pFalseBranch = dynamic_cast<const CEvaluationNode*>(pNode->getChild()->getSibling()->getSibling());
 
               if (pFalseBranch != NULL)
                 {
                   CNormalLogical* pFalseLogical = createLogical(pFalseBranch);
 
                   if (pFalseLogical != NULL)
                     {
                       pResult = new CNormalChoiceLogical();
                       pResult->setCondition(*pLogical);
                       pResult->setTrueExpression(*pTrueLogical);
                       pResult->setFalseExpression(*pFalseLogical);
                       delete pFalseLogical;
                     }
                 }
               else
                 {
                   pResult = new CNormalChoiceLogical();
                   pResult->setCondition(*pLogical);
                   pResult->setTrueExpression(*pTrueLogical);
                   CNormalLogical* pFalseLogical = new CNormalLogical();
                   pResult->setFalseExpression(*pFalseLogical);
                   delete pFalseLogical;
                 }
 
               delete pTrueLogical;
             }
 
           delete pLogical;
         }
     }
 
   return pResult;
 }

CNormalLogicalItem* createLogicalItem ( const CEvaluationNode * pNode )

Definition at line 2582 of file ConvertToCEvaluationNode.cpp.

References CEvaluationNode::CONSTANT, createNormalRepresentation(), CNormalLogicalItem::EQ, CEvaluationNodeLogical::EQ, CNormalLogicalItem::FALSE, CEvaluationNodeConstant::FALSE, CNormalLogicalItem::GE, CEvaluationNodeLogical::GE, CCopasiNode< _Data >::getChild(), CEvaluationNode::getType(), CNormalLogicalItem::GT, CEvaluationNodeLogical::GT, CNormalLogicalItem::LE, CEvaluationNodeLogical::LE, CEvaluationNode::LOGICAL, CNormalLogicalItem::LT, CEvaluationNodeLogical::LT, CNormalLogicalItem::NE, CEvaluationNodeLogical::NE, pResult, CNormalLogicalItem::setLeft(), CNormalLogicalItem::setRight(), CNormalLogicalItem::setType(), CEvaluationNode::subType(), CNormalLogicalItem::TRUE, CEvaluationNodeConstant::TRUE, and CEvaluationNode::type().

Referenced by createLogical().

 {
   CNormalLogicalItem* pResult = NULL;
 
   if (pNode != NULL)
     {
       CEvaluationNode::Type type = pNode->getType();
 
       if (CEvaluationNode::type(type) == CEvaluationNode::CONSTANT)
         {
           if (((CEvaluationNodeConstant::SubType)CEvaluationNode::subType(type)) == CEvaluationNodeConstant::FALSE)
             {
               pResult = new CNormalLogicalItem();
               pResult->setType(CNormalLogicalItem::FALSE);
             }
           else if (((CEvaluationNodeConstant::SubType)CEvaluationNode::subType(type)) == CEvaluationNodeConstant::TRUE)
             {
               pResult = new CNormalLogicalItem();
               pResult->setType(CNormalLogicalItem::TRUE);
             }
         }
 
       if (CEvaluationNode::type(type) == CEvaluationNode::LOGICAL)
         {
           CNormalFraction* pFrac1 = NULL;
           CNormalFraction* pFrac2 = NULL;
 
           switch ((CEvaluationNodeLogical::SubType)CEvaluationNode::subType(type))
             {
               case CEvaluationNodeLogical::EQ:
                 pResult = new CNormalLogicalItem();
                 pResult->setType(CNormalLogicalItem::EQ);
                 pFrac1 = createNormalRepresentation(dynamic_cast<const CEvaluationNode*>(pNode->getChild()));
 
                 if (pFrac1 != NULL)
                   {
                     pFrac2 = createNormalRepresentation(dynamic_cast<const CEvaluationNode*>(pNode->getChild()->getSibling()));
 
                     if (pFrac2 != NULL)
                       {
                         pResult->setLeft(*pFrac1);
                         pResult->setRight(*pFrac2);
                         delete pFrac1;
                         delete pFrac2;
                       }
                     else
                       {
                         delete pFrac1;
                         pFrac1 = NULL;
                       }
                   }
 
                 break;
               case CEvaluationNodeLogical::NE:
                 pResult = new CNormalLogicalItem();
                 pResult->setType(CNormalLogicalItem::NE);
                 pFrac1 = createNormalRepresentation(dynamic_cast<const CEvaluationNode*>(pNode->getChild()));
 
                 if (pFrac1 != NULL)
                   {
                     pFrac2 = createNormalRepresentation(dynamic_cast<const CEvaluationNode*>(pNode->getChild()->getSibling()));
 
                     if (pFrac2 != NULL)
                       {
                         pResult->setLeft(*pFrac1);
                         pResult->setRight(*pFrac2);
                         delete pFrac1;
                         delete pFrac2;
                       }
                     else
                       {
                         delete pFrac1;
                         pFrac1 = NULL;
                       }
                   }
 
                 break;
               case CEvaluationNodeLogical::GT:
                 pResult = new CNormalLogicalItem();
                 pResult->setType(CNormalLogicalItem::GT);
                 pFrac1 = createNormalRepresentation(dynamic_cast<const CEvaluationNode*>(pNode->getChild()));
 
                 if (pFrac1 != NULL)
                   {
                     pFrac2 = createNormalRepresentation(dynamic_cast<const CEvaluationNode*>(pNode->getChild()->getSibling()));
 
                     if (pFrac2 != NULL)
                       {
                         pResult->setLeft(*pFrac1);
                         pResult->setRight(*pFrac2);
                         delete pFrac1;
                         delete pFrac2;
                       }
                     else
                       {
                         delete pFrac1;
                         pFrac1 = NULL;
                       }
                   }
 
                 break;
               case CEvaluationNodeLogical::LT:
                 pResult = new CNormalLogicalItem();
                 pResult->setType(CNormalLogicalItem::LT);
                 pFrac1 = createNormalRepresentation(dynamic_cast<const CEvaluationNode*>(pNode->getChild()));
 
                 if (pFrac1 != NULL)
                   {
                     pFrac2 = createNormalRepresentation(dynamic_cast<const CEvaluationNode*>(pNode->getChild()->getSibling()));
 
                     if (pFrac2 != NULL)
                       {
                         pResult->setLeft(*pFrac1);
                         pResult->setRight(*pFrac2);
                         delete pFrac1;
                         delete pFrac2;
                       }
                     else
                       {
                         delete pFrac1;
                         pFrac1 = NULL;
                       }
                   }
 
                 break;
               case CEvaluationNodeLogical::GE:
                 pResult = new CNormalLogicalItem();
                 pResult->setType(CNormalLogicalItem::GE);
                 pFrac1 = createNormalRepresentation(dynamic_cast<const CEvaluationNode*>(pNode->getChild()));
 
                 if (pFrac1 != NULL)
                   {
                     pFrac2 = createNormalRepresentation(dynamic_cast<const CEvaluationNode*>(pNode->getChild()->getSibling()));
 
                     if (pFrac2 != NULL)
                       {
                         pResult->setLeft(*pFrac1);
                         pResult->setRight(*pFrac2);
                         delete pFrac1;
                         delete pFrac2;
                       }
                     else
                       {
                         delete pFrac1;
                         pFrac1 = NULL;
                       }
                   }
 
                 break;
               case CEvaluationNodeLogical::LE:
                 pResult = new CNormalLogicalItem();
                 pResult->setType(CNormalLogicalItem::LE);
                 pFrac1 = createNormalRepresentation(dynamic_cast<const CEvaluationNode*>(pNode->getChild()));
 
                 if (pFrac1 != NULL)
                   {
                     pFrac2 = createNormalRepresentation(dynamic_cast<const CEvaluationNode*>(pNode->getChild()->getSibling()));
 
                     if (pFrac2 != NULL)
                       {
                         pResult->setLeft(*pFrac1);
                         pResult->setRight(*pFrac2);
                         delete pFrac1;
                         delete pFrac2;
                       }
                     else
                       {
                         delete pFrac1;
                         pFrac1 = NULL;
                       }
                   }
 
                 break;
               default:
                 break;
             }
         }
     }
 
   return pResult;
 }

CNormalFraction* createNormalRepresentation ( const CEvaluationNode * node )

Create an item power from an evaluation node. Must be a CEvaluationNodeChoice.

Returns: CNormalChoice*, pointer to newly created choice node. CNormalChoice * createChoice(const CEvaluationNode* node); Create an item power from an evaluation node.; CNormalBase*, pointer to newly created normal base object power.

Create an item power from an evaluation node. Must be a CEvaluationNodeChoice.

Returns: CNormalChoice*, pointer to newly created choice node. CNormalChoice * createChoice(const CEvaluationNode* node) { there are two kinds of choice nodes, one where mpTrue and mpFalse are fractions and one where they are choice or logical if the child node directly is a logical or a comparison, we create the Logical, otherwise we create the fraction. CNormalChoice* pChoice=NULL; const CEvaluationNodeChoice* pChoiceNode=dynamic_cast<const CEvaluationNodeChoice*>(node); if(pChoiceNode!=NULL) { pChoice=new CNormalChoice(); const CEvaluationNode* pTmpNode=dynamic_cast<const CEvaluationNode*>(pChoiceNode->getChild()); if(pTmpNode!=NULL) { check if this subtree is a logical expression CNormalLogical* pLogic=createLogical(pTmpNode); if(pLogicc) { pChoice->setCondition(<em>pLogic); delete pLogic; pTmpNode=dynamic_cast<const CEvaluationNode>(pTmpNode->getSibling()); if(pTmpNode!=NULL) { if(isLogical(pTmpNode)) { pLogic=createLogical(pTmpNode); if(pLogic!=NULL) { pChoice->setTrueExpression(pLogic); delete pLogic; pTmpNode=dynamic_cast<const CEvaluationNode*>(pTmpNode->getSibling()); if(pTmpNode!=NULL && isLogical(pTmpNode)) { pLogical=createLogical(pTmpNode); if(pLogical!=NULL) { pChoice->setFalseExpression(<em>pLogic); delete pLogic; } else { delete pChoice; pChoice=NULL; } } else { delete pChoice; pChoice=NULL; } } else { delete pChoice; pChoice=NULL; delete pLogic; } } else { CNormalFraction pFrac=createNormalRepresentation(pTmpNode); if(pFrac) { pChoice->setTrueExpression(<em>pFrac); delete pFrac; pTmpNode=dynamic_cast<const CEvaluationNode>(pTmpNode->getSibling()); if(pTmpNode!=NULL) { CNormalFraction* pFrac=createNormalRepresentation(pTmpNode); if(pFrac) { pChoice->setFalseExpression(*pFrac); delete pFrac; } else { delete pChoice; pChoice=NULL; } } else { delete pChoice; pChoice=NULL; } } else { delete pChoice; pChoice=NULL; } } } else { delete pChoice; pChoice=NULL; } } else { delete pChoice; pChoice=NULL; } } else { delete pChoice; pChoice=NULL; } } return pChoice; } Create an item power from an evaluation node.; CNormalBase*, pointer to newly created normal base object power.

Definition at line 1038 of file ConvertToCEvaluationNode.cpp.

References CEvaluationNode::CALL, CEvaluationNode::CHOICE, CEvaluationNode::CONSTANT, CEvaluationNode::copyBranch(), createFraction(), CEvaluationNode::DELAY, CNormalTranslation::expandProducts(), CEvaluationNode::FUNCTION, CEvaluationNode::getType(), CEvaluationNode::LOGICAL, CNormalTranslation::newEvaluateNumbers(), CEvaluationNode::NUMBER, CEvaluationNode::OPERATOR, CEvaluationNode::type(), and CEvaluationNode::VARIABLE.

Referenced by create_normalform(), createCall(), createChoice(), createFunction(), createGeneralPower(), createLogicalItem(), CNormalTranslation::elementaryEliminationDivide(), CNormalTranslation::elementaryEliminationMinus(), CNormalTranslation::elementaryEliminationModulus(), CNormalTranslation::matchPowerBases(), CNormalTranslation::matchSummands(), and CNormalTranslation::normAndSimplify().

 {
   CNormalFraction* pFrac = NULL;
   CEvaluationNode* pTmp2 = node->copyBranch();
   CEvaluationNode* pTmp = CNormalTranslation::expandProducts(pTmp2);
 
   if (pTmp != NULL)
     {
       delete pTmp2;
     }
   else
     {
       pTmp = pTmp2;
     }
 
   pTmp2 = CNormalTranslation::newEvaluateNumbers(pTmp);
 
   if (pTmp2 != NULL)
     {
       delete pTmp;
     }
   else
     {
       pTmp2 = pTmp;
     }
 
   CEvaluationNode::Type type = CEvaluationNode::type(pTmp2->getType());
 
   switch (type)
     {
       case CEvaluationNode::NUMBER:
       case CEvaluationNode::OPERATOR:
       case CEvaluationNode::CONSTANT:
       case CEvaluationNode::VARIABLE:
       case CEvaluationNode::CHOICE:
       case CEvaluationNode::LOGICAL:
       case CEvaluationNode::FUNCTION:
       case CEvaluationNode::CALL:
       case CEvaluationNode::DELAY:
         pFrac = createFraction(pTmp2);
         break;
       default:
         assert(false);
         break;
     }
 
   if (pTmp2 != NULL)
     {
       delete pTmp2;
     }
 
   return pFrac;
 }

CNormalProduct* createProduct ( const CEvaluationNode * node )

Create a product from an evaluation node that is not necessarily a multiply operator.

Returns: CNormalProduct*, pointer to newly created product.

Definition at line 585 of file ConvertToCEvaluationNode.cpp.

References CCopasiNode< _Data >::addChild(), CEvaluationNode::copyBranch(), CNormalTranslation::createChain(), createItemPower(), CEvaluationNodeOperator::DIVIDE, CEvaluationNode::getType(), CEvaluationNode::getValue(), CEvaluationNodeOperator::MULTIPLY, CNormalProduct::multiply(), CEvaluationNode::NUMBER, CNormalTranslation::ONE_NODE, CEvaluationNode::OPERATOR, CNormalProduct::setFactor(), CNormalTranslation::splitProduct(), CEvaluationNode::subType(), CNormalTranslation::TIMES_NODE, and CEvaluationNode::type().

Referenced by createSum().

 {
   CNormalProduct * pProduct = new CNormalProduct();
 
   if (CEvaluationNode::type(node->getType()) == CEvaluationNode::OPERATOR && (CEvaluationNodeOperator::SubType)CEvaluationNode::subType(node->getType()) == CEvaluationNodeOperator::MULTIPLY)
     {
       // find the product chain, if there are divisions, we have to create a
       // general power, else we create a product with as many items as the
       // multiplication vector has entries
       std::vector<const CEvaluationNode*> multiplications, divisions;
       CNormalTranslation::splitProduct(node, multiplications, divisions, false);
       double factor = 1.0;
 
       if (divisions.empty())
         {
           std::vector<const CEvaluationNode*>::const_iterator it = multiplications.begin(), endit = multiplications.end();
           CNormalItemPower* pItemPower = NULL;
 
           while (it != endit)
             {
               // check if the node is a pure number
               // if so, use it to update the factor
               // instead of creating an item for it
               if (CEvaluationNode::type((*it)->getType()) == CEvaluationNode::NUMBER)
                 {
                   factor *= dynamic_cast<const CEvaluationNodeNumber*>(*it)->getValue();
                 }
               else
                 {
                   pItemPower = createItemPower(*it);
                   assert(pItemPower != NULL);
                   pProduct->multiply(*pItemPower);
                   delete pItemPower;
                 }
 
               ++it;
             }
         }
       else
         {
           CEvaluationNodeOperator* pTmpOperator = new CEvaluationNodeOperator(CEvaluationNodeOperator::DIVIDE, "/");
           std::vector<CEvaluationNode*> tmp;
           std::vector<const CEvaluationNode*>::const_iterator it = multiplications.begin(), endit = multiplications.end();
           // check if the multiplications and divisions contain only numbers
           // in that case, done create a general item
           bool empty = false;
 
           while (it != endit)
             {
               // check if the node is a pure number
               // if so, use it to update the factor
               // instead of adding it to the factor
               if (CEvaluationNode::type((*it)->getType()) == CEvaluationNode::NUMBER)
                 {
                   factor *= dynamic_cast<const CEvaluationNodeNumber*>(*it)->getValue();
                 }
               else
                 {
                   tmp.push_back((*it)->copyBranch());
                 }
 
               ++it;
             }
 
           empty = tmp.empty();
           CEvaluationNode* pTmpNode1 = NULL;
 
           // if we have non number nodes, we need to combine those again into a
           // multiplication chain
           // otherwise the numerator is 1
           if (!empty)
             {
               pTmpNode1 = CNormalTranslation::createChain(&CNormalTranslation::TIMES_NODE, &CNormalTranslation::ONE_NODE, tmp);
             }
           else
             {
               pTmpNode1 = CNormalTranslation::ONE_NODE.copyBranch();
             }
 
           pTmpOperator->addChild(pTmpNode1);
           tmp.clear();
           it = divisions.begin();
           endit = divisions.end();
 
           while (it != endit)
             {
               // check if the node is a pure number
               // if so, use it to update the factor
               // instead of adding it to the vector
               if (CEvaluationNode::type((*it)->getType()) == CEvaluationNode::NUMBER)
                 {
                   factor /= dynamic_cast<const CEvaluationNodeNumber*>(*it)->getValue();
                 }
               else
                 {
                   tmp.push_back((*it)->copyBranch());
                 }
 
               ++it;
             }
 
           // if tmp was empty both times, empty must be true
           empty = (empty & tmp.empty());
 
           if (!empty)
             {
               CNormalItemPower* pItemPower = NULL;
 
               if (!tmp.empty())
                 {
                   pTmpNode1 = CNormalTranslation::createChain(&CNormalTranslation::TIMES_NODE, &CNormalTranslation::ONE_NODE, tmp);
                   pTmpOperator->addChild(pTmpNode1);
                   pItemPower = createItemPower(pTmpOperator);
                 }
               else
                 {
                   pItemPower = createItemPower(pTmpNode1);
                 }
 
               assert(pItemPower != NULL);
               pProduct->multiply(*pItemPower);
               delete pItemPower;
             }
 
           delete pTmpOperator;
         }
 
       pProduct->setFactor(factor);
     }
   else if (CEvaluationNode::type(node->getType()) == CEvaluationNode::NUMBER)
     {
       double factor = dynamic_cast<const CEvaluationNodeNumber*>(node)->getValue();
       // set the factor
       pProduct->setFactor(factor);
     }
   else
     {
       // we create the appropriate item
       CNormalItemPower* pItemPower = createItemPower(node);
       assert(pItemPower != NULL);
       pProduct->multiply(*pItemPower);
       delete pItemPower;
     }
 
   return pProduct;
 }

CNormalSum* createSum ( const CEvaluationNode * node )

Create a sum from an evaluation node -node does not need to be a PLUS operator!

Returns: CNormalSum*, pointer to newly created sum

THIS CODE should be obsolete, RG (06/04/6.2008) Create a product from an evaluation node that is not necessarily a multiply operator.

Returns: CNormalProduct*, pointer to newly created product. Create a sum from an evaluation node -node does not need to be a PLUS operator!; CNormalSum*, pointer to newly created sum

Definition at line 842 of file ConvertToCEvaluationNode.cpp.

References CNormalSum::add(), createProduct(), CNormalProduct::getFactor(), CEvaluationNode::getType(), CEvaluationNodeOperator::MINUS, CEvaluationNode::OPERATOR, CEvaluationNodeOperator::PLUS, CNormalProduct::setFactor(), CNormalTranslation::splitSum(), CEvaluationNode::subType(), and CEvaluationNode::type().

Referenced by createFraction().

 {
   CNormalSum* pSum = new CNormalSum();
 
   if (CEvaluationNode::type(node->getType()) == CEvaluationNode::OPERATOR && ((CEvaluationNodeOperator::SubType)CEvaluationNode::subType(node->getType()) == CEvaluationNodeOperator::MINUS || (CEvaluationNodeOperator::SubType)CEvaluationNode::subType(node->getType()) == CEvaluationNodeOperator::PLUS))
     {
       // find a summation chain and create a product node for each addition and
       // subtraction node
       std::vector<const CEvaluationNode*> additions, subtractions;
       CNormalTranslation::splitSum(node, additions, subtractions, false);
       std::vector<const CEvaluationNode*>::iterator it = additions.begin(), endit = additions.end();
       CNormalProduct* pProduct = NULL;
 
       while (it != endit)
         {
           pProduct = createProduct(*it);
           assert(pProduct != NULL);
           pSum->add(*pProduct);
           delete pProduct;
           ++it;
         }
 
       it = subtractions.begin(), endit = subtractions.end();
 
       while (it != endit)
         {
           pProduct = createProduct(*it);
           assert(pProduct != NULL);
           // since these are subtractions, we need to set the factor the -1.0
           // times the old factor
           pProduct->setFactor(-1.0 * pProduct->getFactor());
           pSum->add(*pProduct);
           delete pProduct;
           ++it;
         }
     }
   else
     {
       // create a sum and add the product made from the given node
       CNormalProduct* pProduct = createProduct(node);
       assert(pProduct != NULL);
       pSum->add(*pProduct);
       delete pProduct;
       /* This code should be obsolete, RG (4.6.2008)
       if (node->getData() == "/")
         {
           CNormalFraction* fraction = createFraction(node);
           sum->add(*fraction);
           delete fraction;
           return sum;
         }
       else
         {
           CNormalProduct* product = createProduct(node);
           if (product != NULL)
             {
               sum->add(*product);
               delete product;
             }
           else
             {
               CNormalFraction* pFraction = createFraction(node);
               sum->add(*pFraction);
               delete pFraction;
             }
           return sum;
         }
         */
     }
 
   return pSum;
 }

bool isLogical ( const CEvaluationNode * pNode )

Create an logical node from an evaluation node. Must be a CEvaluationNodeLogical or the not function.

Returns: CNormalLogical*, pointer to newly created logical node. CNormalLogical* createLogical(const CEvaluationNode* node);

Definition at line 1633 of file ConvertToCEvaluationNode.cpp.

References CEvaluationNode::CONSTANT, CEvaluationNodeConstant::FALSE, CEvaluationNode::FUNCTION, CCopasiNode< _Data >::getChild(), CCopasiNode< _Data >::getSibling(), CEvaluationNode::getType(), isLogical(), CEvaluationNode::LOGICAL, CEvaluationNodeFunction::NOT, CEvaluationNode::subType(), CEvaluationNodeConstant::TRUE, and CEvaluationNode::type().

Referenced by isLogical().

 {
   bool result = false;
 
   // go through the tree until one of the following is encountered:
   // CEvaluationNodeLogical, CEvaluationNodeFunction::NOT or
   // CEvaluationNodeConstant::(TRUE|FALSE)
   if ((CEvaluationNode::type(pNode->getType()) == CEvaluationNode::LOGICAL) ||
       ((CEvaluationNode::type(pNode->getType()) == CEvaluationNode::FUNCTION && ((CEvaluationNodeFunction::SubType)CEvaluationNode::subType(pNode->getType())) == CEvaluationNodeFunction::NOT)) ||
       (CEvaluationNode::type(pNode->getType()) == CEvaluationNode::CONSTANT && ((((CEvaluationNodeConstant::SubType)CEvaluationNode::subType(pNode->getType())) == CEvaluationNodeConstant::TRUE) || (((CEvaluationNodeConstant::SubType)CEvaluationNode::subType(pNode->getType())) == CEvaluationNodeConstant::FALSE)))
      )
 
     {
       result = true;
     }
 
   const CEvaluationNode* pChild = dynamic_cast<const CEvaluationNode*>(pNode->getChild());
 
   while (result == false && pChild != NULL)
     {
       result = isLogical(pChild);
       pChild = dynamic_cast<const CEvaluationNode*>(pChild->getSibling());
     }
 
   return result;
 }

Variable Documentation

std::map<std::string, const CEvaluationNode*> str2eval

Definition at line 64 of file ConvertToCEvaluationNode.cpp.

Functions

Variables

Function Documentation

Variable Documentation