bilinear kernel is wrong?

[jszgz]

        g_lt = (1 + (q_lt[..., :N].type_as(p) - p[..., :N])) * (1 + (q_lt[..., N:].type_as(p) - p[..., N:]))
        g_rb = (1 - (q_rb[..., :N].type_as(p) - p[..., :N])) * (1 - (q_rb[..., N:].type_as(p) - p[..., N:]))
        g_lb = (1 + (q_lb[..., :N].type_as(p) - p[..., :N])) * (1 - (q_lb[..., N:].type_as(p) - p[..., N:]))
        g_rt = (1 - (q_rt[..., :N].type_as(p) - p[..., :N])) * (1 + (q_rt[..., N:].type_as(p) - p[..., N:]))

This computation is different from standard bilinear algorithm, is it wrong or another algorithm?

[jszgz]

[jszgz]

For example, when you compute g_lt , 1 + (q_lt[..., :N].type_as(p) - p[..., :N]) is right, but 1 + (q_lt[..., N:].type_as(p) - p[..., N:]) should be 1 + ( p[..., N:] - q_lt[..., N:].type_as(p)) , which is P12(X0,y1) in my figure

[xiab3369]

It's right