<br><font size=2 face="sans-serif">Totally my stupidity - is it better
now. </font>
<br>
<br><font size=2 face="sans-serif">I have added your patches (they broke
CBC but I fixed that)</font>
<br>
<br><font size=2 face="sans-serif">John</font>
<br>
<br>
<br>
<table width=100%>
<tr valign=top>
<td width=40%><font size=1 face="sans-serif"><b>François Galea <francois.galea@prism.uvsq.fr></b>
</font>
<br><font size=1 face="sans-serif">Sent by: coin-lpsolver-bounces@list.coin-or.org</font>
<p><font size=1 face="sans-serif">06/30/2005 07:01 AM</font>
<td width=59%>
<table width=100%>
<tr>
<td>
<div align=right><font size=1 face="sans-serif">To</font></div>
<td valign=top><font size=1 face="sans-serif">coin-lpsolver@list.coin-or.org</font>
<tr>
<td>
<div align=right><font size=1 face="sans-serif">cc</font></div>
<td valign=top>
<tr>
<td>
<div align=right><font size=1 face="sans-serif">Subject</font></div>
<td valign=top><font size=1 face="sans-serif">[Coin-lpsolver] ClpSolve
problem</font></table>
<br>
<table>
<tr valign=top>
<td>
<td></table>
<br></table>
<br>
<br>
<br><font size=2><tt>Hello,<br>
<br>
I have met several segmentation fault errors, when I use the <br>
ClpSimplex::initialSolve( ClpSolve & ) method for starting a new solve.<br>
<br>
The example attached program doesn't work, if I choose to switch <br>
presolve off. If I switch it on, everything works (provided the Clp <br>
sources have been patched with the patch I sent before, otherwise the <br>
names changes won't work). Some other pieces of code I've got crash even
<br>
if I set presolve on.<br>
<br>
I haven't understood enough the internals of ClpSimplex and ClpModel to
<br>
be sure about what is wrong.<br>
<br>
Could somebody tell me what's wrong ?<br>
<br>
Thanks,<br>
<br>
François<br>
<br>
<br>
-- <br>
François Galea<br>
Equipe OPALE - Laboratoire PRiSM<br>
Université de Versailles-Saint Quentin en Yvelines<br>
45 av Etats-Unis F-78035 Versailles CEDEX<br>
Tél. : +33 1 39 25 40 50<br>
<br>
#include <cstdio><br>
<br>
#include "ClpSimplex.hpp"<br>
#include "CoinBuild.hpp"<br>
<br>
int main( int argc, char * argv[] )<br>
{<br>
const static struct { double a, b, c, d; } p [] =<br>
{<br>
{ 5.0, 0.0, 7.0, 3.0 },<br>
{ 8.0, 0.0, 7.0, 6.0 },<br>
{ 0.0, 10.0, 3.0, 9.0 },<br>
{ 1.0, 8.0, 5.0, 9.0 },<br>
{ 2.0, 0.0, 0.0, 7.0 },<br>
{ 0.0, 4.0, 5.0, 0.0 },<br>
{ 0.0, 1.0, 6.0, 2.0 },<br>
{ 2.0, 0.0, 8.0, 3.0 }<br>
};<br>
const static char * columnNames [2] = { "x", "y"
};<br>
double k, obj;<br>
char buf[ 10 ];<br>
char * rowname = buf;<br>
double val[ 2 ];<br>
int ind[ 2 ];<br>
unsigned int i;<br>
const double * sol;<br>
ClpSimplex lp;<br>
CoinBuild * buildObject;<br>
ClpSolve * solveOptions;<br>
<br>
lp.setStrParam( ClpProbName, "foo" );<br>
<br>
buildObject = new CoinBuild( 1 );<br>
buildObject->addCol( 0, NULL, NULL, 0.0, +COIN_DBL_MAX, 4 );<br>
buildObject->addCol( 0, NULL, NULL, 0.0, +COIN_DBL_MAX, 3 );<br>
lp.addColumns( *buildObject );<br>
delete buildObject;<br>
lp.copyColumnNames( columnNames, 0, 2 );<br>
<br>
buildObject = new CoinBuild( 0 );<br>
ind[ 0 ] = 0;<br>
ind[ 1 ] = 1;<br>
for( i = 0; i < (sizeof( p ) / sizeof( p[0] )); i++ )<br>
{<br>
k = p[i].a * p[i].d - p[i].c * p[i].b;<br>
val[0] = (p[i].d-p[i].b)/k;<br>
val[1] = (p[i].a-p[i].c)/k;<br>
<br>
if( ((1.0-val[0]*4.0)/val[1]) > 5.0 )<br>
buildObject->addRow( 2, ind, val, 0.0, 1.0 );<br>
else<br>
buildObject->addRow( 2, ind, val, 1.0, +COIN_DBL_MAX
);<br>
}<br>
lp.setOptimizationDirection( -1.0 );<br>
lp.addRows( *buildObject );<br>
for( i = 0; i < (sizeof( p ) / sizeof( p[0] )); i++ )<br>
{<br>
sprintf( buf, "c%d", i );<br>
lp.copyRowNames( &rowname, i, i+1 );<br>
}<br>
delete buildObject;<br>
<br>
solveOptions = new ClpSolve();<br>
solveOptions->setSolveType( ClpSolve::useDual );<br>
solveOptions->setPresolveType( ClpSolve::presolveOff, 0 );<br>
// solveOptions->setPresolveType( ClpSolve::presolveOn, 5 );<br>
lp.setLogLevel( 1 );<br>
lp.initialSolve( *solveOptions );<br>
delete solveOptions;<br>
<br>
lp.writeMps( "foo.mps" );<br>
<br>
obj = -lp.rawObjectiveValue();<br>
sol = lp.getColSolution();<br>
<br>
printf( "obj = %f\nx = %f\ny = %f\n", obj, sol[ 0 ], sol[
1 ] );<br>
<br>
return 0;<br>
}<br>
_______________________________________________<br>
Coin-lpsolver mailing list<br>
Coin-lpsolver@list.coin-or.org<br>
http://list.coin-or.org/mailman/listinfo/coin-lpsolver<br>
</tt></font>
<br>