fortran | 寒山烟雨
现在位置: 首页 > fortran
2014年08月19日 滴水穿石 ⁄ 共 12621字 评论 7 条 ⁄ 阅读 7,627 views 次
vasp的编译安装可以说是千奇百怪,我采用的是ivf的编译器和mpich2进行最简单的编译,系统是ubuntu14.04,不得不说,就是因为这个系统,当时浪费了我一个星期的时间进行编译。 首先,安装ifort,选择ifort的一大原因是因为ifort有免费适用版,可以申请到一年的适用权。先用邮箱申请一个lic,Intel会把ID和下载地址通过邮箱告诉你,下载下来正常安装即可(教育网速度尚可,能有10M/s的速度,业界良心啊。) 其次,安装mpich2,ubu...
阅读全文
2014年08月12日 小试身手 ⁄ 共 1948字 评论 3 条 ⁄ 阅读 3,667 views 次
逾渗理论作为一个经典的理论,在宏观、微观领域都有极其重要的应用,如重整化群、电阻网络、相变等等。二维格子是逾渗理论最基本的一种模型,其中逾渗阈值则是其中尤为重要的一个概念。 逾渗模拟最大的难点在于团簇分组,机器不比人,寻常递归分组过于繁琐,我这用的是最笨的二次分组法。 1,从左上角开始向右走,若有格点,向左和向上看最近临点,如果有一个点格点,分为改为改组组号,二者都有格点,选择较小者。 2,依次走遍...
阅读全文
2014年07月08日 滴水穿石 ⁄ 共 962字 暂无评论 ⁄ 阅读 2,275 views 次
郎之万顺磁理论模拟,利用马尔科夫链,并加以模拟退火的方法,绘制出了顺磁体在外加磁场的作用下,感应磁场的大小。理论上该模型有严格解,从模拟的结果来看,还是比较符合理论的。 基本算法:(1)随机给定系统初态;(2)再次给定每个粒子态;(3)计算每个粒子该次与上次态的能量差dE;(4)dE<0,跃迁至该态,dE>0,按照麦克斯韦分布跃迁;(5)重复2~4步足够多的次数,系统总E达到一个稳定值,获得此时磁畴方向。...
阅读全文
2014年06月29日 小试身手 ⁄ 共 489字 评论 2 条 ⁄ 阅读 3,594 views 次
利用最常规的模拟方法——欧拉折线法绘制四电荷电场线,fortran写的,下一个版本会有python的Vpython动画版。  Fortran |   copy code |?01module conts02integer :: N=18,Num=403real(4) :: q(4)04data q/1.0,-1.0,1.0,-1.0/05complex :: p(4)06data p/(1.0,1.0),(-1.0,1.0),(-1.0,-1.0),(1.0,-1.0)/07real(4) :: pai=3.141592608r...
阅读全文
2012年06月22日 滴水穿石 ⁄ 共 934字 评论 3 条 ⁄ 阅读 2,603 views 次
不知不觉接触Fortran已经一年多了,在体会到Fortran在做数值计算上的强大之后,也逐渐开始认识到这门语言的不足。曾经叱咤计算风云的Fortran是否还能风光依旧? 在2012年6月份的编程语言排行榜上,Fortran已经沦落到29位。说起Fortran很多编程新手或许都闻所未闻,然而在编程语言的初级阶段,Fortran是第一个真正意义上的高级语言,被广泛应用于科学计算领域。从1951年IBM着手开发Fortran开始,Fortran已经走过60多个年头。版本...
阅读全文
2012年06月13日 小试身手 ⁄ 共 1101字 评论 4 条 ⁄ 阅读 3,315 views 次
光线在均匀介质中沿直线传播,在非均匀介质中的传播,则可利用折射定律进行解析求解,而对于较为复杂的情况,可以利用数值方法模拟求解。 具体可以将连续介质离散化,每一小块区间内的折射率近似看作不变,由起始位置开始逐步求出下一点位置,直到终点。当小区间很小时,散点之间的连线近似可看作光线的路径。 此处,对折射率一维线性变化的介质做了一次模拟,折射率仅与y方向上成线性关系。对于其他折射率一维变化的介质只需修...
阅读全文
2012年06月03日 小试身手 ⁄ 共 727字 评论 8 条 ⁄ 阅读 6,369 views 次
氢原子镜像波函数方程(化为原子单位制)为: \left( -\frac{d^2}{2dr^2}+\frac{l(l+1)}{2r^2}-\frac{1}{r} \right) P(r)=EP(r) 其边界条件为:P(0)=0,P(∞)=0,l为轨道量子数,l=0,1,2,3…分别对应s轨道,p轨道,d轨道,f轨道。。。 本意使用Nomerov方法,数值求解径向波函数,但是对于r=0处有奇点,此处采用差分法求解这个二阶线性非齐次微分方程的本征值问题。Nomerov有五阶精度,而差分只有二阶精度,出现的问题稍后再说。 由...
阅读全文
2012年04月19日 小试身手 ⁄ 共 1987字 评论 4 条 ⁄ 阅读 3,137 views 次
四阶龙格库塔法解含时薛定谔方程,一维线性谐振子附加含时驱动力。  Fortran |   copy code |?001external f002parameter(n=100001,m=11)003Dimension Y(m),D(m),z(m,n),b(m),c(m,n)004complex(8) y,d,z,b,h005real(8) t,a006t=0.0d0</p>007do j=1,m008Y(j)=(0.0,0.0)009enddo010Y(1)=(1.0d0,0...
阅读全文
×