Implementation of Hellan-Hermann_Johnson mixed method with triangular decomposition.