Fortran Coder

查看: 155|回复: 3

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

[复制链接]

12

帖子

5

主题

0

精华

入门

F 币
82 元
贡献
46 点
发表于 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 币
82 元
贡献
46 点
 楼主| 发表于 2020-7-9 09:12:55 | 显示全部楼层
原来是最后打多了两个",",打扰大家了

1538

帖子

12

主题

5

精华

论坛跑堂

臭石头雪球

F 币
1417 元
贡献
936 点

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

发表于 2020-7-9 18:28:24 | 显示全部楼层
(/   /)
在目前,可以写成
[]
更美观一些。

12

帖子

5

主题

0

精华

入门

F 币
82 元
贡献
46 点
 楼主| 发表于 2020-7-17 11:16:06 | 显示全部楼层
fcode 发表于 2020-7-9 18:28
(/   /)
在目前,可以写成
[]

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

本版积分规则

QQ|捐赠本站|Archiver|关于我们 About Us|群聊|Fcode

GMT+8, 2020-8-14 21:31

Powered by Discuz! X3.2

© 2001-2017 Comsenz Inc.

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