Numerical hartree-fock methods for diatomic molecules: A partial-wave expansion approach