//output directory: gen/si3x3-sin-randomMag-sweep-B/si3x3-sin-randomMag.004.044.out/
Msat = 860e3
Aex = 2 * 13e-12
alpha = 0.01
Msat_r1 := 860e3 * (1.0 - 0.0)
Msat.SetRegion(1, Msat_r1)
B := 0.05
f := 100000000.0
periods := 100
OutputFormat = OVF1_TEXT
SetMesh(64, 64, 1, 10e-09, 10e-09, 15e-09, 0, 0, 0)
island1 := rect(150e-9, 40e-9).transl(-210e-9, -105e-9, 0)
island4 := rect(150e-9, 40e-9).transl(-210e-9, 105e-9, 0)
island2 := rect(150e-9, 40e-9).transl(0, -105e-9, 0)
island5 := rect(150e-9, 40e-9).transl(0, 105e-9, 0)
island3 := rect(150e-9, 40e-9).transl(210e-9, -105e-9, 0)
island6 := rect(150e-9, 40e-9).transl(210e-9, 105e-9, 0)
island7 := rect(40e-9, 150e-9).transl(-105e-9, -210e-9, 0)
island8 := rect(40e-9, 150e-9).transl(105e-9, -210e-9, 0)
island9 := rect(40e-9, 150e-9).transl(-105e-9, 0, 0)
island10 := rect(40e-9, 150e-9).transl(105e-9, 0, 0)
island11 := rect(40e-9, 150e-9).transl(-105e-9, 210e-9, 0)
island12 := rect(40e-9, 150e-9).transl(105e-9, 210e-9, 0)
defregion(1, island1)
defregion(2, island2)
defregion(3, island3)
defregion(4, island4)
defregion(5, island5)
defregion(6, island6)
defregion(7, island7)
defregion(8, island8)
defregion(9, island9)
defregion(10, island10)
defregion(11, island11)
defregion(12, island12)
islands := island1
islands = islands.add(island2)
islands = islands.add(island3)
islands = islands.add(island4)
islands = islands.add(island5)
islands = islands.add(island6)
islands = islands.add(island7)
islands = islands.add(island8)
islands = islands.add(island9)
islands = islands.add(island10)
islands = islands.add(island11)
islands = islands.add(island12)
setGeom(islands)
tableAdd(m.Region(1))
tableAdd(m.Region(2))
tableAdd(m.Region(3))
tableAdd(m.Region(4))
tableAdd(m.Region(5))
tableAdd(m.Region(6))
tableAdd(m.Region(7))
tableAdd(m.Region(8))
tableAdd(m.Region(9))
tableAdd(m.Region(10))
tableAdd(m.Region(11))
tableAdd(m.Region(12))
tableAdd(B_ext)
tableAdd(m_full)
tableAdd(m_full.Region(1))
tableAdd(m_full.Region(2))
tableAdd(m_full.Region(3))
tableAdd(m_full.Region(4))
tableAdd(m_full.Region(5))
tableAdd(m_full.Region(6))
tableAdd(m_full.Region(7))
tableAdd(m_full.Region(8))
tableAdd(m_full.Region(9))
tableAdd(m_full.Region(10))
tableAdd(m_full.Region(11))
tableAdd(m_full.Region(12))
m = randomMagSeed(now().nanosecond())
time_start := now()
relax()
snapshot(m)
save(m)
tableautosave((1 / f) / 1000)
B_ext = vector(B*cos(2*pi*f*t), B*sin(2*pi*f*t), 0)
run(periods / f)
snapshot(m)
save(m)
time_end := now()
elapsed := time_end.Sub(time_start)
print("Elapsed time: ", elapsed)
//Elapsed time:  42m37.958018628s