#include <bits/stdc++.h>
using namespace std;
#define all(v) (v).begin(),(v).end()
#define pb(a) push_back(a)
#define rep(i, n) for(int i=0;i<n;i++)
#define foa(e, v) for(auto&& e : v)
using ll = long long;
const ll MOD7 = 1000000007, MOD998 = 998244353, INF = (1LL << 60);
#define dout(a) cout<<fixed<<setprecision(10)<<a<<endl;

bool check(vector<vector<int>> &v) {
    // int n = v.size();
    // int m = v[0].size();
    // set<int> st;
    // rep(i, n) {
    //     rep(j, m - 1) {
    //         if(st.count(v[i][j] + v[i][j + 1])) return false;
    //         st.insert(v[i][j] + v[i][j + 1]);
    //     }
    // }
    // rep(i, n - 1) rep(j, m) {
    //     if(st.count(v[i][j] + v[i + 1][j])) return false;
    //     st.insert(v[i][j] + v[i + 1][j]);
    // } 
    return true;
}

void solve() {
    int n, m; cin >> n >> m;
    vector<vector<int>> v(n, vector(m, 0));
    if(m & 1) {
        rep(i, n) rep(j, m) {
            v[i][(i + j) % m] = i * m + j + 1;
        }
    } else {
        rep(i, n) rep(j, m) v[i][j] = i * m + j + 1;
    }
    if(check(v)) {
        cout << "YES" << endl;
        rep(i, n) {
            rep(j, m) cout << v[i][j] << " ";
            cout << endl;
        }
    } else {
        cout << "NO" << endl;
    }
}
int main() {
    cin.tie(0);
    ios::sync_with_stdio(false);
    int t;
    cin >> t;
    while(t --) {
        solve();
    }
    return 0;
}