最適化問題に対する超高速&安定計算

大規模最適化問題、グラフ探索、機械学習やデジタルツインなどの研究のお話が中心

Lovász number

2014年03月08日 00時04分15秒 | Weblog
日本OR学会の春季発表会でグラフ彩色問題に関する発表があったので、以下のグラフに対して Lovász number を計算してみました。objValDual が Lovász number になります。

今回用いたグラフデータ
http://mat.gsia.cmu.edu/COLOR03/INSTANCES/

◯DSJC125.9.col :
SDPA start at [Fri Mar 7 10:05:12 2014]
param is ./param.sdpa
data is /home/fujisawa/src/makepro/DSJC125.9.col.dat-s : sparse
out is out.DSJC125.9
NumThreads is set as 40
Schur computation : DENSE
mu thetaP thetaD objP objD alphaP alphaD beta
0 1.0e+04 1.0e+00 1.0e+00 -0.00e+00 +1.25e+04 1.0e+00 8.5e-01 2.00e-01
1 2.0e+03 0.0e+00 1.5e-01 +1.42e+02 +1.12e+04 7.5e-01 7.5e-01 2.00e-01
2 6.6e+02 1.1e-18 3.7e-02 +1.81e+02 +1.14e+03 7.8e-01 7.8e-01 2.00e-01
3 1.9e+02 1.6e-18 8.0e-03 +2.35e+02 +1.14e+02 8.9e-01 8.9e-01 2.00e-01
4 3.0e+01 2.2e-18 8.9e-04 +3.13e+02 +1.15e+01 9.8e-01 9.8e-01 2.00e-01
5 3.8e+00 2.2e-18 2.1e-05 +3.82e+02 +1.76e+00 9.3e-01 1.0e+00 2.00e-01
6 1.0e+00 3.3e-18 2.7e-20 +1.30e+02 +2.10e+00 6.9e-01 1.2e+00 1.00e-01
7 3.4e-01 1.8e-17 8.9e-21 +5.92e+01 +1.73e+01 7.2e-01 8.4e-01 1.00e-01
8 1.0e-01 1.8e-17 1.8e-20 +4.32e+01 +3.04e+01 8.3e-01 7.7e-01 1.00e-01
9 2.9e-02 3.5e-17 8.8e-22 +3.90e+01 +3.54e+01 8.7e-01 8.5e-01 1.00e-01
10 6.6e-03 4.4e-17 1.8e-20 +3.80e+01 +3.72e+01 9.1e-01 9.0e-01 1.00e-01
11 1.2e-03 4.4e-17 1.8e-20 +3.78e+01 +3.76e+01 9.2e-01 9.5e-01 1.00e-01
12 1.9e-04 5.3e-17 5.3e-20 +3.78e+01 +3.77e+01 9.3e-01 9.7e-01 1.00e-01
13 2.5e-05 6.2e-17 3.6e-20 +3.78e+01 +3.78e+01 9.4e-01 9.4e-01 1.00e-01
14 3.8e-06 7.0e-17 2.0e-19 +3.78e+01 +3.78e+01 9.6e-01 9.8e-01 1.00e-01
15 4.7e-07 7.0e-17 6.3e-19 +3.78e+01 +3.78e+01 9.5e-01 1.0e+00 1.00e-01
16 5.2e-08 7.0e-17 1.3e-18 +3.78e+01 +3.78e+01 9.5e-01 1.0e+00 1.00e-01
17 5.4e-09 8.8e-17 1.4e-19 +3.78e+01 +3.78e+01 9.5e-01 1.0e+00 1.00e-01

phase.value = pdOPT
Iteration = 17
mu = +5.4098824491255467e-09
relative gap = +1.7905075659021317e-08
gap = +6.7623518873460853e-07
digits = +7.7470238395096018e+00
objValPrimal = +3.7767793203817845e+01
objValDual = +3.7767792527582657e+01
p.feas.error = +8.8817841970012523e-15
d.feas.error = +1.7763568394002505e-15
total time = 0.317898


◯ DSJC250.9.col
SDPA start at [Fri Mar 7 10:05:34 2014]
param is ./param.sdpa
data is /home/fujisawa/src/makepro/DSJC250.9.col.dat-s : sparse
out is out.DSJC250.9
NumThreads is set as 40
Schur computation : DENSE
mu thetaP thetaD objP objD alphaP alphaD beta
0 1.0e+04 1.0e+00 1.0e+00 -0.00e+00 +2.50e+04 6.9e-01 7.6e-01 2.00e-01
1 2.9e+03 3.1e-01 2.4e-01 +9.89e+01 +8.21e+04 1.0e+00 7.1e-01 2.00e-01
2 1.2e+03 2.2e-18 6.9e-02 +1.80e+02 +9.10e+03 7.2e-01 7.2e-01 2.00e-01
3 4.4e+02 4.4e-18 1.9e-02 +2.28e+02 +9.21e+02 7.9e-01 7.9e-01 2.00e-01
4 1.2e+02 4.4e-18 4.2e-03 +2.98e+02 +9.51e+01 9.1e-01 9.1e-01 2.00e-01
5 1.7e+01 6.6e-18 3.9e-04 +4.00e+02 +1.29e+01 9.9e-01 9.9e-01 2.00e-01
6 2.2e+00 8.8e-18 5.0e-06 +4.81e+02 +1.68e+00 7.3e-01 1.0e+00 2.00e-01
7 8.5e-01 8.8e-18 4.4e-21 +2.15e+02 +3.26e+00 7.0e-01 1.3e+00 1.00e-01
8 2.7e-01 1.8e-17 8.9e-21 +9.38e+01 +2.52e+01 7.4e-01 9.2e-01 1.00e-01
9 7.7e-02 3.5e-17 2.2e-22 +6.33e+01 +4.41e+01 8.4e-01 7.9e-01 1.00e-01
10 2.1e-02 7.0e-17 1.8e-20 +5.69e+01 +5.17e+01 8.7e-01 8.8e-01 1.00e-01
11 4.3e-03 7.0e-17 2.2e-20 +5.55e+01 +5.44e+01 9.1e-01 9.4e-01 1.00e-01
12 7.1e-04 7.0e-17 4.4e-21 +5.52e+01 +5.50e+01 9.1e-01 9.5e-01 1.00e-01
13 1.1e-04 7.0e-17 1.8e-20 +5.52e+01 +5.51e+01 9.2e-01 9.6e-01 1.00e-01
14 1.6e-05 7.0e-17 7.1e-20 +5.52e+01 +5.51e+01 9.2e-01 9.4e-01 1.00e-01
15 2.5e-06 7.0e-17 5.3e-20 +5.52e+01 +5.52e+01 9.4e-01 9.7e-01 1.00e-01
16 3.3e-07 7.0e-17 4.9e-20 +5.52e+01 +5.52e+01 9.5e-01 1.0e+00 1.00e-01
17 3.6e-08 7.0e-17 8.3e-19 +5.52e+01 +5.52e+01 9.6e-01 9.1e-01 1.00e-01
18 6.1e-09 8.8e-17 2.4e-18 +5.52e+01 +5.52e+01 9.6e-01 9.1e-01 1.00e-01

phase.value = pdOPT
Iteration = 18
mu = +6.1178430630359770e-09
relative gap = +2.7731202555920717e-08
gap = +1.5294516870767438e-06
digits = +7.5570312969538191e+00
objValPrimal = +5.5152736532046553e+01
objValDual = +5.5152735002594866e+01
p.feas.error = +8.8817841970012523e-15
d.feas.error = +6.0618177144533547e-14
total time = 3.607793


◯DSJC500.9.col
SDPA start at [Fri Mar 7 10:23:39 2014]
param is ./param.sdpa
data is /home/fujisawa/src/makepro/DSJC500.9.col.dat-s : sparse
out is out.DSJC500.9.col
NumThreads is set as 40
Schur computation : DENSE
mu thetaP thetaD objP objD alphaP alphaD beta
0 1.0e+04 1.0e+00 1.0e+00 -0.00e+00 +5.00e+04 4.3e-01 5.9e-01 2.00e-01
1 4.5e+03 5.7e-01 4.1e-01 +6.26e+01 +5.27e+05 1.0e+00 6.7e-01 2.00e-01
2 2.4e+03 1.8e-17 1.4e-01 +1.84e+02 +5.70e+04 7.3e-01 7.3e-01 2.00e-01
3 8.6e+02 3.5e-17 3.7e-02 +2.32e+02 +5.93e+03 7.4e-01 7.4e-01 2.00e-01
4 2.9e+02 3.5e-17 9.9e-03 +2.98e+02 +5.94e+02 8.6e-01 8.6e-01 2.00e-01
5 5.7e+01 3.5e-17 1.4e-03 +3.98e+02 +6.70e+01 9.1e-01 9.1e-01 2.00e-01
6 8.0e+00 3.5e-17 1.3e-04 +5.31e+02 +9.43e+00 1.0e+00 1.0e+00 2.00e-01
7 1.2e+00 3.5e-17 5.3e-20 +5.99e+02 +2.28e+00 5.9e-01 4.8e+00 1.00e-01
8 5.1e-01 3.5e-17 2.6e-19 +2.85e+02 +3.27e+01 7.9e-01 1.6e+00 1.00e-01
9 1.2e-01 3.5e-17 1.1e-19 +1.19e+02 +5.79e+01 7.7e-01 8.0e-01 1.00e-01
10 3.6e-02 3.5e-17 3.6e-20 +9.09e+01 +7.27e+01 8.4e-01 8.3e-01 1.00e-01
11 9.1e-03 7.0e-17 1.8e-20 +8.55e+01 +8.09e+01 8.8e-01 8.9e-01 1.00e-01
12 1.8e-03 7.0e-17 3.6e-20 +8.43e+01 +8.34e+01 9.1e-01 9.3e-01 1.00e-01
13 3.1e-04 7.0e-17 3.8e-20 +8.41e+01 +8.39e+01 9.4e-01 9.4e-01 1.00e-01
14 4.7e-05 7.0e-17 4.4e-21 +8.40e+01 +8.40e+01 9.6e-01 9.4e-01 1.00e-01
15 7.3e-06 8.8e-17 6.7e-20 +8.40e+01 +8.40e+01 9.6e-01 9.4e-01 1.00e-01
16 1.1e-06 8.8e-17 9.8e-20 +8.40e+01 +8.40e+01 9.7e-01 9.3e-01 1.00e-01
17 1.7e-07 1.1e-16 1.3e-20 +8.40e+01 +8.40e+01 1.0e+00 9.5e-01 1.00e-01
18 2.4e-08 1.1e-16 6.7e-21 +8.40e+01 +8.40e+01 1.0e+00 9.8e-01 1.00e-01
19 2.7e-09 1.1e-16 2.2e-20 +8.40e+01 +8.40e+01 1.0e+00 9.8e-01 1.00e-01

phase.value = pdOPT
Iteration = 19
mu = +2.6693292243251675e-09
relative gap = +1.5880926174760332e-08
gap = +1.3346641480893595e-06
digits = +7.7991241731395036e+00
objValPrimal = +8.4041959769852852e+01
objValDual = +8.4041958435188704e+01
p.feas.error = +1.0658141036401503e-14
d.feas.error = +1.1102230246251565e-15
total time = 88.775853


◯DSJC1000.9.col
SDPA start at [Fri Mar 7 10:26:47 2014]
param is ./param.sdpa
data is /home/fujisawa/data/DSJC1000.9.col.dat-s : sparse
out is out.DSJC1000.9.col
NumNodes is set as 16
NumThreads is set as 3
Schur computation : DENSE
mu thetaP thetaD objP objD alphaP alphaD beta
0 1.0e+04 1.0e+00 1.0e+00 -0.00e+00 +1.00e+05 7.3e-01 3.6e-01 2.00e-01
1 6.9e+03 2.7e-01 6.4e-01 +1.10e+02 +2.59e+06 1.0e+00 6.2e-01 2.00e-01
2 4.1e+03 8.8e-18 2.5e-01 +1.92e+02 +6.52e+05 7.8e-01 7.8e-01 2.00e-01
3 1.2e+03 1.8e-17 5.4e-02 +2.37e+02 +7.06e+04 7.7e-01 7.7e-01 2.00e-01
4 3.7e+02 2.6e-17 1.2e-02 +3.03e+02 +7.56e+03 7.7e-01 7.7e-01 2.00e-01
5 1.1e+02 3.5e-17 2.9e-03 +3.91e+02 +7.74e+02 8.0e-01 8.0e-01 2.00e-01
6 3.0e+01 5.6e-16 5.9e-04 +5.13e+02 +7.85e+01 8.7e-01 8.7e-01 2.00e-01
7 5.6e+00 1.1e-15 7.4e-05 +6.75e+02 +9.94e+00 1.0e+00 1.0e+00 2.00e-01
8 7.7e-01 1.1e-15 4.0e-20 +7.77e+02 +3.09e+00 6.4e-01 3.4e+00 1.00e-01
9 2.9e-01 1.1e-15 1.4e-19 +3.39e+02 +4.86e+01 7.5e-01 1.4e+00 1.00e-01
10 7.6e-02 1.1e-15 4.4e-21 +1.63e+02 +8.77e+01 7.8e-01 7.9e-01 1.00e-01
11 2.2e-02 1.1e-15 8.9e-21 +1.31e+02 +1.09e+02 8.6e-01 8.4e-01 1.00e-01
12 5.3e-03 1.1e-15 1.6e-20 +1.24e+02 +1.19e+02 8.9e-01 9.0e-01 1.00e-01
13 1.0e-03 1.1e-15 2.7e-20 +1.23e+02 +1.22e+02 9.2e-01 9.3e-01 1.00e-01
14 1.7e-04 1.1e-15 1.2e-20 +1.23e+02 +1.23e+02 9.4e-01 9.4e-01 1.00e-01
15 2.7e-05 1.1e-15 5.1e-20 +1.23e+02 +1.23e+02 9.6e-01 9.4e-01 1.00e-01
16 4.2e-06 1.1e-15 7.4e-22 +1.23e+02 +1.23e+02 9.6e-01 9.4e-01 1.00e-01
17 6.4e-07 1.1e-15 2.2e-21 +1.23e+02 +1.23e+02 9.5e-01 9.3e-01 1.00e-01
18 1.0e-07 1.1e-15 2.7e-20 +1.23e+02 +1.23e+02 9.3e-01 9.3e-01 1.00e-01
19 1.6e-08 1.1e-15 1.8e-20 +1.23e+02 +1.23e+02 9.5e-01 9.5e-01 1.00e-01
20 2.4e-09 1.1e-15 1.4e-19 +1.23e+02 +1.23e+02 9.8e-01 9.8e-01 1.00e-01
21 2.7e-10 1.1e-15 1.2e-19 +1.23e+02 +1.23e+02 9.8e-01 9.8e-01 1.00e-01

phase.value = pdOPT
Iteration = 21
mu = +2.7154975473786538e-10
relative gap = +2.2136155263879887e-09
gap = +2.7154123927175533e-07
digits = +8.6548978076975249e+00
objValPrimal = +1.2266865511888591e+02
objValDual = +1.2266865484734467e+02
p.feas.error = +1.1368683772161603e-13
d.feas.error = +1.1990408665951691e-14
total time = 996.540333
コメント
  • X
  • Facebookでシェアする
  • はてなブックマークに追加する
  • LINEでシェアする