scipy inline 问题

shihyu

1. from scipy.weave import inline
   
2. 
   
3. vec = 100
   
4. code = """
   
5. double sum = 0.0;
   
6. int npts = Nvec[0];
   
7. for (int i=0;i<npts;++i) {
   
8.     sum += vec[i];
   
9. }
   
10. return_val = sum;
    
11. """
    
12. print inline(code,['vec'])
    
13. print vec.sum()
    
14. 
    

复制代码

/home/shihyu/.python26_compiled/sc_d6f219b85e9586396248b98781b97fd25.cpp: In function ‘PyObject* compiled_func(PyObject*, PyObject*)’:
/home/shihyu/.python26_compiled/sc_d6f219b85e9586396248b98781b97fd25.cpp:667: error: ‘Nvec’ was not declared in this scope
/home/shihyu/.python26_compiled/sc_d6f219b85e9586396248b98781b97fd25.cpp:669: error: invalid types ‘int[int]’ for array subscript
/home/shihyu/.python26_compiled/sc_d6f219b85e9586396248b98781b97fd25.cpp: In function ‘PyObject* compiled_func(PyObject*, PyObject*)’:
/home/shihyu/.python26_compiled/sc_d6f219b85e9586396248b98781b97fd25.cpp:667: error: ‘Nvec’ was not declared in this scope
/home/shihyu/.python26_compiled/sc_d6f219b85e9586396248b98781b97fd25.cpp:669: error: invalid types ‘int[int]’ for array subscript
Traceback (most recent call last):
  File "test.py", line 12, in <module>
    print inline(code,['vec'])
  File "/usr/lib/python2.6/dist-packages/scipy/weave/inline_tools.py", line 335, in inline
    **kw)
  File "/usr/lib/python2.6/dist-packages/scipy/weave/inline_tools.py", line 462, in compile_function
    verbose=verbose, **kw)
  File "/usr/lib/python2.6/dist-packages/scipy/weave/ext_tools.py", line 365, in compile
    verbose = verbose, **kw)
  File "/usr/lib/python2.6/dist-packages/scipy/weave/build_tools.py", line 272, in build_extension
    setup(name = module_name, ext_modules = [ext],verbose=verb)
  File "/usr/lib/python2.6/dist-packages/numpy/distutils/core.py", line 184, in setup
    return old_setup(**new_attr)
  File "/usr/lib/python2.6/distutils/core.py", line 169, in setup
    raise SystemExit, "error: " + str(msg)
scipy.weave.build_tools.CompileError: error: Command "g++ -pthread -fno-strict-aliasing -DNDEBUG -g -fwrapv -O2 -fPIC -I/usr/lib/python2.6/dist-packages/scipy/weave -I/usr/lib/python2.6/dist-packages/scipy/weave/scxx -I/usr/lib/python2.6/dist-packages/numpy/core/include -I/usr/include/python2.6 -c /home/shihyu/.python26_compiled/sc_d6f219b85e9586396248b98781b97fd25.cpp -o /tmp/shihyu/python26_intermediate/compiler_2da6387b1d12110fba46fe47fea9326a/home/shihyu/.python26_compiled/sc_d6f219b85e9586396248b98781b97fd25.o" failed with exit status 1


請問是什麼原因造成?
謝謝

jacuro

没用过scipy.但从错误消息看，应该是你内嵌的那段c/c++代码包含错误(比如Nvec没有声明……)。
你可以试着直接写一段c代码，单独编译，看看有没问题，没问题了再继续inline.

tianfeng0418

还没有学到这里。