An algorithm for constructing orthogonal armlet multi-wavelets with multiplicity r and dilation factor a