import sys
input = sys.stdin.readline

fact = [1]
MOD = 10**9+7

for i in range(1, 4*10**5+10):
    fact.append(fact[-1]*i%MOD)

inv = lambda x: pow(x, MOD-2, MOD)

for _ in range(int(input())):
    H, W = map(int, input().split())
    
    if H==1 or W==1:
        print(0)
    else:
        print(2*(H-1)*fact[H+W-3]%MOD*inv(fact[H-1])%MOD*inv(fact[W-2])%MOD)