Cartan developed a structure theory of such infinite pseudogroups. However this structure theory is difficult to apply to physical PDEs because it requires the analysis of associated overdetermined (pseudogroup) determining systems which are often highly nonlinear and can be very large (e.g. the pseudogroup defining system for the KP equation occupies 30 megabytes when generated automatically by computer). In comparison Lie's method yields infinitesimal determining systems which are linear and generally much smaller.
My first effort  towards a practical algorithm for calculating Cartan structure of infinite Lie pseudogroups of PDEs was a hybrid Lie-Cartan algorithm. To apply the method the infinitesimal determining system is reduced to standard form. This standard form is used to predict a model of the pseudogroup defining system. By requiring the model system to admit the infinitesimal symmetries we obtain a system of PDEs for the unspecified functions in the model. Cartan's algorithm is then run on the model system subject to this differential characterization. Our success in calculating Cartan structure in examples using this method inspired Lisle to conjecture an infinitesimal form of the method. This infinitesimal method was rigorously justified [45,44] for Lie pseudogroups of transitive type using the theory of Singer and Sternberg .
In summary we obtained an effective algorithm which from the infinitesimal determining system can determine whether a pseudogroup of symmetries is isomorphic to a transitive pseudogroup. If so, then the algorithm can determine the Cartan structure of the pseudogroup. We applied a preliminary implementation of this algorithm to physical PDEs with infinite pseudogroups, including the KP and Liouville equations.