def METHOD_NAME(hkl, a, b, c, alpha, beta, gamma, **kws):
    h, k, l = hkl[:, 0], hkl[:, 1], hkl[:, 2]
    alpha, beta, gamma = DEG2RAD*alpha, DEG2RAD*beta, DEG2RAD*gamma
    cos_a, cos_b, cos_g = cos(alpha), cos(beta), cos(gamma)
    x = 1-cos_a**2 - cos_b**2 - cos_g**2 + 2*cos_a*cos_b*cos_g
    y = ((h*sin(alpha)/a)**2 + 2*k*l*(cos_b*cos_g-cos_a)/(b*c) +
         (k*sin(beta)/b)**2  + 2*l*h*(cos_g*cos_a-cos_b)/(c*a) +
         (l*sin(gamma)/c)**2 + 2*h*k*(cos_a*cos_b-cos_g)/(a*b))
    d = np.sqrt(x/y)
    return d