Memory layout (C-order vs Fortran-order) in NumPy - Hands-On Comparison

Memory layout (C-order vs Fortran-order)

📖 Scenario: Imagine you have a small table of numbers representing sales data for 3 products over 4 days. You want to store this data in a way that computers can understand and process efficiently.

🎯 Goal: You will create a 2D array with sales data, then create two copies of this array with different memory layouts: C-order (row-major) and Fortran-order (column-major). Finally, you will print these arrays to see how the data is stored differently.

📋 What You'll Learn

Create a 2D numpy array with exact values

Create a copy of the array with C-order memory layout

Create a copy of the array with Fortran-order memory layout

Print both arrays to compare

💡 Why This Matters

🌍 Real World

Understanding memory layout helps when working with large datasets or optimizing code for speed, especially in scientific computing and machine learning.

💼 Career

Data scientists and engineers often need to optimize data storage and processing. Knowing memory layout concepts helps write efficient numpy code.

1

Create the sales data array

Create a numpy array called sales with these exact values: [[10, 20, 30, 40], [50, 60, 70, 80], [90, 100, 110, 120]]

NumPy

import numpy as np
# Create the numpy array called sales with the given values
# Your code here

Hint

Use np.array() and pass the list of lists exactly as shown.

2

Create C-order and Fortran-order copies

Create two new arrays: sales_c as a copy of sales with order='C', and sales_f as a copy of sales with order='F'.

NumPy

import numpy as np
sales = np.array([[10, 20, 30, 40], [50, 60, 70, 80], [90, 100, 110, 120]])
# Create sales_c as a C-order copy of sales
# Create sales_f as a Fortran-order copy of sales
# Your code here

Hint

Use np.array() with the order parameter set to 'C' or 'F'.

3

Check the memory layout flags

Print the flags['C_CONTIGUOUS'] and flags['F_CONTIGUOUS'] properties for both sales_c and sales_f to see their memory layouts.

NumPy

import numpy as np
sales = np.array([[10, 20, 30, 40], [50, 60, 70, 80], [90, 100, 110, 120]])
sales_c = np.array(sales, order='C')
sales_f = np.array(sales, order='F')
# Print the C_CONTIGUOUS and F_CONTIGUOUS flags for sales_c and sales_f
# Your code here

Hint

Use array.flags['C_CONTIGUOUS'] and array.flags['F_CONTIGUOUS'] to check memory layout.

4

Print the arrays to compare

Print the arrays sales_c and sales_f to see their contents.

NumPy

import numpy as np
sales = np.array([[10, 20, 30, 40], [50, 60, 70, 80], [90, 100, 110, 120]])
sales_c = np.array(sales, order='C')
sales_f = np.array(sales, order='F')
print('sales_c C_CONTIGUOUS:', sales_c.flags['C_CONTIGUOUS'])
print('sales_c F_CONTIGUOUS:', sales_c.flags['F_CONTIGUOUS'])
print('sales_f C_CONTIGUOUS:', sales_f.flags['C_CONTIGUOUS'])
print('sales_f F_CONTIGUOUS:', sales_f.flags['F_CONTIGUOUS'])
# Print sales_c and sales_f arrays
# Your code here

Hint

Use print() to display the arrays.

Practice

1. What does C-order memory layout mean in NumPy arrays?

easy

A. Rows are stored one after another in memory.

B. Columns are stored one after another in memory.

C. Data is stored randomly in memory.

D. Memory layout depends on the array shape only.

Solution

1. Step 1: Understand C-order layout

   C-order means the array stores data row by row, so all elements of the first row come first, then the second row, and so on.

2. Step 2: Compare with other layouts

   Fortran-order stores columns one after another, not rows. Random or shape-based storage is incorrect.

3. Final Answer:

   Rows are stored one after another in memory. -> Option A

4. Quick Check:

   C-order = rows stored sequentially [OK]

Hint: C-order means row-wise storage in memory [OK]

Common Mistakes:

• Confusing C-order with Fortran-order
• Thinking memory layout depends only on shape
• Assuming random storage order

2. Which of the following is the correct way to create a NumPy array with Fortran-order memory layout?

easy

A. np.array([[1, 2], [3, 4]], order='C')

B. np.array([[1, 2], [3, 4]], order='F')

C. np.array([[1, 2], [3, 4]], layout='F')

D. np.array([[1, 2], [3, 4]], memory='Fortran')

Solution

1. Step 1: Recall NumPy array creation syntax

   NumPy uses the parameter 'order' to specify memory layout, with 'C' for C-order and 'F' for Fortran-order.

2. Step 2: Identify correct parameter and value

   np.array([[1, 2], [3, 4]], order='F') uses order='F', which is the correct way to request Fortran-order layout.

3. Final Answer:

   np.array([[1, 2], [3, 4]], order='F') -> Option B

4. Quick Check:

   Use order='F' for Fortran-order arrays [OK]

Hint: Use order='F' to create Fortran-order arrays [OK]

Common Mistakes:

• Using wrong parameter name like 'layout' or 'memory'
• Confusing order='C' with Fortran-order
• Omitting the order parameter

3. What will be the output of the following code?

import numpy as np
arr = np.array([[1, 2], [3, 4]], order='F')
print(arr.flags.f_contiguous)

medium

A. True

B. False

C. Raises an error

D. None

Solution

1. Step 1: Understand the array creation

   The array is created with order='F', so it uses Fortran-order memory layout.

2. Step 2: Check the f_contiguous flag

   The f_contiguous flag is True if the array is stored in Fortran-contiguous order, which matches order='F'.

3. Final Answer:

   True -> Option A

4. Quick Check:

   Fortran-order array has f_contiguous = True [OK]

Hint: f_contiguous is True for Fortran-order arrays [OK]

Common Mistakes:

• Confusing f_contiguous with c_contiguous
• Expecting False for order='F'
• Thinking flags attribute raises error

4. The following code is intended to create a Fortran-order array, but it does not. What is the error?

import numpy as np
arr = np.array([[1, 2], [3, 4]])
arr = np.asfortranarray(arr, order='C')

medium

A. np.array must have order='F' to create Fortran array.

B. order='C' forces C-order, not Fortran-order.

C. np.asfortranarray does not accept 'order' parameter.

D. arr is immutable and cannot change order.

Solution

1. Step 1: Understand np.asfortranarray usage

   np.asfortranarray converts an array to Fortran-order. It does not take an 'order' parameter.

2. Step 2: Analyze the incorrect parameter

   Passing order='C' is an unexpected keyword argument, causing a TypeError that prevents creating the Fortran-order array.

3. Final Answer:

   np.asfortranarray does not accept 'order' parameter. -> Option C

4. Quick Check:

   np.asfortranarray rejects 'order' parameter [OK]

Hint: np.asfortranarray does not accept 'order' parameter [OK]

Common Mistakes:

• Assuming np.asfortranarray accepts order parameter
• Thinking order='C' overrides Fortran layout
• Thinking arrays are immutable in this context

5. You have a large 2D NumPy array representing image data. You want to speed up column-wise operations. Which memory layout should you choose and why?

hard

A. C-order, because rows are stored consecutively, speeding up column access.

B. Either layout works the same for column-wise operations.

C. Use random memory layout for better cache usage.

D. Fortran-order, because columns are stored consecutively, speeding up column access.

Solution

1. Step 1: Understand memory layout impact on access

   For column-wise operations, accessing consecutive memory locations improves speed due to cache efficiency.

2. Step 2: Identify layout matching column-wise access

   Fortran-order stores columns consecutively, so column-wise operations are faster with Fortran-order.

3. Step 3: Evaluate other options

   C-order stores rows consecutively, so column access is slower. Random layout is not a valid option.

4. Final Answer:

   Fortran-order, because columns are stored consecutively, speeding up column access. -> Option D

5. Quick Check:

   For column speed, use Fortran-order [OK]

Hint: For column operations, choose Fortran-order layout [OK]

Common Mistakes:

• Choosing C-order for column-wise speed
• Thinking layout does not affect speed
• Believing random layout improves cache