# -*- coding: utf-8 -*-
"""
Created on Fri Sep 30 08:20:44 2022

@author: Chris
Lab 9 Q2
"""

import numpy as np
from scipy import sparse as sp
from scipy.sparse import linalg as sla

#a)
A = (np.random.randn(20,20)-0.5)**5


A[abs(A)<0.1] = 0.

print(A)

A2 = sp.csc_matrix(A)

print(A2)

#b)
B = (np.random.randn(5000,5000)-0.5)**5
B[abs(B)<0.1] = 0.

print(B)

B3 = sp.csc_matrix(B)

print(B3)


print(B.nbytes/2**20)
print(B3.size/2**20)