24小时热门版块排行榜    

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

maomao1210

金虫 (正式写手)

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

帖主寄言


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

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

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

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


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

felicity6056

至尊木虫 (知名作家)


小木虫(金币+0.5):给个红包,谢谢回帖交流
再请教个问题
block data
        implicit real*8 (a-h,o-z)
        parameter(id=250,jd=250,nfd=10,nfp3=13,maxij=250)
        logical lstop,lsolve,lprint,lblk,ldata
        character*10 title
        common/name/title(nfp3)
        common/cntl/lstop,lsolid
        common/sorc/smax,ssum
        common/error/ier,jer
        common/coef/flow,diff,acof
        common/const/d1,d2,pr,ra,ra1,thita,ar,xlen,lgrid,mgrid
        common/indx/nf,nfmax,np,nrho,ngam,l1,l2,l3,m1,m2,m3,ist,jst,
     1            iter,last,lsolve(nfd),lprint(nfp3),lblk(nfd),mode,
     2            ldata,ntimes(nfd),ipref,jpref,dx,dy
        common/indx2/relax(nfp3),time,dt,xl,yl,rhocon,yox,nozf
        data nfmax,np,nrho,ngam/10,11,12,13/
        data lblk/10*.true./
        data relax,ntimes/nfp3*1.,nfd*1/
        data dt,ipref,jpref,rhocon/1.e10,1,1,1.29/
        data lstop,lsolve,lprint,ldata/25*.false./
        data mode,last/2,100/
        data title(1)/'   u      '/
        data title(2)/'   v      '/
        data title(3)/'   sf     '/
        data title(4)/'   fuel   '/
        data title(5)/'   oxy.   '/
        data title(6)/'   temp   '/
        data title(11)/'   pres.  '/
        end
这一段表示什么意思?最好能给每行解释的详细点,虫子的牛人们!

[ Last edited by felicity6056 on 2009-8-10 at 17:03 ]
努力的过程讲究方法!
24楼2009-08-10 17:02:09
已阅   回复此楼   关注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的回帖
普通表情 高级回复(可上传附件)
信息提示
请填处理意见