def METHOD_NAME(H):
    omega = np.dot(H, H.T)
    x = np.zeros(omega.shape[0]-1)
    y = np.zeros(omega.shape[0])
    sm = np.zeros(omega.shape[0])
    ps = np.zeros(omega.shape[0])
    @jit
    def calc_sm(omega, x, y, sm):  # pragma: no cover
        m = omega.shape[0]
        x_shape = x.shape[0]
        y_shape = y.shape[0]
        count0 = 0
        for k in range(m):
            for i in range(m):
                count0 = 0
                for j in range(m):
                    if j != i:
                        x_ix = j if j < i else j-1
                        x[x_ix] = np.abs(omega[i, k]+omega[j, k])
                        if not x[x_ix]:
                            count0 += 1
                x.sort()
                y[i] = x[count0 + (x_shape - count0 + 1)//2 - 1]
            y.sort()
            sm[k] = y[(y_shape + 1)//2 - 1] * 1.1926
        return sm
    sm = calc_sm(omega, x, y, sm)
    for i in range(omega.shape[0]):
        ps[i] = np.max(np.abs(omega[i, :])/sm)
    return ps