Numerical 3+1 General Relativistic Magnetohydrodynamics: A Local Characteristic Approach