def METHOD_NAME(mapdl, cleared, cube_geom_and_mesh):
    for icoord in range(6):
        mapdl.csys(icoord)
        mapdl.dsys(icoord)
        assert np.allclose(
            mapdl.mesh.nodes_in_current_CS, mapdl.nlist().to_array()[:, 1:4], atol=1e-3
        )  # nlist is not as accurate as 'nodes_in_current_CS'