#include<iostream>
using namespace std;
const int o=1e6+10;
int T,n,m,a[o],l[o],cnt;long long s1,s2;bool flg;
inline int hsh(int x,int y){return (x-1)*m+y;}
int main(){
	ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
	for(cin>>T;T--;cnt=flg=0){
		cin>>n>>m;
		if(n*(n+1ll)*m*(m+1)%8){cout<<"No\n";continue;}
		if(m*(m+1)%4) swap(n,m),flg=1;
		s1=m;s2=m*(m+1ll)/4;
		for(int i=m;i;--i) for(;i<=s1&&i*(i+1ll)/2+s1-i<=s2;s1-=i,s2-=i*(i+1ll)/2) l[++cnt]=i;
		for(int i=1,t=0;i<=cnt;++i) for(int j=1;j<=l[i];++j) a[++t]=(i&1);
		for(int i=2;i<=n;++i) for(int j=1;j<=m;++j) a[hsh(i,j)]=a[j];
		cout<<"Yes\n";
		if(flg) for(int i=1;i<=m;++i,cout<<"\n") for(int j=1;j<=n;++j) cout<<a[hsh(j,i)]<<" ";
		else for(int i=1;i<=n;++i,cout<<"\n") for(int j=1;j<=m;++j) cout<<a[hsh(i,j)]<<" ";
	}
	return 0;
}