code: def _nominal_kernel(x, y, out):
      
      for i in range(x.shape[0]):
          for j in range(y.shape[0]):
              out[i, j] += (x[i, :] == y[j, :]).sum()
  
      return out
docs: Number of features that match exactly