Each magnetic island (macrospin) in artificial spin ice usually have a form anisotropy (for example length, width, and thickness of 220 nm \(\times \)  80 nm \(\times \)  25 nm in [1]} or 63 nm \(\times \)  26 nm \(\times \)  6 nm in [2]} or even 1000 nm \(\times \)  75 nm \(\times \)  10 nm in [3]}), which limits the degree of freedom of the island magnetic moment and causes it to align along the long axis, at the same time remaining single-domain. In this work, the artificial spin ice model with dipole-dipole interaction described in the previous study [4]}, [5]} also in [6]}, [7]} is investigated. According to this model, each macrospin is represented by a point dipole, the magnetic moment of which have the Ising-like behavior. Detailed study of the usage of dipole interaction between two-dimensional magnetic particles was carried out by Politi and Pini in [8]} and Leon and Pozo in [9]}. Thus, in this approach the kagome artificial spin ice model of hexagonal shape (for symmetry) with free boundary conditions is investigated (see Fig. REF a). Each macrospin represented by point dipole (denoted by dots) with Ising-like magnetic moment (denoted by arrows) which equal \(\sigma _i \vec{m}_i\) , where \(\sigma _i=\pm 1\)  and \(\vec{m}_i=\lbrace 0,-1\rbrace \)  for \(S_1\) , \(\vec{m}_i=\lbrace \sqrt{3/4},0.5\rbrace \)  for \(S_2\)  and \(\vec{m}_i=\lbrace \sqrt{3/4},-0.5\rbrace \)  for \(S_3\)  on the Fig.REF a example. The system size changing was occurred by successively adding layers to the hexagonal form to preserve the shape symmetry. The dipole-dipole interaction energy for two \(ij\)  spins was
\(E^{ij}_{dip}=D \left(\frac{({\vec{m}}_i{\vec{m}}_j)}{\vert {\vec{r}}_{ij}\vert ^3}-3\frac{({\vec{m}}_i{\vec{r}}_{ij})({\vec{m}}_j{\vec{r}}_{ij})}{\vert {\vec{r}}_{ij}\vert ^5}\right),\) 
