Hermite-type interpolation of squared-magnitude functions