fix DefRegionCell: x,y,z was still swapped (z,y,x)
add square shape
relax() bugfix: should not save output
cleanup cubic anisotropy test