On the Numerical Computation of Second Order Control Barrier Functions