使用weave.inline时出现分段错误

Question

I am new to embedding C++ code into Python. I am testing weave.inline. However I get a segmentation fault when running my code. Could someone tell me what I am doing wrong? Here is my code:

from scipy import weave

def cpp_call(np_points):

assert(type(np_points) == type(1))

code = """
double z[np_points+1][np_points+1];

for (int x = 0; x <= np_points; x++)
{
    for (int y = 0; y <= np_points; y++)
    {
        z[x][y] = x*y;
    }
}
"""

return weave.inline(code,'np_points')

Answer 1

I see two issues here:

1) indention - your python function needs to be indented past the def line

2) your argument to the weave.inline() should be a list. See here for details

So the corrected code should look like:

from scipy import weave

def cpp_call(np_points):
  assert(type(np_points) == type(1))  
  code = """
  double z[np_points+1][np_points+1];

  for (int x = 0; x <= np_points; x++)
  {
    for (int y = 0; y <= np_points; y++)
    {
        z[x][y] = x*y;
    }
  }
  """ 
  return weave.inline(code,['np_points']) #Note the very important change here

This code runs fine for me.