Why custom ufuncs matter in NumPy - Performance Analysis

Time Complexity: Why custom ufuncs matter

O(n)

Understanding Time Complexity

We want to see how using custom ufuncs affects the time it takes to run numpy operations.

How does the speed change when we replace normal Python loops with custom ufuncs?

Scenario Under Consideration

Analyze the time complexity of the following code snippet.

import numpy as np

def square(x):
    return x * x

vectorized_square = np.frompyfunc(square, 1, 1)
arr = np.arange(1_000_000)
result = vectorized_square(arr)

This code creates a custom ufunc to square each element in a large array.

Identify Repeating Operations

Look at what repeats in this code.

• Primary operation: Applying the square function to each element.
• How many times: Once for every element in the array (1,000,000 times).

How Execution Grows With Input

As the array gets bigger, the time to run grows in a clear way.

Input Size (n)	Approx. Operations
10	10
100	100
1000	1000

Pattern observation: The number of operations grows directly with the number of elements.

Final Time Complexity

Time Complexity: O(n)

This means the time to complete the operation grows in a straight line with the input size.

Common Mistake

[X] Wrong: "Custom ufuncs make the operation run in constant time regardless of input size."

[OK] Correct: Even with custom ufuncs, each element must be processed, so time grows with input size.

Interview Connect

Understanding how custom ufuncs scale helps you explain performance improvements clearly and confidently.

Self-Check

"What if we used a built-in numpy ufunc instead of a custom one? How would the time complexity change?"

Practice

1. Why do custom ufuncs matter in numpy?

easy

A. They let you apply your own functions element-wise on arrays easily.

B. They are used to create plots from arrays.

C. They replace the need for Python loops in all cases.

D. They automatically fix errors in array data.

Solution

1. Step 1: Understand the purpose of custom ufuncs

   Custom ufuncs allow you to apply your own math or logic functions element-wise on numpy arrays, making operations simpler.

2. Step 2: Compare with other options

   The other options describe unrelated or incorrect uses of custom ufuncs.

3. Final Answer:

   They let you apply your own functions element-wise on arrays easily. -> Option A

4. Quick Check:

   Custom ufuncs = element-wise custom functions [OK]

Hint: Custom ufuncs = your function applied to each array item [OK]

Common Mistakes:

• Thinking custom ufuncs create plots
• Believing they fix data errors automatically
• Confusing them with general Python loops

2. Which of the following is the correct way to create a custom ufunc using np.vectorize?

easy

A. ufunc = np.vectorize(my_function())

B. ufunc = np.vectorize[my_function]

C. ufunc = np.vectorize(my_function)

D. ufunc = np.vectorize{my_function}

Solution

1. Step 1: Recall the syntax of np.vectorize

   The correct syntax to create a vectorized function is np.vectorize(function_name) without brackets after the function name.

2. Step 2: Analyze each option

   ufunc = np.vectorize(my_function) uses parentheses correctly. ufunc = np.vectorize[my_function] uses square brackets which is invalid. ufunc = np.vectorize{my_function} uses curly braces which is invalid. ufunc = np.vectorize(my_function()) calls the function immediately, which is incorrect.

3. Final Answer:

   ufunc = np.vectorize(my_function) -> Option C

4. Quick Check:

   np.vectorize(function) uses parentheses [OK]

Hint: Use parentheses with np.vectorize, not brackets or braces [OK]

Common Mistakes:

• Using square brackets instead of parentheses
• Calling the function inside np.vectorize
• Using curly braces instead of parentheses

3. What is the output of this code?

import numpy as np

def square_plus_one(x):
    return x**2 + 1

vec_func = np.vectorize(square_plus_one)
arr = np.array([1, 2, 3])
print(vec_func(arr))

medium

A. [2 5 10]

B. [1 4 9]

C. [1 2 3]

D. [3 5 7]

Solution

1. Step 1: Understand the function applied

   The function squares the input and adds 1: x**2 + 1.

2. Step 2: Apply function to each array element

   For arr = [1, 2, 3]: 1**2+1=2, 2**2+1=5, 3**2+1=10.

3. Final Answer:

   [2 5 10] -> Option A

4. Quick Check:

   Square plus one applied element-wise = [2 5 10] [OK]

Hint: Apply function to each element: square then add one [OK]

Common Mistakes:

• Forgetting to add 1 after squaring
• Confusing output with original array
• Mixing up element-wise operation

4. Identify the error in this code that tries to create a custom ufunc:

import numpy as np

def add_three(x):
    return x + 3

ufunc = np.vectorize(add_three())
arr = np.array([1, 2, 3])
print(ufunc(arr))

medium

A. Missing import statement for numpy.

B. Calling add_three() inside np.vectorize instead of passing the function.

C. Using np.vectorize instead of np.frompyfunc.

D. Array arr is not defined.

Solution

1. Step 1: Check how np.vectorize is called

   np.vectorize expects a function, but add_three() calls the function immediately, passing its result instead.

2. Step 2: Understand the impact of the error

   Calling add_three() returns an integer, not a function, causing np.vectorize to fail or behave incorrectly.

3. Final Answer:

   Calling add_three() inside np.vectorize instead of passing the function. -> Option B

4. Quick Check:

   Pass function name, don't call it inside np.vectorize [OK]

Hint: Pass function name to np.vectorize, don't call it [OK]

Common Mistakes:

• Calling function inside np.vectorize
• Confusing function call with function reference
• Ignoring error messages about wrong argument type

5. You want to create a custom ufunc that returns the sign of each element in an array: -1 if negative, 0 if zero, and 1 if positive. Which code correctly creates and applies this ufunc to np.array([-2, 0, 3])?

hard

A. def sign_func(x):
if x <= 0:
return 1
else:
return -1
sign_ufunc = np.vectorize(sign_func)
print(sign_ufunc(np.array([-2, 0, 3])))

B. def sign_func(x):
return x > 0
sign_ufunc = np.vectorize(sign_func)
print(sign_ufunc(np.array([-2, 0, 3])))

C. def sign_func(x):
return x
sign_ufunc = np.vectorize(sign_func)
print(sign_ufunc(np.array([-2, 0, 3])))

D. def sign_func(x):
if x < 0:
return -1
elif x == 0:
return 0
else:
return 1
sign_ufunc = np.vectorize(sign_func)
print(sign_ufunc(np.array([-2, 0, 3])))

Solution

1. Step 1: Understand the sign function requirements

   The function must return -1 for negative, 0 for zero, and 1 for positive values.

2. Step 2: Check each option's logic

   def sign_func(x): if x < 0: return -1 elif x == 0: return 0 else: return 1 sign_ufunc = np.vectorize(sign_func) print(sign_ufunc(np.array([-2, 0, 3]))) correctly implements the conditions. def sign_func(x): return x > 0 sign_ufunc = np.vectorize(sign_func) print(sign_ufunc(np.array([-2, 0, 3]))) returns True/False, not -1/0/1. def sign_func(x): return x sign_ufunc = np.vectorize(sign_func) print(sign_ufunc(np.array([-2, 0, 3]))) returns the value itself, not the sign. def sign_func(x): if x <= 0: return 1 else: return -1 sign_ufunc = np.vectorize(sign_func) print(sign_ufunc(np.array([-2, 0, 3]))) reverses sign logic and misses zero case.

3. Final Answer:

   def sign_func(x): if x < 0: return -1 elif x == 0: return 0 else: return 1 sign_ufunc = np.vectorize(sign_func) print(sign_ufunc(np.array([-2, 0, 3]))) -> Option D

4. Quick Check:

   Correct sign logic with if-elif-else = def sign_func(x): if x < 0: return -1 elif x == 0: return 0 else: return 1 sign_ufunc = np.vectorize(sign_func) print(sign_ufunc(np.array([-2, 0, 3]))) [OK]

Hint: Use if-elif-else to cover negative, zero, positive cases [OK]

Common Mistakes:

• Returning boolean instead of -1,0,1
• Ignoring zero case
• Reversing sign conditions