#include<iostream>
#include<algorithm>
#include<vector>
#include<map>
#include<random>
#include<string>
#include<time.h>
#include<cstring>
#include<set>
#include<cmath>
#include<queue>
#define int long long
#define endl '\n'
#define F first
#define S second
#define sz size()
#define B begin()
#define pb push_back
#define D end()
#define pii pair<int,int> 
#define no cout<<"No"<<'\n'
#define yes cout<<"Yes"<<'\n'
#define ln cout<<'\n'
using namespace std;
typedef long long ll;
const int N=2e6+10;
const int P=1331;
const int mod=1e9+7;
vector<int> s;
void solve()
{
int n,m;cin>>n>>m;
	int cnt=n*(n+1)/2*(m+1)*m/2;
	if(cnt&1)
	{
		no;
		return ;
	}
	yes;
	if((m*(m+1)/2)%2==0)
	{
		int is=0;
		int res=(m+1)*m/2;
		res/=2;
		res-=m;
		string a;
		while(res)
		{
			auto it=upper_bound(s.B+1,s.D,res)-s.B;
			res-=s[it-1];
			for(int i=1;i<=it-1;i++) a.pb('0'+is);
			is++;is%=2;
		}
		for(int i=a.sz;i<m;i++)
		{
			a.pb('0'+is);
			is++;is%=2;
		}
		vector<int>ans[n+1];
		for(int i=1;i<=n;i++)
		{
			for(auto x:a)
			{
				ans[i].pb(x-'0');
			}
		}
		for(int i=1;i<=n;i++)
		{
			for(int j=0;j<m;j++)
			{
				cout<<ans[i][j]<<' ';
			}ln;
		}
	}else
	{
//	cout<<"afgaga"<<'\n';
		int is=0;
		int res=(n+1)*n/2;
		res/=2;
		res-=n;
		//cout<<res<<'\n';
		string a;
		while(res)
		{
			auto it=upper_bound(s.B+1,s.D,res)-s.B;
			res-=s[it-1];
			for(int i=1;i<=it-1;i++) a.pb('0'+is);
			is++;is%=2;
		}
		for(int i=a.sz;i<n;i++)
		{
			a.pb('0'+is);
			is++;is%=2;
		}
	//	cout<<a<<'\n';
		vector<vector<int>>ans(n+1,vector<int>(m+1));
		for(int i=1;i<=m;i++)
		{
			for(int j=1;j<=n;j++)
			{
				ans[j][i]=a[j-1]-'0';
			}
		}
		for(int i=1;i<=n;i++)
		{
			for(int j=1;j<=m;j++)
			{
				cout<<ans[i][j]<<' ';
			}ln;
		}
	}
	
}
signed main()
{
	ios::sync_with_stdio(false);
	cin.tie(0);
	int tc=1;
	cin>>tc;
	s.pb(0);s.pb(0);
	for(int i=2;i<=1e6;i++)
	{
		int x=s[i-1]+i-1;
		s.pb(x);
		if(s[i]>500000000000000) break;
	}
	//for(int i=1;i<s.sz;i++) cout<<s[i]<<'\n';
	while(tc--)
	{
		solve();
	}
	return 0;
}