# -*- coding: utf-8 -*-
"""
Created on Fri Sep 30 10:33:26 2022

@author: 102194
"""


import numpy as np
import matplotlib.pyplot as plt
from scipy import sparse as sp
from scipy.sparse import linalg as sla
from scipy import linalg as la

AD = (np.random.randn(500,500)-1)**5
AD[abs(AD)<0.1] = 0.


A = sp.csc_matrix(AD)

lu = sla.splu(A)

#print(lu.L)
#plt.spy(lu.L)
plt.spy(lu.U)