Advertisement
Guest User

Untitled

a guest
Apr 28th, 2017
57
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.46 KB | None | 0 0
  1. # coding: utf-8
  2. # Here your code !
  3.  
  4. # 疎行列とベクトルの積
  5.  
  6. import numpy as np
  7. from scipy.sparse import lil_matrix, csr_matrix
  8.  
  9. def csr_SpMV(A, b):
  10. return A.dot(b)
  11.  
  12. #全体行列の大きさ
  13. n = 10 # 10×10 行列
  14.  
  15. A = lil_matrix( (n, n) )
  16.  
  17. for i in range(1, n - 1):
  18. for j in range(2):
  19. A[i, i + j] = i + j
  20.  
  21. A[0, 0] = 1
  22.  
  23. A[n - 1, n - 1] = n - 1
  24.  
  25. A = A.tocsr()
  26.  
  27. # 密行列に変換して内容確認
  28. print(A.todense())
  29. print(A.dot(A).todense())
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement