I want confusion matrix for this training process

for epoch in range(num_epochs):
for i ,(images, labels) in enumerate(train_data):
# Run the forward pass
images = images.unsqueeze(1).float()
outputs = model(images)
loss = criterion(outputs, labels)
loss_list.append(loss.item())

    # Backprop and perform Adam optimisation
    optimizer.zero_grad()
    loss.backward()
    optimizer.step()
#optimizer update the weight parameters to minimize the loss function

    # Track the accuracy
    total = labels.size(0)
    _, predicted = torch.max(outputs.data, 1)
    correct = (predicted == labels).sum().item()
    acc_list.append(correct / total)

    if (i + 1) % 35 == 0:
        print('Epoch [{}/{}], Step [{}/{}], Loss: {:.4f}, Accuracy: {:.2f}%'
              .format(epoch + 1, num_epochs, i + 1, total_step, loss.item(),
                      (correct / total) * 100))

Please refer

https://github.com/tgi25/pytorch/blob/master/Lesson10%20-%20LeNet.ipynb - [4 Predictions]