FeynHiggsFast: a program for a fast calculation of masses and mixing   angles in the Higgs Sector of the MSSM
  FeynHiggsFast is a Fortran code for the calculation of the masses and the mixing angle of the neutral CP-even Higgs bosons in the MSSM up to two-loop order. It is based on a compact analytical approximation formula of the complete diagrammatic one-loop and the dominant two-loop contributions. At the one-loop level a leading logarithmic result is used, taking into account all sectors of the MSSM. At the two-loop level at O(alpha alpha_s) the leading logarithmic and non-logarithmic contributions are taken into account. The approximation formula is valid for arbitrary choices of the parameters in the Higgs sector of the model. Comparing its quality to the full diagrammatic result, we find agreement better than 2 GeV for most parts of the MSSM parameter space. 
S. Heinemeyer, W. Hollik, G. Weiglein