#include<bits/stdc++.h>
using namespace std;
#define int long long
#define endl '\n'
const int N=1e6+10,mod=998244353;

int gcd(int a,int b){return b?gcd(b,a%b):a;}
typedef pair<int,int> PII;

vector<int> num;
int path[N];

bool dfs(int u,int n1,int n2,int op){
	if(n2==u*(u+1)/4&&n1==u){
		return true;
	}
	
	if(u*(u+1)/4-n2==u-n1){
		for(int p=n1+1;p<=u;p++) path[p]=op,op=1-op;
		return true;
	}
	
	int sheng=u*(u+1)/4-n2;
	auto it=upper_bound(num.begin(),num.end(),sheng)-num.begin();
	it--;
	
	for(int i=it;i>=1;i--){
		int st=1;
		for(int j=min((u-n1)/i,(u*(u+1)/4-n2)/num[i]);j>=1;j--){
			if(i*j+n1>u) break;
			if(num[i]*j+n2>(u+1)*u/4) break;
			if(dfs(u,n1+j*i,n2+num[i]*j,(op+j)%2)){ 
				int nums=0;
				for(int p=n1+1;p<=n1+i*j;p++){
					 path[p]=op;
					 nums++;
					 if(nums%i==0) op=1-op;
				}
				return true;
			}
		} 
	}
	return false;
}

void slove(){
	int n,m;cin>>n>>m;
	
	vector<vector<int>>g(n+4,vector<int>(m+5,0));
	if(m%4==3||m%4==0){
		dfs(m,0,0,0);
		for(int i=1;i<=n;i++){
			for(int j=1;j<=m;j++) g[i][j]=path[j];
		}
	}
	else if(n%4==0|n%4==3){
		dfs(n,0,0,0);
		for(int i=1;i<=m;i++){
			for(int j=1;j<=n;j++) g[j][i]=path[j];
		}
	}
	else {
		cout<<"No"<<endl;
		return ;
	}
	cout<<"Yes"<<endl;
	for(int i=1;i<=n;i++){	
		for(int j=1;j<=m;j++){
		 	cout<<g[i][j]<<" ";
		}
		cout<<endl;
	}
}
signed main(){
	for(int i=0;i<=1e4;i++) num.push_back(i*(i+1)/2);
	ios::sync_with_stdio(false);cin.tie(nullptr);cout.tie(nullptr);
	int T=1;
	cin>>T;
	while(T--) slove();
}