Reliable approximation of the solution of time-harmonic electromagnetic boundary value problems involving bianisotropic media for numerical simulator testing