#include <iostream>
#include <vector>
#include <cmath>
#include <algorithm>
#include <map>
#include <unordered_map>
#include <unordered_set>
#include <set>
#include <queue>
#include <deque>
#include <numeric>

#define ll              long long 
#define all(v)          v.begin(),v.end()
#define ld              long double
#define pll             std::pair<ll,ll>
#define pi              std::pair<int,int>
#define vi              std::vector<int>
#define vll             std::vector<ll>
#define len(x)          (int)x.size()
#define vec(T)          std::vector<T>

// overload << for pair
template <typename T>
std::ostream& operator<<(std::ostream& os, const std::pair<T, T>& pair) {
    std::cout << "(" << pair.first << ", " << pair.second << ")";
    return os;
}


// overload << for vec
template <typename T>
std::ostream& operator<<(std::ostream& os, const std::vector<T>& vec) {
    os << "[";
    for (size_t i = 0; i < vec.size(); ++i) {
        os << vec[i] << (i == vec.size() - 1 ? "\0" : ", ");
    }
    os << "]";
    return os;
}


// overload << for map and unordered_map Output
template<typename Map>
void printMap(const Map& m) {
    std::cout << '{';
    for (auto it = m.begin(); it != m.end(); it++) {
        std::cout << it -> first << ": " << it -> second << (std::next(it) == m.end() ? "\0" : ", ");
    }
    std::cout << "}";
}
template <typename K, typename V>
std::ostream& operator<<(std::ostream& os, const std::map<K, V>& m) {
    printMap(m);
    return os;
}
template <typename K, typename V>
std::ostream& operator<<(std::ostream& os, const std::unordered_map<K, V>& m) {
    printMap(m);
    return os;
}

// overload << for set and unordered_set Output
template<typename Set>
void printSet(const Set& s) {
    std::cout << '{';
    for (auto it = s.begin(); it != s.end(); it++) {
        std::cout << *it << (std::next(it) == s.end() ? "\0" : ", ");
    }
    std::cout << "}";
}
template<typename T>
std::ostream& operator<<(std::ostream& os, const std::set<T>& s) {
    printSet(s);
    return os;
}
template<typename T>
std::ostream& operator<<(std::ostream& os, const std::unordered_set<T>& s) {
    printSet(s);
    return os;
}

// General print function
template <typename T>
void print(const T& val) {
    std::cout << val << std::endl;
}

// Variadic template print function
template<typename T, typename... Args>
void print(const T& t, const Args&... args) {
    std::cout << t << ' ';
    print(args...);
    if (sizeof...(args) == 1) std::cout << '\n';

}

void solve(){  
    int n, m;
    std::cin >> n >> m;
    if (m % 2 == 0) {
        std::cout << "YES\n";
        for (int i = 0; i < n; i++) {
            for (int j = 0; j < m; j++) {
                std::cout << i * m + j + 1 << " ";
            }
            std::cout << "\n";
            
        }
        return;
    }
    std::cout << "YES\n";
    for (int i = 0; i < n; i++) {
        if (i % 2 == 1) {
            for (int j = m - 1; j >= 0; j--) {
                std::cout << i * m + j + 1 << " ";
            }
            std::cout << "\n";
        }
        else {
            std::cout << i * m + 1 << " ";
            for (int j = m - 1; j > 0; j--) {
                std::cout << i * m + j + 1 << " ";
            }
            std::cout << "\n";
        }
    }
}



    
int main(){
    std::ios::sync_with_stdio(false);
    std::cin.tie(nullptr);
    int t = 1;
    std::cin >> t;
    while (t--) solve();
}