Relativistically extended Blanchard recurrence relation for hydrogenic matrix elements

General recurrence relations for arbitrary non-diagonal, radial hydrogenic matrix elements are derived in Dirac relativistic quantum mechanics. Our approach is based on a generalization of the second hypervirial method previously employed in the non-relativistic Schr\"odinger case. A relativistic version of the Pasternack-Sternheimer relation is thence obtained in the diagonal (i.e. total angular momentum and parity the same) case, from such relation an expression for the relativistic virial theorem is deduced. To contribute to the utility of the relations, explicit expressions for the radial matrix elements of functions of the form $r^\lambda$ and $\beta r^\lambda$ ---where $\beta$ is a Dirac matrix--- are presented.