/*
	Micromagnetic standard problem 4,
	with 3D discretization and minimize instead of relax
*/

// geometry
	setgridsize(128, 32, 2)
	setcellsize(500e-9/128, 125e-9/32, 3e-9/2)

// material
	Msat = 800e3
	Aex = 13e-12
	m  = uniform(1, .1, 0)

// minimize
	minimize()
	TOL := 1e-3
	expectv("m", m.average(), vector(0.9669952, 0.12521563, 0), TOL)

