We present a novel method for optimizing parameters of the Ignition and Growth (I&G) reactive flow model for high explosives. The I&G model predicts the shock initiation response of explosives subjected to dynamic loading. However, calibrating the model is a time-consuming process because it has several parameters, typically requires multiple data sets, and there are no direct correlations between parameters and data. In this study, we couple a scalable optimization algorithm to simulations of shock initiation experiments in the multi-physics code ALE3D. We develop four I&G model parameter sets for the HMX (1,3,5,7-tetranitro-1,3,5,7-tetrazocane)-based explosive LX-07 [90 wt.% HMX, 10 wt.% Viton A] based on minimizing the difference between calculations and measurements at the embedded pressure gauges in 1D gun shock experiments. Our study shows that a cost function based on both shock time of arrival and pressure pulse shape demonstrates the best agreement with experimental data. Published by AIP Publishing.