#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int MAXN=303;
const int MAXV=1e6;
const ll MAX=0x3f3f3f3f3f3f3f3f;
const int MAX_I=0x3f3f3f3f;
int n,m,b[MAXN][MAXN],g[MAXN<<1][MAXN<<1],cnt[MAXN<<1];
ll dis[MAXN<<1];
queue<int> q;
bitset<MAXN<<1> vis;
bool spfa()
{
    fill_n(dis,n+m+1,MAX);
    dis[0]=cnt[0]=0;
    vis.reset(),vis.set(0);
    while(!q.empty()) q.pop();
    q.push(0);
    while(!q.empty())
    {
        int u=q.front();
        q.pop(),vis.reset(u);
        for(int v=0;v<=n+m;v++)
            if(dis[v]>dis[u]+g[u][v])
            {
                dis[v]=dis[u]+g[u][v],cnt[v]=cnt[u]+1;
                if(cnt[v]>=n+m+1) return 0;
                if(!vis[v]) q.push(v),vis.set(v);
            }
    }
    return 1;
}
int main()
{
    int T;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d%d",&n,&m);
        for(int i=1;i<n;i++)
            for(int j=1;j<m;j++)
                scanf("%d",&b[i][j]);
        for(int i=1;i<=n+1;i++) b[i][m]=b[i][m+1]=0;
        for(int i=1;i<=m+1;i++) b[n][i]=b[n+1][i]=0;
        for(int i=n;i>=1;i--)
            for(int j=m;j>=1;j--)
                b[i][j]-=b[i+1][j+1]+b[i+1][j]+b[i][j+1];
        // memset(g,0x3f,sizeof(g));
        for(int i=0;i<=n+m;i++)
            fill_n(g[i],n+m+1,MAX_I);
        for(int i=1;i<=n;i++)
            for(int j=1;j<=m;j++)
            {
                if((i+j)&1) g[j+n][i]=MAXV-b[i][j],g[i][j+n]=b[i][j];
                else g[i][j+n]=MAXV-b[i][j],g[j+n][i]=b[i][j];
            }
        for(int i=1;i<=n+m;i++)
            g[0][i]=0;
        if(!spfa())
        {
            printf("NO\n");
            continue;
        }
        printf("YES\n");
        for(int i=1;i<=n;i++)
            for(int j=1;j<=m;j++)
            {
                if((i+j)&1) b[i][j]+=dis[i]-dis[j+n];
                else b[i][j]+=dis[j+n]-dis[i];
                printf("%d%c",b[i][j]," \n"[j==m]);
            }
    }
    return 0;
}