版块导航
正在加载中...
客户端APP下载
登录
注册
帖子
帖子
用户
本版
应《网络安全法》要求,自2017年10月1日起,未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用,请尽快对帐号进行手机号验证,感谢您的理解与支持!
24小时热门版块排行榜
>
休闲灌水
(1279)
>
论坛更新日志
(578)
>
基金申请
(250)
>
有奖起名
(162)
>
虫友互识
(130)
>
健康生活
(87)
>
论文投稿
(37)
>
留学DIY
(32)
>
导师招生
(29)
>
硕博家园
(27)
>
考博
(25)
>
文献求助
(20)
>
考研
(11)
>
教师之家
(10)
>
博后之家
(8)
>
公派出国
(8)
小木虫论坛-学术科研互动平台
»
计算模拟区
»
程序语言
»
Fortran
»
【求助】此FORTRAN90程序没有错误,但是运行不出来。
5
1/1
返回列表
查看: 1216 | 回复: 4
只看楼主
@他人
存档
新回复提醒
(忽略)
收藏
在APP中查看
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖
zyj8119
木虫
(著名写手)
应助: 65
(初中生)
贵宾: 0.003
金币: 915.1
散金: 1440
红花: 35
帖子: 2936
在线: 1329.4小时
虫号: 664177
注册: 2008-11-29
性别: GG
专业: 理论和计算化学
[交流]
【求助】此FORTRAN90程序没有错误,但是运行不出来。
CODE:
MODULE convexity
TYPE point
COMPLEX::z
END TYPE point
CONTAINS
FUNCTION orientation(p,i) RESULT(orient)
TYPE(point),DIMENSION(:),INTENT(in)::p
INTEGER::i,n
REAL::orient,f
COMPLEX::z1,z2,z3
f(z1,z2,z3)=REAL(z2-z1)*AIMAG(z3-z2)-REAL(z3-z2)*AIMAG(z2-z1)
n=SIZE(p)
IF(i==n-1)THEN
orient=f(p(n-1).z,p(n).z,p(1).z)
ELSE IF(i==n)THEN
orient=f(p(n).z,p(1).z,p(2).z)
ELSE
orient=f(p(i).z,p(i+1).z,p(i+2).z)
END IF
END FUNCTION orientation
!
SUBROUTINE convex_polygon(p,convex)
LOGICAL::convex
TYPE(point),DIMENSION(:),INTENT(in)::p
INTEGER::i,n,anticlockwise
convex=.TRUE.;n=SIZE(P)
IF(orientation(p,1)>0)THEN
anticlockwise=1
ELSE
anticlocakwise=-1
END IF
DO i=2,n
IF(anticlockwise*orientation(p,i)<0)THEN
convex=.FALSE.;EXIT
END IF
END DO
END SUBROUTINE convex_polygon
END MODULE convexity
!
PROGRAM test
USE convexity
INTEGER,parameter::n=6
TYPE(point),DIMENSION(1:n)::polygon
LOGICAL::convex
INTEGER::i
OPEN(1,FILE='convex')
DO I=1,2
READ(1,*)polygon
CALL convex_polygon(polygon,convex)
IF(convex==.TRUE.)PRINT*,'polygon is convex'
IF(convex/=.TRUE.)PRINT*,'polygon is not convex'
END DO
END PROGRAM test
回复此楼
» 猜你喜欢
售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急
已经有4人回复
售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急
已经有4人回复
售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急
已经有4人回复
售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗?
已经有20人回复
今天务委会开完了,明天出结果吗
已经有22人回复
在坚冰还盖着北海的时候,我看到了怒放的梅花。
已经有9人回复
明天应该可查了!?
已经有6人回复
有没有大神帮我看看基金代码
已经有10人回复
放榜前的不淡定
已经有19人回复
高级回复
好好学习,天天向上。
1楼
2010-08-29 09:54:02
已阅
回复此楼
关注TA
给TA发消息
送TA红花
TA的回帖
janecool
铁杆木虫
(著名写手)
应助: 54
(初中生)
贵宾: 0.006
金币: 5996.2
散金: 342
红花: 12
帖子: 2070
在线: 626.6小时
虫号: 830148
注册: 2009-08-18
性别: GG
专业: 粒子物理学和场论
★
zyj8119(金币+4):先尝试下看看。 2010-08-29 15:16:16
resonant(金币+1):感谢参与应助:-) 2010-08-29 19:54:03
引用回帖:
Originally posted by
zyj8119
at 2010-08-29 10:09:13:
什么意思?
你可以先在do循环前给N赋予一个数值,看程序是否能运行,如果还不能,我现在也不知道是为什么。
赞
一下
(1人)
回复此楼
高级回复
相信会有属于你的精彩
4楼
2010-08-29 10:15:29
已阅
回复此楼
关注TA
给TA发消息
送TA红花
TA的回帖
查看全部 5 个回答
janecool
铁杆木虫
(著名写手)
应助: 54
(初中生)
贵宾: 0.006
金币: 5996.2
散金: 342
红花: 12
帖子: 2070
在线: 626.6小时
虫号: 830148
注册: 2009-08-18
性别: GG
专业: 粒子物理学和场论
★
zyj8119(金币
+1
):谢谢参与
问题可能出自第一个DO循环里面,里面的N应该是一个已知量
赞
一下
(1人)
回复此楼
相信会有属于你的精彩
2楼
2010-08-29 10:03:47
已阅
回复此楼
关注TA
给TA发消息
送TA红花
TA的回帖
zyj8119
木虫
(著名写手)
应助: 65
(初中生)
贵宾: 0.003
金币: 915.1
散金: 1440
红花: 35
帖子: 2936
在线: 1329.4小时
虫号: 664177
注册: 2008-11-29
性别: GG
专业: 理论和计算化学
引用回帖:
Originally posted by
janecool
at 2010-08-29 10:03:47:
问题可能出自第一个DO循环里面,里面的N应该是一个已知量
什么意思?
回复此楼
好好学习,天天向上。
3楼
2010-08-29 10:09:13
已阅
回复此楼
关注TA
给TA发消息
送TA红花
TA的回帖
ifmc1234
木虫
(著名写手)
应助: 18
(小学生)
金币: 2797.9
散金: 6001
红花: 14
帖子: 1761
在线: 572.6小时
虫号: 825500
注册: 2009-08-10
专业: 理论和计算化学
★ ★
zyj8119(金币
+1
):谢谢参与
nono2009(金币+1):鼓励应助。 2010-09-09 08:02:23
引用回帖:
Originally posted by
zyj8119
at 2010-08-29 10:09:13:
什么意思?
你没把n定义为全局变量。所以在子模块里n值是0
赞
一下
(2人)
回复此楼
5楼
2010-09-08 18:48:11
已阅
回复此楼
关注TA
给TA发消息
送TA红花
TA的回帖
查看全部 5 个回答
如果回帖内容含有宣传信息,请如实选中。否则帐号将被全论坛禁言
普通表情
龙
兔
虎
猫
高级回复
(可上传附件)
百度网盘
|
360云盘
|
千易网盘
|
华为网盘
在新窗口页面中打开自己喜欢的网盘网站,将文件上传后,然后将下载链接复制到帖子内容中就可以了。
最具人气热帖推荐
[查看全部]
作者
回/看
最后发表
[
考博
]
售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急
+3
Qj98d49dHFb9
2026-08-25
4/200
2026-08-26 03:08
by
cNXvBfCpiZOM
[
考研
]
售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急
+3
7K1CJE38xLG4
2026-08-25
4/200
2026-08-26 02:04
by
cNXvBfCpiZOM
[
基金申请
]
2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗?
+17
木水思豆
2026-08-25
20/1000
2026-08-26 01:26
by
chengyan1220
[
基金申请
]
今天务委会开完了,明天出结果吗
+18
angus9576
2026-08-25
22/1100
2026-08-26 00:41
by
merchancy
[
基金申请
]
在坚冰还盖着北海的时候,我看到了怒放的梅花。
(金币+10)
+6
ziyangfang
2026-08-25
9/450
2026-08-25 20:26
by
huagongfeihu
[
基金申请
]
有没有大神帮我看看基金代码
28
+4
1234567wang
2026-08-24
10/500
2026-08-25 19:15
by
lfy8008
[
基金申请
]
放榜前的不淡定
20
+4
snowwithsea
2026-08-19
19/950
2026-08-25 17:57
by
chick875
[
基金申请
]
某些机构,以效率低为荣,以效率低作为存在感
+9
yuleib84
2026-08-25
10/500
2026-08-25 17:14
by
alexon
[
基金申请
]
今日不放榜?网传国自然预计 8 月 27 日可查结果
+17
医学老男孩
2026-08-20
22/1100
2026-08-25 15:36
by
医学老男孩
[
基金申请
]
如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》
+8
scalable
2026-08-24
8/400
2026-08-25 12:52
by
jnhyjjm
[
基金申请
]
没有任何消息-是不是就凉了
+9
图啦图啦
2026-08-24
10/500
2026-08-25 11:59
by
南海小哥
[
基金申请
]
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
+5
process2012
2026-08-23
7/350
2026-08-25 09:42
by
huixian257
[
基金申请
]
今天基金会出结果吗?20260819
+17
kkkl_v
2026-08-19
18/900
2026-08-25 09:41
by
windflowerwy
[
基金申请
]
范进中举一文的中心思想
+7
炎黄贵胄
2026-08-22
8/400
2026-08-25 08:48
by
ZJTJZ
[
基金申请
]
建议基金发布提前给出明确的时间点
+13
kulium
2026-08-21
16/800
2026-08-24 16:27
by
superceng
[
基金申请
]
让我中一个面上吧!
+13
大萍1987
2026-08-20
16/800
2026-08-24 10:23
by
太傻了
[
基金申请
]
时间戳今天,20号变了
+5
archvillain
2026-08-20
5/250
2026-08-22 06:12
by
hui_daxiao
[
基金申请
]
时间戳又变了
+13
wuchongjun
2026-08-20
19/950
2026-08-21 17:21
by
紫杉醇
[
基金申请
]
今天放榜没戏了吧
+9
yuleib84
2026-08-19
11/550
2026-08-21 10:06
by
gltch
[
基金申请
]
基金啊基金
+4
longfie172
2026-08-20
4/200
2026-08-21 08:58
by
mark mao
信息提示
关闭
请填处理意见
关闭
确定