#include <bits/stdc++.h>
using namespace std;
int n, m, v[1001][1001];
char a[1001][1001];
bool check(int x, int y) {
  if (x + 2 < n && y + 2 < m) {
    if (a[x + 1][y] == '#' && a[x + 2][y] == '#' && a[x][y + 1] == '#' &&
        a[x][y + 2] == '#' && a[x + 1][y + 2] == '#' &&
        a[x + 2][y + 1] == '#' && a[x + 2][y + 2] == '#') {
      v[x + 1][y] = v[x + 2][y] = v[x][y + 1] = v[x][y + 2] = v[x + 1][y + 2] =
          v[x + 2][y + 1] = v[x + 2][y + 2] = 1;
      return 1;
    }
  }
  if (0 <= x - 1 && x + 1 < n && y + 2 < m) {
    if (a[x - 1][y] == '#' && a[x + 1][y] == '#' && a[x - 1][y + 1] == '#' &&
        a[x + 1][y + 1] == '#' && a[x - 1][y + 2] == '#' &&
        a[x][y + 2] == '#' && a[x + 1][y + 2] == '#') {
      v[x + 1][y] = v[x + 2][y] = v[x][y + 1] = v[x][y + 2] = v[x + 1][y + 2] =
          v[x + 2][y + 1] = v[x + 2][y + 2] = 1;
      return 1;
    }
  }
  if (0 <= x - 2 && y + 2 < m) {
    if (a[x - 1][y] == '#' && a[x - 2][y] == '#' && a[x][y + 1] == '#' &&
        a[x][y + 2] == '#' && a[x - 2][y + 1] == '#' &&
        a[x - 2][y + 2] == '#' && a[x - 1][y + 2] == '#') {
      v[x - 1][y] = v[x - 2][y] = v[x][y + 1] = v[x][y + 2] = v[x - 2][y + 2] =
          v[x - 2][y + 2] = v[x - 1][y + 2] = 1;
      return 1;
    }
  }
  if (x + 2 < n && 0 <= y - 1 && y + 1 < m) {
    if (a[x][y - 1] == '#' && a[x][y + 1] == '#' && a[x + 1][y - 1] == '#' &&
        a[x + 1][y + 1] == '#' && a[x + 2][y - 1] == '#' &&
        a[x + 2][y] == '#' && a[x + 2][y + 1] == '#') {
      v[x][y - 1] = v[x][y + 1] = v[x + 1][y - 1] = v[x + 1][y + 1] =
          v[x + 2][y - 1] = v[x + 2][y] = v[x + 2][y + 1] = 1;
      return 1;
    }
  }
  if (0 <= x - 2 && 0 <= y - 1 && y + 1 < m) {
    if (a[x][y - 1] == '#' && a[x][y - 1] == '#' && a[x - 1][y - 1] == '#' &&
        a[x - 1][y + 1] == '#' && a[x - 2][y - 1] == '#' &&
        a[x - 2][y] == '#' && a[x - 2][y + 1] == '#') {
      v[x][y - 1] = v[x][y - 1] = v[x - 1][y - 1] = v[x - 1][y + 1] =
          v[x - 2][y - 1] = v[x - 2][y] = v[x - 2][y + 1] = 1;
      return 1;
    }
  }
  if (x + 2 < n && 0 <= y - 2) {
    if (a[x][y - 1] == '#' && a[x][y - 2] == '#' && a[x + 1][y] == '#' &&
        a[x + 1][y - 2] == '#' && a[x + 2][y] == '#' &&
        a[x + 2][y - 1] == '#' && a[x + 2][y - 2] == '#') {
      v[x][y - 1] = v[x][y - 1] = v[x + 1][y] = v[x + 1][y - 2] = v[x + 2][y] =
          v[x + 2][y - 1] = v[x + 2][y - 2] = 1;
      return 1;
    }
  }
  if (0 <= x - 1 && x + 1 < n && 0 <= y - 2) {
    if (a[x - 1][y] == '#' && a[x - 1][y - 1] == '#' &&
        a[x - 1][y - 2] == '#' && a[x][y - 2] == '#' && a[x + 1][y] == '#' &&
        a[x + 1][y - 1] == '#' && a[x + 1][y - 2] == '#') {
      v[x - 1][y] = v[x - 1][y - 1] = v[x - 1][y - 2] = v[x][y - 2] =
          v[x + 1][y] = v[x + 1][y - 1] = v[x + 1][y - 2] = 1;
      return 1;
    }
  }
  if (0 <= x - 2 && 0 <= y - 2) {
    if (a[x][y - 1] == '#' && a[x][y - 2] == '#' && a[x - 1][y] == '#' &&
        a[x - 1][y - 2] == '#' && a[x - 2][y] == '#' &&
        a[x - 2][y - 1] == '#' && a[x - 2][y - 1] == '#') {
      v[x][y - 1] = v[x][y - 1] = v[x - 1][y] = v[x - 1][y - 2] = v[x - 2][y] =
          v[x - 2][y - 1] = v[x - 2][y - 1] = 1;
      return 1;
    }
  }
  return 0;
}
int main() {
  ios_base::sync_with_stdio(false);
  cin.tie(NULL);
  int t, i, j;
  memset(v, 0, sizeof(v));
  cin >> n >> m;
  for (i = 0; i < n; i++) cin >> a[i];
  for (i = 0; i < n; i++) {
    for (j = 0; j < m; j++) {
      if (a[i][j] == '#' && !v[i][j] && !check(i, j)) {
        cout << "NO\n";
        return 0;
      }
    }
  }
  cout << "YES\n";
  return 0;
}
