On the hp-adaptive solution of complete electrode model forward problems of electrical impedance tomography