N = 10**6
fac = [0]*N
inv = [0]*N
fac[0] = inv[0] = 1
mod = 10**9+7

for i in range(1,N):
    fac[i] = (fac[i-1]*i)%mod
    inv[i] = (inv[i-1]*pow(i,mod-2,mod))%mod

T = int(input())
for i in range(T):
    H, W = map(int, input().split())
    if H==1 or W==1:
        print(0)
    else:
        ncr = (fac[H-2+W-2]*inv[H-2]*inv[W-2])%mod
        print((2*(H+W-3)*ncr)%mod)