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
32
33
34
|
/*
Test quantity components
*/
setgridsize(100, 50, 2)
setcellsize(5e-9, 4e-9, 3e-9)
setgeom(xrange(0, inf))
Msat = 800e3
Aex = 13e-12
m = uniform(1, 2, 0)
// Magnetization components should be automatically normalized.
tol := 1e-6
expect("mx", m.comp(0).average(), 1/sqrt(5), tol)
expect("my", m.comp(1).average(), 2/sqrt(5), tol)
expect("mz", m.comp(2).average(), 0/sqrt(5), tol)
// Test components of external field
B_ext = vector(1, 2, 3)
expect("B_ext", B_ext.comp(0).average(), 1., tol)
expect("B_ext", B_ext.comp(1).average(), 2., tol)
expect("B_ext", B_ext.comp(2).average(), 3., tol)
// Test that we can save components
save(m)
save(m.comp(0))
save(m.comp(1))
save(m.comp(2))
save(B_eff.Comp(0))
save(B_exch.Comp(0))
save(B_demag.Comp(0))
save(anisU.Comp(0))
|
