LINUX.ORG.RU

История изменений

Исправление dikiy, (текущая версия) :

ты рукожоп, с таким кодом тебе никакой blas с numpy не поможет :D

этот код бессмысленный. Так как на выходе всегда ноль ибо сумма по a_m a_n a_l вычисляется (при твоих вводных) в ноль.

for i in xrange(el):
    for m in xrange(iInt):
        for n in xrange(iInt):
            for l in xrange(iInt):
                koef = a[m]*a[n]*a[l]
                Bt = np.matmul(D, B)
                dK = np.matmul(koef*B, Bt)
                K += dK


И еще раз повторяю: используй einsum.

Исходная версия dikiy, :

ты рукожоп, с таким кодом тебе никакой blas с numpy не поможет :D

этот код бессмысленный. Так как на выходе всегда ибо сумма по a_m a_n a_l вычисляется (при твоих вводных) в ноль.

for i in xrange(el):
    for m in xrange(iInt):
        for n in xrange(iInt):
            for l in xrange(iInt):
                koef = a[m]*a[n]*a[l]
                Bt = np.matmul(D, B)
                dK = np.matmul(koef*B, Bt)
                K += dK


И еще раз повторяю: используй einsum.