where \(\gamma >0\)  controls the contribution of the cluster loss, \(L_c\)  and is typically set at \(\gamma = 0.1\)  [1]}. Here it is important to also keep the reconstruction loss \(L_r\)  to enable good representations of the images [2]}. The cluster loss is often defined using the Kullback-Leibler (KL) divergence [3]}:
\(L_c=KL(P||Q)=\sum _{i=1}^N \sum _{k=1}^K p_{ik} log \frac{p_{ik}}{q_{ik}}\) 
