A robust and efficient harmonic balance (HB) using direct solution of HB Jacobian