We develop a numerical approach for quantifying entanglement in mixed quantum states by convex-roof entanglement measures, based on the optimal entanglement witness operator and the minimax optimization method. Our approach is applicable to general e
ntanglement measures and states and is an efficient alternative to the conventional approach based on the optimal pure-state decomposition. Compared with the conventional one, it has two important merits: (i) that the global optimality of the solution is quantitatively verifiable, and (ii) that the optimization is considerably simplified by exploiting the common symmetry of the target state and measure. To demonstrate the merits, we quantify Greenberger-Horne-Zeilinger (GHZ) entanglement in a class of three-qubit full-rank mixed states composed of the GHZ state, the W state, and the white noise, the simplest mixtures of states with different genuine multipartite entanglement, which have not been quantified before this work. We discuss some general properties of the form of the optimal witness operator and of the convex structure of mixed states, which are related to the symmetry and the rank of states.
