摘要:
In mathematics, the Courant–Friedrichs–Lewy condition (CFL condition) isa necessary condition for convergence while solving certain partial differential equations (usually hyperbolic PDEs) numerically by the method of finite differences.[1]虽然在Citcom解地幔对流中,每一个时间步是有限元法,但是在时间演化上,仍然是有限差分。我现在,还没见到过有限元用在时 阅读全文
随笔档案-2012年11月
(转)u检验、t检验、F检验、X2检验
2012-11-22 00:59 by H_OU, 5047 阅读, 收藏, 编辑
摘要:
常用显著性检验 1.t检验 适用于计量资料、正态分布、方差具有齐性的两组间小样本比较。包括配对资料间、样本与均数间、两样本均数间比较三种,三者的计算公式不能混淆。 2.t'检验 应用条件与t检验大致相同,但t′检验用于两组间方差不齐时,t′检验的计算公式实际上是方差不齐时t检验的校正公式。 3.U检验 应用条件与t检验基本一致,只是当大样本时用U检验,而小样本时则用t检验,t检验可以代替U检验。 4.方差分析 用于正态分布、方差齐性的多组间计量比较。常见的有单因素分组的多样本均数比较及双因素分组的多个样本均数的比较,方差分析首先是比较各组间总的差异,如总差异有显著性,再进行组间... 阅读全文
contains
2012-11-20 11:16 by H_OU, 298 阅读, 收藏, 编辑
摘要:
prigram main...... ......call p1()call p2().......... containssubroutine p1() ............end subroutine p1subroutine p2() ............end subroutine p2end! 子程序中的变量无需重新说明,变量值无需通过哑实结合传递。elemental代表的是你可以传数组给这个函数,而这个函数会一个一个地处理这个数组的,并且是相互独立的(可并行)external代表这后面声名的子程序是在别处定义的,不属于当前子程序的内部子程序contains代表下... 阅读全文
[转]在linux下如何使用Makefile对fortran程序进行编译
2012-11-14 20:38 by H_OU, 6566 阅读, 收藏, 编辑
摘要:
-在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, 1008 阅读, 收藏, 编辑
摘要:
一些解决变态数学公式的算法网址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, 338 阅读, 收藏, 编辑
摘要:
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, 502 阅读, 收藏, 编辑
摘要:
分类的关键是计算相关性。我们首先对两个文本计算出它们的内容词,或者说实词的向量,然后求这两个向量的夹角。当这两个向量夹角为零时,新闻就相关;当它们垂直或者说正交时,新闻则无关。当然,夹角的余弦等同于向量的内积。从理论上讲,这种算法非常好。但是计算时间特别长。通常,我们要处理的文章的数量都很大,至少在百万篇以上,二次回标有非常长,比如说有五十万个词(包括人名地名产品名称等等)。如果想通过对一百万篇文章两篇两篇地成对比较,来找出所有共同主题的文章,就要比较五千亿对文章。现在的计算机一秒钟最多可以比较一千对文章,完成这一百万篇文章相关性比较就需要十五年时间。注意,要真正完成文章的分类还要反复重复上述 阅读全文
【转】grads遇到的一些问题
2012-11-07 15:43 by H_OU, 2475 阅读, 收藏, 编辑
摘要:
1.【转】 grads中同时使用多个地图叠加的方法使用GrADS作图时,有时候会用到不同分辨率的地图叠加在一起,比如说,世界地图+省地图+区县地图,以前笨笨总是多次设置,然后画地图,今天偶然发现可以一次设置多个地图,然后一起画的。'set mpdset cnworld jiangsu jiangsu_q''draw map'可以将地图放在grads-1.9b4/data目录下2.Error opening stroke character data set 错误需要加两个系统变量 GADDIR=d:\pcgrads\dat GASCRP=d:\pcgrads\l 阅读全文