Describir: Efficient computation of the Nagaoka–Hayashi bound for multiparameter estimation with separable measurements