Numerical analysis of the ultrarelativistic and magnetized Bondi--Hoyle problem