Fast Computation of Posterior Mode in Multi-Level Hierarchical Models