The standard approach to solve the above problem is by using the projected mirror gradient descent (MGD), where the projection onto \( \mathcal {C}_{\rm fuzzy}\)  is computed according to the Kullback-Leibler (\(\rm KL\) ) metric [1]}. Adding the entropic regularisation \(\epsilon H(P)\)  to the problem yields the following iterative algorithm
\(P_{t+1} = \mathcal {P}_{\mathcal {C}_{\rm fuzzy}}^{\rm KL}\big (P_{t} \, \odot \, e^{- \alpha q_t}\big ),\) 
