Higher order optimization and adaptive numerical solution for optimal control of monodomain equations in cardiac electrophysiology