The high accuracy and reliability of our method, as demonstrated here, position it for routine crystal structure prediction in drug formulation. In the future, we plan to extend our method to support more complex and/or multi-component systems, such as co-crystals, solvates, hydrates and salts. We also plan to integrate our method with other computational tools, such as solubility, permeability, mechanical properties and crystal morphology predictions, to enable a comprehensive analysis of polymorphs in terms of their structure, stability, function, and performance.