Although it is not best solution, but it will work.
scipy will work without changing torch tensor to numpy as eg.
import scipy.linalg
m = torch.tensor([[.5,.5],[.7,.9]])
print(scipy.linalg.fractional_matrix_power(m, (-1/2)))

array([[ 2.69776664, -1.10907208],
       [-1.55270001,  1.81051025]])