Fortran Coder

查看: 1204|回复: 2

[输入输出] 利用format 输出矩阵时如何将format里面的10直接变成N

[复制链接]

22

帖子

11

主题

0

精华

入门

F 币
97 元
贡献
54 点
发表于 2024-8-19 16:21:26 | 显示全部楼层 |阅读模式
[Fortran] 纯文本查看 复制代码
  Module m_wik
  Contains
  Subroutine Wik_()
  Implicit none
  Integer,parameter::N=10
  Integer::a=0,b=0,i,k,j,c,d
  real*8::G(N),M(N),Wik(N,N),x(N),y(N),z(N),H,pi=acos(-1.0),L=16

  Do i=1,N
    x(i)=(L/2)*(1-cos(((i-1.0)/(N-1.0))*pi))
  END Do
  open(unit=16,file='X.txt')
  write(16,105)x
105 format(T4,'x=',10(/F21.8))
  close(16)
  y=x
  z=x
  DO
100 a=a+1
    Do
      b=1+b
      if (a==b)then
        i=a
        k=b
        Wik(i,k)=sum(1/(x(i)-y),y/=x(i))
        if(a==N.and.b==N)then
          exit
        end if
      end if
      if(a/=b.and.b<=N)then
        i=a
        k=b
        G(i)=product(x(i)-y,y/=x(i))
        M(k)=product(z(k)-y,y/=z(k))
        H=1/(x(i)-x(k))
        Wik(i,k)=H*G(i)/M(k)
      end if
      if (b==N)then
        b=0
        Goto 100
      end if
    end do
    exit
  end do
  open(unit=15,file='Wik.txt')
  write(15,104)((Wik(c,d),d=1,N),c=1,N)!写入Wik
104 format(10(F20.10,2x))
  close(15)
  end Subroutine Wik_
  end Module m_wik

这个format里面的10改成N

这个format里面的10改成N

174

帖子

2

主题

1

精华

大师

Vim

F 币
1061 元
贡献
497 点

规矩勋章

发表于 2024-8-19 16:30:17 | 显示全部楼层
ifort可以用扩展
[Fortran] 纯文本查看 复制代码
format( <N>(F20.10,2x))


gfortran目前没什么直接的语法,需要用字符串自己造一个

22

帖子

11

主题

0

精华

入门

F 币
97 元
贡献
54 点
 楼主| 发表于 2024-8-19 16:45:13 | 显示全部楼层

对了,就是这个,非常感谢,解决了
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

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

GMT+8, 2025-4-24 09:04

Powered by Tencent X3.4

© 2013-2025 Tencent

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