The traditional methodology to computationally predict elasticity is to apply different strain tensors on optimized atomic structures. However, an applied tensor can allow the system to escape from one local, metastable minimum to another, nearby minimum, which introduces deviations from a systematic exploration of structures in the elasticity calculations. These structural changes that occur during the calculation of a property for a given materials structure are unfortunately difficult to be captured correctly, using only energyminimized structures. Unless a global geometry optimization scheme is employed for each alloy composition, the accuracy of elasticity predictions for a given structure is hard to improve even if more data was included, because it is difficult to tell just from atomic coordinates alone whether a certain structure is metastable or not.