#include<bits/stdc++.h>
#define int long long
#define MAXN 1000000
using namespace std;
int T;
int N,M;
int ans[MAXN+5],tp=0;
int st[MAXN+5],newp;
int B[MAXN+5];
int flag;
void DFS(int n,int num){
    if(flag)return;
    if(num==0 && n==0){
        newp=0;
        for(int i=1;i<=tp;i++){
            st[++newp]=ans[i];
        }
        flag=1;
        return;
    }
    if(n*2>num)return;
    if(n*2==num){
        for(int i=1;i*2<=num;i++){
            ans[++tp]=1;
        }
        DFS(0,0);
    }
    int k=upper_bound(B+1,B+1+MAXN,num)-B;
    for(int i=k-1;i>=1;i--){
        if(n*(i+1)<num)break;
        ans[++tp]=i;
        DFS(n-i,num-(i*(i+1)));
        tp--;
    }
}
signed main()
{
    for(int i=1;i<=MAXN;i++){
        B[i]=i*(i+1);
    }
    cin>>T;
    while(T--){
        cin>>N>>M;
        int c=0;
        if(M%4==0 || (M+1)%4==0){
            swap(N,M);
            c=1;
        }
        if(N%4!=0 &&(N+1)%4!=0){
            cout<<"No\n";
            continue;
        }
        else{
            cout<<"Yes\n";
        }
        vector<vector<int> > a(N+1,vector<int>(M+1));
        tp=0;
        flag=0;
        DFS(N,N*(N+1)/2);
        /*for(int i=1;i<=newp;i++){
            cout<<st[i]<<" ";
        }
        cout<<'\n';*/

        int op=0,k=1,tmp=1;
        while(tmp<=N){
            for(int i=1;i<=st[k];i++){
                for(int j=1;j<=M;j++){
                    a[tmp+i-1][j]=op;
                }
            }
            k++;
            tmp+=st[k];
            op^=1;
        }
        if(!c){
            for(int i=1;i<=N;i++){
                for(int j=1;j<=M;j++){
                    cout<<a[i][j]<<" ";
                }
                cout<<'\n';
            }
        }
        else{
            for(int j=1;j<=M;j++){
                for(int i=1;i<=N;i++){
                    cout<<a[i][j]<<" ";
                }
                cout<<'\n';
            }
        }
    }

    return 0;
}