Full vectorial simulation of multilayer anisotropic waveguides with an accurate and automated finite-element program