回覆列表
  • 1 # 藍風24

    如果已經安裝了IMSL,在你編寫的程式中使用use 語句呼叫你要用的函式所在的庫就可以了.

    下面是抄了的一個例子:

    !lin_sol_gen_int 為IMSL mod檔案,在include目錄下

    use lin_sol_gen_int

    use rand_gen_int

    use error_option_packet

    implicit none

    ! This is Example 1 for LIN_SOL_GEN.

    integer, parameter :: n=32

    real(kind(1e0)), parameter :: one=1e0

    real(kind(1e0)) err

    real(kind(1e0)) A(n,n), b(n,n), x(n,n), res(n,n), y(n**2)

    ! Generate a random matrix.

    call rand_gen(y)

    A = reshape(y,(/n,n/))

    ! Generate random right-hand sides.

    call rand_gen(y)

    b = reshape(y,(/n,n/))

    ! Compute the solution matrix of Ax=b.

    call lin_sol_gen(A, b, x)

    ! Check the results for small residuals.

    res = b - matmul(A,x)

    err = maxval(abs(res))/sum(abs(A)+abs(b))

    if (err

    write (*,*) "Example 1 for LIN_SOL_GEN is correct."

    end if

    end

  • 中秋節和大豐收的關聯?
  • 半導體行業會崛起嗎?