numpy - 是否可以在 scipy.weave.inline 中使用 blitz++ 索引和 blitz 函数

Question

scipy 文档给出了使用 Blitz++ 样式操作weave.blitz()和使用weave.inline(). 还支持Blitz weave.inline()++ 风格的索引和缩减。那将非常方便。如果weave.inline()确实允许 Blitz++ 样式索引，你能告诉我如何从weave.inline()代码中的 numpy 数组中获取 Blitz 数组。非常感激。

Answer 1

这里是一个例子，在调用 weave.inline() 时设置 type_converter = weave.converters.blitz

# -*- coding: utf-8 -*-
import scipy.weave as weave
import numpy as np
import time

def my_sum(a):
    n=int(len(a))
    code="""
    int i;

    double counter;
    counter =0;
    for(i=0;i<n;i++){
        counter=counter+a(i);
    }
    return_val=counter;
    """

    err=weave.inline(  
        code,   
        ['a','n'],   
        type_converters=weave.converters.blitz, 
        compiler="gcc"  
    )
    return err

a = np.arange(0, 10000000, 1.0)
print my_sum(a)