An appropriate choice of contact forces for anthropomorphic robotic grasping devices is important for achieving a balanced grasp. Too little applied force may cause an object to slip or be dropped, and too much applied force may cause damage to delicate objects. Prior methods of grasping force optimization (GFO) in the literature can be difficult to compare due to variability in the parameters, such as the type of grasping device, the object grasped, and the contact model, among other factors. Additionally, methods are typically tested on a very small number of scenarios and may not be as robust in other settings. This paper presents a detailed analysis of three optimization approaches based on the literature, comparing them on the basis of accuracy and computational efficiency. Numerical examples are provided for three types of grasp commonly performed by the human hand (cylindrical grasp, tip grasp, and tripod grasp) using both soft finger (SF) contact and hard finger (HF) contact friction models. For each method and grasping example, an external force is applied to the object in eighteen different directions to provide a more complete picture of the methods' performance. Contact points between the hand and the object are predetermined (given). A comparison of the results showed that the nonlinear and linear matrix inequality (LMI) approaches perform best in terms of accuracy, while the computational efficiency of the linear method is stronger unless the number of contact points and segments becomes too large. In this case, the nonlinear method performs more quickly. Future work will extend the problem of GFO to real-time implementation, and a related work (briefly addressed here) examines the sensitivity of optimization methods to variability in the contact locations.