mel-force-dmzdz-dmydz.mx3

Nx := 16
Ny := 16
Nz := 1024

csX := 1e-9
csY := 2e-9
csZ := 0.5e-9

setgridsize(Nx, Ny, Nz)
setcellsize(csX, csY, csZ)

Msat = 1000e3
B1 = 1000e3
B2 = 1000e3
enabledemag = false
m = uniform(1.0, 1.0, 1.0)

mask := newVectorMask(Nx, Ny, Nz)

Period := 32e-9
kx := 2 * pi / Period

pre := 1.0 / (2.0 * B1.Average() * kx)
pre2 := 1.0 / (B2.Average() * kx)

for ii := 0; ii < Nx; ii++ {
	for jj := 0; jj < Ny; jj++ {
		for kk := 0; kk < Nz; kk++ {
			r := index2coord(ii, jj, kk)
			z := r.Z()
			my := cos(kx * z)
			mz := sin(kx * z)
			mask.setVector(ii, jj, kk, vector(0.0, my, mz))
		}
	}
}

m.setarray(mask)

save(m)
save(F_mel)

Fmel := F_mel.HostCopy()

errx := -inf
errz := -inf

for ii := 0; ii < Nx; ii++ {
	for jj := 0; jj < Ny; jj++ {
		for kk := 0; kk < Nz; kk++ {
			r := index2coord(ii, jj, kk)
			z := r.Z()

			ref := sin(kx*z) * cos(kx*z)
			ref2 := (cos(kx*z)*cos(kx*z) - sin(kx*z)*sin(kx*z))

			val := Fmel.get(2, ii, jj, kk)
			val2 := Fmel.get(1, ii, jj, kk)

			ez := abs(val*pre - ref)
			ex := abs(val2*pre2 - ref2)

			if ex > errx {
				errx = ex
			}
			if ez > errz {
				errz = ez
			}
		}
	}
}

ii := Nx / 2
jj := Ny / 2
kk := Nz / 2

r := index2coord(ii, jj, kk)
z := r.Z()

ref := sin(kx*z) * cos(kx*z)
ref2 := (cos(kx*z)*cos(kx*z) - sin(kx*z)*sin(kx*z))

val := Fmel.get(2, ii, jj, kk)
val2 := Fmel.get(1, ii, jj, kk)

ez := abs(val*pre - ref)
ex := abs(val2*pre2 - ref2)

ERRMAX := 0.004
ERRMIN := 3e-6
expect("max(Fmel,x)", errx, 0., ERRMAX)
expect("max(Fmel,z)", errz, 0., ERRMAX)
expect("min(Fmel,x)", ex, 0., ERRMIN)
expect("min(Fmel,z)", ez, 0., ERRMIN)