Grid-characteristic numerical method for medical ultrasound