#include<bits/stdc++.h>
#define ll long long
#define ld long double
using namespace std;

const int N=1e6+10,inf=1e9;
vector<ll>e;
int a[N],b[N];
int n,m;

inline void solve(){
    
    cin >> n >> m;
    
    ll s1=1ll*(n+1)*n/2;
    ll s2=1ll*(m+1)*m/2;
    if ((s1&1)&&(s2&1)){
        cout << "No\n";
        return;
    }
    int bj=0;
    if (s2&1){
        swap(s1,s2);
        bj=1;
    }
    s2/=2;
    ll x=1;
    int r=0,o=0,mm;
    if (!bj) mm=m; else mm=n;
    
    s2-=mm;
    while(s2){
        
        int p;
        p = std::upper_bound(e.begin(), e.end(), s2) - e.begin();

        b[++r]=p;

        s2-=e[p-1];
    }
    // if (r!=mm){
    //     cout << "No\n";
    //     return;
    // }

    int t=0;
    for(int i=1;i<=r;++i){
        for(int j=t+1;j<=t+b[i]+1;++j)
            a[j]=o;
        o^=1; t+=b[i]+1;
    }
    for(int i=t+1;i<=mm;++i){
        a[i]=o;
        o^=1;
    }

    cout << "Yes\n";
    if (bj){
        for(int i=1;i<=n;++i){
            for(int j=1;j<=m;++j)
                cout << a[i] << ' ';
            cout << '\n';
        }
    }
    else{
        for(int i=1;i<=n;++i){
            for(int j=1;j<=m;++j)
                cout << a[j] << ' ';
            cout << '\n';
        }
    }
}
int main(){
    ios::sync_with_stdio(false);
    cin.tie(0),cout.tie(0);

    ll x=1;
    while(x<=N-10){
        e.push_back(x*(x+1)/2);
        x++;
    }

    int T;
    cin >> T;
    while(T--){
        solve();
    }
    return 0;
}