/*
	Test for m.inRegion() quantity
*/

setgridsize(500, 250, 2)
setcellsize(1e-9, 1e-9, 1e-9)

m = vortex(1, 1)

defregion(1, circle(5e-9) )
defregion(2, circle(5e-9).transl(200e-9, 0, 0 ))

tol := 1e-3

expect("mcore", m.region(1).average()[2],  0.42637932, tol)
expect("mdomain", m.region(2).average()[1],  0.999984, tol)
