[Clp] Wrong problemStatus with setupForRepeatedUse(2)

[Sebastian Theophil]

Hi John, 


the code pasted below produces the wrong problemStatus in resolve().
After the modifications, the problem should be dual infeasible, but is
reported as being primal infeasible. The created LP file will clearly
show that the problem is not primal infeasible.

			OsiClpSolverInterface simplex;
			simplex.readMps("C:\\initial.mps");
			simplex.setDblParam( OsiPrimalObjectiveLimit,
-1.0 ); 
			simplex.setupForRepeatedUse(2);
			simplex.initialSolve();

			double fObjective = simplex.getObjValue(); //
objective is ~ 0.54

			// replace objective variable with its value in
binding constraints
			simplex.setRowLower( 12,
simplex.getRowLower()[12] - fObjective );
			simplex.modifyCoefficient( 12, 7, 0.0 );

			simplex.setRowLower( 14,
simplex.getRowLower()[14] - fObjective );
			simplex.modifyCoefficient( 14, 7, 0.0 );

			simplex.resolve();
			fObjective = simplex.getObjValue(); // is ~ 0.72

			simplex.setRowLower( 13,
simplex.getRowLower()[13] - fObjective );
			simplex.modifyCoefficient( 13, 7, 0.0 );

			simplex.setRowLower( 10,
simplex.getRowLower()[10] - fObjective );
			simplex.modifyCoefficient( 10, 7, 0.0 );

			simplex.resolve(); 
			simplex.writeLp("C:\\pretty");
			_ASSERT(simplex.isProvenPrimalInfeasible()); //
should be dual infeasible


I would very much appreciate if you could take a look at the problem. 

Regards, Sebastian

--
Sebastian Theophil . stheophil at think-cell.com
Software Engineer
 
think-cell Software GmbH . Invalidenstr. 34 . 10115 Berlin, Germany 
http://www.think-cell.com . phone +49-30-666473-10 . toll-free (US) +1-800-891-8091
Directors: Dr. Markus Hannebauer, Dr. Arno Schoedl . Amtsgericht Berlin-Charlottenburg, HRB 85229

-------------- next part --------------
A non-text attachment was scrubbed...
Name: initial.mps
Type: application/octet-stream
Size: 1878 bytes
Desc: initial.mps
URL: <http://list.coin-or.org/pipermail/clp/attachments/20091106/21278339/attachment.obj>