Efficient computation of Jacobian matrices for entropy stable summation by parts schemes