import torch
class myrelu(torch.autograd.Function):#自定義子類
    # 通過建立torch.autograd的子類來實(shí)現(xiàn)自定義的autograd函數(shù)，并完成張量的正向和反向傳播
    @staticmethod
    def forward(ctx, x):
        # 在正向傳播中，接受到一個(gè)上下文對(duì)象和一個(gè)包含輸入的張量，必須返回一個(gè)包含輸出的張量，可以使用上下文對(duì)象來緩存對(duì)象，以便在反向傳播中使用
        ctx.save_for_backward(x)
        return x.clamp(min=0)
    @staticmethod
    def backward(ctx, grad_output):
        """
        在反向傳播中，我們接收到上下文對(duì)象和一個(gè)張量，
        其包含了相對(duì)于正向傳播過程中產(chǎn)生的輸出的損失的梯度。
        我們可以從上下文對(duì)象中檢索緩存的數(shù)據(jù)，
        并且必須計(jì)算并返回與正向傳播的輸入相關(guān)的損失的梯度。
        """
        x, = ctx.saved_tensors
        grad_x = grad_output.clone()
        grad_x[x < 0] = 0
        return grad_x

#%%
device=torch.device('cuda' if torch.cuda.is_available() else 'cpu')
# n是批量大小，d_in是輸入維度
# h是隱藏的維度，d_out是輸出維度
n,d_in,h,d_out=64,1000,100,10
# 創(chuàng)建隨機(jī)輸入和輸出數(shù)據(jù),requires_grad默認(rèn)設(shè)置為False,表示不需要后期微分操作
x=torch.randn(n,d_in,device=device)
y=torch.randn(n,d_out,device=device)
# 隨機(jī)初始化權(quán)重，requires_grad默認(rèn)設(shè)置為True,表示想要計(jì)算其微分
w1=torch.randn(d_in,h,device=device,requires_grad=True)
w2=torch.randn(h,d_out,device=device,requires_grad=True)

learning_rate=1e-6
for i in range(500):
    #前向傳播，使用tensor上的操作計(jì)算預(yù)測值y
  #調(diào)用自定義的myrelu.apply函數(shù)
    y_pred=myrelu.apply(x.mm(w1)).mm(w2)
       
    
    #使用tensor中的操作計(jì)算損失值，loss.item()得到loss這個(gè)張量對(duì)應(yīng)的數(shù)值
    loss=(y_pred-y).pow(2).sum()
    print(i,loss.item())
    
    #使用autograd計(jì)算反向傳播，這個(gè)調(diào)用將計(jì)算loss對(duì)所有的requires_grad=True的tensor梯度，
    #調(diào)用之后，w1.grad和w2.grad將分別是loss對(duì)w1和w2的梯度張量
    loss.backward()
    #使用梯度下降更新權(quán)重，只想對(duì)w1和w2的值進(jìn)行原地改變：不想更新構(gòu)建計(jì)算圖
    #所以使用torch.no_grad()阻止pytorch更新構(gòu)建計(jì)算圖
    with torch.no_grad():
        w1-=learning_rate*w1.grad
        w2-=learning_rate*w2.grad
        #反向傳播后手動(dòng)將梯度置零
        w1.grad.zero_()
        w2.grad.zero_() 

#%%使用TensorFlow
import tensorflow.compat.v1 as tf   #為了用placeholder不惜一切代價(jià)
tf.disable_v2_behavior()
import numpy as np
#%%
# 建立計(jì)算圖

# n是批量大小，d_in是輸入維度
# h是隱藏的維度，d_out是輸出維度
n,d_in,h,d_out=64,1000,100,10

# 為輸入和目標(biāo)數(shù)據(jù)創(chuàng)建placeholder，在執(zhí)行計(jì)算圖時(shí)，他們將會(huì)被真實(shí)的數(shù)據(jù)填充
x=tf.placeholder(tf.float32,shape=(None,d_in))
y=tf.placeholder(tf.float32,shape=(None,d_out))

# 為權(quán)重創(chuàng)建variable并用隨機(jī)數(shù)據(jù)初始化,TensorFlow的variable在執(zhí)行計(jì)算圖時(shí)不會(huì)改變
w1 = tf.Variable(tf.random_normal((d_in,h)))
w2=tf.Variable(tf.random_normal((h,d_out)))
# 前向傳播：使用TensorFlow的張量運(yùn)算計(jì)算預(yù)測值y（這段代碼不執(zhí)行任何數(shù)值運(yùn)算，只是建立了稍后要執(zhí)行的計(jì)算圖）
h=tf.matmul(x,w1)
h_relu=tf.maximum(h,tf.zeros(1))
y_pred=tf.matmul(h_relu,w2)
# 使用TensorFlow的張量運(yùn)算損失loss
loss=tf.reduce_sum((y-y_pred)**2.0)
# 計(jì)算loss對(duì)于權(quán)重w1和w2的導(dǎo)數(shù)
grad_w1,grad_w2=tf.gradients(loss,[w1,w2])
# 使用梯度下降更新權(quán)重，為了實(shí)際更新權(quán)重，我們需要在執(zhí)行計(jì)算圖時(shí)計(jì)算new_w1和new_w2
# 注：在TensorFlow中，更新權(quán)重值得行為是計(jì)算圖的一部分，但在Pytorch中發(fā)生在計(jì)算圖形之外
learning_rate=1e-6
new_w1=w1.assign(w1-learning_rate*grad_w1)
new_w2=w2.assign(w2-learning_rate*grad_w2)

# 現(xiàn)在搭建好了計(jì)算圖，開始一個(gè)TensorFlow回話來執(zhí)行計(jì)算圖
with tf.Session() as sess:
    # 運(yùn)算一次計(jì)算圖來出事話variable w1和w2
    sess.run(tf.global_variables_initializer())
    # 創(chuàng)建numpy數(shù)組存儲(chǔ)輸入x和目標(biāo)y的實(shí)際數(shù)據(jù)
    x_value=np.random.randn(n,d_in)
    y_value=np.random.randn(n,d_out)
    for i in range(500):
        # 多次運(yùn)行計(jì)算圖，每次執(zhí)行時(shí)，都有feed_dict參數(shù)，
        # 將x_value綁定到x，將y_value綁定到y(tǒng).每次執(zhí)行計(jì)算圖都要計(jì)算損失，
        # new_w1和new_w2，這些張量的值以numpy數(shù)組的形式返回
        loss_value,i,i=sess.run([loss,new_w1,new_w2],
                                feed_dict={x:x_value,y:y_value})
        print(loss_value)