A feasibility study of optimal control techniques for cardiac defibrillation on anatomical three spatial dimensional rabbit ventricle geometry in the presence of bilateral control constraints is presented. The work addresses the numerical treatment of multi-scale and multi-domain simulations of the bidomain equations and is based on the primal-dual active set method to solve the optimality system for this large scale optimization problem. Numerical results are presented for a successful defibrillation study. Robustness of the optimization algorithm w.r.t to variations in the model parameters is demonstrated. A feasibility study for multiple small boundary control support is included as well. Finally, the numerical convergence of the optimization algorithm and the parallel efficiency is demonstrated.