In order to show and prove the proposed approach, systems of ferromagnetic sphere, regular box (a = b = c) and tetragonal box (a = b, c = 2a) were simulated with the use of different values of the scaling factor n. The chosen calculation procedure was the same as the one presented in [22,24,25]. The system parameters were as follows: number of system nodes 40 × 40 × 40 = 64,000, r = 0.28 nm, Jij = 1 × 10−2 eV, Ki = 0 (perfect soft magnet), Si = 1, kBT = 1 × 10−5 eV, D = 2.18 eVnm3 and n = 1, 10, 50 and 100. The other parameters appearing in the presented algorithm are Nrelax = Navr = 400, Pcl = 0.001 and θ = π/100. The initial 3D systems (for n = 1) and their magnetic states (H = 0) for n equal to 10, 50 and 100 are depicted in Figure 3, Figure 4 and Figure 5, respectively. The arrows represent magnetic moments assigned to the nodes, and color depends on the arrow direction. For all cases with n = 1, the magnetic moment alignment is ferromagnetic because the contribution of the dipolar energy is marginal. With the development of the system size, the appearance of different magnetic structures dependent on the n parameter and the shape of the objects can be observed. In the case of n = 50 and for all analyzed objects, the depicted configurations of magnetic moments have pure vortex (spheres) and vortex-like (boxes) characters. For the higher value of n, more complex magnetic structures were detected. Such magnetic structures are expected and experimentally observed in, for example, magnetically soft amorphous and nanocrystalline iron-based alloys, such as the so-called fingerprint domains.