Fortran Coder

查看: 9327|回复: 4
打印 上一主题 下一主题

[数值问题] 运行时出错,access violation.

[复制链接]

15

帖子

5

主题

0

精华

熟手

F 币
129 元
贡献
65 点

规矩勋章

跳转到指定楼层
楼主
发表于 2020-11-18 00:33:02 | 只看该作者 |只看大图 回帖奖励 |倒序浏览 |阅读模式
real function F(x,y,m,n)
    implicit none
    real::x,y,pi
    integer :: m,n
    PI=4.d0*ATAN(1.d0)
    F = exp((0.d0,-1.d0)*(2.d0*pi*x/m+2.d0*pi*y/n))
    end function

    program test
    include 'link_fnl_shared.h'
    implicit none
    real:: F


    real, parameter :: ERRABS = 0.001
    real, parameter :: ERRREL = 0.001
    integer, parameter :: IRULE = 6

    real:: G,H
    real :: ANS
    real :: ERREST
    integer :: i,j
    real::x1,y1
    real :: A,B,k
    A = -0.5
    B = 0.5

    do i=1,10
    do j=1,10
    k=F(x1,y1,i,j)
    call TWODQ (k, A, B, G, H, ERRABS, ERRREL, IRULE, ANS, ERREST)
    write(*,*) ANS
    end do
    end do
    end program test

    real function G(X)
    implicit none
    real X
    G = -0.5
    return
    end function

    real function H(X)
    implicit none
    real X
    H = 0.5
    return
    end function

目的是为了求e^-i(2pi*x/n+2pi*y/m)的二重积分,m和n是变的,但是运行会报错

分享到:  微信微信
收藏收藏 点赞点赞 点踩点踩

15

帖子

5

主题

0

精华

熟手

F 币
129 元
贡献
65 点

规矩勋章

沙发
 楼主| 发表于 2020-11-18 00:34:47 | 只看该作者
错误如图所示

XT{G4HJ_R9L`]3L`)MO$JC7.png (14.11 KB, 下载次数: 248)

XT{G4HJ_R9L`]3L`)MO$JC7.png

213

帖子

2

主题

0

精华

宗师

F 币
2132 元
贡献
875 点

规矩勋章

板凳
发表于 2020-11-18 04:53:04 | 只看该作者
x1, y1在调用 F 前没有赋值。

15

帖子

5

主题

0

精华

熟手

F 币
129 元
贡献
65 点

规矩勋章

地板
 楼主| 发表于 2020-11-18 09:35:58 | 只看该作者
风平老涡 发表于 2020-11-18 04:53
x1, y1在调用 F 前没有赋值。

能具体讲一下么

213

帖子

2

主题

0

精华

宗师

F 币
2132 元
贡献
875 点

规矩勋章

5#
发表于 2020-11-19 01:15:53 | 只看该作者
Botton 发表于 2020-11-18 09:35
能具体讲一下么

k=F(x1,y1,i,j) 在调用这个函数时,i和J都有赋值,X1和Y1没有赋值。
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

捐赠本站|Archiver|关于我们 About Us|小黑屋|Fcode ( 京ICP备18005632-2号 )

GMT+8, 2024-11-24 02:32

Powered by Tencent X3.4

© 2013-2024 Tencent

快速回复 返回顶部 返回列表