Fast numerical algorithms for fitting multiresolution hybrid shape models to brain MRI