[已收纳]python科学计算之np.r_和np.c_

15195775117

np.r_和np.c_感觉是比较实用的功能，一是可以新建更灵活的索引序列，二是合并数组

一、索引序列

np.r_比arange更灵活些，支持多段序列合并:

print(np.r_[1:15:3,-1:-15:-3])

相当于np.arange(1,15,3)和np.arange(-1,-15,-3)的连接

结果：[  1   4   7  10  13  -1  -4  -7 -10 -13]

应用实例，求π:
n=10000000
pi=np.sum(4.0/np.r_[1:n:4,-3:-n:-4])
print(pi)

结果：3.1415924535898356

二、合并数组
a=[[1,2,3],[4,5,6]]
b=[[-1,-2,-3],[-4,-5,-6]]
print(np.r_[a,b])#上下合并，相当于np.vstack((a,b))
print(np.c_[a,b])#左右合并，相当于np.hstack((a,b))

结果：

[[ 1  2  3]

[ 4  5  6]

[-1 -2 -3]

[-4 -5 -6]]

[[ 1  2  3 -1 -2 -3]

[ 4  5  6 -4 -5 -6]]