I am wondering how to integrate
$$ I(a, b) = \int_{-\infty}^{\infty}dx \frac{\sin(x^2 + a)}{x^2 + a} \exp(i[x-b]^2) $$
Thus far I tried integration by parts and contour integration but could not find a solution. This integral arises in connection with the phase-matching condition in nonlinear crystals. I suspect there is no solution in closed form but an expression for $I(a,b)$ that is suitable for efficient numerical computation would be most appreciated.