Files
TinyTorch/modules/source/05_dense/dense_dev.py
Vijay Janapa Reddi 4fced96023 Fix module test execution issues
- Fixed test functions to only run when modules executed directly
- Added proper __name__ == '__main__' guards to all test calls
- Fixed syntax errors from incorrect replacements in Module 13 and 15
- Modules now import properly without executing tests
- ProductionBenchmarkingProfiler (Module 14) and ProductionMLSystemProfiler (Module 16) fully working
- Other profiler classes present but require full numpy environment to test completely
2025-09-16 00:17:32 -04:00

1468 lines
53 KiB
Python

# ---
# jupyter:
# jupytext:
# text_representation:
# extension: .py
# format_name: percent
# format_version: '1.3'
# jupytext_version: 1.17.1
# ---
# %% [markdown]
"""
# Networks - Neural Network Architectures
Welcome to the Networks module! This is where we compose layers into complete neural network architectures.
## Learning Goals
- Understand networks as function composition: `f(x) = layer_n(...layer_2(layer_1(x)))`
- Build the Sequential network architecture for composing layers
- Create common network patterns like MLPs (Multi-Layer Perceptrons)
- Visualize network architectures and understand their capabilities
- Master forward pass inference through complete networks
## Build → Use → Reflect
1. **Build**: Sequential networks that compose layers into complete architectures
2. **Use**: Create different network patterns and run inference
3. **Reflect**: How architecture design affects network behavior and capability
## What You'll Learn
By the end of this module, you'll understand:
- How simple layers combine to create complex behaviors
- The fundamental Sequential architecture pattern
- How to build MLPs with any number of layers
- Different network architectures (shallow, deep, wide)
- How neural networks approximate complex functions
"""
# %% nbgrader={"grade": false, "grade_id": "networks-imports", "locked": false, "schema_version": 3, "solution": false, "task": false}
#| default_exp core.dense
#| export
import numpy as np
import sys
import os
from typing import List, Union, Optional, Callable
import matplotlib.pyplot as plt
# Import all the building blocks we need - try package first, then local modules
try:
from tinytorch.core.tensor import Tensor
from tinytorch.core.layers import Dense
from tinytorch.core.activations import ReLU, Sigmoid, Tanh, Softmax
except ImportError:
# For development, import from local modules
sys.path.append(os.path.join(os.path.dirname(__file__), '..', '01_tensor'))
sys.path.append(os.path.join(os.path.dirname(__file__), '..', '02_activations'))
sys.path.append(os.path.join(os.path.dirname(__file__), '..', '03_layers'))
from tensor_dev import Tensor
from activations_dev import ReLU, Sigmoid, Tanh, Softmax
from layers_dev import Dense
# %% nbgrader={"grade": false, "grade_id": "networks-welcome", "locked": false, "schema_version": 3, "solution": false, "task": false}
print("🔥 TinyTorch Networks Module")
print(f"NumPy version: {np.__version__}")
print(f"Python version: {sys.version_info.major}.{sys.version_info.minor}")
print("Ready to build neural network architectures!")
# %% [markdown]
"""
## 📦 Where This Code Lives in the Final Package
**Learning Side:** You work in `modules/source/04_networks/networks_dev.py`
**Building Side:** Code exports to `tinytorch.core.networks`
```python
# Final package structure:
from tinytorch.core.networks import Sequential, create_mlp # Network architectures!
from tinytorch.core.layers import Dense, Conv2D # Building blocks
from tinytorch.core.activations import ReLU, Sigmoid, Tanh # Nonlinearity
from tinytorch.core.tensor import Tensor # Foundation
```
**Why this matters:**
- **Learning:** Focused modules for deep understanding
- **Production:** Proper organization like PyTorch's `torch.nn.Sequential`
- **Consistency:** All network architectures live together in `core.networks`
- **Integration:** Works seamlessly with layers, activations, and tensors
"""
# %% [markdown]
"""
## 🔧 DEVELOPMENT
"""
# %% [markdown]
"""
## Step 1: Understanding Neural Networks as Function Composition
### What is a Neural Network?
A neural network is simply **function composition** - chaining simple functions together to create complex behaviors:
```
f(x) = f_n(f_{n-1}(...f_2(f_1(x))))
```
### Real-World Analogy: Assembly Line
Think of an assembly line in a factory:
- **Input:** Raw materials (data)
- **Stations:** Each worker (layer) transforms the product
- **Output:** Final product (predictions)
### The Power of Composition
```python
# Simple functions
def add_one(x): return x + 1
def multiply_two(x): return x * 2
def square(x): return x * x
# Composed function
def complex_function(x):
return square(multiply_two(add_one(x)))
# This is what neural networks do!
```
### Why This Matters
- **Universal Approximation:** MLPs can approximate any continuous function
- **Hierarchical Learning:** Early layers learn simple features, later layers learn complex patterns
- **Composability:** Mix and match layers to create custom architectures
- **Scalability:** Add more layers or make them wider as needed
### From Modules We've Built
- **Tensors:** The data containers that flow through networks
- **Activations:** The nonlinear transformations that enable complex behaviors
- **Layers:** The building blocks that transform data
Now let's build our first network architecture!
"""
# %% [markdown]
"""
## Step 2: Building the Sequential Network
### What is Sequential?
**Sequential** is the most fundamental network architecture - it applies layers in order:
```
Sequential([layer1, layer2, layer3])
→ f(x) = layer3(layer2(layer1(x)))
```
### Why Sequential Matters
- **Foundation:** Every neural network library has this pattern
- **Simplicity:** Easy to understand and implement
- **Flexibility:** Can compose any layers in any order
- **Building Block:** Foundation for more complex architectures
### The Sequential Pattern
```python
# PyTorch style
model = nn.Sequential(
nn.Linear(784, 128),
nn.ReLU(),
nn.Linear(128, 10)
)
# Our TinyTorch style
model = Sequential([
Dense(784, 128),
ReLU(),
Dense(128, 10)
])
```
Let's implement this fundamental architecture!
"""
# %% nbgrader={"grade": false, "grade_id": "sequential-class", "locked": false, "schema_version": 3, "solution": true, "task": false}
#| export
class Sequential:
"""
Sequential Network: Composes layers in sequence
The most fundamental network architecture.
Applies layers in order: f(x) = layer_n(...layer_2(layer_1(x)))
"""
def __init__(self, layers: Optional[List] = None):
"""
Initialize Sequential network with layers.
Args:
layers: List of layers to compose in order (optional, defaults to empty list)
TODO: Store the layers and implement forward pass
APPROACH:
1. Store the layers list as an instance variable
2. Initialize empty list if no layers provided
3. Prepare for forward pass implementation
EXAMPLE:
Sequential([Dense(3,4), ReLU(), Dense(4,2)])
creates a 3-layer network: Dense → ReLU → Dense
HINTS:
- Use self.layers to store the layers
- Handle empty initialization case
"""
### BEGIN SOLUTION
self.layers = layers if layers is not None else []
### END SOLUTION
def forward(self, x: Tensor) -> Tensor:
"""
Forward pass through all layers in sequence.
Args:
x: Input tensor
Returns:
Output tensor after passing through all layers
TODO: Implement sequential forward pass through all layers
APPROACH:
1. Start with the input tensor
2. Apply each layer in sequence
3. Each layer's output becomes the next layer's input
4. Return the final output
EXAMPLE:
Input: Tensor([[1, 2, 3]])
Layer1 (Dense): Tensor([[1.4, 2.8]])
Layer2 (ReLU): Tensor([[1.4, 2.8]])
Layer3 (Dense): Tensor([[0.7]])
Output: Tensor([[0.7]])
HINTS:
- Use a for loop: for layer in self.layers:
- Apply each layer: x = layer(x)
- The output of one layer becomes input to the next
- Return the final result
"""
### BEGIN SOLUTION
# Apply each layer in sequence
for layer in self.layers:
x = layer(x)
return x
### END SOLUTION
def __call__(self, x: Tensor) -> Tensor:
"""Make the network callable: sequential(x) instead of sequential.forward(x)"""
return self.forward(x)
def add(self, layer):
"""Add a layer to the network."""
self.layers.append(layer)
# %% [markdown]
"""
### 🧪 Unit Test: Sequential Network
Let's test your Sequential network implementation! This is the foundation of all neural network architectures.
**This is a unit test** - it tests one specific class (Sequential network) in isolation.
"""
# %% nbgrader={"grade": true, "grade_id": "test-sequential-immediate", "locked": true, "points": 10, "schema_version": 3, "solution": false, "task": false}
# Test Sequential network immediately after implementation
print("🔬 Unit Test: Sequential Network...")
# Create a simple 2-layer network: 3 → 4 → 2
try:
network = Sequential([
Dense(input_size=3, output_size=4),
ReLU(),
Dense(input_size=4, output_size=2),
Sigmoid()
])
print(f"Network created with {len(network.layers)} layers")
print("✅ Sequential network creation successful")
# Test with sample data
x = Tensor([[1.0, 2.0, 3.0]])
print(f"Input: {x}")
# Forward pass
y = network(x)
print(f"Output: {y}")
print(f"Output shape: {y.shape}")
# Verify the network works
assert y.shape == (1, 2), f"Expected shape (1, 2), got {y.shape}"
print("✅ Sequential network produces correct output shape")
# Test that sigmoid output is in valid range
assert np.all(y.data >= 0) and np.all(y.data <= 1), "Sigmoid output should be between 0 and 1"
print("✅ Sequential network output is in valid range")
# Test that layers are stored correctly
assert len(network.layers) == 4, f"Expected 4 layers, got {len(network.layers)}"
print("✅ Sequential network stores layers correctly")
# Test batch processing
x_batch = Tensor([[1.0, 2.0, 3.0], [4.0, 5.0, 6.0]])
y_batch = network(x_batch)
assert y_batch.shape == (2, 2), f"Expected batch shape (2, 2), got {y_batch.shape}"
print("✅ Sequential network handles batch processing")
except Exception as e:
print(f"❌ Sequential network test failed: {e}")
raise
# Show the network architecture
print("🎯 Sequential network behavior:")
print(" Applies layers in sequence: f(g(h(x)))")
print(" Input flows through each layer in order")
print(" Output of layer i becomes input of layer i+1")
print("📈 Progress: Sequential network ✓")
# %% [markdown]
"""
## Step 3: Building Multi-Layer Perceptrons (MLPs)
### What is an MLP?
A **Multi-Layer Perceptron** is the classic neural network architecture:
```
Input → Dense → Activation → Dense → Activation → ... → Dense → Output
```
### Why MLPs are Important
- **Universal approximation**: Can approximate any continuous function
- **Foundation**: Basis for understanding all neural networks
- **Versatile**: Works for classification, regression, and more
- **Simple**: Easy to understand and implement
### MLP Architecture Pattern
```
create_mlp(3, [4, 2], 1) creates:
Dense(3→4) → ReLU → Dense(4→2) → ReLU → Dense(2→1) → Sigmoid
```
### Real-World Applications
- **Tabular data**: Customer analytics, financial modeling
- **Feature learning**: Learning representations from raw data
- **Classification**: Spam detection, medical diagnosis
- **Regression**: Price prediction, time series forecasting
### The MLP Factory Pattern
Instead of manually creating each layer, we'll build a function that creates MLPs automatically!
"""
# %% nbgrader={"grade": false, "grade_id": "create-mlp", "locked": false, "schema_version": 3, "solution": true, "task": false}
#| export
def create_mlp(input_size: int, hidden_sizes: List[int], output_size: int,
activation=ReLU, output_activation=Sigmoid) -> Sequential:
"""
Create a Multi-Layer Perceptron (MLP) network.
Args:
input_size: Number of input features
hidden_sizes: List of hidden layer sizes
output_size: Number of output features
activation: Activation function for hidden layers (default: ReLU)
output_activation: Activation function for output layer (default: Sigmoid)
Returns:
Sequential network with MLP architecture
TODO: Implement MLP creation with alternating Dense and activation layers.
APPROACH:
1. Start with an empty list of layers
2. Add layers in this pattern:
- Dense(input_size → first_hidden_size)
- Activation()
- Dense(first_hidden_size → second_hidden_size)
- Activation()
- ...
- Dense(last_hidden_size → output_size)
- Output_activation()
3. Return Sequential(layers)
EXAMPLE:
create_mlp(3, [4, 2], 1) creates:
Dense(3→4) → ReLU → Dense(4→2) → ReLU → Dense(2→1) → Sigmoid
HINTS:
- Start with layers = []
- Track current_size starting with input_size
- For each hidden_size: add Dense(current_size, hidden_size), then activation
- Finally add Dense(last_hidden_size, output_size), then output_activation
- Return Sequential(layers)
"""
layers = []
current_size = input_size
# Add hidden layers with activations
for hidden_size in hidden_sizes:
layers.append(Dense(current_size, hidden_size))
layers.append(activation())
current_size = hidden_size
# Add output layer with output activation
layers.append(Dense(current_size, output_size))
layers.append(output_activation())
return Sequential(layers)
# %% [markdown]
"""
### 🧪 Unit Test: MLP Creation
Let's test your MLP creation function! This builds complete neural networks with a single function call.
**This is a unit test** - it tests one specific function (create_mlp) in isolation.
"""
# %% nbgrader={"grade": true, "grade_id": "test-mlp-immediate", "locked": true, "points": 10, "schema_version": 3, "solution": false, "task": false}
# Test MLP creation immediately after implementation
print("🔬 Unit Test: MLP Creation...")
# Create a simple MLP: 3 → 4 → 2 → 1
try:
mlp = create_mlp(input_size=3, hidden_sizes=[4, 2], output_size=1)
print(f"MLP created with {len(mlp.layers)} layers")
print("✅ MLP creation successful")
# Test the structure - should have 6 layers: Dense, ReLU, Dense, ReLU, Dense, Sigmoid
expected_layers = 6 # 3 Dense + 2 ReLU + 1 Sigmoid
assert len(mlp.layers) == expected_layers, f"Expected {expected_layers} layers, got {len(mlp.layers)}"
print("✅ MLP has correct number of layers")
# Test layer types
layer_types = [type(layer).__name__ for layer in mlp.layers]
expected_pattern = ['Dense', 'ReLU', 'Dense', 'ReLU', 'Dense', 'Sigmoid']
assert layer_types == expected_pattern, f"Expected pattern {expected_pattern}, got {layer_types}"
print("✅ MLP follows correct layer pattern")
# Test with sample data
x = Tensor([[1.0, 2.0, 3.0]])
y = mlp(x)
print(f"MLP input: {x}")
print(f"MLP output: {y}")
print(f"MLP output shape: {y.shape}")
# Verify the output
assert y.shape == (1, 1), f"Expected shape (1, 1), got {y.shape}"
print("✅ MLP produces correct output shape")
# Test that sigmoid output is in valid range
assert np.all(y.data >= 0) and np.all(y.data <= 1), "Sigmoid output should be between 0 and 1"
print("✅ MLP output is in valid range")
except Exception as e:
print(f"❌ MLP creation test failed: {e}")
raise
# Test different architectures
try:
# Test shallow network
shallow_net = create_mlp(input_size=3, hidden_sizes=[4], output_size=1)
assert len(shallow_net.layers) == 4, f"Shallow network should have 4 layers, got {len(shallow_net.layers)}"
# Test deep network
deep_net = create_mlp(input_size=3, hidden_sizes=[4, 4, 4], output_size=1)
assert len(deep_net.layers) == 8, f"Deep network should have 8 layers, got {len(deep_net.layers)}"
# Test wide network
wide_net = create_mlp(input_size=3, hidden_sizes=[10], output_size=1)
assert len(wide_net.layers) == 4, f"Wide network should have 4 layers, got {len(wide_net.layers)}"
print("✅ Different MLP architectures work correctly")
except Exception as e:
print(f"❌ MLP architecture test failed: {e}")
raise
# Show the MLP pattern
print("🎯 MLP creation pattern:")
print(" Input → Dense → Activation → Dense → Activation → ... → Dense → Output_Activation")
print(" Automatically creates the complete architecture")
print(" Handles any number of hidden layers")
print("📈 Progress: Sequential network ✓, MLP creation ✓")
# %% [markdown]
"""
## Step 4: Understanding Network Architectures
### Architecture Patterns
Different network architectures solve different problems:
#### **Shallow vs Deep Networks**
```python
# Shallow: 1 hidden layer
shallow = create_mlp(10, [20], 1)
# Deep: Many hidden layers
deep = create_mlp(10, [20, 20, 20], 1)
```
#### **Narrow vs Wide Networks**
```python
# Narrow: Few neurons per layer
narrow = create_mlp(10, [5, 5], 1)
# Wide: Many neurons per layer
wide = create_mlp(10, [50], 1)
```
### Why Architecture Matters
- **Capacity:** More parameters can learn more complex patterns
- **Depth:** Enables hierarchical feature learning
- **Width:** Allows parallel processing of features
- **Efficiency:** Balance between performance and computation
### Different Activation Functions
```python
# ReLU networks (most common)
relu_net = create_mlp(10, [20], 1, activation=ReLU)
# Tanh networks (centered around 0)
tanh_net = create_mlp(10, [20], 1, activation=Tanh)
# Multi-class classification
classifier = create_mlp(10, [20], 3, output_activation=Softmax)
```
Let's test different architectures!
"""
# %% [markdown]
"""
### 🧪 Unit Test: Architecture Variations
Let's test different network architectures to understand their behavior.
**This is a unit test** - it tests architectural variations in isolation.
"""
# %% [markdown]
"""
### 📊 Visualization: Network Architecture Comparison
This function creates and visualizes different neural network architectures to demonstrate how activation functions and layer configurations affect network behavior and output characteristics.
"""
# %% nbgrader={"grade": true, "grade_id": "test-architectures", "locked": true, "points": 10, "schema_version": 3, "solution": false, "task": false}
def plot_network_architectures():
"""Visualize different network architectures."""
# Create different architectures
relu_net = create_mlp(input_size=3, hidden_sizes=[4], output_size=1, activation=ReLU)
tanh_net = create_mlp(input_size=3, hidden_sizes=[4], output_size=1, activation=Tanh)
classifier = create_mlp(input_size=3, hidden_sizes=[4], output_size=3, output_activation=Softmax)
# Create input data
x = Tensor([[1.0, 2.0, 3.0]])
# Get outputs
y_relu = relu_net(x)
y_tanh = tanh_net(x)
y_multi = classifier(x)
# Plot the results
fig, axs = plt.subplots(1, 3, figsize=(15, 4))
axs[0].set_title("ReLU Network Output")
axs[0].bar(['Output'], [y_relu.data[0][0]], color='skyblue')
axs[1].set_title("Tanh Network Output")
axs[1].bar(['Output'], [y_tanh.data[0][0]], color='salmon')
axs[2].set_title("Softmax Classifier Output")
axs[2].bar([f"Class {i}" for i in range(3)], y_multi.data[0], color='lightgreen')
plt.tight_layout()
plt.show()
plot_network_architectures()
# %% [markdown]
"""
### 🧪 Unit Test: Network Architecture Variations
This test validates different neural network architectures created with various activation functions. It ensures that networks with ReLU, Tanh, and Softmax activations work correctly, and tests both shallow and deep network configurations for comprehensive architecture validation.
"""
# %%
def test_unit_network_architectures():
"""Unit test for different network architectures."""
# Test different architectures
print("🔬 Unit Test: Network Architecture Variations...")
try:
# Test different activation functions
relu_net = create_mlp(input_size=3, hidden_sizes=[4], output_size=1, activation=ReLU)
tanh_net = create_mlp(input_size=3, hidden_sizes=[4], output_size=1, activation=Tanh)
# Test different output activations
classifier = create_mlp(input_size=3, hidden_sizes=[4], output_size=3, output_activation=Softmax)
# Test with sample data
x = Tensor([[1.0, 2.0, 3.0]])
# Test ReLU network
y_relu = relu_net(x)
assert y_relu.shape == (1, 1), "ReLU network should work"
print("✅ ReLU network works correctly")
# Test Tanh network
y_tanh = tanh_net(x)
assert y_tanh.shape == (1, 1), "Tanh network should work"
print("✅ Tanh network works correctly")
# Test multi-class classifier
y_multi = classifier(x)
assert y_multi.shape == (1, 3), "Multi-class classifier should work"
# Check softmax properties
assert abs(np.sum(y_multi.data) - 1.0) < 1e-6, "Softmax outputs should sum to 1"
print("✅ Multi-class classifier with Softmax works correctly")
# Test different architectures
shallow = create_mlp(input_size=4, hidden_sizes=[5], output_size=1)
deep = create_mlp(input_size=4, hidden_sizes=[5, 5, 5], output_size=1)
wide = create_mlp(input_size=4, hidden_sizes=[20], output_size=1)
x_test = Tensor([[1.0, 2.0, 3.0, 4.0]])
# Test all architectures
for name, net in [("Shallow", shallow), ("Deep", deep), ("Wide", wide)]:
y = net(x_test)
assert y.shape == (1, 1), f"{name} network should produce correct shape"
print(f"{name} network works correctly")
print("✅ All network architectures work correctly")
except Exception as e:
print(f"❌ Architecture test failed: {e}")
raise
print("🎯 Architecture insights:")
print(" Different activations create different behaviors")
print(" Softmax enables multi-class classification")
print(" Architecture affects network capacity and learning")
print("📈 Progress: Sequential ✓, MLP creation ✓, Architecture variations ✓")
# %% [markdown]
"""
### 📊 Visualization Demo: Network Architectures
Let's visualize the different network architectures for educational purposes:
"""
# %% [markdown]
"""
## Step 5: Comprehensive Test - Complete Network Applications
### Real-World Network Applications
Let's test our networks on realistic scenarios:
#### **Classification Problem**
```python
# 4 features → 2 classes (binary classification)
classifier = create_mlp(4, [8, 4], 2, output_activation=Softmax)
```
#### **Regression Problem**
```python
# 3 features → 1 continuous output
regressor = create_mlp(3, [10, 5], 1, output_activation=lambda: Dense(0, 0)) # Linear output
```
#### **Deep Learning Pattern**
```python
# Complex feature learning
deep_net = create_mlp(10, [64, 32, 16], 1)
```
This comprehensive test ensures our networks work for real ML applications!
"""
# %% nbgrader={"grade": true, "grade_id": "test-integration", "locked": true, "points": 15, "schema_version": 3, "solution": false, "task": false}
# Comprehensive test - complete network applications
print("🔬 Comprehensive Test: Complete Network Applications...")
try:
# Test 1: Multi-class Classification (Iris-like dataset)
print("\n1. Multi-class Classification Test:")
iris_classifier = create_mlp(input_size=4, hidden_sizes=[8, 6], output_size=3, output_activation=Softmax)
# Simulate iris features: [sepal_length, sepal_width, petal_length, petal_width]
iris_samples = Tensor([
[5.1, 3.5, 1.4, 0.2], # Setosa
[7.0, 3.2, 4.7, 1.4], # Versicolor
[6.3, 3.3, 6.0, 2.5] # Virginica
])
iris_predictions = iris_classifier(iris_samples)
assert iris_predictions.shape == (3, 3), "Iris classifier should output 3 classes for 3 samples"
# Check softmax properties
row_sums = np.sum(iris_predictions.data, axis=1)
assert np.allclose(row_sums, 1.0), "Each prediction should sum to 1"
print("✅ Multi-class classification works correctly")
# Test 2: Regression Task (Housing prices)
print("\n2. Regression Task Test:")
# Create a regressor without final activation (linear output)
class Identity:
def __call__(self, x): return x
housing_regressor = create_mlp(input_size=3, hidden_sizes=[10, 5], output_size=1, output_activation=Identity)
# Simulate housing features: [size, bedrooms, location_score]
housing_samples = Tensor([
[2000, 3, 8.5], # Large house, good location
[1200, 2, 6.0], # Medium house, ok location
[800, 1, 4.0] # Small house, poor location
])
housing_predictions = housing_regressor(housing_samples)
assert housing_predictions.shape == (3, 1), "Housing regressor should output 1 value per sample"
print("✅ Regression task works correctly")
# Test 3: Deep Network Performance
print("\n3. Deep Network Test:")
deep_network = create_mlp(input_size=10, hidden_sizes=[20, 15, 10, 5], output_size=1)
# Test with realistic batch size
batch_data = Tensor(np.random.randn(32, 10)) # 32 samples, 10 features
deep_predictions = deep_network(batch_data)
assert deep_predictions.shape == (32, 1), "Deep network should handle batch processing"
assert not np.any(np.isnan(deep_predictions.data)), "Deep network should not produce NaN"
print("✅ Deep network handles batch processing correctly")
# Test 4: Network Composition
print("\n4. Network Composition Test:")
# Create a feature extractor and classifier separately
feature_extractor = Sequential([
Dense(input_size=10, output_size=5),
ReLU(),
Dense(input_size=5, output_size=3),
ReLU()
])
classifier_head = Sequential([
Dense(input_size=3, output_size=2),
Softmax()
])
# Test composition
raw_data = Tensor(np.random.randn(5, 10))
features = feature_extractor(raw_data)
final_predictions = classifier_head(features)
assert features.shape == (5, 3), "Feature extractor should output 3 features"
assert final_predictions.shape == (5, 2), "Classifier should output 2 classes"
row_sums = np.sum(final_predictions.data, axis=1)
assert np.allclose(row_sums, 1.0), "Composed network predictions should be valid"
print("✅ Network composition works correctly")
print("\n🎉 Comprehensive test passed! Your networks work correctly for:")
print(" • Multi-class classification (Iris flowers)")
print(" • Regression tasks (housing prices)")
print(" • Deep learning architectures")
print(" • Network composition and feature extraction")
except Exception as e:
print(f"❌ Comprehensive test failed: {e}")
print("📈 Final Progress: Complete network architectures ready for real ML applications!")
# Run the test
test_unit_network_architectures()
# %% [markdown]
"""
### 🏗️ Class: MLP (Multi-Layer Perceptron)
This class provides a convenient wrapper around Sequential networks specifically designed for standard MLP architectures. It maintains parameter information and provides a clean interface for creating and managing multi-layer perceptrons with consistent structure.
"""
# %% nbgrader={"grade": false, "grade_id": "networks-compatibility", "locked": false, "schema_version": 3, "solution": false, "task": false}
#| export
class MLP:
"""
Multi-Layer Perceptron (MLP) class.
A convenient wrapper around Sequential networks for standard MLP architectures.
Maintains parameter information and provides a clean interface.
Args:
input_size: Number of input features
hidden_size: Size of the single hidden layer
output_size: Number of output features
activation: Activation function for hidden layer (default: ReLU)
output_activation: Activation function for output layer (default: Sigmoid)
"""
def __init__(self, input_size: int, hidden_size: int, output_size: int,
activation=ReLU, output_activation=None):
self.input_size = input_size
self.hidden_size = hidden_size
self.output_size = output_size
# Build the network layers
layers = []
# Input to hidden layer
layers.append(Dense(input_size, hidden_size))
layers.append(activation())
# Hidden to output layer
layers.append(Dense(hidden_size, output_size))
if output_activation is not None:
layers.append(output_activation())
self.network = Sequential(layers)
def forward(self, x):
"""Forward pass through the MLP network."""
return self.network.forward(x)
def __call__(self, x):
"""Make the MLP callable."""
return self.forward(x)
# %% [markdown]
"""
### 🧪 Unit Test: Sequential Network Implementation
This test validates the Sequential network class functionality, ensuring proper layer composition, forward pass execution, and network architecture validation for multi-layer neural networks.
"""
# %%
def test_unit_sequential_networks():
"""Unit test for the Sequential network implementation."""
print("🔬 Unit Test: Sequential Networks...")
# Test basic Sequential network
net = Sequential([
Dense(input_size=3, output_size=4),
ReLU(),
Dense(input_size=4, output_size=2),
Sigmoid()
])
x = Tensor([[1.0, 2.0, 3.0]])
y = net(x)
assert y.shape == (1, 2), "Sequential network should produce correct output shape"
assert np.all(y.data > 0), "Sigmoid output should be positive"
assert np.all(y.data < 1), "Sigmoid output should be less than 1"
print("✅ Sequential networks work correctly")
# Run the test
test_unit_sequential_networks()
# %% [markdown]
"""
### 🧪 Unit Test: MLP Creation Function
This test validates the `create_mlp` function, ensuring it correctly constructs Multi-Layer Perceptrons with various architectures, activation functions, and layer configurations for different machine learning tasks.
"""
# %%
def test_unit_mlp_creation():
"""Unit test for the MLP creation function."""
print("🔬 Unit Test: MLP Creation...")
# Test different MLP architectures
shallow = create_mlp(input_size=4, hidden_sizes=[5], output_size=1)
deep = create_mlp(input_size=4, hidden_sizes=[8, 6, 4], output_size=2)
x = Tensor([[1.0, 2.0, 3.0, 4.0]])
# Test shallow network
y_shallow = shallow(x)
assert y_shallow.shape == (1, 1), "Shallow MLP should work"
# Test deep network
y_deep = deep(x)
assert y_deep.shape == (1, 2), "Deep MLP should work"
print("✅ MLP creation works correctly")
# Run the test
test_unit_mlp_creation()
# %% [markdown]
"""
### 🧪 Unit Test: Network Applications in Real ML Scenarios
This comprehensive test validates network performance on real machine learning tasks including classification and regression, ensuring the implementations work correctly with actual datasets and practical applications.
"""
# %%
def test_unit_network_applications():
"""Comprehensive unit test for network applications in real ML scenarios."""
print("🔬 Comprehensive Test: Network Applications...")
# Test multi-class classification
iris_classifier = create_mlp(input_size=4, hidden_sizes=[8, 6], output_size=3, output_activation=Softmax)
iris_samples = Tensor([[5.1, 3.5, 1.4, 0.2], [7.0, 3.2, 4.7, 1.4], [6.3, 3.3, 6.0, 2.5]])
iris_predictions = iris_classifier(iris_samples)
assert iris_predictions.shape == (3, 3), "Iris classifier should work"
row_sums = np.sum(iris_predictions.data, axis=1)
assert np.allclose(row_sums, 1.0), "Predictions should sum to 1"
# Run the test
test_unit_network_applications()
# %% [markdown]
"""
## 🧪 Module Testing
Time to test your implementation! This section uses TinyTorch's standardized testing framework to ensure your implementation works correctly.
**This testing section is locked** - it provides consistent feedback across all modules and cannot be modified.
"""
# %% nbgrader={"grade": false, "grade_id": "standardized-testing", "locked": true, "schema_version": 3, "solution": false, "task": false}
# =============================================================================
# STANDARDIZED MODULE TESTING - DO NOT MODIFY
# This cell is locked to ensure consistent testing across all TinyTorch modules
# =============================================================================
# %% [markdown]
"""
## 🔬 Integration Test: End-to-End Network Forward Pass
"""
# %%
def test_module_full_network_forward_pass():
"""
Integration test for a complete forward pass through a multi-layer network.
Tests a complete forward pass through a multi-layer network,
integrating Tensors, Dense layers, Activations, and the Sequential container.
"""
print("🔬 Running Integration Test: Full Network Forward Pass...")
# 1. Define a simple 2-layer MLP
# Input (3) -> Dense(4) -> ReLU -> Dense(2) -> Output
model = Sequential([
Dense(3, 4),
ReLU(),
Dense(4, 2)
])
# 2. Create a batch of input Tensors
# Batch of 5 samples, each with 3 features
input_tensor = Tensor(np.random.randn(5, 3))
# 3. Perform a forward pass through the entire network
output_tensor = model(input_tensor)
# 4. Assert the final output is correct
assert isinstance(output_tensor, Tensor), "Network output must be a Tensor"
assert output_tensor.shape == (5, 2), f"Expected output shape (5, 2), but got {output_tensor.shape}"
print("✅ Integration Test Passed: Full network forward pass is successful.")
# Run the integration test
test_module_full_network_forward_pass()
# %% [markdown]
"""
## 🔧 ML Systems: Network Stability & Error Handling
Now that you have complete neural networks, let's develop **production robustness skills**. This section teaches you to identify and fix stability issues that can break training in production systems.
### **Learning Outcome**: *"I understand why numerical stability matters in production and can detect/fix stability issues"*
---
## Network Stability Monitor (Medium Guided Implementation)
As an ML systems engineer, you need to ensure networks remain stable during training. Let's build tools to detect numerical instability and understand gradient flow issues.
"""
# %%
import time
import numpy as np
class NetworkStabilityMonitor:
"""
Stability monitoring toolkit for neural networks.
Helps ML engineers detect numerical instability, gradient problems,
and other issues that can break training in production systems.
"""
def __init__(self):
self.stability_history = []
self.warning_threshold = 1e6
self.error_threshold = 1e10
def check_tensor_stability(self, tensor, tensor_name="tensor"):
"""
Check if a tensor has numerical stability issues.
TODO: Implement tensor stability checking.
STEP-BY-STEP IMPLEMENTATION:
1. Check for NaN values using np.isnan()
2. Check for infinite values using np.isinf()
3. Check for extremely large values (> 1e6)
4. Calculate value statistics (min, max, mean, std)
5. Return stability report with warnings
EXAMPLE:
monitor = NetworkStabilityMonitor()
tensor = Tensor([1.0, 2.0, np.inf])
report = monitor.check_tensor_stability(tensor, "weights")
print(f"Stable: {report['is_stable']}")
print(f"Issues: {report['issues']}")
HINTS:
- Use tensor.data to get numpy array
- Check: np.any(np.isnan(tensor.data))
- Check: np.any(np.isinf(tensor.data))
- Check: np.any(np.abs(tensor.data) > self.warning_threshold)
- Return dict with analysis
"""
### BEGIN SOLUTION
data = tensor.data
# Check for numerical issues
has_nan = np.any(np.isnan(data))
has_inf = np.any(np.isinf(data))
has_large = np.any(np.abs(data) > self.warning_threshold)
has_extreme = np.any(np.abs(data) > self.error_threshold)
# Calculate statistics (avoiding issues if all values are problematic)
finite_mask = np.isfinite(data)
if np.any(finite_mask):
finite_data = data[finite_mask]
stats = {
'min': np.min(finite_data),
'max': np.max(finite_data),
'mean': np.mean(finite_data),
'std': np.std(finite_data),
'finite_count': np.sum(finite_mask),
'total_count': data.size
}
else:
stats = {
'min': np.nan,
'max': np.nan,
'mean': np.nan,
'std': np.nan,
'finite_count': 0,
'total_count': data.size
}
# Compile issues
issues = []
if has_nan:
issues.append("Contains NaN values")
if has_inf:
issues.append("Contains infinite values")
if has_extreme:
issues.append(f"Contains extremely large values (>{self.error_threshold:.0e})")
elif has_large:
issues.append(f"Contains large values (>{self.warning_threshold:.0e})")
is_stable = len(issues) == 0
return {
'tensor_name': tensor_name,
'is_stable': is_stable,
'issues': issues,
'has_nan': has_nan,
'has_inf': has_inf,
'has_large_values': has_large,
'statistics': stats
}
### END SOLUTION
def analyze_gradient_flow(self, network, input_tensor, target_output):
"""
Analyze gradient flow through a network to detect vanishing/exploding gradients.
TODO: Implement gradient flow analysis.
STEP-BY-STEP IMPLEMENTATION:
1. Perform forward pass through network
2. Simulate simple loss calculation (MSE)
3. Estimate gradient magnitudes using finite differences
4. Check for vanishing gradients (very small)
5. Check for exploding gradients (very large)
6. Return gradient flow analysis
EXAMPLE:
monitor = NetworkStabilityMonitor()
analysis = monitor.analyze_gradient_flow(network, input_data, target)
print(f"Gradient health: {analysis['gradient_status']}")
HINTS:
- Forward pass: output = network(input_tensor)
- Simple loss: 0.5 * np.sum((output.data - target_output.data)**2)
- Use small perturbations to estimate gradients
- Vanishing: gradients < 1e-6, Exploding: gradients > 1e3
"""
### BEGIN SOLUTION
# Forward pass
output = network(input_tensor)
# Calculate simple MSE loss
loss = 0.5 * np.sum((output.data - target_output.data)**2)
# Estimate gradient magnitudes using finite differences
# This is a simplified approach - real backprop would be more accurate
epsilon = 1e-5
gradient_estimates = []
# Check first layer weights (simplified analysis)
if hasattr(network, 'layers') and len(network.layers) > 0:
first_layer = network.layers[0]
if hasattr(first_layer, 'weights'):
# Perturb a small sample of weights to estimate gradients
original_weight = first_layer.weights.data[0, 0]
# Forward pass with small perturbation
first_layer.weights.data[0, 0] = original_weight + epsilon
output_plus = network(input_tensor)
loss_plus = 0.5 * np.sum((output_plus.data - target_output.data)**2)
# Estimate gradient
grad_estimate = (loss_plus - loss) / epsilon
gradient_estimates.append(abs(grad_estimate))
# Restore original weight
first_layer.weights.data[0, 0] = original_weight
# Analyze gradient magnitudes
if gradient_estimates:
avg_grad = np.mean(gradient_estimates)
max_grad = np.max(gradient_estimates)
if avg_grad < 1e-8:
gradient_status = "Vanishing gradients detected"
elif max_grad > 1e3:
gradient_status = "Exploding gradients detected"
elif avg_grad < 1e-6:
gradient_status = "Potentially vanishing gradients"
elif max_grad > 100:
gradient_status = "Potentially exploding gradients"
else:
gradient_status = "Healthy gradient flow"
else:
gradient_status = "Unable to analyze gradients"
return {
'loss': loss,
'gradient_estimates': gradient_estimates,
'avg_gradient': np.mean(gradient_estimates) if gradient_estimates else 0,
'max_gradient': np.max(gradient_estimates) if gradient_estimates else 0,
'gradient_status': gradient_status
}
### END SOLUTION
def comprehensive_stability_check(self, network, input_tensor, target_output):
"""
Perform comprehensive stability analysis of a neural network.
This function is PROVIDED to demonstrate complete stability monitoring.
Students use it to understand production stability requirements.
"""
print("🔧 COMPREHENSIVE NETWORK STABILITY CHECK")
print("=" * 50)
stability_report = {
'overall_status': 'STABLE',
'issues_found': [],
'recommendations': []
}
# Check input stability
input_check = self.check_tensor_stability(input_tensor, "input")
if not input_check['is_stable']:
stability_report['overall_status'] = 'UNSTABLE'
stability_report['issues_found'].extend([f"Input: {issue}" for issue in input_check['issues']])
stability_report['recommendations'].append("Normalize or clip input data")
print(f"📊 Input Check: {'✅ STABLE' if input_check['is_stable'] else '❌ UNSTABLE'}")
if input_check['issues']:
for issue in input_check['issues']:
print(f" - {issue}")
# Check each layer's weights and outputs
if hasattr(network, 'layers'):
for i, layer in enumerate(network.layers):
if hasattr(layer, 'weights'):
weight_check = self.check_tensor_stability(layer.weights, f"layer_{i}_weights")
if not weight_check['is_stable']:
stability_report['overall_status'] = 'UNSTABLE'
stability_report['issues_found'].extend([f"Layer {i}: {issue}" for issue in weight_check['issues']])
stability_report['recommendations'].append(f"Re-initialize layer {i} weights")
print(f"🔗 Layer {i} Weights: {'✅ STABLE' if weight_check['is_stable'] else '❌ UNSTABLE'}")
if weight_check['issues']:
for issue in weight_check['issues']:
print(f" - {issue}")
# Check network output
try:
output = network(input_tensor)
output_check = self.check_tensor_stability(output, "network_output")
if not output_check['is_stable']:
stability_report['overall_status'] = 'UNSTABLE'
stability_report['issues_found'].extend([f"Output: {issue}" for issue in output_check['issues']])
stability_report['recommendations'].append("Check activation functions and weight initialization")
print(f"📤 Output Check: {'✅ STABLE' if output_check['is_stable'] else '❌ UNSTABLE'}")
if output_check['issues']:
for issue in output_check['issues']:
print(f" - {issue}")
except Exception as e:
stability_report['overall_status'] = 'CRITICAL'
stability_report['issues_found'].append(f"Network forward pass failed: {str(e)}")
stability_report['recommendations'].append("Check network architecture and input compatibility")
print(f"📤 Output Check: ❌ CRITICAL - Forward pass failed")
# Gradient flow analysis
try:
gradient_analysis = self.analyze_gradient_flow(network, input_tensor, target_output)
print(f"🌊 Gradient Flow: {gradient_analysis['gradient_status']}")
if "exploding" in gradient_analysis['gradient_status'].lower():
stability_report['overall_status'] = 'UNSTABLE'
stability_report['recommendations'].append("Use gradient clipping or reduce learning rate")
elif "vanishing" in gradient_analysis['gradient_status'].lower():
stability_report['overall_status'] = 'UNSTABLE'
stability_report['recommendations'].append("Use ReLU activations or residual connections")
except Exception as e:
print(f"🌊 Gradient Flow: ❌ Analysis failed - {str(e)}")
print(f"\n🎯 OVERALL STATUS: {stability_report['overall_status']}")
if stability_report['recommendations']:
print(f"\n💡 RECOMMENDATIONS:")
for rec in stability_report['recommendations']:
print(f" - {rec}")
return stability_report
def create_unstable_network_demo():
"""
Create networks with known stability issues for demonstration.
This function is PROVIDED to show common stability problems.
Students use it to practice detecting and fixing issues.
"""
print("⚠️ STABILITY ISSUES DEMONSTRATION")
print("=" * 50)
# Create networks with different stability issues
demo_networks = {}
# 1. Network with exploding weights
print("\n1. 🔥 Exploding Weights Network:")
exploding_net = Sequential([
Dense(10, 5),
ReLU(),
Dense(5, 2)
])
# Manually set large weights to simulate training instability
exploding_net.layers[0].weights.data *= 100 # Very large weights
demo_networks['exploding'] = exploding_net
print(" Created network with artificially large weights")
# 2. Network with NaN weights (simulate numerical overflow)
print("\n2. 💀 NaN Weights Network:")
nan_net = Sequential([
Dense(10, 5),
ReLU(),
Dense(5, 2)
])
# Inject NaN values
nan_net.layers[0].weights.data[0, 0] = np.nan
demo_networks['nan'] = nan_net
print(" Created network with NaN values in weights")
# 3. Healthy network for comparison
print("\n3. ✅ Healthy Network:")
healthy_net = Sequential([
Dense(10, 5),
ReLU(),
Dense(5, 2)
])
demo_networks['healthy'] = healthy_net
print(" Created properly initialized network")
return demo_networks
# %% [markdown]
"""
### 🎯 Learning Activity 1: Stability Detection Practice (Medium Guided Implementation)
**Goal**: Learn to detect numerical instability issues that can break neural network training in production.
Complete the missing implementations in the `NetworkStabilityMonitor` class above, then use your monitor to detect stability issues.
"""
# %%
# Initialize the network stability monitor
monitor = NetworkStabilityMonitor()
print("🔧 NETWORK STABILITY MONITORING")
print("=" * 50)
# Create test networks with different stability characteristics
demo_networks = create_unstable_network_demo()
# Create test data
input_data = Tensor(np.random.randn(3, 10)) # Batch of 3 samples
target_data = Tensor(np.random.randn(3, 2)) # Target outputs
print(f"\n🔍 STABILITY ANALYSIS RESULTS:")
print(f"=" * 40)
# Test each network
for network_name, network in demo_networks.items():
print(f"\n📊 Testing {network_name.upper()} Network:")
# Students use their implemented stability checker
stability_report = monitor.comprehensive_stability_check(network, input_data, target_data)
# Show what this means for production
if stability_report['overall_status'] == 'STABLE':
print(f" 🎯 Production Impact: Safe to deploy")
elif stability_report['overall_status'] == 'UNSTABLE':
print(f" ⚠️ Production Impact: May cause training failures")
else:
print(f" 💀 Production Impact: Would crash in production")
print(f"\n💡 STABILITY ENGINEERING INSIGHTS:")
print(f" - NaN values spread through entire network (one bad value ruins everything)")
print(f" - Large weights cause exponential growth through layers")
print(f" - Stability monitoring prevents silent training failures")
print(f" - Early detection saves compute resources and time")
# %% [markdown]
"""
### 🎯 Learning Activity 2: Production Stability Patterns (Review & Understand)
**Goal**: Understand common stability issues in production ML systems and learn industry best practices for preventing them.
"""
# %%
print("🏭 PRODUCTION STABILITY PATTERNS")
print("=" * 50)
# Test different input scenarios that cause instability
print("\n🔍 Input Data Stability Scenarios:")
stability_scenarios = [
("Normal Data", np.random.randn(5, 10)),
("Large Values", np.random.randn(5, 10) * 1000),
("Extreme Values", np.random.randn(5, 10) * 1e8),
("Mixed with NaN", np.random.randn(5, 10)),
("All Zeros", np.zeros((5, 10))),
("All Ones", np.ones((5, 10)) * 1e6)
]
# Inject NaN in mixed scenario
stability_scenarios[3] = ("Mixed with NaN", np.random.randn(5, 10))
scenario_data = stability_scenarios[3][1].copy()
scenario_data[0, 0] = np.nan
stability_scenarios[3] = ("Mixed with NaN", scenario_data)
# Test each scenario
healthy_network = demo_networks['healthy']
for scenario_name, test_data in stability_scenarios:
print(f"\n📊 {scenario_name}:")
try:
input_tensor = Tensor(test_data)
input_check = monitor.check_tensor_stability(input_tensor, scenario_name)
print(f" Input Status: {'✅ STABLE' if input_check['is_stable'] else '❌ UNSTABLE'}")
if input_check['issues']:
print(f" Issues: {', '.join(input_check['issues'])}")
# Try network forward pass
try:
output = healthy_network(input_tensor)
output_check = monitor.check_tensor_stability(output, f"{scenario_name}_output")
print(f" Output Status: {'✅ STABLE' if output_check['is_stable'] else '❌ UNSTABLE'}")
if output_check['issues']:
print(f" Output Issues: {', '.join(output_check['issues'])}")
except Exception as e:
print(f" ❌ Forward pass failed: {str(e)}")
except Exception as e:
print(f" ❌ Could not create tensor: {str(e)}")
print(f"\n🎯 PRODUCTION STABILITY LESSONS:")
print(f"=" * 40)
print(f"\n1. 🛡️ INPUT VALIDATION:")
print(f" - Always validate input data before processing")
print(f" - Clip extreme values to reasonable ranges")
print(f" - Check for NaN/inf values in data pipelines")
print(f"\n2. 🔧 MONITORING STRATEGY:")
print(f" - Monitor weight magnitudes during training")
print(f" - Track gradient norms to detect vanishing/exploding")
print(f" - Log activation statistics to catch distribution shift")
print(f"\n3. 🚨 EARLY WARNING SYSTEM:")
print(f" - Set thresholds for weight magnitudes")
print(f" - Alert when gradients become too large/small")
print(f" - Automatically stop training on stability issues")
print(f"\n4. 🛠️ PREVENTIVE MEASURES:")
print(f" - Proper weight initialization (Xavier/He)")
print(f" - Gradient clipping for exploding gradients")
print(f" - Batch normalization for internal stability")
print(f" - Learning rate scheduling to prevent instability")
print(f"\n💡 SYSTEMS ENGINEERING INSIGHT:")
print(f"Stability monitoring is like production health checks:")
print(f"- Prevent silent failures that waste compute resources")
print(f"- Enable automatic recovery strategies (restart training)")
print(f"- Provide debugging information for model developers")
print(f"- Critical for unattended training jobs in production")
# %% [markdown]
"""
## 🎯 MODULE SUMMARY: Neural Network Architectures
Congratulations! You've successfully implemented complete neural network architectures:
### What You've Accomplished
✅ **Sequential Networks**: Chained layers for complex transformations
✅ **MLP Creation**: Multi-layer perceptrons with flexible architectures
✅ **Network Architectures**: Different activation patterns and output types
✅ **Integration**: Real-world applications like classification and regression
### Key Concepts You've Learned
- **Sequential Processing**: How layers chain together for complex functions
- **MLP Design**: Multi-layer perceptrons as universal function approximators
- **Architecture Choices**: How depth, width, and activations affect learning
- **Real Applications**: Classification, regression, and feature extraction
### Next Steps
1. **Export your code**: `tito package nbdev --export 04_networks`
2. **Test your implementation**: `tito test 04_networks`
3. **Build complete models**: Combine with training for full ML pipelines
4. **Move to Module 5**: Add convolutional layers for image processing!
**Ready for CNNs?** Your network foundations are now ready for specialized architectures!
"""