#include <bits/stdc++.h>
#pragma GCC optimize("Ofast")
#pragma GCC target("avx,avx2,fma")
using namespace std;
const long long int N = 200043;
const long long int MOD = 1000000007;
long long int add(long long int x, long long int y) {
  x = ((x % MOD) + (y % MOD)) % MOD;
  while (x >= MOD) x -= MOD;
  while (x < 0) x += MOD;
  return x;
}
long long int mul(long long int x, long long int y) {
  return ((x % MOD) * (y % MOD)) % MOD;
}
long long int binpow(long long int x, long long int y) {
  long long int z = 1;
  while (y) {
    if (y & 1) z = mul(z, x);
    x = mul(x, x);
    y >>= 1;
  }
  return z;
}
long long int inv(long long int x) { return binpow(x, MOD - 2); }
long long int divide(long long int x, long long int y) {
  return mul(x, inv(y));
}
long long int n, m;
long long int ar[1005][1005];
long long int dp1[1005][1005];
long long int dp2[1005][1005];
long long int dp3[1005][1005];
long long int dp4[1005][1005];
int main() {
  ios_base::sync_with_stdio(false);
  cin.tie(NULL);
  long long int t, i, j;
  t = 1;
  while (t--) {
    cin >> n >> m;
    for (i = 1; i <= n; i++) {
      for (j = 1; j <= m; j++) cin >> ar[i][j];
    }
    for (i = 1; i <= n; i++) {
      for (j = 1; j <= m; j++)
        dp1[i][j] = ar[i][j] + max(dp1[i - 1][j], dp1[i][j - 1]);
    }
    for (i = n; i >= 1; i--) {
      for (j = m; j >= 1; j--) {
        dp2[i][j] = ar[i][j] + max(dp2[i][j + 1], dp2[i + 1][j]);
      }
    }
    for (i = n; i >= 1; i--) {
      for (j = 1; j <= m; j++)
        dp3[i][j] = ar[i][j] + max(dp3[i][j - 1], dp3[i + 1][j]);
    }
    for (i = 1; i <= n; i++) {
      for (j = m; j >= 1; j--)
        dp4[i][j] = ar[i][j] + max(dp4[i - 1][j], dp4[i][j + 1]);
    }
    long long int ans = 0;
    for (i = 2; i < n; i++) {
      for (j = 2; j < m; j++) {
        ans = max(dp1[i - 1][j] + dp2[i + 1][j] + dp3[i][j - 1] + dp4[i][j + 1],
                  ans);
        ans = max(dp1[i][j - 1] + dp2[i][j + 1] + dp3[i + 1][j] + dp4[i - 1][j],
                  ans);
      }
    }
    cout << ans << endl;
  }
}
