The article proposes a nonlinear optimal control method for the dynamic model of the parallel double inverted pendulum. Three different forms of this system are considered: (i) two poles mounted on the same cart resulting into a model with four state variables and one control input, (ii) two poles mounted on the same cart resulting into a model with six state variables and one control input, (iii) two different cart-pole systems connected through an elastic link, resulting into a model with eight state variables and two control inputs. To implement the proposed nonlinear optimal control method the dynamic model of the parallel double pendulum undergoes approximate linearization around a temporary operating point which is updated at each sampling instance. This operating point is defined within each sampling period by the present value of the state vector of the parallel double inverted pendulum and by the last sampled value of the control inputs vector. The linearization process is based on first-order Taylor series expansion and on the computation of the associated Jacobian matrices. For the approximately linearized model of the parallel double inverted pendulum a stabilizing H-infinity feedback controller is designed. To compute the feedback gains of this controller an algebraic Riccati equation is also being solved at each time-step of the control algorithm. The global stability properties of the control scheme are proven through Lyapunov analysis. First it is demonstrated that the H-infinity tracking performance criterion holds, which signifies robustness of the control method under model uncertainty and external perturbations. Moreover, under moderate conditions it is proven that the control loop of the parallel double inverted pendulum is globally asymptotically stable. To implement state estimation-based control the H-infinity Kalman Filter is used as a robust observer. The nonlinear optimal control scheme for the parallel double inverted pendulum achieves fast and accurate tracking of reference setpoints under moderate variations of the control inputs.