Numerical simulation of tensorial impedance metasurfaces