Ex03 做计算常用的批量处理方法（三）
前面练习中我们在0.01的文件夹基础上，通过一个命令，复制得到了从0.02到0.09的文件夹。但是， 所有文件夹中的输入文件都是一样的，我们还需要把INCAR中的SIGMA参数值 SIGMA = 0.01 改成与文件夹对应的数值。 首先我们可以逐个进行编辑，但太浪费时间，这也不是大师兄的风格。Ex03练习分为2小节：
- 新命令 sed 的学习
- for + sed 组合
最终我们会结合for循环和sed命令，来学会批量处理输入文本的另一个方法。还是要强调一下：大家要主动，多去网上找资料，并系统性的学习linux下面的基本命令。光指望着本书中的这么一点，是很难提高的。
复习上一节的功课
1 | iciq-lq@ln3:/THFS/home/iciq-lq/LVASPTHW$ ls |
上面的操作依次为：
新建ex03文件夹
将ex02中的0.01 复制到ex03中
- 进入ex03/0.01的目录中
- 使用cat 查看0.01中INCAR的内容
- 使用cat -n 查看INCAR中的内容，该选项在输出中给每一行标出了行数，方便我们查看。
注意:
1）Linux 下面的命令都有很多的选项，用以丰富我们的不同需求，比如上面的 cat -n，可以使用 man cat 这个命令查看 cat 的其他选项。使用这个命令后，如果想退出，敲 q 键即可;
2）另外，我们也可以使用 cat –help 来查看，效果与man cat 一样。
sed 命令修改INCAR
前面我们提到，可以使用vim打开INCAR然后修改SIGMA的参数。除了vim当然还有文本编辑器等其他的工具。但这些工具都有个缺点，就是得把文件打开后才能修改。下面我们使用sed命令，不打开文本，直接对里面的内容进行替换操作。
1 | iciq-lq@ln3:/THFS/home/iciq-lq/LVASPTHW/ex03/0.01$ sed '3s/0.01/0.02/g' INCAR |
精解：
1) 单引号中是我们的操作， 3s 表示的是选择第三行，因为我们知道 0.01 在第三行中出现，s 是substitute 的缩写
2) 3s 后面跟一个斜杠 / 用来和后面被替换的内容分开，这里0.01 表示选择第三行的0.01;
3) 0.01后面再用一个斜杠，将其和替换后的数字分开(0.01 0.02 0.03 等)，表示将0.01替换为斜杠后面的内容;
4) 再加一个斜杠，后面的g 代表 global ，意思是全部替换。
5) 输入完毕后，我们选择要执行该命令的对象(要替换的文件)，也就是当前目录下INCAR 文件。
6) 命令的意思就是：我们用sed命令，将INCAR中的第三行的0.01全部替换成0.02。
从上面实例中最后的cat INCAR命令结果不难发现，实际上我们并没有将INCAR文件中的0.01替换成0.02。也就是说这个命令只是输出了替换后的结果，但没有更新INCAR文件。那怎么样才可以更新INCAR文件呢？ 我们可以这样做：
1 | iciq-lq@ln3:/THFS/home/iciq-lq/LVASPTHW/ex03/0.01$ sed '3s/0.01/0.02/g' INCAR > INCAR_new |
箭头（>）的意思是：我们将命令的输出存到一个新的文件INCAR_new中，
通过mv命令之前的INCAR替换掉。
但，这样做也太麻烦了，更简单一点，如下：
前面例子的INCAR中SIGMA的值已经不是0.01了，我们先从ex02/0.01中复制一个过来。
1 | iciq-lq@ln3:/THFS/home/iciq-lq/LVASPTHW/ex03/0.01$ cp ../../ex02/0.01/INCAR . |
详解：
sed –i 是sed 的命令和其附加选项， -i 表示直接对源文件进行编辑，也就是说编辑之后源文件被新文件替换掉。因此，使用这个参数的时候要小心，小心，再小心。要格外小心！！！
- 最保险的做法就是运行前，先对操作的对象进行备份：如下：
1 | iciq-lq@ln3:/THFS/home/iciq-lq/LVASPTHW/ex03/0.01$ cp INCAR INCAR_back |
- 其次是，先不加 -i 运行下sed命令，确保输出的是正确结果后，然后再加上 -i 运行.
1 | iciq-lq@ln3:/THFS/home/iciq-lq/LVASPTHW/ex03/0.01$ cp INCAR_back INCAR |
小结:
sed 是一个非常强大的命令，对于做计算的我们来说，熟练正确地使用sed可以极大的提高我们的工作效率，大家务必硬着头皮掌握这个命令。这个网站列举了一些基本的用法：
大家参照着进行练习，也可以百度里面搜索一些其他的 sed 使用技巧，如果你有认为很好的sed 技巧，也可以发邮件分享给大师兄（lqcata@gmail.com）。
本节我们学习2个批量操作：
- 通过sed单个命令进行批量操作
- 以及sed + for循环的批量操作
1 sed 批量将0.01到0.09中所有INCAR中的0.01替换成0.05。到现在为止，相信大家都可以看懂下面的命令操作。就不再啰嗦解释了。有一点需要注意的是grep 命令中的星号，检查输入输出的时候用*非常方便。
1 | iciq-lq@ln3:/THFS/home/iciq-lq/LVASPTHW/ex03$ ls |
2 for 循环结合sed
前面我们使用sed命令，将文件夹中所有的0.01替换成了0.05。但我们的目标是，每个文件夹中的SIGMA值与文件夹相同。既然我们知道了sed可以对单个文件进行操作，那么我们也可以结合for循环，来实现一个批量操作的目的。命令如下：
1 | iciq-lq@ln3:/THFS/home/iciq-lq/LVASPTHW/ex03$ ls |
注意：
- 这里我们用的是双引号 “ “ ，sed 命令中你会见到大部分都用单引号 ‘ ‘ 。但如果这里使用单引号，则所有的 0.01 都会被替换成 $i (单引号中的$i 是纯字符)，因为单引号中的所有内容都会被当做字符来处理，也就是里面是什么就输出什么。使用双引号，则可以读取变量 $i 的值，下面的例子大家一看就知道怎么回事了:
1 | iciq-lq@ln3:/THFS/home/iciq-lq/LVASPTHW/ex03$ abc=bigbro |
这里单引号中的内容被原封不动地打印出来了。而双引号的话，则可以顺利地把变量调用起来。
- for i in ; 这里的 指的是当前目录下所有的文件以及文件夹，本例中没有文件，只有从0.01， 0.02， 0.03 到 0.09 的文件夹；所以： for i in * = for i in 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09
1.5.2 使用sed 命令将INCAR中的 0.05 (所有文件夹中的都是0.05)替换成文件夹的数字;
1 | sed -i "3s/0.05/$i/g" $i/INCAR |
当$i的值为0.01的时候，我们就把0.01/INCAR中的0.05替换为0.01；为0.02的时候，就把0.02/INCAR中的0.05替换为0.02，依次类推，直至for循环完成。
1.6 该命令瞬间运行完成，我们使用grep SIGMA 来查看下这些文件夹中的INCAR参数 ，圆满完成任务!
总结
在ex03的练习中，我们并不着急去提交任务，反正已经走了计算的路，以后提交任务的机会还多着呢，也不差这一两天的学习时间。先从简单到复杂些，逐渐掌握Linux的一节操作命令，对以后的学习帮助很大。
1) 学会 man command 或者 command –help 查看命令的具体参数;
2) 大量使用sed 命令进行操作练习;
3) 知道 * 在for循环中代表的含义
3) 熟知单引号和双引号的使用，以及区别
4) 主动搜索相关的Linux命令的相关知识，积极操练。
5) 学会从单一操作，通过for循环转化成批量操作的思路。
- for i in XXX; do YYY; done
- XXX就是我们要操作的范围或者对象
- YYY就是单一的一个操作。