From Scratch to PyTorch

Hallo. I’m just starting with PyTorch.

I have built the following program in python and I want to translate it to PyTorch so I can compare the result but I still haven’t found how to do it. I’m currently attenting a course for PyTorch but it’s really slow and I need some help. Can anyone tell me how to write this code using PyTorch?
It’s the code in the paper Deep Learning: An Introduction for Applied
Mathematicians - Catherine F. Higham, Desmond J. Higham but I have translated it from MATLAB to Python

P.S. Thank you everyone in advance for your time and sorry if it’s against the rules of the Forum. Please delete it if it is.
Also if you have any suggestion about Courses or YouTube videos/series that you found helpful while you were getting starting with PyTorch please let me know.

import numpy as np
from numpy.linalg import norm as norm

def activate(x,W,b): #Evaluates sigmoid function. x is the input vector, y is the output vector W contains the weights, b contains the shifts. The i-th component of y is activate((Wx+b)_i) where activate(z) = 1/(1+exp(-z))
    z=W@x+breturn 1/(1+np.exp(-z))

def cost(W2,W3,W4,b2,b3,b4):
	costvec = np.zeros(10)
	for i in range(10):
		costvec[i] = norm(y[:,i]-a4,2)
	costval = norm(costvec,2)**2
	return costval

x1 = np.array([0.1,0.3,0.1,0.6,0.4,0.6,0.5,0.9,0.4,0.7])
x2 = np.array([0.1,0.4,0.5,0.9,0.2,0.3,0.6,0.2,0.4,0.6])
y=np.array([ [1,1,1,1,1,0,0,0,0,0],[0,0,0,0,0,1,1,1,1,1] ])

#forward and back propagate
eta = 0.05 #learning rate
Niter = 10**6 #number of SG iterations
savecost = np.zeros(Niter) #Value of cost at each iteration
for counter in range(Niter):
	k=np.random.randint(0,10) #choose a training point at random
	#forward pass
	a2 = activate(x,W2,b2)
	a3 = activate(a2,W3,b3)
	a4 = activate(a3,W4,b4)
	#backward pass
	delta4 = a4*(1-a4)*(a4-y[:,k])
	delta3 = a3*(1-a3)*(W4.T@delta4)
	delta2 = a2*(1-a2)*(W3.T@delta3)
	#gradient step
	W2 = W2 - eta*delta2@x
	W3 = W3 - eta*delta3.reshape(3,1)@a2.reshape(1,2)
	W4 = W4 - eta*delta4.reshape(2,1)@a3.reshape(1,3)
	b2 = b2 - eta*delta2
	b3 = b3 - eta*delta3
	b4 = b4 - eta*delta4
	newcost = cost(W2,W3,W4,b2,b3,b4)
	print(newcost , 'i=',counter)

import matplotlib.pyplot as plt
iterr=[i for i in range(counter+1)]

I think the best approach would be if you could explain what you have done so far, what is not working, and where you are stuck.
Once we know this information we might be able to guide you through it. Especially if you are taking a course I would recommend to try to debug it yourself first and avoid looking at perfect solutions. :wink: