/// <summary>
/// This funnction returns all possible workflows in an incidence matrix format, given the foundation matrix
/// as well as the vIndep vector which represent the independent variables.
/// This function is the top level function which performs the incidence matrix method.
/// </summary>
/// <param name="IMf"></param>
/// <param name="independent"></param>
/// <returns></returns>
private static IncidenceMatrix[] AllWorkflows(IncidenceMatrix IMf, int[] independent)
{
IncidenceMatrix IMi = IncidenceMatrix.Build.DenseOfMatrix(IMf);
// Number of outputs per row (affect model r)
IMi.ReplaceAlltoOne();
// valrf & valcf in [1]
double[] Noutvalr = new double[IMf.RowCount];
for (int model = 0; model < IMf.RowCount; model++)
{
int NOutputs = IMf.NOutputs(model);
int NVariables = IMi.NVariables(model); // Number of variables affecting model r
int NInputs = (NVariables - NOutputs);
Noutvalr[model] = ValRow(NOutputs, NInputs);
}
double[] Noutvalc = new double[IMf.ColumnCount];
for (int c = 0; c < IMf.ColumnCount; c++)
{
int NModels = IMi.Column(c).FindValue(ANY); // Number of models related to variable c
int NInputs = (NModels - 1);
Noutvalc[c] = ValColumn(NInputs);
}
// Assign user defined inputs/outputs
for (int variable = 0; variable < IMi.ColumnCount; variable++)
{
if (independent[variable] == MARKED_IN)
{
IMi.SetInput(variable); // Asign variable as input in IMi
}
else if (independent[variable] == MARKED_OUT)
{
IMi.SetOutput(variable, IMf);
}
}
List <int> notAssigned = FindNotAssignedModels(IMi, IMf);
if (notAssigned.Count == 0)
{
return(new IncidenceMatrix[1] {
IMf
});
}
IncidenceMatrixFirstStep(IMi, Noutvalr, Noutvalc, true);
FillAsInput(IMi, Noutvalr, Noutvalc);
IncidenceMatrixFirstStep(IMi, Noutvalr, Noutvalc, true);
////////////////////////
List <int> reversed = FindReversedModels(IMi, IMf);
bool reversedModelsInSCC = false;
List <int> modelsInScc = ModelsInSCCFromIMM(IMi);
foreach (int model in modelsInScc)
{
if (FindValueInArray(reversed, model)) // if ith model is in SCC and reversed then break
{
reversedModelsInSCC = true;
break;
}
}
if (!reversedModelsInSCC) //i.e. no models in SCC are modified/reversed
{
if (IMi.FindLocations(ANY).Count != 0)
{
// copy default SCC inputs and outputs to new incidence matrix, imMatrixi
IMi.CopyDefaultInputs(IMf, modelsInScc);
}
return(new IncidenceMatrix[1] {
IMi
});
}
else
{
// When clicking 'edit' then get incidence before it has been modified or if creating new object, then get default incidence
DesignStructureMatrix DSM = IMMtoDSM(IMf);
modelsInScc.Clear();
List <Cluster> clusters = Decompose(DSM);
foreach (Cluster cluster in clusters)
{
if (cluster.Count > 1)
{
modelsInScc.AddRange(cluster);
}
}
reversedModelsInSCC = false;
foreach (int model in modelsInScc)
{
// foreach model in SCC is any model...
if (FindValueInArray(reversed, model)) // if ith model is in SCC and reversed then break
{
reversedModelsInSCC = true;
break;
}
}
if (!reversedModelsInSCC) //i.e. copy existing SCC configurations
{
IMi.CopyDefaultInputs(IMf, modelsInScc);
}
}
////////////////////
if (IMi.FindLocations(ANY).Count > 0)
{
return(IncidenceExplore(IMi, IMf, Noutvalr, Noutvalc));
//incmset=uniqueincmset(incmset);Unque object has to be added here....
}
else
{
return(new IncidenceMatrix[1] {
IMi
});
}
throw new ArgumentException(NoIMsFoundErrorMessage);
}