我想将1x1矩阵乘以8x8矩阵。显然我知道你不能将这两种形状的矩阵相乘。我的问题是如何从1x1矩阵中“提取”该值,使其相当于将8x8矩阵乘以标量。换句话说,有没有办法将1x1矩阵转换为标量?
到目前为止,这是我的代码,其中n
是我的1x1矩阵,flux
是我的8x8矩阵:
n=0
for i in range(delta_E.shape[0]):
n+= 100/(210*(Sig_f_cell[i])*flux[i]*delta_E[i]*(1.6022e-13)*V_core)
flux = (np.linalg.inv(L))*G
目标:将通量乘以n的值
看起来n是标量但是当我乘以它时我得到了这个错误:
ValueError Traceback (most recent call last)
<ipython-input-26-0df98fb5a138> in <module>()
----> 1 Design_Data (1.34,.037,90)
<ipython-input-25-5ef77d3433bc> in Design_Data(pitch, Pu_fraction, FE_length)
201 print('Number of fuel elements : ',N_FE)
202
--> 203 return n*flux
204
C:\Users\Katey\Anaconda3\lib\site-packages\numpy\matrixlib\defmatrix.py in __mul__(self, other)
341 if isinstance(other, (N.ndarray, list, tuple)) :
342 # This promotes 1-D vectors to row vectors
--> 343 return N.dot(self, asmatrix(other))
344 if isscalar(other) or not hasattr(other, '__rmul__') :
345 return N.dot(self, other)
ValueError: shapes (1,1) and (8,1) not aligned: 1 (dim 1) != 8 (dim 0)
我也试过加倍n[0]*flux
,我得到同样的错误。
答案 0 :(得分:2)
您可以使用numpy.multiply
功能。给定1x1数组和8x8数组,此函数将通过1x1数组中的元素将8x8数组中的每个元素复用。如果我理解你的问题,这就是你正在寻找的。这是文档中的使用示例
>>> x1 = np.arange(9.0).reshape((3, 3))
>>> x2 = np.arange(3.0)
>>> np.multiply(x1, x2)
array([[ 0., 1., 4.],
[ 0., 4., 10.],
[ 0., 7., 16.]])
您可以在https://docs.scipy.org/doc/numpy/reference/generated/numpy.multiply.html找到此功能的文档。