IDoubleArray F = doublesNew.fromFile(arg.getArgument("pathways", 0));
double[] Q = io.readDoubleColumn(arg.getArgument("pathways", 1), 0);
int[] A = str.toIntArray(arg.getArgument("pathways", 2));
int[] B = str.toIntArray(arg.getArgument("pathways", 3));
PathwayDecomposition decomp = new PathwayDecomposition(F, Q, A, B);
List<int[]> allPaths = new ArrayList<int[]>();
List<Double> allFluxes = new ArrayList<Double>();
int[] path = null;
double totalFlux = 0;
do
{
path = decomp.nextPathway();
if (path != null)
{
allPaths.add(path);
allFluxes.add(decomp.getCurrentFlux().doubleValue());
totalFlux += decomp.getCurrentFlux().doubleValue();
}
}
while (path != null);
System.out.println("path\tflux\tcum. flux\tcum. fraction");