Why linear algebra is the foundation of scientific computing in SciPy - Challenge Your Understanding

Challenge - 5 Problems

❓ Predict Output

intermediate

Matrix multiplication output

What is the output of the following code that multiplies two matrices using SciPy?

SciPy

import numpy as np
from scipy.linalg import blas

A = np.array([[1, 2], [3, 4]])
B = np.array([[5, 6], [7, 8]])

result = blas.dgemm(alpha=1.0, a=A, b=B)
print(result)

A

[[5. 12.]
 [21. 32.]]

B

[[12. 14.]
 [30. 34.]]

C

[[17. 20.]
 [39. 46.]]

D

[[19. 22.]
 [43. 50.]]

💡 Hint

Recall how matrix multiplication sums products of rows and columns.

Matrix multiplication multiplies rows of the first matrix by columns of the second. For example, element (0,0) is 1*5 + 2*7 = 19.

❓ data_output

intermediate

Shape of matrix after transpose and multiplication

Given a matrix A of shape (3, 4), what is the shape of the matrix resulting from A.T @ A?

SciPy

import numpy as np
A = np.random.rand(3, 4)
result = A.T @ A
print(result.shape)

A (4, 4)

B (3, 3)

C (3, 4)

D (4, 3)

💡 Hint

Transpose flips rows and columns. Matrix multiplication shape depends on inner dimensions.

A.T has shape (4, 3). Multiplying (4, 3) by (3, 4) results in (4, 4).

❓ visualization

advanced

Visualizing eigenvectors of a matrix

Which option correctly plots the eigenvectors of matrix M as arrows starting from the origin?

SciPy

import numpy as np
import matplotlib.pyplot as plt

M = np.array([[2, 1], [1, 2]])
eigenvalues, eigenvectors = np.linalg.eig(M)

plt.figure()
plt.axhline(0, color='grey')
plt.axvline(0, color='grey')
for vec in eigenvectors.T:
    plt.arrow(0, 0, vec[0], vec[1], head_width=0.05, head_length=0.1, fc='blue', ec='blue')
plt.xlim(-1, 1)
plt.ylim(-1, 1)
plt.title('Eigenvectors of M')
plt.show()

A Plots eigenvalues as points on x-axis

B Plots two arrows from origin pointing along eigenvectors of M

C Plots matrix M as heatmap

D Plots random vectors unrelated to M

💡 Hint

Eigenvectors show directions that stretch or shrink under M.

The code plots arrows from origin along eigenvectors, showing directions preserved by M.

🧠 Conceptual

advanced

Why is linear algebra essential in scientific computing?

Which statement best explains why linear algebra is the foundation of scientific computing?

A It provides tools to represent and solve systems of linear equations common in simulations and data analysis.

B It focuses only on scalar arithmetic which is rarely used in computing.

C It is mainly about calculus and derivatives, unrelated to matrices.

D It deals with string manipulation and text processing.

💡 Hint

Think about how many problems can be expressed as linear systems.

Most scientific problems can be modeled as linear systems or transformations, making linear algebra essential.

🔧 Debug

expert

Identify the error in matrix inversion code

What error does the following code raise when trying to invert matrix A?

SciPy

import numpy as np
A = np.array([[1, 2], [2, 4]])
inv_A = np.linalg.inv(A)
print(inv_A)

A ValueError: shapes not aligned

B TypeError: unsupported operand type(s)

C LinAlgError: Singular matrix

D No error, prints inverse matrix

💡 Hint

Check if matrix A is invertible by its determinant.

Matrix A is singular because its rows are linearly dependent, so it cannot be inverted.