1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
|
/*
Micromagnetic standard problem 4 (a) according to
http://www.ctcms.nist.gov/~rdm/mumag.org.html
*/
// geometry
setgridsize(128, 32, 1)
setcellsize(500e-9/128, 125e-9/32, 3e-9)
// material
Msat = 1600e3
Aex = 13e-12
E_total.get() // try to trigger bad lex2
Msat = 800e3
alpha = 0.02
m = uniform(1, .1, 0)
// relax
relax()
save(m)
TOL := 1e-5
expectv("m", m.average(), vector(0.9669684171676636, 0.1252732127904892, 0), TOL)
// reversal
tableautosave(10e-12)
autosave(m, 100e-12)
autosnapshot(m, 50e-12)
B_ext = vector(-24.6E-3, 4.3E-3, 0)
run(1e-9)
expectv("m", m.average(), vector(-0.9846124053001404, 0.12604089081287384, 0.04327124357223511), TOL)
|
