fortran 90 将切割矩阵传递给函数并且出现额外值

Question

我正在尝试将一些函数从原始 lapack 转换为 GO gonum 库，这样做很常见

         a(1,1)=199
         a(1,2)=2
         a(1,3)=3
         a(1,4)=4
         
         a(2,1)=5
         a(2,2)=6
         a(2,3)=7
         a(2,4)=8
         
         a(3,1)=9
         a(3,2)=10
         a(3,3)=11
         a(3,4)=12
         
         a(4,1)=13
         a(4,2)=14
         a(4,3)=15
         a(4,4)=16
         
         Call nagf_file_print_matrix_real_gen('General', ' ', 4, 4, a, 4, 'a pre matrix ', ifail)
         
         CALL dlarf2( 'Left', m-i+1, n-i, a( i, i ), 1, tau( i ),a( i, i+1 ), lda, work )

然后在 dlarf2 我做

SUBROUTINE dlarf2( SIDE, M, N, V, INCV, TAU, C, LDC, WORK )
   CHARACTER          SIDE
   INTEGER            INCV, LDC, M, N
   DOUBLE PRECISION   TAU
   DOUBLE PRECISION   C( LDC, * ), V( * ), WORK( * )
   DOUBLE PRECISION   ONE, ZERO
   parameter( one = 1.0d+0, zero = 0.0d+0 )
   LOGICAL            APPLYLEFT
   INTEGER            I, LASTV, LASTC
   EXTERNAL           dgemv, dger
   LOGICAL            LSAME
   INTEGER            ILADLR, ILADLC
   EXTERNAL           lsame, iladlr, iladlc
   Call nagf_file_print_matrix_real_gen('General', ' ', 4, 4, c, 4, 'c pre iladc ', ifail)

我得到了这个

其中矩阵 a 是 4x4 的双精度矩阵，为什么会出现 1？