JN182-优快云博客

原创 gfortran和gcc的降级安装中的那些坑

Linux环境是win10中的Linux子系统（WSL)，WSL搭建完毕后我迅速下载了gcc。安装之后发现是gcc-9.3.0版本的，无法运行老代码。于是准备降级安装低版本的gcc编译器。这个过程中我踩了很多坑，先说成功的一次安装吧。GCC降级安装>>ls /usr/bin/gcc*/usr/bin/gcc /usr/bin/gcc-ar /usr/bin/gcc-nm /usr/bin/gcc-ranlib/usr/bin/gcc-4.8 /usr/bi

2020-08-04 09:31:05 5342 2

原创联想ThinkPad E14开箱测评跑分

特别感谢黑科技小王子治文兄推荐的机型，教我测评，把我从小白带进了笔记本发烧友的行列中。抓住7月的尾巴，2020年7月31日夜里22点58分下单联想ThinkPad E14。今天建军节，中午刚好到货，迫不及待开箱，京东货源，售价￥3999，可以说是ThinkPad的最低价系列之一了。非常适合网页程序办公开发人员，把钱花在CPU和内存上。虽然显卡很鸡肋，但是对于不打游戏不玩设计的人来说，不太影响使用体验。开箱过程我录了一个视频，过两天上传到b站上。硬件参数处理器处理器 AMD Ryzen 5 4.

2020-08-01 22:28:56 13650 2

原创 LSQR算法的程序设计（带阻尼与无阻尼）

LSQR算法是计算大型稀疏线性方程组的算法，由Paige和Saunders于1982年提出。该方法主要为求解以下线性方程组（A为m*n的矩阵，m>n），且保持二阶残差范数最小：Ax=bMin∥Ax−b∥2Ax=b\\ Min\parallel Ax-b \parallel_{2}Ax=bMin∥Ax−b∥2对于传统的最小二乘问题，可以直接求出方程的解，但是对于大型稀疏矩阵，往往是病态问题（超定方程组或欠定方程组），因此，LSQR方法不是直接求解x，而是遵守最小二乘的原则，在法方程限定的空间中

2020-07-13 16:14:20 4615 1

原创大型稀疏矩阵的一维存储及其乘法运算

大型稀疏矩阵对存储和CPU要求高，很多集成化软件对于特大型矩阵的运算束手无策，因此需要一种优良的数据结构筛除稀疏矩阵中的0，从而实现数据的高效存储和运算。对于m行n列的大型稀疏矩阵A（large sparse matrix），它大多数元素都为0，则可以使用三个一维数组存储该矩阵：A(k),JA(k),NA(n).每个数据的含义：A(k)：若A共有m*n个元素，其中有k个非零元素，每个非零元素按行存储，A(i)为第i个非零元素。JA(k)：与数组A一样，有k个元素；JA(i)代表第i个非零元素在矩阵A的

2020-07-01 23:52:01 1446

原创基于Lanczos方法的矩阵双对角化或三对角化

LSQR算法是计算大型稀疏线性方程组的算法，由Paige和Saunders于1982年提出。该方法主要为求解以下线性方程组（A为m*n的矩阵，m>n），且保持二阶残差范数最小：Ax=bMin∥Ax−b∥2Ax=b\\ Min\parallel Ax-b \parallel_{2}Ax=bMin∥Ax−b∥2对于传统的最小二乘问题，可以直接求出方程的解，但是对于大型稀疏矩阵，往往是病态问题（超定方程组或欠定方程组），因此，LSQR方法不是直接求解x，而是遵守最小二乘的原则，在法方程限定的空间中

2020-06-26 22:17:20 6044 2

原创 fortran特定位置文件读写（open、write和read命令）

fortran在科学计算中常常需要从文件中读出和写入大量数据，因此指定读写位置就尤其重要，本文介绍三个命令的主要选项。打开命令open的选项当我们使用open命令时，默认的打开选项为：open(UNIT=number, FILE=‘filename’, FORM=‘FORMATTED’, STATUS=‘UNKNOWN’, ACCESS=‘SEQUENTIAL’, RECL=length, ERR=label, IOSTAT=iostat, BLANK=‘NULL’, POSITION=‘ASIS’,

2020-06-18 10:40:57 35746 2

原创如何编译含有多个文件的fortran程序包（Linux多文件编译）

如何编译含有多个文件的fortran程序包（Linux多文件编译）初学Linux，拿到的程序都是压缩包，解压之后发现各个程序互相嵌套，而且往往不是用同一种语言写的，因此在编译的过程，需要厘清每个文件之间的逻辑关系。下面就谈一谈我初学多文件编译的心得吧：对于有Makefile的程序包这样的程序包往往关系复杂，但是编译起来却很简单：不需要按照文件的结构依次编译每个文件，只需要在shell中输入make即可全部编译完成。如果非常关心各个程序之间的相互关系，只需要打开Makefile查看即可。关于Make

2020-06-15 21:54:02 5016

原创 fk算法绘制层状介质理论地震图

fk算法绘制层状介质理论地震图经过了一个学期理论地震学推公式写程序的折磨，作为小白，从对Linux系统一窍不通到多文件编译运行，能正演出了第一张理论地震图还是非常开心的！朱露培老师的fk3.2程序包下载：http://www.eas.slu.edu/People/LZhu/home.html首先在shell中将程序解压，编译，得到的主要文件有：fk.f 计算层状介质中爆炸源/单力/力偶的格林函数（.grn文件）syn.c 用已知的震源机制和fk.f得到的格林函数合成理论地震图（.z

2020-06-15 21:11:46 4467 13

weixin_44498266的博客