Computation of nonlinear ultrasound fields using a linearized contrast source method