// Just added 5/18/04
void ReactionsWidget1::slotBtnDelete()
{
assert(CCopasiRootContainer::getDatamodelList()->size() > 0);
CCopasiDataModel* pDataModel = (*CCopasiRootContainer::getDatamodelList())[0];
assert(pDataModel != NULL);
CModel * pModel = pDataModel->getModel();
if (pModel == NULL)
return;
CReaction * pReaction =
dynamic_cast< CReaction * >(CCopasiRootContainer::getKeyFactory()->get(mKey));
if (pReaction == NULL) return;
QMessageBox::StandardButton choice =
CQMessageBox::confirmDelete(NULL, "reaction",
FROM_UTF8(pReaction->getObjectName()),
pReaction->getDeletedObjects());
switch (choice)
{
case QMessageBox::Ok: // Yes or Enter
{
pDataModel->getModel()->removeReaction(mKey);
mpRi->setFunctionWithEmptyMapping("");
protectedNotify(ListViews::REACTION, ListViews::DELETE, mKey);
protectedNotify(ListViews::REACTION, ListViews::DELETE, "");//Refresh all as there may be dependencies.
break;
}
default: // No or Escape
break;
}
}
void CFixLocalReactionParameters::changeModel()
{
CCopasiParameter * pParameter = NULL;
CModelValue * pModelValue = NULL;
CReaction * pReaction = NULL;
std::stringstream NameStream;
std::stringstream Message;
std::string OldCN;
std::string NewCNBase;
std::string NewCN;
std::string Infix;
std::string::size_type Start;
// Loop through all changes.
std::multimap< CCopasiParameter *, const CExpression * >::const_iterator itChanges = mChanges.begin();
std::multimap< CCopasiParameter *, const CExpression * >::const_iterator endChanges = mChanges.end();
for (; itChanges != endChanges; ++itChanges)
{
if (pParameter != itChanges->first)
{
// We have a new parameter
pParameter = itChanges->first;
OldCN = "<" + pParameter->getCN() + ",Reference=";
// Create a global quantity of type FIXED.
std::string Name = pParameter->getObjectName();
pReaction = static_cast< CReaction * >(pParameter->getObjectAncestor("Reaction"));
Name += "{" + pReaction->getObjectName() + "}";
pModelValue = mpModel->createModelValue(Name, pParameter->getValue< C_FLOAT64 >());
// In case the created name is not unique we append _n with increasing n
// until we succeed;
C_INT32 index = 0;
while (pModelValue == NULL)
{
NameStream.str("");
NameStream << Name << "_" << index++;
pModelValue = mpModel->createModelValue(NameStream.str(), pParameter->getValue< C_FLOAT64 >());
}
NewCNBase = "<" + pModelValue->getCN() + ",Reference=";
// If the parameter is actually used in the reaction
// it is changed to the global quantity.
if (pReaction->isLocalParameter(pParameter->getObjectName()))
pReaction->setParameterMapping(pParameter->getObjectName(), pModelValue->getKey());
Message << " " << pParameter->getObjectName() << " in " << pReaction->getObjectName()
<< " is replaced by " << pModelValue->getObjectName() << std::endl;
}
// We need to distinguish between initial and other expressions.
if (itChanges->second->getObjectName().compare(0, 7, "Initial") == 0)
NewCN = NewCNBase + "Initial";
else
NewCN = NewCNBase;
// Replace the OldCN of the parameter with the NewCN of global quantity in all expressions.
Infix = itChanges->second->getInfix();
// There may be more than one occurrence.
Start = 0;
while ((Start = Infix.find(OldCN), Start) != std::string::npos)
Infix.replace(Start, OldCN.length(), NewCN);
const_cast< CExpression * >(itChanges->second)->setInfix(Infix);
}
CCopasiMessage(CCopasiMessage::WARNING, MCXML + 14, Message.str().c_str());
}
bool CModelAdd::addReactions(std::string name)
{
bool info = false;
//create copies of the relevant reactions
size_t i, imax = mmModel->getReactions().size();
size_t ic, icmax = mmModel->getCompartments().size();
for (ic = 0; ic < icmax; ++ic)
{
const CCompartment* sourceComp = &mmModel->getCompartments()[ic];
if (!sourceComp) return info;
for (i = 0; i < imax; ++i)
{
CReaction * sourceReac = &mmModel->getReactions()[i];
if (reactionInvolvesCompartment(sourceReac, sourceComp))
{
std::string newName = sourceReac->getObjectName() + "_" + name;
CReaction* newReac = mpModel->createReaction(newName);
if (!newReac) return info;
//copy the chemical equation. If the involved metabs are among those that
//were copied with the compartment, replace them. Otherwise keep the original metab
newReac->setReversible(sourceReac->isReversible());
std::map<std::string, std::string>::const_iterator mapIt;
std::string targetKey;
size_t j, jmax = sourceReac->getChemEq().getSubstrates().size();
for (j = 0; j < jmax; ++j)
{
const CChemEqElement * sourceElement = &sourceReac->getChemEq().getSubstrates()[j];
//check if the metab is in the map. If yes, translate it, otherwise not.
mapIt = keyMap.find(sourceElement->getMetaboliteKey());
if (mapIt == keyMap.end())
{
targetKey = sourceElement->getMetaboliteKey();
}
else
targetKey = mapIt->second;
newReac->addSubstrate(targetKey, sourceElement->getMultiplicity());
}
jmax = sourceReac->getChemEq().getProducts().size();
for (j = 0; j < jmax; ++j)
{
const CChemEqElement * sourceElement = &sourceReac->getChemEq().getProducts()[j];
//check if the metab is in the map. If yes, translate it, otherwise not.
mapIt = keyMap.find(sourceElement->getMetaboliteKey());
if (mapIt == keyMap.end())
{
targetKey = sourceElement->getMetaboliteKey();
}
else
targetKey = mapIt->second;
newReac->addProduct(targetKey, sourceElement->getMultiplicity());
}
jmax = sourceReac->getChemEq().getModifiers().size();
for (j = 0; j < jmax; ++j)
{
const CChemEqElement * sourceElement = &sourceReac->getChemEq().getModifiers()[j];
//check if the metab is in the map. If yes, translate it, otherwise not.
mapIt = keyMap.find(sourceElement->getMetaboliteKey());
if (mapIt == keyMap.end())
{
targetKey = sourceElement->getMetaboliteKey();
}
else
targetKey = mapIt->second;
newReac->addModifier(targetKey);
}
//set the kinetic function
newReac->setFunction(const_cast<CFunction*>(sourceReac->getFunction()));
//mapping and local parameters
for (j = 0; j < newReac->getFunctionParameters().size(); ++j)
{
switch (newReac->getFunctionParameters()[j]->getUsage())
{
case CFunctionParameter::SUBSTRATE:
case CFunctionParameter::PRODUCT:
case CFunctionParameter::MODIFIER:
//translate the metab keys
{
bool isVector = (newReac->getFunctionParameters()[j]->getType() == CFunctionParameter::VFLOAT64);
//we assume that only SUBSTRATE, PRODUCT, MODIFIER can be vectors
if (isVector)
newReac->clearParameterMapping(j);
size_t k;
for (k = 0; k < sourceReac->getParameterMappings()[j].size(); ++k)
{
mapIt = keyMap.find(sourceReac->getParameterMappings()[j][k]);
if (mapIt == keyMap.end())
{
targetKey = sourceReac->getParameterMappings()[j][k];
}
else
targetKey = mapIt->second;
if (isVector)
newReac->addParameterMapping(j, targetKey);
else
newReac->setParameterMapping(j, targetKey);
}
}
break;
case CFunctionParameter::TIME:
//just copy the key
{
mapIt = keyMap.find(sourceReac->getParameterMappings()[j][0]);
if (mapIt == keyMap.end())
{
targetKey = sourceReac->getParameterMappings()[j][0];
}
else
targetKey = mapIt->second;
newReac->setParameterMapping(j, targetKey);
}
break;
case CFunctionParameter::VOLUME:
//translate the compartment key if necessary
if (sourceReac->getParameterMappings()[j][0] == sourceComp->getKey())
newReac->setParameterMapping(j, keyMap[sourceComp->getKey()]);
else
{
mapIt = keyMap.find(sourceReac->getParameterMappings()[j][0]);
if (mapIt == keyMap.end())
{
targetKey = sourceReac->getParameterMappings()[j][0];
}
else
targetKey = mapIt->second;
newReac->setParameterMapping(j, targetKey);
}
//TODO: this needs to be adapted when sets of compartments will be copied
break;
case CFunctionParameter::PARAMETER:
if (sourceReac->isLocalParameter(j))
newReac->setParameterValue(newReac->getFunctionParameters()[j]->getObjectName(),
sourceReac->getParameterValue(newReac->getFunctionParameters()[j]->getObjectName()));
else
{
mapIt = keyMap.find(sourceReac->getParameterMappings()[j][0]);
if (mapIt == keyMap.end())
{
targetKey = sourceReac->getParameterMappings()[j][0];
}
else
targetKey = mapIt->second;
newReac->setParameterMapping(j, targetKey);
}
break;
default:
return info;
break;
}
}
}
}
}
return true;
}
