def test_binary_spacegroup( self ):
if ( not has_ase_with_ce ):
self.skipTest("ASE version does not have CE")
return
bs, db_name = get_bulkspacegroup_binary()
cf = CorrFunction( bs )
cf_vals = cf.get_cf( bs.atoms )
ecis = {name:1.0 for name in cf_vals.keys()}
calc = CE( bs, ecis )
bs.atoms.set_calculator( calc )
for i in range(25):
if ( bs.atoms[i].symbol == "Al" ):
new_symb = "Mg"
old_symb = "Al"
else:
new_symb = "Al"
old_symb = "Mg"
calc.calculate( bs.atoms, ["energy"], [(i,old_symb,new_symb)] )
updated_cf = calc.get_cf()
brute_force = cf.get_cf_by_cluster_names( bs.atoms, updated_cf.keys() )
for key,value in brute_force.items():
self.assertAlmostEqual( value, updated_cf[key] )
def test_double_swaps_ternary( self ):
if ( not has_ase_with_ce ): # Disable this test
self.skipTest("ASE version has not cluster expansion")
return
db_name = "test_db_ternary.db"
conc_args = {
"conc_ratio_min_1":[[4,0,0]],
"conc_ratio_max_1":[[0,4,0]],
"conc_ratio_min_1":[[2,2,0]],
"conc_ratio_max_2":[[1,1,2]]
}
max_dia = get_max_cluster_dia_name()
size_arg = {max_dia:4.05}
ceBulk = CEBulk( crystalstructure="fcc", a=4.05, size=[4,4,4], basis_elements=[["Al","Mg","Si"]], \
conc_args=conc_args, db_name=db_name, max_cluster_size=3, **size_arg)
ceBulk.reconfigure_settings()
corr_func = CorrFunction( ceBulk )
cf = corr_func.get_cf( ceBulk.atoms )
#prefixes = [name.rpartition("_")[0] for name in cf.keys()]
#prefixes.remove("")
eci = {name:1.0 for name in cf.keys()}
calc = CE( ceBulk, eci )
n_tests = 10
# Insert 25 Mg atoms and 25 Si atoms
n = 18
for i in range(n):
calc.calculate( ceBulk.atoms, ["energy"], [(i,"Al","Mg")])
calc.calculate( ceBulk.atoms, ["energy"], [(i+n,"Al","Si")])
updated_cf = calc.get_cf()
brute_force = corr_func.get_cf_by_cluster_names( ceBulk.atoms, updated_cf.keys() )
for key in updated_cf.keys():
self.assertAlmostEqual( brute_force[key], updated_cf[key])
# Swap atoms
for i in range(n_tests):
indx1 = np.random.randint(low=0,high=len(ceBulk.atoms))
symb1 = ceBulk.atoms[indx1].symbol
indx2 = indx1
symb2 = symb1
while( symb2 == symb1 ):
indx2 = np.random.randint( low=0, high=len(ceBulk.atoms) )
symb2 = ceBulk.atoms[indx2].symbol
calc.calculate( ceBulk.atoms, ["energy"], [(indx1,symb1,symb2),(indx2,symb2,symb1)])
update_cf = calc.get_cf()
brute_force = corr_func.get_cf_by_cluster_names( ceBulk.atoms, update_cf.keys() )
for key,value in brute_force.items():
self.assertAlmostEqual( value, update_cf[key])