A finite-element-based vectorial optical mode solver, furnished with Bayliss-Gunzburger-Turkel-like transparent boundary conditions, is used to rigorously analyze photonic crystal fibers (PCFs). Both the real and imaginary part of the modal indices can be computed in a relatively small computational domain. The leakage loss, the dispersion properties, the vectorial character, as well as the degeneracy of modes of the fibers can be studied through the finite element results. Results for PCFs with either circular or non-circular microstructured holes, solidor air-core will be presented, including the air-core air-silica Bragg fiber. Using the mode solver, the single-modeness of a commercial endlessly single-mode PCF was also investigated