Lennard-Jones 38

SystemLennard-Jones
Atoms38
Runs200
Cap10000 force calls

The main metric is Force calls. Scoring rules.

Charts

Success rate on Lennard-Jones 38, in rank order. The rate is the share of runs that finished. 9 ranked rows did not record a success rate and are not drawn.
Success rate on Lennard-Jones 38, in rank order. The rate is the share of runs that finished. 9 ranked rows did not record a success rate and are not drawn. PDF.
Dolan-Moré profile for Minimum mode on Lennard-Jones 38. The curve is the fraction of starts whose cost is at most tau times the best finite cost in this family. Tau is 1, 2, 4, 8, and 16. An unsolved start stays in the denominator and does not cross a finite tau. Dolan, E. D. and Moré, J. J., Math. Program. 91:201-213, 2002. Colored lines are the six largest areas. The other curves are gray.
Dolan-Moré profile for Minimum mode on Lennard-Jones 38. The curve is the fraction of starts whose cost is at most tau times the best finite cost in this family. Tau is 1, 2, 4, 8, and 16. An unsolved start stays in the denominator and does not cross a finite tau. Dolan, E. D. and Moré, J. J., Math. Program. 91:201-213, 2002. Colored lines are the six largest areas. The other curves are gray. PDF.
Force calls on each recorded start for Lennard-Jones 38. The box is the quartiles of those counts. A blue point is a success and a vermillion cross is a recorded failure (464). The horizontal position is the recorded count. The vertical offset only keeps marks from covering one another.
Force calls on each recorded start for Lennard-Jones 38. The box is the quartiles of those counts. A blue point is a success and a vermillion cross is a recorded failure (464). The horizontal position is the recorded count. The vertical offset only keeps marks from covering one another. PDF.
TSASE on Lennard-Jones 38. Each point is one recorded version of one method. Versions run left to right in numeric order. PAR-2 is omitted where that row has none. No line is drawn across a missing value.
TSASE on Lennard-Jones 38. Each point is one recorded version of one method. Versions run left to right in numeric order. PAR-2 is omitted where that row has none. No line is drawn across a missing value. PDF.
eOn on Lennard-Jones 38. Each point is one recorded version of one method. Versions run left to right in numeric order. PAR-2 is omitted where that row has none. No line is drawn across a missing value.
eOn on Lennard-Jones 38. Each point is one recorded version of one method. Versions run left to right in numeric order. PAR-2 is omitted where that row has none. No line is drawn across a missing value. PDF.
Performance profile for Lennard-Jones 3800.51124816Cost ratioFraction of startseOn (3.5.0) · Dimer L-BFGSEon svn (r2549) · Dimer L-BFGSeOn (3.5.0) · Dimer CGeOn (3.6.0) · DimerEon svn (r2549) · Dimer CGrgsaddle (0.2.2) · DimerpARTn (1.1.0) · ARTn Lanczospysisyphus (0.1.dev1) · LanczosTSASE (1.0) · SSDimerSoftSaddle (2019.04) · LanczosSoftSaddle (2019.04) · DavidsonSella (0.1.1) · lowest eigenvectorrgsaddle (0.2.2) · LanczoseOn (3.5.0) · Lanczos CGeOn (3.6.0) · LanczoseOn (develop f6817b53) (3.6.0) · LanczosEon svn (r2549) · Lanczos CG
Performance profile. Higher means more starts solved within that cost ratio. Colored lines are the six largest areas. Thin lines are the other 11.
  • Sella (0.1.1) · lowest eigenvector
  • rgsaddle (0.2.2) · Lanczos
  • eOn (3.5.0) · Lanczos CG
  • eOn (3.6.0) · Lanczos
  • eOn (develop f6817b53) (3.6.0) · Lanczos
  • Eon svn (r2549) · Lanczos CG
Cost distribution for Lennard-Jones 38081617.8 to 56.8667: 2156.8667 to 181.675: 0181.675 to 580.41: 0580.41 to 1854.27: 01854.27 to 5923.94: 35923.94 to 18925.6: 217.818925.6
Cost distribution of 26 ranked rows. Logarithmic bins, from 17.8 to 18925.6.

Performance profile SVG Cost and success SVG

Performance profile for Lennard-Jones 38
Performance profile. Select the plot for the full SVG.
Cost and success for Lennard-Jones 38
Cost and success. Select the plot for the full SVG.
Best cost for lj38lj3810¹10²10³10⁴10⁵cost (force calls, log scale)Sella: 17.8Sella17.8OPTIM: 25OPTIM25pele: 25pele25TSASE: 25TSASE25rgsaddle: 26.12rgsaddle26.12eOn: 26.385eOn26.385Eon svn: 26.385Eon svn26.385pARTn: 2619.84pARTn2619.84pysisyphus: 4017.14pysisyphus4017.14SoftSaddle: 18700.6SoftSaddle18700.6
A shorter bar is a lower cost on a log scale shared by the codes of this problem.

Leaderboard

A shorter bar is a lower cost on the log scale for that problem. Select a column heading to sort.
RankCodeBestMedianSuccessPAR-2
1 Sella 17.8 17.8 100.0% 17.8
2 OPTIM 25 25 100.0% 25
3 pele 25 25 100.0% 25
4 TSASE 25 28 100.0% 25
5 rgsaddle 26.12 28.76 100.0% 26.12
6 eOn 26.385 27.3575 100.0% 26.385
7 Eon svn 26.385 28.165 100.0% 26.385
8 pARTn 2619.84 2619.84 87.0% 2619.84
9 pysisyphus 4017.14 4017.14 80.0% 4017.14
10 SoftSaddle 18700.6 18813.1 6.5% 18700.6

The tables below list every ranked row, by method family and overall.

A run short of the full set, or a run that is not the stock method, stays off the ranked table.

Minimum mode

Minimum mode. One row is one code, one version, and one method. A lower cost is better. PAR-2 counts a finished run as the calls it used and a failed run as twice the cap. An archived row shows the average it published, and its PAR-2 applies that rule to the failures the entry records. Hessians is the mean analytic Hessian count. FD calls is that count times one force call per coordinate.
RankRowCostUnitHessiansFD callsRuns
1Sella (0.1.1) · lowest eigenvector17.8force calls00200
2TSASE (410) · tsase-lanczos25force calls00200
3rgsaddle (0.2.2) · Lanczos26.12force calls00200
4Eon svn (r2549) · Lanczos CG26.385force calls00200
5eOn (3.5.0) · Lanczos CG26.385force calls00200
6eOn (3.6.0) · Lanczos26.385force calls00200
7eOn (develop f6817b53) (3.6.0) · Lanczos26.385force calls00200
8TSASE (408) · tsase-dimer-bfgs27force calls00200
9eOn (r2156) · eon-dimer-lbfgs28force calls00200
10eOn (r2156) · eon-lanczos28force calls00200
11Eon svn (r2549) · Dimer L-BFGS28.33force calls00200
12eOn (3.5.0) · Dimer L-BFGS28.33force calls00200
13TSASE (408) · tsase-dimer-cg29force calls00200
14eOn (r2156) · eon-dimer-cg30force calls00200
15Eon svn (r2549) · Dimer CG30.355force calls00200
16eOn (3.5.0) · Dimer CG30.355force calls00200
17eOn (3.6.0) · Dimer30.355force calls00200
18rgsaddle (0.2.2) · Dimer31.4force calls00200
19TSASE (1.0) · SSDimer2420.61force calls00200
20pARTn (1.1.0) · ARTn Lanczos2619.84force calls00200
21pysisyphus (0.1.dev1) · Lanczos4017.14force calls00200
22SoftSaddle (2019.04) · Lanczos18700.6force calls00200
23SoftSaddle (2019.04) · Davidson18925.6force calls00200

Single-ended saddle

Single-ended saddle. One row is one code, one version, and one method. A lower cost is better. PAR-2 counts a finished run as the calls it used and a failed run as twice the cap. An archived row shows the average it published, and its PAR-2 applies that rule to the failures the entry records. Hessians is the mean analytic Hessian count. FD calls is that count times one force call per coordinate.
RankRowCostUnitHessiansFD callsRuns
1Sella (0.1.1) · sella17.8force calls00200

Connect / double-ended

Connect / double-ended. One row is one code, one version, and one method. A lower cost is better. PAR-2 counts a finished run as the calls it used and a failed run as twice the cap. An archived row shows the average it published, and its PAR-2 applies that rule to the failures the entry records. Hessians is the mean analytic Hessian count. FD calls is that count times one force call per coordinate.
RankRowCostUnitHessiansFD callsRuns
1OPTIM · hybrid eigenvector following25force calls00200
2pele · hybrid eigenvector following25force calls00200

Overall

One row is one code, one version, and one method. A lower cost is better. PAR-2 counts a finished run as the calls it used and a failed run as twice the cap. An archived row shows the average it published, and its PAR-2 applies that rule to the failures the entry records. Hessians is the mean analytic Hessian count. FD calls is that count times one force call per coordinate.
RowFamilyFamily rankUnitPAR-2PublishedMedian callsMean callsHessiansFD callsFailuresRuns
Sella (0.1.1) · lowest eigenvectorMinimum mode1force calls17.81717.800crash 0, divergence 0, cap 0, wrong_stationary 0200
Sella (0.1.1) · sellaSingle-ended saddle1force calls17.817.817.8000 failed200
OPTIM · hybrid eigenvector followingConnect / double-ended1force calls25252325000 failed200
TSASE (410) · tsase-lanczosMinimum mode2force calls252525000 failed200
pele · hybrid eigenvector followingConnect / double-ended2force calls252525000 failed200
rgsaddle (0.2.2) · LanczosMinimum mode3force calls26.122626.1200crash 0, divergence 0, cap 0, wrong_stationary 0200
Eon svn (r2549) · Lanczos CGMinimum mode4force calls26.3852526.38500crash 0, divergence 0, cap 0, wrong_stationary 0200
eOn (3.5.0) · Lanczos CGMinimum mode5force calls26.3852526.38500crash 0, divergence 0, cap 0, wrong_stationary 0200
eOn (3.6.0) · LanczosMinimum mode6force calls26.3852526.38500crash 0, divergence 0, cap 0, wrong_stationary 0200
eOn (develop f6817b53) (3.6.0) · LanczosMinimum mode7force calls26.3852526.38500crash 0, divergence 0, cap 0, wrong_stationary 0200
TSASE (408) · tsase-dimer-bfgsMinimum mode8force calls272727000 failed200
eOn (r2156) · eon-dimer-lbfgsMinimum mode9force calls282828000 failed200
eOn (r2156) · eon-lanczosMinimum mode10force calls282828000 failed200
Eon svn (r2549) · Dimer L-BFGSMinimum mode11force calls28.332728.3300crash 0, divergence 0, cap 0, wrong_stationary 0200
eOn (3.5.0) · Dimer L-BFGSMinimum mode12force calls28.332728.3300crash 0, divergence 0, cap 0, wrong_stationary 0200
TSASE (408) · tsase-dimer-cgMinimum mode13force calls292929000 failed200
eOn (r2156) · eon-dimer-cgMinimum mode14force calls303030000 failed200
Eon svn (r2549) · Dimer CGMinimum mode15force calls30.3552930.35500crash 0, divergence 0, cap 0, wrong_stationary 0200
eOn (3.5.0) · Dimer CGMinimum mode16force calls30.3552930.35500crash 0, divergence 0, cap 0, wrong_stationary 0200
eOn (3.6.0) · DimerMinimum mode17force calls30.3552930.35500crash 0, divergence 0, cap 0, wrong_stationary 0200
rgsaddle (0.2.2) · DimerMinimum mode18force calls31.42331.400crash 0, divergence 0, cap 0, wrong_stationary 0200
TSASE (1.0) · SSDimerMinimum mode19force calls2420.61233267.94500crash 0, divergence 22, cap 0, wrong_stationary 0200
pARTn (1.1.0) · ARTn LanczosMinimum mode20force calls2619.842323.9900crash 0, divergence 26, cap 0, wrong_stationary 0200
pysisyphus (0.1.dev1) · LanczosMinimum mode21force calls4017.142222.19500crash 0, divergence 40, cap 0, wrong_stationary 0200
SoftSaddle (2019.04) · LanczosMinimum mode22force calls18700.61918.9100crash 0, divergence 187, cap 0, wrong_stationary 0200
SoftSaddle (2019.04) · DavidsonMinimum mode23force calls18925.6695696.47500crash 0, divergence 189, cap 0, wrong_stationary 0200

Success on this problem: harness calls under the cap and absolute overlap of at least 0.99 with the reference lowest mode.