Exact evaluation of bit- and symbol-error rates in a 2-D constellation is a fundamental problem of digital communications, which only for particular modulations and/or bits-to-symbol mapping has closed-form solutions. Here, we propose a general, numerically efficient algorithmic method, which yields the exact results for arbitrary modulation symbols' set and arbitrary bits-to-symbol mapping, and which deals with the case of nonuniform signaling. These three conditions define any digital transmission using memoryless modulation, so the proposed method is a general tool solving all problems tackled in the literature under constraints imposed on one or more of the parameters defining the modulation or signaling type. Such an evaluation tool is of practical importance during the design of the modulation. Our analysis and numerical simulations show the advantages offered by the new method when compared with the bounding techniques