import numpy as np from SimPEG.Utils import sdiag # Creating a numpy array x = np.array([1, 2, 3]) # Creating a sparse diagonal matrix D = sdiag(x) print("Sparse Diagonal Matrix:") print(D.toarray())
Sparse Diagonal Matrix: [[1 0 0] [0 2 0] [0 0 3]]
import numpy as np from SimPEG.Utils import sdiag # Creating a numpy array x = np.array([1, 2, 3]) # Creating a sparse diagonal matrix D = sdiag(x) # Defining a vector y = np.array([4, 5, 6]) # Applying the sparse diagonal z = D.dot(y) print("Resulting Vector with Sparse Diagonal Applied:") print(z)
Resulting Vector with Sparse Diagonal Applied: [ 4 10 18]The package library for SimPEG Utils is "Utils" within the broader SimPEG package for geophysical simulation.