第七色在线视频,2021少妇久久久久久久久久,亚洲欧洲精品成人久久av18,亚洲国产精品特色大片观看完整版,孙宇晨将参加特朗普的晚宴

為了賬號(hào)安全,請(qǐng)及時(shí)綁定郵箱和手機(jī)立即綁定
已解決430363個(gè)問(wèn)題,去搜搜看,總會(huì)有你想問(wèn)的

從兩個(gè)列表快速創(chuàng)建對(duì)角矩陣

從兩個(gè)列表快速創(chuàng)建對(duì)角矩陣

冉冉說(shuō) 2021-08-24 18:00:14
給定 2 個(gè)列表,我想創(chuàng)建一個(gè)Diagonal matrix.一個(gè)列表將填充diagonal-constant,另一個(gè)將填充矩陣。例如:fast_matrix([1,2], [6,7,8])應(yīng)該輸出2個(gè)矩陣:[2] # unused[1, 6, 7, 8][6, 1, 7, 8][6, 7, 1, 8][6, 7, 8, 1]和[1] # unused[2, 6, 7, 8][6, 2, 7, 8][6, 7, 2, 8][6, 7, 8, 2]我的代碼在 2.5 秒內(nèi)在我的電腦上進(jìn)行了 10000 次轉(zhuǎn)換。from pprint import pprintimport timeitdef not_so_fast_matrix(A, B):    rt_obj = []    for i,_ in enumerate(A):        for z in range(len(B) + 1):            new_from = A.copy()            new_from.remove(A[i])            new_list = B.copy()            new_list.insert(z, A[i])            rt_obj.append({'remain': new_from, 'to_list': new_list})    return rt_obj# pprint(not_so_fast_matrix([1,2], [6,7,8]))A = ([1,2,3,4,5,6,7,8,9,10])B = ([60,70,80,90,100,200,300])t = timeit.Timer(lambda: not_so_fast_toeplitz(A, B))print("not_so_fast_matrix took: {:.3f}secs for 10000 iterations".format(t.timeit(number=10000)))我想知道使用另一種方法是否可以更快。Circulantfromscipy.linalg看起來(lái)像我想要的但沒(méi)有滾動(dòng):from scipy.linalg import circulantprint(circulant([1, 8,7,6])) # <- Should be invertedoutputs:[[1 6 7 8] [8 1 6 7] [7 8 1 6] [6 7 8 1]]元素被向右移動(dòng)(推)。
查看完整描述

1 回答

?
精慕HU

TA貢獻(xiàn)1845條經(jīng)驗(yàn) 獲得超8個(gè)贊

轉(zhuǎn)置和展平矩陣具有重復(fù)來(lái)自 B 的元素的結(jié)構(gòu)。這種方法使用該屬性來(lái)創(chuàng)建對(duì)角線上具有錯(cuò)誤值的矩陣藍(lán)圖,然后用正確的值填充對(duì)角線。


import numpy


def create_matrices(A, B):

    # Create blueprint of result matrix

    n = len(B) + 1

    b = numpy.empty((n * n,), dtype=numpy.int32)

    b[:-1] = numpy.repeat(B, n + 1)

    b = b.reshape((n, n)).T  # <- added transposition

    # Change diagonal elements

    for a in A:

        m = b.copy()

        numpy.fill_diagonal(m, a)

        print(m)


create_matrices([1, 2], [6, 7, 8])


查看完整回答
反對(duì) 回復(fù) 2021-08-24
  • 1 回答
  • 0 關(guān)注
  • 186 瀏覽
慕課專(zhuān)欄
更多

添加回答

舉報(bào)

0/150
提交
取消
微信客服

購(gòu)課補(bǔ)貼
聯(lián)系客服咨詢(xún)優(yōu)惠詳情

幫助反饋 APP下載

慕課網(wǎng)APP
您的移動(dòng)學(xué)習(xí)伙伴

公眾號(hào)

掃描二維碼
關(guān)注慕課網(wǎng)微信公眾號(hào)