Fortran
文章平均质量分 89
weixin_38594380
这个作者很懒,什么都没留下…
展开
-
Fortran 稀疏矩阵
Fortran 处理稀疏矩阵稀疏矩阵Ax=b在Fortran里面使用稀疏矩阵最基础的是用BLAS(Basic Linear Algebra Subprograms),但是在后来的MKL库中有集成BLAS。如果是解Ax=b的线性方程组,建议使用Pardiso,同样在MKL库中有集成,可以去官网查找资料。这里给出两篇博文,介绍了如何在Fortran里面解线性方程组以及BLAS库函数的示例https://blog.csdn.net/chd_lkl/article/details/83011186http原创 2021-04-26 14:36:35 · 2277 阅读 · 6 评论 -
Fortran 编译注意事项
Fortran 编译注意事项高精度运算如果要做精度要求很高的计算,一般在声明实数型变量的时候就要指定是单精度还是双精度,用real(kind=4)便表示单精度,相应地还可以改为8,10,12分别表示双精度及更高的精度。在学习Fortran的过程中应该是注意到,算式2/3和2.0/3.0计算结果是不一样的,前者是以整型数据计算后者是按浮点数计算。但是!后者的计算精度有时候还是不够,比如在计算1.0/3.0并赋值给一个双精度变量的时候,计算结果并不是0.3333333333333333,这里要改写成1.0原创 2020-09-14 16:19:36 · 1140 阅读 · 0 评论 -
《Fortran 95 程序设计》阅读笔记三
Chapter 8 函数子程序 SUBROUTINE把经常需要用到的程序封装成为一个子程序,通过 call 进行调用 call sub1(para1, para2)subroutine sub1() ...... returnend subroutine sub1 ! 这里的sub1作为子函数名也可以省去其中最后一行一般是 return 指返回到先前调用子程序的位置,不写的话算作默认添加,当然return也可以放在任何地方。其次,关于子程序的变量和行代码是与主程序相互独立的,利用行代码进原创 2020-09-03 16:28:29 · 565 阅读 · 0 评论 -
《Fortran 95 程序设计》阅读笔记二
Chapter 5 流程控制与逻辑运算If 语句规范的表达式写作:if (逻辑判断式) then ......else ......end if如果只有单行命令的话还可以简写为 if (逻辑判断式) write(*,*) "Hello World!"特别注意在Fortran 90里面关于等于和不等于的判断表达式分别为 == 和 \=逻辑判断语句的连接:连接语句运算规则.and.交集.or.并集.not.逻辑反向,若后面表达式不成立则整个表达式原创 2020-09-02 18:11:14 · 621 阅读 · 0 评论 -
《Fortran 95 程序设计》阅读笔记一
Chapter 3 Fortran程序设计基础Fortran 代码分为固定格式fixed format和自由格式free format,前者基本已经被淘汰。后者在使用的时候需要注意:感叹号表示注释每行可以编写132个字符行号放在每行程序最前面一行程序的最后如果是 & 或者第一个字符是 & 则表示接着上一行Fortran不区分大小写!!!Fortran的数据类型整型 INTEGER:又分为长整型10910^9109 量级和短整型10510^5105量级浮点数 REAL:原创 2020-09-01 19:40:53 · 1388 阅读 · 0 评论