<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<br>
The problem may be related to the fact that there are more equality constraints<br>
(i.e. 96) than unknowns (i.e. 93). But it seems that only 81 of the 96 rows of<br>
the equality matrix are numerically linearly independent :<br>
<br>
-->rank(eqMat)<br>
 ans  =<br>
    81.<br>
<br>
There are no such rank defficient tests neither for karmarkar, nor for linpro.<br>
Actually, there were no unit test at all for linpro until June 2011... But I<br>
expect that linpro is more robust in general, so I am not surprised by the<br>
result. But karmarkar does not seem to work here, and returns with the exitflag<br>
equals to -1, meaning that the initial guess x0 could not be found. More testing<br>
should be done to create a simpler rank defficient linear problem with, say 4-5<br>
unknowns and see what happens. Do you have such a simple example and its exact<br>
solution ?<br></blockquote><div>This kind of rank deficiency is very frequent in this kind of problem (biochemical networks) and has never been a problem using the simplex method. Since the only thing I know about the karmarkar algorithm is that it is an interior point method, I have no clue about what the error message about rank deficiency means. I have tried eliminating the redundancies but the reesult is the same so the error may refer to a completely different issue. <br>
</div><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<br>
When you say that you have problems when you use linpro through the Java<br>
interface, what do you mean exactly ? Can you provide a script ?<br></blockquote><div> </div><div>If I  put this problem in an .sce file and execute it, I get the result you mention. No problem so far. If I call linpro from the java interface with the Siclab.exec() method it also works fine as long as I use a 32 bits machine. But calling linpro through  the exec() method in any linux 64 bits machine crashes the java virtual machine. This same problem happened whenever I called the javasci engine with scilab 4.2 but in 5.3 it is restricted to the linpro function, any other operation I've tried so far (including karmarkar) seem to work fine. <br>
<br>Alberto<br></div><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<br>
Regards,<br>
<font color="#888888"><br>
Michaël Baudin<br>
<br>
<br>
</font></blockquote></div><br>