
Convergence analysis of a nested iterative scheme proposed by Bank,Welfert
and Yserentant (BWY) ([Numer. Math., 666: 645666, 1990]) for solving saddle
point system is presented. It is shown that this scheme converges under weaker
conditions: the contraction rate for solving the $(1,1)$ block matrix is bound
by $(\sqrt{5}1)/2$. Similar convergence result is also obtained for a class of
inexact Uzawa method with even weaker contraction bound $\sqrt{2}/2$.
Preconditioned generalized minimal residual method using BWY method as a
preconditioner is shown to converge with realistic assumptions.

Finite element exterior calculus (FEEC) has been developed as a systematical
framework for constructing and analyzing stable and accurate numerical method
for partial differential equations by employing differential complexes. This
paper is devoted to analyze convergence of adaptive mixed finite element
methods for Hodge Laplacian equations based on FEEC without considering
harmonic forms. More precisely, a residual type posteriori error estimates is
obtained by using the Hodge decomposition, the regular decomposition and
bounded commuting quasiinterpolants. An additional marking strategy is added
to ensure the quasiorthogonality. Using this quasiorthogonality, the uniform
convergence of adaptive mixed finite element methods is obtained without
assuming the initial mesh size is small enough.

Due to the indefiniteness and poor spectral properties, the discretized
linear algebraic system of the vector Laplacian by mixed finite element methods
is hard to solve. A block diagonal preconditioner has been developed and shown
to be an effective preconditioner by Arnold, Falk, and Winther [Acta Numerica,
15:1155, 2006]. The purpose of this paper is to propose alternative and
effective block diagonal and block triangular preconditioners for solving this
saddle point system. A variable Vcycle multigrid method with the standard
pointwise GaussSeidel smoother is proved to be a good preconditioner for a
discrete vector Laplacian operator. This multigrid solver will be further used
to build preconditioners for the saddle point systems of the vector Laplacian
and the Maxwell equations with divergent free constraint. The major benefit of
our approach is that the pointwise GaussSeidel smoother is more algebraic and
can be easily implemented as a blackbox smoother.

Superconvergence and a posteriori error estimators of recovery type are
analyzed for the 4node hybrid stress quadrilateral finite element method
proposed by Pian and Sumihara (Int. J. Numer. Meth. Engrg., 1984, 20:
16851695) for linear elasticity problems. Uniform superconvergence of order
$O(h^{1+\min\{\alpha,1\}})$ with respect to the Lam\'{e} constant $\lambda$ is
established for both the recovered gradients of the displacement vector and the
stress tensor under a mesh assumption, where $\alpha>0$ is a parameter
characterizing the distortion of meshes from parallelograms to quadrilaterals.
A posteriori error estimators based on the recovered quantities are shown to be
asymptotically exact. Numerical experiments confirm the theoretical results.