numpy의 숫자 통합

Question

나는 아주 간단하게하고 싶지만 numpy의 깊이에서 찾을 수 없습니다. 나는 값으로 주어진 함수를 수치 적으로 그리고 지속적으로 통합하고 싶다. 즉, 단순히 입력 배열의 시작 부분을 합한 배열을 원한다는 의미입니다. 예 :

입력 :

[ 4, 3, 5, 8 ] 

출력 :

[ 4, 7, 12, 20 ] # [ sum(i[0:1]), sum(i[0:2]), sum(i[0:3]), sum(i[0:4]) ] 

꽤 똑바로 앞으로 소리, 그래서 이것이 내가 현재 찾을 수있어 일부 numpy 기능을 쉽게해야합니다 희망입니다.

나는 scipy.integrate.quad()과 같은 물건을 찾았지만 주어진 범위 (a에서 b까지)와 하나의 값을 반환하는 것처럼 보입니다. 출력으로 배열이 필요합니다.

Answer 1

당신은 numpy.cumsum 찾고 :

>>> numpy.cumsum([ 4, 3, 5, 8 ]) 
    array([ 4, 7, 12, 20]) 

Answer 2

간단히 numpy.cumsum()이 필요합니다.

import numpy as np 
a = np.array([ 4, 3, 5, 8 ]) 
print np.cumsum(a) # prints [ 4 7 12 20] 

Answer 3

당신은 당신이 quadpy (pip install quadpy)는 scipy.integrate.quad() 반대로 내 프로젝트는 compution를 벡터화 않을 수 있습니다. 스틱은 많은 간격을 제공하고,이 간격을 통해 모든 정수 값을 반환합니다.

import numpy 
import quadpy 

a = 0.0 
b = 3.0 
h = 1.0e-2 
n = int((b-a)/h) 

x0 = numpy.linspace(a, b, num=n, endpoint=False) 
x1 = x0 + h 
intervals = numpy.stack([x0, x1]) 

vals = quadpy.line_segment.integrate(
     lambda x: numpy.sin(x), 
     intervals, 
     quadpy.line_segment.GaussLegendre(5) 
     ) 

res = numpy.cumsum(vals) 

import matplotlib.pyplot as plt 
plt.plot(x1, numpy.sin(x1), label='f') 
plt.plot(x1, res, label='F') 
plt.legend() 
plt.show()