FIX: re-initiate R/T in every layer and add delay outside the macro "if" (#161)

Author: Dengda98

pygrt/C_extension/src/common/kernel_template.c_

@@ -75,11 +75,11 @@ void __GRT_MATRIX_FUNC__(GRT_MODEL1D *mod1d, const real_t k)
     // 自由表面
     RT_MATRIX *M_top = &mod1d->M_top;
 
-    // 定义物理层内的反射透射系数矩阵
-    grt_init_RT_matrix(M);
-
     // 从顶到底进行矩阵递推, 公式(5.5.3)
     for(size_t iy=1; iy<mod1d->n; ++iy){ // 因为n>=3, 故一定会进入该循环
+        
+        // 定义物理层内的反射透射系数矩阵
+        grt_init_RT_matrix(M);
 
         // 只有动态解才可以使用这个 if 简化，
         // 对于静态解，即使是震源层、接收层这种虚拟层位也需要显式计算R/T矩阵
@@ -90,13 +90,13 @@ void __GRT_MATRIX_FUNC__(GRT_MODEL1D *mod1d, const real_t k)
             __grt_RT_matrix_PSV(mod1d, iy, M);
             __grt_RT_matrix_SH(mod1d, iy, M);
             if(M->stats==GRT_INVERSE_FAILURE)  goto BEFORE_RETURN;
-
-            // 加入时间延迟因子(第iy-1界面与第iy界面之间)
-            __grt_delay_RT_matrix(mod1d, iy, M);
         #ifdef __DYNAMIC_KERNEL__
         }
         #endif
 
+        // 加入时间延迟因子(第iy-1界面与第iy界面之间)
+        __grt_delay_RT_matrix(mod1d, iy, M);
+
         // FA
         if(iy <= imin){ 
             if(iy == 1){ // 初始化FA