-
Notifications
You must be signed in to change notification settings - Fork 1
/
magnetoelastic.go
50 lines (42 loc) · 1.6 KB
/
magnetoelastic.go
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
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
package cuda
import (
"github.com/kuchkin/mumax3-gneb/data"
"github.com/kuchkin/mumax3-gneb/util"
)
// Add magneto-elasticit coupling field to the effective field.
// see magnetoelasticfield.cu
func AddMagnetoelasticField(Beff, m *data.Slice, exx, eyy, ezz, exy, exz, eyz, B1, B2, Msat MSlice) {
util.Argument(Beff.Size() == m.Size())
util.Argument(Beff.Size() == exx.Size())
util.Argument(Beff.Size() == eyy.Size())
util.Argument(Beff.Size() == ezz.Size())
util.Argument(Beff.Size() == exy.Size())
util.Argument(Beff.Size() == exz.Size())
util.Argument(Beff.Size() == eyz.Size())
N := Beff.Len()
cfg := make1DConf(N)
k_addmagnetoelasticfield_async(Beff.DevPtr(X), Beff.DevPtr(Y), Beff.DevPtr(Z),
m.DevPtr(X), m.DevPtr(Y), m.DevPtr(Z),
exx.DevPtr(0), exx.Mul(0), eyy.DevPtr(0), eyy.Mul(0), ezz.DevPtr(0), ezz.Mul(0),
exy.DevPtr(0), exy.Mul(0), exz.DevPtr(0), exz.Mul(0), eyz.DevPtr(0), eyz.Mul(0),
B1.DevPtr(0), B1.Mul(0), B2.DevPtr(0), B2.Mul(0),
Msat.DevPtr(0), Msat.Mul(0),
N, cfg)
}
// Calculate magneto-elasticit force density
// see magnetoelasticforce.cu
func GetMagnetoelasticForceDensity(out, m *data.Slice, B1, B2 MSlice, mesh *data.Mesh) {
util.Argument(out.Size() == m.Size())
cellsize := mesh.CellSize()
N := mesh.Size()
cfg := make3DConf(N)
rcsx := float32(1.0 / cellsize[X])
rcsy := float32(1.0 / cellsize[Y])
rcsz := float32(1.0 / cellsize[Z])
k_getmagnetoelasticforce_async(out.DevPtr(X), out.DevPtr(Y), out.DevPtr(Z),
m.DevPtr(X), m.DevPtr(Y), m.DevPtr(Z),
B1.DevPtr(0), B1.Mul(0), B2.DevPtr(0), B2.Mul(0),
rcsx, rcsy, rcsz,
N[X], N[Y], N[Z],
mesh.PBC_code(), cfg)
}