#include <bits/stdc++.h>
using namespace std;
const int INF = 0x3f3f3f3f;
const double eps = 1e-8;
const double PI = acos(-1.0);
const int MAXN = 1e3 + 5;
const int MAXM = 2e6 + 5;
const int MOD = 1e9 + 7;
int sgn(double x) {
  if (fabs(x) < eps) return 0;
  if (x < 0)
    return -1;
  else
    return 1;
}
long long pw(long long a, long long n, long long mod) {
  long long ret = 1;
  while (n) {
    if (n & 1) ret = ret * a % mod;
    a = a * a % mod;
    n >>= 1;
  }
  return ret;
}
long long gcd(long long a, long long b) { return b == 0 ? a : gcd(b, a % b); }
long long lcm(long long a, long long b) { return a / gcd(a, b) * b; }
void inc(long long &a, long long b) {
  a += b;
  if (a >= MOD) a -= MOD;
}
void dec(long long &a, long long b) {
  a -= b;
  if (a < 0) a += MOD;
}
int mul(int a, int b) {
  long long c = 1ll * a * b;
  return c - c / MOD * MOD;
}
int n, m;
char G[MAXN][MAXN];
bool flag;
bool jud(int d, int u, int l, int r) {
  if (d < 1 || u > n || l < 1 || r > m) return false;
  for (int i = d; i < u + 1; i++)
    for (int j = l; j < r + 1; j++) {
      if (i == d + 1 && j == l + 1) continue;
      if (G[i][j] != '#') return false;
    }
  return true;
}
bool ok(int x, int y) {
  if (jud(x, x + 2, y, y + 2)) return true;
  if (jud(x, x + 2, y - 1, y + 1)) return true;
  if (jud(x, x + 2, y - 2, y)) return true;
  if (jud(x - 1, x + 1, y, y + 2)) return true;
  if (jud(x - 1, x + 1, y - 2, y)) return true;
  if (jud(x - 2, x, y, y + 2)) return true;
  if (jud(x - 2, x, y - 1, y + 1)) return true;
  if (jud(x - 2, x, y - 2, y)) return true;
  return false;
}
int main() {
  scanf("%d%d", &n, &m);
  for (int i = 1; i < n + 1; i++) scanf("%s", G[i] + 1);
  flag = true;
  for (int i = 1; i < n + 1; i++)
    for (int j = 1; j < m + 1; j++) {
      if (G[i][j] == '#' && !ok(i, j)) {
        flag = false;
      }
    }
  if (flag)
    puts("YES");
  else
    puts("NO");
}
