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

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

Mac OS 10.15.5 上での SDPA

2020年07月06日 21時30分42秒 | Weblog
Mac OS 10.15.5 上での SDPA 7.4.2 の動作を確認出来ました。

gcc version 10.1.0 (Homebrew GCC 10.1.0)
SDPA 7.4.2
MUMPS 5.3.3

export CC=gcc-10
export CXX=g++-10
export F77=gfortran-10
export CFLAGS="-funroll-all-loops -O2 -m64 -fopenmp -DINTSIZE64 -DDIMACS_PRINT"
export CXXFLAGS="-funroll-all-loops -O2 -m64 -fopenmp -DINTSIZE64 -DDIMACS_PRINT"
export FFLAGS="-funroll-all-loops -O2 -m64 -fopenmp -fdefault-integer-8 -fallow-argument-mismatch"
./configure --with-blas="-L./lib -lopenblas -lgomp" --with-lapack="-L./lib -lopenblas -lgomp"
make

$ time ./sdpa ../data/DSJC500.9.col.dat-s out
SDPA start at [Sat Jul 4 21:24:59 2020]
param is ./param.sdpa
data is ../data/DSJC500.9.col.dat-s : sparse
out is out
NumThreads is set as 8
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 2.6e-17 3.7e-02 +2.32e+02 +5.93e+03 7.4e-01 7.4e-01 2.00e-01
4 2.9e+02 5.6e-16 9.9e-03 +2.98e+02 +5.94e+02 8.6e-01 8.6e-01 2.00e-01
5 5.7e+01 5.6e-16 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 4.0e-20 +5.99e+02 +2.28e+00 5.9e-01 4.8e+00 1.00e-01
8 5.1e-01 3.5e-17 1.6e-19 +2.85e+02 +3.27e+01 7.9e-01 1.6e+00 1.00e-01
9 1.2e-01 3.5e-17 5.6e-20 +1.19e+02 +5.79e+01 7.7e-01 8.0e-01 1.00e-01
10 3.6e-02 3.5e-17 3.1e-20 +9.09e+01 +7.27e+01 8.4e-01 8.3e-01 1.00e-01
11 9.1e-03 5.3e-17 1.1e-20 +8.55e+01 +8.09e+01 8.8e-01 8.9e-01 1.00e-01
12 1.8e-03 7.0e-17 2.2e-20 +8.43e+01 +8.34e+01 9.1e-01 9.3e-01 1.00e-01
13 3.1e-04 7.0e-17 2.2e-21 +8.41e+01 +8.39e+01 9.4e-01 9.4e-01 1.00e-01
14 4.7e-05 7.0e-17 1.1e-19 +8.40e+01 +8.40e+01 9.6e-01 9.4e-01 1.00e-01
15 7.3e-06 7.0e-17 8.0e-20 +8.40e+01 +8.40e+01 9.6e-01 9.4e-01 1.00e-01
16 1.1e-06 8.8e-17 5.8e-20 +8.40e+01 +8.40e+01 9.7e-01 9.3e-01 1.00e-01
17 1.7e-07 1.1e-16 7.8e-20 +8.40e+01 +8.40e+01 1.0e+00 9.5e-01 1.00e-01
18 2.4e-08 1.1e-16 1.5e-19 +8.40e+01 +8.40e+01 1.0e+00 9.8e-01 1.00e-01
19 2.7e-09 1.1e-16 1.0e-19 +8.40e+01 +8.40e+01 1.0e+00 9.8e-01 1.00e-01

phase.value = pdOPT
Iteration = 19
mu = +2.6693293513069261e-09
relative gap = +1.5880907912783365e-08
gap = +1.3346626133170503e-06
digits = +7.7991246725486878e+00
objValPrimal = +8.4041959769852852e+01
objValDual = +8.4041958435190239e+01
p.feas.error = +1.0658141036401503e-14
d.feas.error = +5.1070259132757201e-15
total time = 55.933365


コメント
  • X
  • Facebookでシェアする
  • はてなブックマークに追加する
  • LINEでシェアする