Forward discrete probability propagation method for device performance characterization under process variations