Optimizing regression weights for NN outputs with PyTorch

So I'm basically trying to fit a regression on the relation of the input and output of a neural network model. Then the idea is, that these estimated regression weights should be optimized to some specific target value (let's say matrix of ones). So the weights are optimized, but have a direct relation to the neural network weights. For a simple problem, I can find a solution on paper, but I fail to implement this. The loss is not decreasing, but oscillating.

# Some arbitrary PyTorch model
my_nn = Net()

# Our target regression weights
target = torch.ones(2, 2)

# Optimizer for NN weights

# Arbitrary loss function
criterion = nn.MSELoss()

# Training loop
my_nn.train()
for epoch in range(500):
print('Running Epoch: ', epoch)

# output of model is 2-dim
y_pred = my_nn(X)

# weights for regression model
# optimizer for regression model
optimizer1 = optim.SGD([w1], lr=0.001, momentum=0.9)
for nest_epoch in range(1000):

y_hat = torch.mm(X, w1)
loss = criterion(y_hat, y_pred)

# Need to add retain_graph, as it doesn't work otherwise
# but it does work, this loss will get close to 0
loss.backward(retain_graph=True)
optimizer1.step()

# this loss does not decrease
loss = criterion(target, w1)
print('Loss:', loss.item())
loss.backward()
optimizer.step()
$$`$$

My solution is to update the variable manually and creating a new graph for autograd. Works perfectly well now!

my_nn.train()
for epoch in range(300):

Y = my_nn(X)
loss_output = criterion(Y, y)
w = torch.zeros_like(target)