24小时热门版块排行榜    

CyRhmU.jpeg
查看: 10747  |  回复: 157
本帖产生 1 个 程序强帖 ,点击这里进行查看
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

maomao1210

金虫 (正式写手)

[交流] 【交流】Fortran语言答疑专帖已有46人参与

帖主寄言


其实语言并不是最重要的,形势的载体而已,fortran擅长工程计算,因为工作需要,偶尔用用fortran。在此开贴目的有二:

第一,希望能和大家交流的同时提高和丰富自己;

第二,认识来自五湖四海的朋友。

资料目前还没有整理,有机会整理上传一些。呵呵。


[ Last edited by nono2009 on 2009-11-18 at 10:34 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

maomao1210

金虫 (正式写手)

★ ★ ★ ★ ★
余泽成(金币+5):谢谢参与应助! 2010-12-06 16:32:26
引用回帖:
Originally posted by astringent at 2010-12-02 04:59:14:
您好,我想请教一个问题,
1       15733
1       18193
1       20251
2       15733
2       20251
3       20824
4       15733
4       16159
4       18193
4       20251
4       20824
5      ...

依据我的理解,写了一个简单的示范程序,可以参考一下,我没有按照你的意思来写,我感觉你或许说的有点问题,呵呵,


      Implicit None
      Integer Max_Parameter,Max_step
      Parameter(Max_Parameter=20,Max_step=20)
      Integer H(Max_Parameter,Max_step)
c        Allocatable H(:,
      Character(Len=60) FileName
        Integer Ihere,Istep,Io,Ibegin
       
c        Write(*,*) 'Please input name of the input files'
       
c        Read(*,*) FileName

c        Open(10,File=Trim(FileName)//'.txt',Status='OlD')
      Open(10,File='aaa.txt',Status='OlD')
      
      Read(10,*) Ihere
      Ibegin=0
      Rewind(10)
14        Read(10,*,End=123) Istep,Io

      If(Istep.eq.Ihere) Then
           Ibegin=Ibegin+1
         H(Ibegin,Istep)=Io
      Else
        Ibegin=1
          Ihere=Istep
          H(Ibegin,Istep)=Io
      EndIF  
      Goto 14
123   Continue
      write(*,*) H(1,1), H(1,2), H(1,3), H(1,4), H(1,5), H(1,6)
        write(*,*) H(2,1), H(2,2), H(2,3), H(2,4), H(2,5), H(2,6)
        write(*,*) H(3,1), H(3,2), H(3,3), H(3,4), H(3,5), H(3,6)
        write(*,*) H(4,1), H(4,2), H(4,3), H(4,4), H(4,5), H(4,6)
        write(*,*) H(5,1), H(5,2), H(5,3), H(5,4), H(5,5), H(5,6)
        End









[ Last edited by maomao1210 on 2010-12-2 at 09:45 ]
88楼2010-12-02 09:42:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 158 个回答

stereochemistry


小木虫(金币+0.5):给个红包,谢谢回帖交流
你好,最近写程序遇见一个问题,一个对称矩阵想线性存储,请问怎么实现呢?
2楼2009-06-01 12:32:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

maomao1210

金虫 (正式写手)

★ ★ ★ ★ ★
wangen994(金币+2,VIP+0):鼓励一下,哈哈 6-4 11:55
gwdavid(金币+3,VIP+0):辛苦了!答疑帖加大奖励力度!;) 6-7 10:32
wangen994(金币+0,VIP+0):请你讲九月份十月份的参与应助的帖子整理附在http://emuch.net/bbs/viewthread.php?tid=1358729&fpage=1后面,以便发放津贴 11-9 19:00
哦。我给你举个例子吧,比如对称矩阵 A[3,3]
                              A(1,1)     A(1,2)     A(1,3)
                              A(2,1)     A(2,2)     A(2,3)
                              A(3,1)     A(3,2)     A(3,3)
只要存储成一个一维数组即可: B(1)=A(1,1),B(2)=A(2,1),B(3)=A(2,2), B(4)=A(3,1),B(5)=A(3,2),B(6)=A(3,3).
还要记住这个: 行和列与存储该值的关系为:II=MAX(I,J)*(MAX(I,J)+1)/2+MIN(I,J), 那么B(II)==A(I,J).
不知道我讲的能听懂与否。如果不懂,继续发问。
3楼2009-06-01 12:41:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

anbb1009

金虫 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
我同意楼上的观点,因为我以前也这样做过,而且效果不错,但顺便问一下:目前fortran中用的最好的求 非线性方程一组实根的方法有哪几种,有没有现成的子程序?谢谢
4楼2009-06-05 14:35:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复(可上传附件)
信息提示
请填处理意见