代码改变世界

随笔档案-2012年11月14日

[转]在linux下如何使用Makefile对fortran程序进行编译

2012-11-14 20:38 by H_OU, 6685 阅读, 收藏,
摘要: -在linux下写过不少的fortran程序,但很少写makefile文件,通常情况下我是手动用ifort或pgf90进行编译,虽然也在大型程序中写过一些代码,但因为都有现成的makefile文件,我也只是修改一下配置参数,先make clean,再make,只要生成可执行程序就算搞定了。今天心血来潮,自己写了几个Makefile,供大家参考。一、不用Makefile的直接编译比如有2个简单程序,主程序为main.F90,子程序为sub.F90. 如果直接编译的话,键入ifort main.F90 sub.F90就可以生成一个a.out文件,即为可执行程序了。你也可以指定一个可执行程序的名字, 阅读全文

【转】一些解决变态数学公式的算法地址

2012-11-14 20:26 by H_OU, 1067 阅读, 收藏,
摘要: 一些解决变态数学公式的算法网址2011-09-26 09:48:16|分类: 默认分类 |标签:矩阵分析 |字号大中小订阅http://www.acm.org/calgo/contents/各种算法的程序。http://www.kinjo.ac.jp/~kawabe/qc/mig.htmlMIG,分子积分生成程序。直接下载,可运行于Windows/UNIX。 http://forcal.go.nease.net/Forcal数学软件,用于数学计算的库函数。http://www.cnvf.com/Fortran 博物馆。http://www.matlab-world.com/MATLAB大观园。 阅读全文

【转】fortran 广义逆 子程序

2012-11-14 16:56 by H_OU, 364 阅读, 收藏,
摘要: SUBROUTINE BGINV(M,N,A,AA,L,EPS,U,V,KA,S,E,WORK)DIMENSION A(M,N),U(M,M),V(N,N),AA(N,M)DIMENSION S(KA),E(KA),WORK(KA)DOUBLE PRECISION A,U,V,AA,S,E,WOEKCALL BMUAV(A,M,N,U,V,L,EPS,KA,S,E,WORK) !一般实矩阵奇异值分解IF (L.EQ.0) THENK=110 IF(A(K,K).NE.0.0) THENK=K+1IF (K.LE.MIN(M,N)) GOTO 10ENDIFK=K-1IF(K.NE.0) THE 阅读全文

【转】SVD分解

2012-11-14 11:13 by H_OU, 534 阅读, 收藏,
摘要: 分类的关键是计算相关性。我们首先对两个文本计算出它们的内容词,或者说实词的向量,然后求这两个向量的夹角。当这两个向量夹角为零时,新闻就相关;当它们垂直或者说正交时,新闻则无关。当然,夹角的余弦等同于向量的内积。从理论上讲,这种算法非常好。但是计算时间特别长。通常,我们要处理的文章的数量都很大,至少在百万篇以上,二次回标有非常长,比如说有五十万个词(包括人名地名产品名称等等)。如果想通过对一百万篇文章两篇两篇地成对比较,来找出所有共同主题的文章,就要比较五千亿对文章。现在的计算机一秒钟最多可以比较一千对文章,完成这一百万篇文章相关性比较就需要十五年时间。注意,要真正完成文章的分类还要反复重复上述 阅读全文