BingqingCheng/ASAP

Broken Examples

Ryan-Rhys opened this issue · 6 comments

Working Examples

  1. B-RSS-train-set
  2. revpbe0-water-train-set
  3. Tungstun-train-set (Tungsten?)

Broken Examples

  1. Ni-icosahedral-np
  2. TiO2-RSS

Fail with the error:

File "/Users/Ryan-Rhys/ml_physics/asap/ASAP/asap/asaplib/kde/ml_density_estimation.py", line 164, in find_optimal_bandwidth assert abs(h_optimal - hmin) > 1e-4, "Lower boundary reached for bandwidth" AssertionError: Lower boundary reached for bandwidth

  1. Lennard-Jones-nucleus
  2. ice-phases

Fail with the error:

File "../../asap/pca.py", line 202, in <module> main(args.fmat, args.fxyz, args.tags, args.colors, args.colorscolumn, args.prefix, args.output, args.peratom, args.keepraw, args.scale, args.d, args.pc1, args.pc2, args.plotatomic, args.adjusttext) File "../../asap/pca.py", line 66, in main if len(desc)==0: raise ValueError('Please supply descriptor in a xyz file or a standlone descriptor matrix')

  1. single-triplet-fission

Fails with the error:

frame_select.py: error: unrecognized arguments: -n 1000

  1. water-nnp-trainset

has no README file.

This particular singlet-triplet-fission error is now fixed on my branch but a new error arises where the variable y_all is undefined in the case that the chosen algorithm is random as opposed to sort:

File "../../asap/frame_select.py", line 75, in <module> main(args.fxyz, args.y, args.prefix, args.n, args.algo, args.mat) File "../../asap/frame_select.py", line 61, in main np.savetxt(prefix+"-"+algorithm+"-n-"+str(nkeep)+'-'+fy, np.asarray(y_all)[sbs], fmt='%4.8f') UnboundLocalError: local variable 'y_all' referenced before assignment

fixed this yall error

  1. water-npp-trainset example is missing a README file.
  2. small molecules example is now broken following the change to the way soap hyperparameters are generated.
  3. TiO2 example currently fails with a bandwidth estimation error in the kernel density estimation routine. This error occurs even using just single dimension.

All other examples working as expected.

For me (on b61e883) the TiO2-RSS example fails with this error:

Traceback (most recent call last):
  File "/home/bonan/miniconda3/envs/asap/bin/pca.py", line 7, in <module>
    exec(compile(f.read(), __file__, 'exec'))
  File "/home/bonan/appdir/ASAP/scripts/pca.py", line 190, in <module>
    main(args.fmat, args.fxyz, args.tags, args.colors, args.colorscolumn, args.prefix, args.output, args.peratom,
  File "/home/bonan/appdir/ASAP/scripts/pca.py", line 94, in main
    plotcolor, colorlabel, colorscale = set_color_function(fcolor, asapxyz, colorscol, len(proj), False)
  File "/home/bonan/appdir/ASAP/asaplib/plot/plot_colors.py", line 18, in set_color_function
    if os.path.isfile(fcolor):
  File "/home/bonan/miniconda3/envs/asap/lib/python3.8/genericpath.py", line 30, in isfile
    st = os.stat(path)
TypeError: stat: path should be string, bytes, os.PathLike or integer, not NoneType

elif os.path.isfile(fcolor):
# load the column=colorscol for color functions
try:
loadcolor = np.genfromtxt(fcolor, dtype=float)
# print(np.shape(loadcolor))
if colorscol > 0 or len(np.shape(loadcolor)) > 1:

change if os.path.isfile(fcolor): to if fcolor and os.path.isfile(fcolor): solves the problem.

But there is another exception raised when saving the plot:

  File "/home/bonan/miniconda3/envs/asap/bin/pca.py", line 7, in <module>
    exec(compile(f.read(), __file__, 'exec'))
  File "/home/bonan/appdir/ASAP/scripts/pca.py", line 190, in <module>
    main(args.fmat, args.fxyz, args.tags, args.colors, args.colorscolumn, args.prefix, args.output, args.peratom,
  File "/home/bonan/appdir/ASAP/scripts/pca.py", line 154, in main
    fig.savefig('PCA_4_' + prefix + '-c-' + fcolor + '.png')
TypeError: can only concatenate str (not "NoneType") to str

@zhubonan Now it's fixed. The error was because the color function was not set. Now I set the default to using index of samples.

I think all the examples @Ryan-Rhys mentioned are working okay now