Massive form factors at $\mathcal{O}({\alpha }_s^3)$

We report on our recent calculation of massive quark form factors using a semi-numerical approach based on series expansions of the master integrals around singular and regular kinematic points and numerical matching. The methods allows to cover the whole kinematic range of negative and positive values of the virtuality $s$ with at least seven significant digits accuracy.