jac_large_angle[:, 0, 0] = A
jac_large_angle[:, 0, 1] = -B
jac_large_angle[:, 1, 0] = B
jac_large_angle[:, 1, 1] = A
jac[large_angle_inds] = jac_large_angle
return jac.squeeze_()
After Change
jac[large_angle_inds])
B = ((1. - c) / angle).unsqueeze_(dim=1).unsqueeze_(
dim=2).expand_as(jac[large_angle_inds]) * \
cls.wedge(phi.__class__([1.]))
jac[large_angle_inds] = A + B