An interior point algorithm for the p*(k) matrix horizontal LCP