Nowadays, the problem of classification of integrable nonlinear partial differential equations and their discrete analogues in 1+1 dimensions is well–studied. Within the framework of the symmetry approach, there was obtained a complete description of integrable representatives of a number of classes of equations that are interesting from the point of view of application, see [17], [34], [26], [2]. The problem of exhaustive classification of integrable equations containing a large number of independent variables remains less studied due to its extreme complexity. The symmetry approach, which has proven to be the most effective tool for classifying equations of dimension 1+1, is not quite suitable for integrable classification of multidimensional equations. As it is noted in [27], in this problem the symmetry approach loses its efficiency due to problems with nonlocalities involved in higher symmetries