A theory of higher rank multiresolution analysis is given in the setting of abelian multiscalings. This theory enables the construction, from a higher rank MRA, of finite wavelet sets whose multidilations have translates forming an orthonormal basis in $L^2(\mathbb R^d)$. While tensor products of uniscaled MRAs provide simple examples we construct many nonseparable higher rank wavelets. In particular we construct \emph{Latin square wavelets} as rank~$2$ variants of Haar wavelets. Also we construct nonseparable scaling functions for rank $2$ variants of Meyer wavelet scaling functions, and we construct the associated nonseparable wavelets with compactly supported Fourier transforms. On the other hand we show that compactly supported scaling functions for biscaled MRAs are necessarily separable.

For a general vector field we exhibit two Hilbert spaces, namely the space of so called \emph{closed functions} and the space of \emph{exact functions} and we calculate the codimension of the space of exact functions inside the larger space of closed functions. In particular we provide a new approach for the known cases: the Glauber field and the second-order Ginzburg--Landau field and for the case of the fourth-order Ginzburg--Landau field.

Eventual positivity problems for real convergent Maclaurin series lead to density questions for sets of harmonic functions. These are solved for large classes of series, and in so doing, asymptotic estimates are obtained for the values of the series near the radius of convergence and for the coefficients of convolution powers.