def METHOD_NAME(self, p):
    x = p[..., 0]
    y = p[..., 1]
    eips = 1
    shape = p.shape + (2, )
    val = np.zeros(shape, dtype=np.float64) 
    val[..., 0, 0] = eips*(1+x**2)
    val[..., 0, 1] = eips*x*y
    val[..., 1, 0] = val[..., 0, 1]
    val[..., 1, 1] = eips*(1+y**2)
    return val 