Fortran Coder

标题: 矩阵算法合集问题续II [打印本页]

作者: PotsyYZhou 时间: 2022-9-30 15:12
标题: 矩阵算法合集问题续II
按提示，变换了一下格式如下：

[Fortran] 纯文本查看 复制代码

subroutine eulerian_inverse ( n, a )

  do j = 1, n

    do i = 1, n

      if ( i .lt. j ) then

        a(i,j) = 0.0D+00

      else if ( i .eq. j ) then

        a(i,j) = 1.0D+00

      else if ( j .lt. i ) then

        dot = 0.0D+00

        do k = j, i - 1

          dot = dot + b(i,k) * a(k,j)

        end do

        a(i,j) = - dot

      end if

    end do

  end do

  return

end

   还是不行，报错如下：
main.f90:14269:8:14269 |       a(i, j) = 0.0D+00    |       1Error: Unclassifiable statement at (1)
main.f90:14281:8:14281 |       a(i, j) = -dot    |       1Error: The function result on the lhs of the assignment at (1) must have the pointer attribute.
？？？！

作者: li913 时间: 2022-10-2 09:39
a需要定义为数组。

欢迎光临 Fortran Coder (http://bbs.fcode.cn/)