<s> import skimage . graph . mcp as mcp <EOL> from numpy . testing import ( assert_array_equal , <EOL> assert_almost_equal , <EOL> ) <EOL> import numpy as np <EOL> a = np . ones ( ( <NUM_LIT:8> , <NUM_LIT:8> ) , dtype = np . float32 ) <EOL> horizontal_ramp = np . array ( [ [ <NUM_LIT:0.> , <NUM_LIT:1.> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , ] , <EOL> [ <NUM_LIT:0.> , <NUM_LIT:1.> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , ] , <EOL> [ <NUM_LIT:0.> , <NUM_LIT:1.> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , ] , <EOL> [ <NUM_LIT:0.> , <NUM_LIT:1.> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , ] , <EOL> [ <NUM_LIT:0.> , <NUM_LIT:1.> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , ] , <EOL> [ <NUM_LIT:0.> , <NUM_LIT:1.> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , ] , <EOL> [ <NUM_LIT:0.> , <NUM_LIT:1.> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , ] , <EOL> [ <NUM_LIT:0.> , <NUM_LIT:1.> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , ] ] ) <EOL> vertical_ramp = np . array ( [ [ <NUM_LIT:0.> , <NUM_LIT:0.> , <NUM_LIT:0.> , <NUM_LIT:0.> , <NUM_LIT:0.> , <NUM_LIT:0.> , <NUM_LIT:0.> , <NUM_LIT:0.> , ] , <EOL> [ <NUM_LIT:1.> , <NUM_LIT:1.> , <NUM_LIT:1.> , <NUM_LIT:1.> , <NUM_LIT:1.> , <NUM_LIT:1.> , <NUM_LIT:1.> , <NUM_LIT:1.> , ] , <EOL> [ <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , ] , <EOL> [ <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , ] , <EOL> [ <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , ] , <EOL> [ <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , ] , <EOL> [ <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , ] , <EOL> [ <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , <NUM_LIT> , ] ] ) <EOL> def test_anisotropy ( ) : <EOL> seeds_for_horizontal = [ ( i , <NUM_LIT:0> ) for i in range ( <NUM_LIT:8> ) ] <EOL> seeds_for_vertcal = [ ( <NUM_LIT:0> , i ) for i in range ( <NUM_LIT:8> ) ] <EOL> for sy in range ( <NUM_LIT:1> , <NUM_LIT:5> ) : <EOL> for sx in range ( <NUM_LIT:1> , <NUM_LIT:5> ) : <EOL> sampling = sy , sx <EOL> m1 = mcp . MCP_Geometric ( a , sampling = sampling , fully_connected = True ) <EOL> costs1 , traceback = m1 . find_costs ( seeds_for_horizontal ) <EOL> m2 = mcp . MCP_Geometric ( a , sampling = sampling , fully_connected = True ) <EOL> costs2 , traceback = m2 . find_costs ( seeds_for_vertcal ) <EOL> assert_array_equal ( costs1 , horizontal_ramp * sx ) <EOL> assert_array_equal ( costs2 , vertical_ramp * sy ) <EOL> if __name__ == "<STR_LIT:__main__>" : <EOL> np . testing . run_module_suite ( ) </s>
