clp.exe may crash.

[yk011309]

clp.exe crashes with the following mps file when primalTolerance is small.

> clp.exe test.mps -primalTolerance 1e-6 -minimize -solve -solution result.txt -> crashed. > clp.exe test.mps -primalTolerance 1e-5 -minimize -solve -solution result.txt -> didn't crash.

I tested with Clp-releases.1.17.8-w64-msvc16-md and Clp-releases.1.17.8-w64-msvc17-md.

test.mps

NAME TestConvexQP ROWS N OBJECT G R1 G R2 G R3 G R4 G R5 G R6 G R7 G R8 G R9 G R10 G R11 G R12 G R13 G R14 COLUMNS X1 OBJECT -6.33174068731535E-16 X1 R1 -1.38237231824552 X1 R2 2.57705493348984 X1 R3 0.893710887963403 X1 R4 -1.40467816181794 X1 R5 -0.74909952135441 X1 R6 -0.0750949408363687 X1 R7 -0.731082739214938 X1 R8 -1.38209298467333 X1 R9 2.57905334394858 X1 R10 1.05764436210498 X1 R11 -1.40458079066066 X1 R12 -0.746246377481727 X1 R13 -0.0780573394061342 X1 R14 -0.726226183869852 X2 OBJECT -5.69390793737391E-16 X2 R1 -0.0163560436138557 X2 R2 1.99468668353512 X2 R3 -2.08061722789719 X2 R4 -1.7452198444874 X2 R5 0.00576562468973074 X2 R6 1.72289997073075 X2 R7 0.00317872161630151 X2 R8 -0.0162427791511776 X2 R9 1.99549700105277 X2 R10 -2.01414531505974 X2 R11 -1.74518036233091 X2 R12 0.00692252038666297 X2 R13 1.72169877432745 X2 R14 0.00514796264702411 X3 OBJECT -6.14525791364784E-16 X3 R1 -0.128877715086185 X3 R2 1.37055762721054 X3 R3 -0.0716728882465328 X3 R4 -1.97720981439181 X3 R5 -0.319475756546151 X3 R6 0.00441989753315607 X3 R7 -0.315769691514621 X3 R8 -0.128659639793125 X3 R9 1.3721177832026 X3 R10 0.0563097245659944 X3 R11 -1.97713379687808 X3 R12 -0.317248311483645 X3 R13 0.00210715749555091 X3 R14 -0.311978186176328 X4 OBJECT -1.15098902631061E-15 X4 R1 -0.0432908552288926 X4 R2 1.07806573967557 X4 R3 0.327339448908376 X4 R4 -0.678501456388584 X4 R5 -1.35945672835786 X4 R6 0.106541328004554 X4 R7 -0.483422818595674 X4 R8 -0.0435646187786731 X4 R9 1.07610717831539 X4 R10 0.16667487337562 X4 R11 -0.678596885926302 X4 R12 -1.3622529794081 X4 R13 0.109444655169813 X4 R14 -0.488182532297562 X5 OBJECT -2.09684700158697E-16 X5 R1 -0.0993212094007268 X5 R2 0.0618017282483161 X5 R3 1.30352157601599 X5 R4 0.0265692044750114 X5 R5 -0.295791120941798 X5 R6 -1.93279732830741 X5 R7 -0.289370365468515 X5 R8 -0.0992885550823477 X5 R9 0.0620353440319846 X5 R10 1.32268553047975 X5 R11 0.0265805872413084 X5 R12 -0.295457586137262 X5 R13 -1.93314363507356 X5 R14 -0.288802630257806 X6 OBJECT 2.92300905702092E-16 X6 R1 -0.160617952702279 X6 R2 1.20755881297912 X6 R3 0.0441285699135943 X6 R4 -0.711979517569525 X6 R5 -0.552017697819496 X6 R6 0.0946476743325082 X6 R7 -1.40296588288707 X6 R8 -0.160254311772758 X6 R9 1.21016037580263 X6 R10 0.257539798818364 X6 R11 -0.711852758233496 X6 R12 -0.548303429313324 X6 R13 0.0907911763073798 X6 R14 -1.3966435411575 RHS RHS1 R1 -37.6429498997557 RHS1 R2 97.2260867768474 RHS1 R3 92.1174211328235 RHS1 R4 -76.1564401334131 RHS1 R5 -78.9968183755308 RHS1 R6 -78.6017251879198 RHS1 R7 -78.9197383853335 RHS1 R8 -37.6327822133182 RHS1 R9 97.2988733287527 RHS1 R10 98.0872492496583 RHS1 R11 -76.1528955227783 RHS1 R12 -78.8929569832673 RHS1 R13 -78.7095974542488 RHS1 R14 -78.7429479316258 BOUNDS FR BOUND X1 FR BOUND X2 FR BOUND X3 FR BOUND X4 FR BOUND X5 FR BOUND X6 QUADOBJ X1 X1 0.14280365295569 X2 X1 0.112326679560112 X3 X1 -0.0123426882592386 X4 X1 0.0130221420022721 X5 X1 -0.00303279324360872 X6 X1 0.0214287229889673 X2 X2 0.287384690609157 X3 X2 -0.0105760572630717 X4 X2 0.0659032217822261 X5 X2 0.0136386481909511 X6 X2 0.0684447360576953 X3 X3 0.0186310316323555 X4 X3 -0.00534137201121943 X5 X3 0.00551850579515465 X6 X3 -0.00717579847817494 X4 X4 0.074974826336554 X5 X4 0.000161393296927416 X6 X4 -0.0160052249623282 X5 X5 0.0191971633222669 X6 X5 -0.0012507650563196 X6 X6 0.0684963196238557 ENDATA

[jjhforrest]

Quadratic has not been kept up to date and "initialSolve" got confused. You will not get any advantages from -solve and I would suggest that you use -primals (or -barrier).

[yk011309]

Using -primalS, -dualS, -either to test.mps crashed clp.exe as did using -solve. When using -barrier clp.exe did not crash and got the optimal solution.

However, when -barrier is used for test2.mps below, "Stopped on iterations or time" is output and the optimal solution cannot be obtained. Using anything other than -barrier gives the optimal solution.

test2.mps

NAME TestConvexQP2 ROWS N OBJECT G R1 G R2 G R3 G R4 G R5 G R6 G R7 COLUMNS X1 OBJECT 2.55795384873636E-13 X1 R1 -1.25631423183625 X1 R2 1.98604848145615 X1 R3 0.99185595596329 X1 R4 -0.954603162502961 X1 R5 -0.52646433314256 X1 R6 -0.0817190489498771 X1 R7 -0.509774744102942 X2 OBJECT 1.61248792096558E-12 X2 R1 -0.0233025790374162 X2 R2 2.02878160251631 X2 R3 -2.01045269630086 X2 R4 -1.77043702879025 X2 R5 -0.00541356354344802 X2 R6 1.7219043013822 X2 R7 -0.00700117872522277 X3 OBJECT 1.45661260830821E-13 X3 R1 -0.112130956868261 X3 R2 1.29419055683126 X3 R3 0.0476644193715896 X3 R4 -1.91800188968916 X3 R5 -0.288366583516715 X3 R6 0.00162596942009061 X3 R7 -0.28353413552517 X4 OBJECT 3.83693077310454E-13 X4 R1 -0.0113112590499856 X4 R2 0.92404334623397 X4 R3 0.149804797298355 X4 R4 -0.563209429909544 X4 R5 -1.30589441980973 X4 R6 0.108505685491132 X4 R7 -0.432678037356059 X5 OBJECT -1.45661260830821E-13 X5 R1 -0.0811470859700689 X5 R2 -0.0234956312216178 X5 R3 1.31319665965662 X5 R4 0.0914822900989671 X5 R5 -0.263757723649542 X5 R6 -1.93367177507406 X5 R7 -0.257583151448648 X6 OBJECT 2.50466314355435E-13 X6 R1 -0.138826951875191 X6 R2 1.10913752160014 X6 R3 0.24633224727497 X6 R4 -0.63519567293963 X6 R5 -0.510861899191156 X6 R6 0.0901673764487839 X6 R7 -1.35976940466292 RHS RHS1 R1 -17.2597424123373 RHS1 R2 14225.9130627521 RHS1 R3 22936.8714349602 RHS1 R4 -8833.86192490507 RHS1 R5 -17192.4491752307 RHS1 R6 -36571.2210227025 RHS1 R7 -13972.2371885318 BOUNDS FR BOUND X1 FR BOUND X2 FR BOUND X3 FR BOUND X4 FR BOUND X5 FR BOUND X6 QUADOBJ X1 X1 0.108246709675597 X2 X1 0.0981788720437024 X3 X1 -0.00845140083983216 X4 X1 0.00561835927314888 X5 X1 -0.00820344962936441 X6 X1 0.0186714268943699 X2 X2 0.386635725120289 X3 X2 -0.0062400187535258 X4 X2 0.0732691582982062 X5 X2 -0.00515108026110568 X6 X2 0.0826523917402175 X3 X3 0.0228650739989253 X4 X3 -0.0031017729527417 X5 X3 0.00656871316535296 X6 X3 -0.0035688309802128 X4 X4 0.0709331318070513 X5 X4 -0.00254740864681019 X6 X4 -0.0124685458835092 X5 X5 0.0192526536066132 X6 X5 -0.00380880744331457 X6 X6 0.0718042446464811 ENDATA

[jjhforrest]

I managed to get first problem to error as described with some builds. Can't get it to error with master - that says - Unscaled problem has dual infeasibilities. With -scaling off all versions work.

As it looks like a scaling/accuracy problem - how are you solving the problems? If you are creating an mps file then it would be better to create it with more accuracy writing values as exact representations of doubles (like formatType==2 in writeMps

[yk011309]

I am generating test cases in C#. I changed ToString() to ToString("R") as you advised.

First, I generated test_R.mps with improved accuracy of test.mps. > clp.exe test.mps -scaling X -primalS -solution test.txt The results are shown in the table below. scaling auto dynamic equi geo off rows test.mps Crash Crash OK Crash OK OK test_R.mps OK OK Crash OK Crash OK In all test cases I generated clp.exe did not crash if -scaling was rows. But it could be a coincidence so I would like to remove the root cause of the crash if possible.

I then generated test2_R.mps, which is an improved version of test2.mps, and a new test case, test3.mps. > clp.exe test.mps -scaling X -barrier -solution test.txt The results are shown in the table below. scaling auto dynamic equi geo off rows test2.mps Stopped Stopped Stopped Stopped Stopped Stopped test2_R.mps OK OK OK OK OK OK test3_R.mps OK Crash OK OK Crash OK Test2_R.mps has improved results. However, I found a test case where clp.exe crashes when -barrier is used.

test_R.mps

NAME TestConvexQPRoundTrip ROWS N OBJECT G R1 G R2 G R3 G R4 G R5 G R6 G R7 G R8 G R9 G R10 G R11 G R12 G R13 G R14 COLUMNS X1 OBJECT -6.3317406873153459E-16 X1 R1 -1.3823723182455239 X1 R2 2.5770549334898445 X1 R3 0.89371088796340281 X1 R4 -1.4046781618179374 X1 R5 -0.74909952135441049 X1 R6 -0.075094940836368709 X1 R7 -0.73108273921493794 X1 R8 -1.3820929846733299 X1 R9 2.5790533439485821 X1 R10 1.0576443621049834 X1 R11 -1.4045807906606562 X1 R12 -0.74624637748172684 X1 R13 -0.078057339406134157 X1 R14 -0.72622618386985249 X2 OBJECT -5.6939079373739127E-16 X2 R1 -0.016356043613855675 X2 R2 1.9946866835351202 X2 R3 -2.0806172278971924 X2 R4 -1.7452198444873976 X2 R5 0.0057656246897307417 X2 R6 1.7228999707307491 X2 R7 0.0031787216163015053 X2 R8 -0.016242779151177555 X2 R9 1.9954970010527744 X2 R10 -2.0141453150597384 X2 R11 -1.7451803623309112 X2 R12 0.006922520386662967 X2 R13 1.7216987743274506 X2 R14 0.0051479626470241069 X3 OBJECT -6.1452579136478391E-16 X3 R1 -0.12887771508618545 X3 R2 1.3705576272105386 X3 R3 -0.071672888246532784 X3 R4 -1.9772098143918055 X3 R5 -0.31947575654615096 X3 R6 0.0044198975331560661 X3 R7 -0.31576969151462142 X3 R8 -0.12865963979312511 X3 R9 1.3721177832025986 X3 R10 0.056309724565994416 X3 R11 -1.9771337968780769 X3 R12 -0.31724831148364541 X3 R13 0.0021071574955509147 X3 R14 -0.31197818617632811 X4 OBJECT -1.1509890263106115E-15 X4 R1 -0.043290855228892589 X4 R2 1.078065739675568 X4 R3 0.32733944890837563 X4 R4 -0.67850145638858383 X4 R5 -1.359456728357856 X4 R6 0.10654132800455392 X4 R7 -0.48342281859567426 X4 R8 -0.0435646187786731 X4 R9 1.0761071783153946 X4 R10 0.16667487337561995 X4 R11 -0.6785968859263023 X4 R12 -1.3622529794081006 X4 R13 0.10944465516981253 X4 R14 -0.48818253229756203 X5 OBJECT -2.0968470015869656E-16 X5 R1 -0.099321209400726834 X5 R2 0.061801728248316132 X5 R3 1.3035215760159888 X5 R4 0.026569204475011418 X5 R5 -0.29579112094179788 X5 R6 -1.9327973283074056 X5 R7 -0.28937036546851452 X5 R8 -0.099288555082347726 X5 R9 0.062035344031984616 X5 R10 1.3226855304797531 X5 R11 0.026580587241308393 X5 R12 -0.2954575861372622 X5 R13 -1.9331436350735625 X5 R14 -0.288802630257806 X6 OBJECT 2.92300905702092E-16 X6 R1 -0.1606179527022788 X6 R2 1.2075588129791246 X6 R3 0.044128569913594304 X6 R4 -0.71197951756952516 X6 R5 -0.55201769781949572 X6 R6 0.094647674332508216 X6 R7 -1.4029658828870666 X6 R8 -0.16025431177275784 X6 R9 1.2101603758026327 X6 R10 0.25753979881836375 X6 R11 -0.71185275823349614 X6 R12 -0.5483034293133241 X6 R13 0.09079117630737979 X6 R14 -1.3966435411574976 RHS RHS1 R1 -37.642949899755664 RHS1 R2 97.226086776847438 RHS1 R3 92.117421132823452 RHS1 R4 -76.156440133413071 RHS1 R5 -78.996818375530836 RHS1 R6 -78.601725187919826 RHS1 R7 -78.919738385333531 RHS1 R8 -37.632782213318244 RHS1 R9 97.298873328752663 RHS1 R10 98.087249249658257 RHS1 R11 -76.152895522778337 RHS1 R12 -78.892956983267325 RHS1 R13 -78.7095974542488 RHS1 R14 -78.742947931625764 BOUNDS FR BOUND X1 FR BOUND X2 FR BOUND X3 FR BOUND X4 FR BOUND X5 FR BOUND X6 QUADOBJ X1 X1 0.14280365295569 X2 X1 0.11232667956011218 X3 X1 -0.012342688259238639 X4 X1 0.013022142002272096 X5 X1 -0.0030327932436087189 X6 X1 0.021428722988967284 X2 X2 0.28738469060915678 X3 X2 -0.010576057263071748 X4 X2 0.065903221782226115 X5 X2 0.013638648190951131 X6 X2 0.068444736057695257 X3 X3 0.018631031632355539 X4 X3 -0.0053413720112194286 X5 X3 0.00551850579515465 X6 X3 -0.0071757984781749447 X4 X4 0.074974826336553987 X5 X4 0.00016139329692741594 X6 X4 -0.016005224962328236 X5 X5 0.019197163322266918 X6 X5 -0.0012507650563196 X6 X6 0.068496319623855662 ENDATA

test3_R.mps

NAME TestConvexQPRoundTripInfeasible ROWS N OBJECT G R1 G R2 G R3 G R4 G R5 G R6 G R7 G R8 G R9 G R10 G R11 G R12 G R13 G R14 G R15 G R16 G R17 G R18 G R19 G R20 G R21 G R22 G R23 G R24 G R25 G R26 G R27 G R28 G R29 G R30 G R31 G R32 G R33 G R34 G R35 G R36 G R37 G R38 G R39 G R40 G R41 G R42 G R43 G R44 G R45 G R46 G R47 G R48 G R49 G R50 G R51 G R52 G R53 G R54 G R55 G R56 G R57 G R58 G R59 G R60 G R61 G R62 G R63 G R64 G R65 G R66 G R67 G R68 G R69 G R70 G R71 G R72 G R73 G R74 G R75 G R76 G R77 G R78 G R79 G R80 G R81 G R82 G R83 G R84 COLUMNS X1 OBJECT 0.71410195383874608 X1 R1 1.7461065890401677 X1 R2 -1.7461065890401677 X1 R3 1 X1 R4 -1 X1 R5 1 X1 R6 -1 X1 R7 1 X1 R8 -1 X1 R9 1 X1 R10 -1 X1 R11 1 X1 R12 -1 X1 R13 -0.74610658904016769 X1 R14 0.74610658904016769 X1 R15 -0.74610658904016769 X1 R16 0.74610658904016769 X1 R17 -0.74610658904016769 X1 R18 0.74610658904016769 X1 R19 -0.74610658904016769 X1 R20 0.74610658904016769 X1 R21 -0.74610658904016769 X1 R22 0.74610658904016769 X1 R43 1.7461065890401677 X1 R44 -1.7461065890401677 X1 R45 1 X1 R46 -1 X1 R47 1 X1 R48 -1 X1 R49 1 X1 R50 -1 X1 R51 1 X1 R52 -1 X1 R53 1 X1 R54 -1 X1 R55 -0.74610658904016769 X1 R56 0.74610658904016769 X1 R57 -0.74610658904016769 X1 R58 0.74610658904016769 X1 R59 -0.74610658904016769 X1 R60 0.74610658904016769 X1 R61 -0.74610658904016769 X1 R62 0.74610658904016769 X1 R63 -0.74610658904016769 X1 R64 0.74610658904016769 X2 OBJECT -0.32533914707899636 X2 R1 0.98332866760501225 X2 R2 -0.98332866760501225 X2 R3 -1 X2 R4 1 X2 R13 -1.9833286676050124 X2 R14 1.9833286676050124 X2 R15 -0.98332866760501225 X2 R16 0.98332866760501225 X2 R17 -0.98332866760501225 X2 R18 0.98332866760501225 X2 R19 -0.98332866760501225 X2 R20 0.98332866760501225 X2 R21 -0.98332866760501225 X2 R22 0.98332866760501225 X2 R23 1 X2 R24 -1 X2 R25 1 X2 R26 -1 X2 R27 1 X2 R28 -1 X2 R29 1 X2 R30 -1 X2 R43 0.98332866760501225 X2 R44 -0.98332866760501225 X2 R45 -1 X2 R46 1 X2 R55 -1.9833286676050124 X2 R56 1.9833286676050124 X2 R57 -0.98332866760501225 X2 R58 0.98332866760501225 X2 R59 -0.98332866760501225 X2 R60 0.98332866760501225 X2 R61 -0.98332866760501225 X2 R62 0.98332866760501225 X2 R63 -0.98332866760501225 X2 R64 0.98332866760501225 X2 R65 1 X2 R66 -1 X2 R67 1 X2 R68 -1 X2 R69 1 X2 R70 -1 X2 R71 1 X2 R72 -1 X3 OBJECT -0.10297267340369098 X3 R1 0.080493799416719816 X3 R2 -0.080493799416719816 X3 R5 -1 X3 R6 1 X3 R13 -0.080493799416719816 X3 R14 0.080493799416719816 X3 R15 -1.0804937994167199 X3 R16 1.0804937994167199 X3 R17 -0.080493799416719816 X3 R18 0.080493799416719816 X3 R19 -0.080493799416719816 X3 R20 0.080493799416719816 X3 R21 -0.080493799416719816 X3 R22 0.080493799416719816 X3 R23 -1 X3 R24 1 X3 R31 1 X3 R32 -1 X3 R33 1 X3 R34 -1 X3 R35 1 X3 R36 -1 X3 R43 0.080493799416719816 X3 R44 -0.080493799416719816 X3 R47 -1 X3 R48 1 X3 R55 -0.080493799416719816 X3 R56 0.080493799416719816 X3 R57 -1.0804937994167199 X3 R58 1.0804937994167199 X3 R59 -0.080493799416719816 X3 R60 0.080493799416719816 X3 R61 -0.080493799416719816 X3 R62 0.080493799416719816 X3 R63 -0.080493799416719816 X3 R64 0.080493799416719816 X3 R65 -1 X3 R66 1 X3 R73 1 X3 R74 -1 X3 R75 1 X3 R76 -1 X3 R77 1 X3 R78 -1 X4 OBJECT -0.12103523796485483 X4 R1 0.52049963119723075 X4 R2 -0.52049963119723075 X4 R7 -1 X4 R8 1 X4 R13 -0.52049963119723075 X4 R14 0.52049963119723075 X4 R15 -0.52049963119723075 X4 R16 0.52049963119723075 X4 R17 -1.5204996311972308 X4 R18 1.5204996311972308 X4 R19 -0.52049963119723075 X4 R20 0.52049963119723075 X4 R21 -0.52049963119723075 X4 R22 0.52049963119723075 X4 R25 -1 X4 R26 1 X4 R31 -1 X4 R32 1 X4 R37 1 X4 R38 -1 X4 R39 1 X4 R40 -1 X4 R43 0.52049963119723075 X4 R44 -0.52049963119723075 X4 R49 -1 X4 R50 1 X4 R55 -0.52049963119723075 X4 R56 0.52049963119723075 X4 R57 -0.52049963119723075 X4 R58 0.52049963119723075 X4 R59 -1.5204996311972308 X4 R60 1.5204996311972308 X4 R61 -0.52049963119723075 X4 R62 0.52049963119723075 X4 R63 -0.52049963119723075 X4 R64 0.52049963119723075 X4 R67 -1 X4 R68 1 X4 R73 -1 X4 R74 1 X4 R79 1 X4 R80 -1 X4 R81 1 X4 R82 -1 X5 OBJECT -0.28366946676513416 X5 R1 0.073229375252495135 X5 R2 -0.073229375252495135 X5 R9 -1 X5 R10 1 X5 R13 -0.073229375252495135 X5 R14 0.073229375252495135 X5 R15 -0.073229375252495135 X5 R16 0.073229375252495135 X5 R17 -0.073229375252495135 X5 R18 0.073229375252495135 X5 R19 -1.0732293752524951 X5 R20 1.0732293752524951 X5 R21 -0.073229375252495135 X5 R22 0.073229375252495135 X5 R27 -1 X5 R28 1 X5 R33 -1 X5 R34 1 X5 R37 -1 X5 R38 1 X5 R41 1 X5 R42 -1 X5 R43 0.073229375252495135 X5 R44 -0.073229375252495135 X5 R51 -1 X5 R52 1 X5 R55 -0.073229375252495135 X5 R56 0.073229375252495135 X5 R57 -0.073229375252495135 X5 R58 0.073229375252495135 X5 R59 -0.073229375252495135 X5 R60 0.073229375252495135 X5 R61 -1.0732293752524951 X5 R62 1.0732293752524951 X5 R63 -0.073229375252495135 X5 R64 0.073229375252495135 X5 R69 -1 X5 R70 1 X5 R75 -1 X5 R76 1 X5 R79 -1 X5 R80 1 X5 R83 1 X5 R84 -1 X6 OBJECT -0.66742026208402516 X6 R1 0.53643699285113067 X6 R2 -0.53643699285113067 X6 R11 -1 X6 R12 1 X6 R13 -0.53643699285113067 X6 R14 0.53643699285113067 X6 R15 -0.53643699285113067 X6 R16 0.53643699285113067 X6 R17 -0.53643699285113067 X6 R18 0.53643699285113067 X6 R19 -0.53643699285113067 X6 R20 0.53643699285113067 X6 R21 -1.5364369928511308 X6 R22 1.5364369928511308 X6 R29 -1 X6 R30 1 X6 R35 -1 X6 R36 1 X6 R39 -1 X6 R40 1 X6 R41 -1 X6 R42 1 X6 R43 0.53643699285113067 X6 R44 -0.53643699285113067 X6 R53 -1 X6 R54 1 X6 R55 -0.53643699285113067 X6 R56 0.53643699285113067 X6 R57 -0.53643699285113067 X6 R58 0.53643699285113067 X6 R59 -0.53643699285113067 X6 R60 0.53643699285113067 X6 R61 -0.53643699285113067 X6 R62 0.53643699285113067 X6 R63 -1.5364369928511308 X6 R64 1.5364369928511308 X6 R71 -1 X6 R72 1 X6 R77 -1 X6 R78 1 X6 R81 -1 X6 R82 1 X6 R83 -1 X6 R84 1 RHS RHS1 R1 -1.9832427778049702 RHS1 R2 -1.9867572221950303 RHS1 R3 0.00070519999999962835 RHS1 R4 -3.9707051999999994 RHS1 R5 -0.0005845000000004319 RHS1 R6 -3.9694154999999984 RHS1 R7 0.0011155999999985511 RHS1 R8 -3.9711155999999974 RHS1 R9 0.00024540000000072837 RHS1 R10 -3.9702453999999996 RHS1 R11 0.0012371000000044319 RHS1 R12 -3.9712371000000033 RHS1 R13 -0.0010520221950309772 RHS1 R14 -3.9689479778049694 RHS1 R15 -0.0023417221950321476 RHS1 R16 -3.9676582778049676 RHS1 R17 -0.00064162219502961193 RHS1 R18 -3.9693583778049693 RHS1 R19 -0.001511822195029211 RHS1 R20 -3.9684881778049714 RHS1 R21 -0.00052012219503083656 RHS1 R22 -3.969479877804968 RHS1 R23 -1.9862896999999997 RHS1 R24 -1.9837103000000003 RHS1 R25 -1.9845895999999996 RHS1 R26 -1.985410400000001 RHS1 R27 -1.9854597999999992 RHS1 R28 -1.9845401999999996 RHS1 R29 -1.9844680999999991 RHS1 R30 -1.9855319000000016 RHS1 R31 -1.9832998999999998 RHS1 R32 -1.9867001000000002 RHS1 R33 -1.9841700999999992 RHS1 R34 -1.9858299000000006 RHS1 R35 -1.9831783999999981 RHS1 R36 -1.9868216000000007 RHS1 R37 -1.9858701999999997 RHS1 R38 -1.9841298000000007 RHS1 R39 -1.9848784999999998 RHS1 R40 -1.9851215000000009 RHS1 R41 -1.9840083 RHS1 R42 -1.9859917000000005 RHS1 R43 -23.773242777804967 RHS1 R44 -16.176757222195029 RHS1 R45 -25.589294799999998 RHS1 R46 -14.3607052 RHS1 R47 -29.3905845 RHS1 R48 -10.559415499999997 RHS1 R49 -33.1888844 RHS1 R50 -6.7611155999999966 RHS1 R51 -36.9897546 RHS1 R52 -2.960245399999998 RHS1 R53 -40.788762899999995 RHS1 R54 0.838762899999999 RHS1 R55 -21.791052022195029 RHS1 R56 -18.158947977804967 RHS1 R57 -25.592341722195027 RHS1 R58 -14.357658277804969 RHS1 R59 -29.39064162219503 RHS1 R60 -10.559358377804966 RHS1 R61 -33.191511822195025 RHS1 R62 -6.758488177804967 RHS1 R63 -36.990520122195029 RHS1 R64 -2.9594798778049665 RHS1 R65 -23.776289699999996 RHS1 R66 -16.1737103 RHS1 R67 -27.5745896 RHS1 R68 -12.375410399999996 RHS1 R69 -31.375459799999994 RHS1 R70 -8.5745402 RHS1 R71 -35.1744681 RHS1 R72 -4.7755318999999972 RHS1 R73 -23.773299899999998 RHS1 R74 -16.176700099999998 RHS1 R75 -27.574170099999996 RHS1 R76 -12.3758299 RHS1 R77 -31.373178399999993 RHS1 R78 -8.5768216 RHS1 R79 -23.775870199999996 RHS1 R80 -16.1741298 RHS1 R81 -27.574878499999997 RHS1 R82 -12.375121499999999 RHS1 R83 -23.7740083 RHS1 R84 -16.175991699999997 BOUNDS FR BOUND X1 FR BOUND X2 FR BOUND X3 FR BOUND X4 FR BOUND X5 FR BOUND X6 QUADOBJ X1 X1 0.12773471343638093 X2 X1 0.10196952961371175 X3 X1 -0.0093292386606955227 X4 X1 0.015355698544653802 X5 X1 0.00179148994111826 X6 X1 0.022234396172699303 X2 X2 0.27053515944157813 X3 X2 -0.0044801782928801 X4 X2 0.071346059989825059 X5 X2 0.025238285238904652 X6 X2 0.072776031202962388 X3 X3 0.025804696537700223 X4 X3 -0.0025202318116454596 X5 X3 0.0061732011452994137 X6 X3 -0.0046901140692241858 X4 X4 0.086962424907745722 X5 X4 0.0052205577718211216 X6 X4 -0.012441625024233263 X5 X5 0.028499112283986966 X6 X5 0.0035736918205115237 X6 X6 0.07904798363568713 ENDATA

[jjhforrest]

Looking at comments in barrier code, quadratic barrier seems a work in progress (or a mess to put it another way). As I said the initial solve part of Clp does not seem to have been updated for quadratic. Possibly the best solution is to keep it simple. There is an example testQP.cpp in Clp/examples which seems to work on all of your problems.

[yk011309]

I understand that the contents of testQP.cpp need to be included in clp.exe to solve quadratic, is my understanding correct? (I'm using a translation system, so I'm worried if I understand correctly. )

It's difficult for me, so I'd like to wait for it to be implemented someday.

[jjhforrest]

I have found a bug and will fix it. Clp primal minimises the true objective plus a large weight times the sum of infeasibilities. The user can set this weight but if the code finds that there is no objective it sets it to 1.0. With quadratic the code first solves the problem just with the linear objective and then starts using the quadratic objective. The bug was that it thought there was no objective (in linear part) and then leaves weight at 1.0. This can lead to problems. readMps as default treats all values <1.0e-14 as zero and in the first problem all values in linear objective are smaller than that. If I change X4 OBJECT -1.15098902631061E-15 to X4 OBJECT -1.15098902631061E-14 the code worked. See if scaling your problem so largest entry in linear objective is >1.0e-14 solves problem.

[jjhforrest]

Just noticed that bug had already been fixed in master.

[yk011309]

Thanks for the additional information. I would like to experiment with my test case. However, the deadline for my main business is approaching, so I plan to start the experiment in about three weeks.