some minor performance improvement and code reduction

- remove the the PD check of Omega_j because Omega_j seems to always stay PD
- vetorize the recursive formula for a (Bolduc eq 15)