Despite the great importance of the RD theory, the analytic expression of the RD function is only known for some typical distributions with specific distortion measures [1]}, [2]}, making it difficult to solve real problems of interest.
To date, the most well-known numerical method for computing the RD function is the Blahut-Arimoto (BA) algorithm [3]}, [4]}, which solves the classical RD function by alternating minimization.
