#include<bits/stdc++.h>
using namespace std;
#define il inline
#define fi first
#define se second
#define mk make_pair
#define eb emplace_back
#define rep(i,l,r) for(int i=(l); i<=(r); ++i)
#define rep_(i,l,r) for(int i=(l); i>=(r); --i)
typedef long long lr;
typedef double db;
typedef pair<int,int> pii;
typedef vector<int> vi;
constexpr int mod1=998244353,mod2=1e9+7;
constexpr db pi=3.141592653589793,eps=1e-9;
constexpr int inf32=0x3f3f3f3f,Inf32=0xc0c0c0c0;
constexpr lr inf64=0x3f3f3f3f3f3f3f3f,Inf64=0xc0c0c0c0c0c0c0c0;
template<typename T>il T Max(T x,T y) { return (x>y)? x:y; }
template<typename T>il T Min(T x,T y) { return (x<y)? x:y; }
template<typename T>il T gcd(T x,T y) { return (!y)? x:gcd(y,x%y); }
template<typename T>il T Abs(T x) { return (x>0)? x:(-x); }
template<typename T>il T Rnd(T l,T r,mt19937_64 &eng)
{
	uniform_int_distribution<T> uid(l,r);
	return uid(eng);
}
mt19937_64 eng(chrono::high_resolution_clock::now().time_since_epoch().count());
constexpr int N=200200;
int n,m;
il void Solve()
{
	cin>>n>>m;
	if((n*m)&3)
	{
		cout<<"NO\n";
		return;
	}
	if(!(n&3))
	{
		cout<<"YES\n";
		int cnt=0;
		rep(i,1,n)
		{
			rep(j,1,m)
				cout<<cnt+j<<' ';
			cout<<'\n';
			if(!(i&3))
				cnt+=m;
		}
		return;
	}
	if(!(m&3))
	{
		cout<<"YES\n";
		int cnt=1;
		rep(i,1,n)
		{
			rep(j,1,m)
			{
				cout<<cnt<<' ';
				if(!(j&3))
					++cnt;
			}
			cout<<'\n';
		}
		return;
	}
	if(n==2&&m==2)
	{
		cout<<"NO\n";
		return;
	}
	if(m==2)
	{
		cout<<"YES\n1 1\n1 2\n1 2\n3 2\n3 2\n3 3\n";
		int cnt=4;
		rep(i,7,n)
		{
			rep(j,1,m)
				cout<<cnt+j-1<<' ';
			cout<<'\n';
			if(!((i-6)&3))
				cnt+=m;
		}
		return;
	}
	cout<<"YES\n";
	int cnt=4;
	cout<<"1 1 1 3 3 3 ";
	rep(i,7,m)
	{
		cout<<cnt<<' ';
		if(!((i-6)&3))
			++cnt;
	}
	cout<<'\n';
	cout<<"1 2 2 2 2 3 ";
	rep(i,7,m)
	{
		cout<<cnt<<' ';
		if(!((i-6)&3))
			++cnt;
	}
	cout<<'\n';
	rep(i,3,n)
	{
		rep(j,1,m)
			cout<<cnt+j-1<<' ';
		cout<<'\n';
		if(!((i-2)&3))
			cnt+=m;
	}
	return;
}
int main()
{
#ifdef LOCAL
	string fpre="test",isuf="in",osuf="out";
	assert(freopen((fpre+"."+isuf).c_str(),"r",stdin));
	assert(freopen((fpre+"."+osuf).c_str(),"w",stdout));
#endif
	ios::sync_with_stdio(0),cin.tie(0),cout.tie(0);
	int T;
	cin>>T;
	while(T--)
		Solve();
	return 0;
}