One investigates the solution set of Sv(x)=0, where S represents a matrix and v(x) represents the flux vector, v(x)>=0 .So the set of all possible solutions defines a convex polyhedral cone.
s=[[-1,1,0,0,0,0,0,0,0,0,0,1],[-1,1,1,-1,1,1,0,0,0,0,0,0],[1,-1,-1,0,0,0,0,0,0,0,0,0],[0,0,1,-1,1,0,0,0,1,-1,1,0],[0,0,0,1,-1,-1,0,0,0,0,0,0],[0,0,0,0,0,1,-1,1,0,0,0,0],[0,0,0,0,0,0,-1,1,1,-1,1,1],[0,0,0,0,0,0,1,-1,-1,0,0,0],[0,0,0,0,0,0,0,0,0,1,-1,-1]]
Should I use that,
polytope > $s=new Matrix<Rational>([[0,-1,1,0,0,0,0,0,0,0,0,0,1],[0,-1,1,1,-1,1,1,0,0,0,0,0,0],[0,1,-1,-1,0,0,0,0,0,0,0,0,0],[0,0,0,1,-1,1,0,0,0,1,-1,1,0],[0,0,0,0,1,-1,-1,0,0,0,0,0,0],[0,0,0,0,0,0,1,-1,1,0,0,0,0],[0,0,0,0,0,0,0,-1,1,1,-1,1,1],[0,0,0,0,0,0,0,1,-1,-1,0,0,0],[0,0,0,0,0,0,0,0,0,0,1,-1,-1]]);
polytope > $p=new Cone<Rational> (INEQUALITIES=>$s);
Now I want to get the the minimal set of generating vectors Ei. Ei are the edges of the cone. Each solution flux vector v(x) can be represented as a linear combination of the Ei with nonnegative coefficients Ji. i is index.
I want use print $p->RAYS; but failed. Could you tell me why? Many thanks!