@@ -303,8 +303,8 @@ using AtomsCalculators: potential_energy
303303
304304function assess_model (model, train_dataset)
305305
306- plot ([- 164 ,- 158 ], [- 164 ,- 158 ]; lc= :black , label= " " )
307-
306+ plot ([- 164 ,- 157 ], [- 164 ,- 157 ]; lc= :black , label= " " )
307+
308308 model_energies = []
309309 model_std = []
310310 for atoms in Si_dataset
@@ -320,8 +320,8 @@ function assess_model(model, train_dataset)
320320 @sprintf (" RMSE (MAE) For Entire Dataset: %.0f (%.0f) meV/atom" , 1000 * rmse, 1000 * mae),
321321 titlefontsize = 8 ,
322322 yerror = model_std,
323- xlabel= " Energy [eV/atom]" , xlims= (- 164 ,- 158 ),
324- ylabel= " Model Energy [eV/atom]" , ylims= (- 164 ,- 158 ),
323+ xlabel= " Energy [eV/atom]" , xlims= (- 164 ,- 157 ),
324+ ylabel= " Model Energy [eV/atom]" , ylims= (- 164 ,- 157 ),
325325 aspect_ratio = :equal , color= 1 )
326326
327327 model_energies = [ustrip (potential_energy (atoms, model)/ length (atoms)) for atoms in train_dataset]
@@ -407,7 +407,7 @@ r_cut = 6.0
407407rdf = ACEpotentials. get_rdf (tial_data, r_cut * u " Å" )
408408plt_TiTi = histogram (rdf[(:Ti , :Ti )], bins= 100 , xlabel = " " , c = 1 ,
409409 ylabel = " RDF - TiTi" , label = " " , yticks = [], xlims = (0 , r_cut) )
410- plt_TiAl = histogram (rdf[(:Ti , :Ti )], bins= 100 , xlabel = " " , c = 2 ,
410+ plt_TiAl = histogram (rdf[(:Ti , :Al )], bins= 100 , xlabel = " " , c = 2 ,
411411 ylabel = " RDF - TiAl" , label = " " , yticks = [], xlims = (0 , r_cut) )
412412plt_AlAl = histogram (rdf[(:Al , :Al )], bins= 100 , xlabel = L " r [\A A]" , c = 3 ,
413413 ylabel = " RDF - AlAl" , label = " " , yticks = [], xlims = (0 , r_cut), )
0 commit comments