Hi,
I am trying to use graph.MCP_Geometric(costs) on a cluster. The script and methods I use work on my personal computer which has 64GB memory. I am now trying to implement the same code on a cluster using IPython parallel. Here's the problem, each node/task only has about 20GB memory so they cannot even complete the
mcp=graph.MCP_Geometric(costs)
command. My cost data has shape (12837, 43345).
I am trying IPython parallel since I think it should allow me to overcome the lack of memory in each node. E.g. If I define a parallel function like
@dview.parallel(block=True) def costsnx(costs): mcp=graph.MCP_Geometric(costs) return mcp
and run it
mcpp=costsnx(costs)
I get a list with mcp objects like <skimage.graph._mcp.MCP_Geometric at 0x2d6fdd0>, one for each process I created. So in principle this would solve the memory problem, but the individual mcp objects do not seem to work e.g. executing
mcpp[30].find_costs([[0,0]])
generates
TypeError: object of type 'NoneType' has no len()
Also, there is no way to join them to recreate the real graph. I tried also including the
mcp.find_costs(location)
command in the parallel function, but it cannot find the location. I was thinking that perhaps there would be a way to use the intermediate steps of the graph.MCP_Geometric function to get the indices or something that could be constructed in parallel and then joined back into a unique element to compute minimum_costs.
Any ideas on how this could be implemented? Or any idea on how to tackle the problem at all?
I appreciate any help or pointers.
Thanks
OOM