#include <iostream>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <queue>
#include <vector>
using namespace std;
const int N=610;
typedef long long int ll;
typedef pair <int,int> PII;
int T,n,m,b[N][N],a[N][N],cnt[N];
bool vis[N];
ll dist[N];
queue <int> q;
vector <PII> g[N];
inline int read()
{
    int x=0,f=1;char c=getchar();
    while(c<'0' ||c>'9') {if(c=='-') f=-1;c=getchar();}
    while(c>='0' && c<='9') {x=x*10+c-48,c=getchar();}
    return x*f;
}
void addedge(int u,int v,int w)
{
    g[u].push_back({v,w});
}
bool spfa()
{
    while(q.empty()==false) q.pop();
    for(int i=1;i<=n+m;i++) dist[i]=cnt[i]=0,vis[i]=true,q.push(i);
    while(q.empty()==false)
    {
        int u=q.front();
        q.pop();
        vis[u]=false;
        for(int i=0;i<g[u].size();i++)
        {
            PII p=g[u][i];
            int v=p.first;
            if(dist[v]>dist[u]+p.second)
            {
                dist[v]=dist[u]+p.second;
                if(vis[v]==false)
                {
                    q.push(v),vis[v]=true,cnt[v]++;
                    if(cnt[v]>n+m+1) return false;
                }
            }
        }
    }
    return true;
}
int main()
{
    T=read();
    while(T--)
    {
        n=read(),m=read();
        for(int i=1;i<=n+m;i++) g[i].clear();
        memset(a,0,sizeof(a));
        for(int i=1;i<=n-1;i++)
            for(int j=1;j<=m-1;j++) b[i][j]=read();
        for(int i=n-1;i>=1;i--)
            for(int j=m-1;j>=1;j--)
                a[i][j]=b[i][j]-a[i][j+1]-a[i+1][j]-a[i+1][j+1];
        for(int i=1;i<=n;i++)
        {
            for(int j=1;j<=m;j++)
            {
                int x=-a[i][j],y=1000000-a[i][j];
                if((i+j)&1) addedge(j+n,i,-x),addedge(i,j+n,y);
                else addedge(j+n,i,y),addedge(i,j+n,-x);
            }
        }
        if(spfa()==false)
        {
            puts("NO");
            continue;
        }
        puts("YES");
        for(int i=1;i<=n;i++)
        {
            for(int j=1;j<=m;j++)
            {
                if((i+j)&1) a[i][j]+=dist[j+n]-dist[i];
                else a[i][j]+=dist[i]-dist[j+n];
            }
        }
        for(int i=1;i<=n;i++)
        {
            for(int j=1;j<=m;j++) printf("%d ",a[i][j]);
            puts("");
        }
    }
    return 0;
}