# Set up lattice
vector_a = [0.02, 2.69, 2.7153]*Angstrom
vector_b = [2.7153, 0.0, 2.7153]*Angstrom
vector_c = [2.7153, 2.7153, 0.01]*Angstrom
lattice = UnitCell(vector_a, vector_b, vector_c)

# Define elements
elements = [Silicon, Silicon]

# Define coordinates
fractional_coordinates = [[ 0.007025069439,  0.006459590918, -0.015020491416],
                          [ 0.248715281446,  0.28949668641 ,  0.238503614316]]

# Set up configuration
bulk_configuration = BulkConfiguration(
    bravais_lattice=lattice,
    elements=elements,
    fractional_coordinates=fractional_coordinates
    )

# The current bulk_configuration has lattice vectors and atomic positions 
# that do not obey FCC symetry. This function will produce a new configuration
# that is perfectly symmetric.
bulk_configuration = symmetrizeConfiguration(bulk_configuration, 0.5*Ang)