- 博客(6)
- 收藏
- 关注
RSS订阅原创 Fortran 稀疏矩阵
Fortran 处理稀疏矩阵稀疏矩阵Ax=b在Fortran里面使用稀疏矩阵最基础的是用BLAS（Basic Linear Algebra Subprograms），但是在后来的MKL库中有集成BLAS。如果是解Ax=b的线性方程组，建议使用Pardiso，同样在MKL库中有集成，可以去官网查找资料。这里给出两篇博文，介绍了如何在Fortran里面解线性方程组以及BLAS库函数的示例https://blog.youkuaiyun.com/chd_lkl/article/details/83011186http
2021-04-26 14:36:35
2674
6
原创 VS下编译Fortran
VS下编译Fortran参考博客https://www.liangzl.com/get-article-detail-161716.html完成Intel Parallel Studio XE 2020的安装，在这个环境下可以编译Fortran如果需要用到数学库，记得勾选MKL，并在之后参考博客http://www.guanjihuan.com/archives/3021完成配置。但是要注意里面的第二条要稍作修改，因为找不到那些文件，我在这改成了：mkl_intel_thread.lib m
2021-02-27 22:01:59
4365
原创 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
1378
原创 《Fortran 95 程序设计》阅读笔记三
Chapter 8 函数子程序 SUBROUTINE把经常需要用到的程序封装成为一个子程序，通过 call 进行调用 call sub1(para1, para2)subroutine sub1() ...... returnend subroutine sub1 ! 这里的sub1作为子函数名也可以省去其中最后一行一般是 return 指返回到先前调用子程序的位置，不写的话算作默认添加，当然return也可以放在任何地方。其次，关于子程序的变量和行代码是与主程序相互独立的，利用行代码进
2020-09-03 16:28:29
638
原创 《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
815
原创 《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
1525
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人