Fortran Coder

查看: 5910|回复: 3
打印 上一主题 下一主题

[数值问题] Error: Syntax error in array constructor

[复制链接]

12

帖子

5

主题

0

精华

入门

F 币
84 元
贡献
47 点
跳转到指定楼层
楼主
发表于 2020-7-9 09:03:20 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式
请问我这样写的数组的code哪里出问题,是(/ /)不能有两段吗?
[Fortran] 纯文本查看 复制代码
    !! metric_p(:) = h(:)
    m_p_k(:, i, j, k) = ( (/ & 
    ( (1._dl - 2._dl * xi) * f_x(1, i, j, k) ** 2 - 2._dl * xi * fields_fk(1, i, j, k) * I44_D2F_XX_FLD(fields_fk, 1, i, j, k) + f_x(2, i, j, k) ** 2 ) , &
    ( (1._dl - 2._dl * xi) * f_y(1, i, j, k) ** 2 - 2._dl * xi * fields_fk(1, i, j, k) * I44_D2F_YY_FLD(fields_fk, 1, i, j, k) + f_y(2, i, j, k) ** 2 ) , &
    ( (1._dl - 2._dl * xi) * f_z(1, i, j, k) ** 2 - 2._dl * xi * fields_fk(1, i, j, k) * I44_D2F_ZZ_FLD(fields_fk, 1, i, j, k) + f_z(2, i, j, k) ** 2 ) , &
    ( (1._dl - 2._dl * xi) * f_y(1, i, j, k) * f_z(1, i, j, k) - 2._dl * xi * fields_fk(1, i, j, k) * I44_D2F_YZ_FLD(fields_fk, 1, i, j, k) + f_y(2, i, j, k) * f_z(2, i, j, k) ) , &
    ( (1._dl - 2._dl * xi) * f_x(1, i, j, k) * f_z(1, i, j, k) - 2._dl * xi * fields_fk(1, i, j, k) * I44_D2F_XZ_FLD(fields_fk, 1, i, j, k) + f_x(2, i, j, k) * f_z(2, i, j, k) ) , &
    ( (1._dl - 2._dl * xi) * f_x(1, i, j, k) * f_y(1, i, j, k) - 2._dl * xi * fields_fk(1, i, j, k) * I44_D2F_XY_FLD(fields_fk, 1, i, j, k) + f_x(2, i, j, k) * f_y(2, i, j, k) ) , &
    /) &
    / (metric%dx) ** 2 / 144._dl * (2._dl / Mplsq) / (1 - xi * fields_fk(1, i, j, k) ** 2 / Mplsq) &
    -(/ &
    (2._dl * h12_xy(i, j, k) + 2._dl * h13_xz(i, j, k) - h11_yy(i, j, k) - h11_zz(i, j, k) - h22_xx(i, j, k) - h33_xx(i, j, k) ) , &
    (2._dl * h12_xy(i, j, k) + 2._dl * h23_yz(i, j, k) - h22_xx(i, j, k) - h22_zz(i, j, k) - h11_yy(i, j, k) - h33_yy(i, j, k)) , &
    (2._dl * h13_xz(i, j, k) + 2._dl * h23_yz(i, j, k) - h33_xx(i, j, k) - h33_yy(i, j, k) - h11_zz(i, j, k) - h22_zz(i, j, k)) , &
    (h12_xz(i, j, k) + h13_xy(i, j, k) - h23_xx(i, j, k) - h11_yz(i, j, k)) , &
    (h12_yz(i, j, k) + h23_xy(i, j, k) - h13_yy(i, j, k) - h22_xz(i, j, k)) , &
    (h13_yz(i, j, k) + h23_xz(i, j, k) - h12_zz(i, j, k) - h33_xy(i, j, k)) , &
    /) &
    / (metric%dx) ** 2 / 144._dl ) / ak ** 2 - 3._dl * Hk * metric_pk(:, i, j, k) + (R / 6._dl * 4._dl - 2._dl * Hk ** 2) * metric_hk(:, i, j, k)

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

12

帖子

5

主题

0

精华

入门

F 币
84 元
贡献
47 点
沙发
 楼主| 发表于 2020-7-9 09:12:55 | 只看该作者
原来是最后打多了两个",",打扰大家了

2022

帖子

12

主题

5

精华

论坛跑堂

臭石头雪球

F 币
1598 元
贡献
689 点

美女勋章热心勋章星光勋章新人勋章贡献勋章管理勋章帅哥勋章爱心勋章规矩勋章元老勋章水王勋章

板凳
发表于 2020-7-9 18:28:24 | 只看该作者
(/   /)
在目前,可以写成
[]
更美观一些。

12

帖子

5

主题

0

精华

入门

F 币
84 元
贡献
47 点
地板
 楼主| 发表于 2020-7-17 11:16:06 | 只看该作者
fcode 发表于 2020-7-9 18:28
(/   /)
在目前,可以写成
[]

谢谢提醒!
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

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

GMT+8, 2024-11-23 13:13

Powered by Tencent X3.4

© 2013-2024 Tencent

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