“…In a recent work, Asgari et al [26] showed that probabilistic CA based on the Glauber algorithm [27] is a fast and reliable simulation method for obtaining the critical point of the two-layer Ising and Potts model in the isotropic case (K x = K y = K z ), where K x and K y are the nearest-neighbor interactions within each layer in the x and y directions, respectively and K z is the inter-layer coupling. They have also shown that this approach is useful for the case of different inter-layer coupling of the twolayer Potts model (K x = K y = K z ) and constructed a critical curve for this model [28]. Although most of the works that have been done until now are for qualitative descriptions or for introducing fast methods for solving various Ising models, we have shown in our previous works that the probabilistic CA increases the calculations precision [26,28].…”