摘要： 之前写过一个随笔，描述怎么用　gnuplot 绘制球谐函数图：https://www.cnblogs.com/luyi07/p/14713231.html 其中提到，在画球谐函数这事上，python的缺点是图片不能旋转，图片小不够清楚华丽，代码细节多（其实也还好，多一点点）。 现在，真香定律显现，我 阅读全文

摘要： 1. 注释 1.1 单行注释 以 # 开头的即注释。第一行注释如下写的话，有特别意义，代表用什么 shell 执行这个脚本： #!/bin/bash 以 # 开头的其他行都是不影响执行的注释。 1.2 多行注释 如下编写多行注释 :<<! 这里是注释，爱写多少行写多少行 ! 2. 变量 可以如下定义 阅读全文

摘要： 1. 定义 范德蒙德行列式定义为： \[ V(x_1, x_2, \cdots, x_n) \equiv \left| \begin{matrix} 1 & 1 & \cdots & 1 \\ x_1 & x_2 & \cdots & x_n \\ x^2_1 & x^2_2 & \cdots & 阅读全文

摘要： 1. 投影框架下的本征态 如果变分得到的能量最低波函数为 \(|PC\rangle\)，这里 PC 是 pair condensate 的缩写，表示我们最近做的对凝聚组态，但下面的讨论并不局限于这一种组态。 如果要得到角动量为 \((J,M)\) 的近似本征态，我们先从 \(|PC\rangle\) 阅读全文

摘要： 1. 参考来源 参考 lapack 官网上的源码，其中有 dgesv 函数的使用说明：http://www.netlib.org/lapack/lapack-3.1.1/html/dgesv.f.html 其算法就是教科书式的 LU 分解+换行。 2. 测试 然后写一个小代码测试了一下，做一个最简单 阅读全文

摘要： 之前写了投影技术的理论基础，即转动算符、D矩阵、投影算符等理论知识：https://www.cnblogs.com/luyi07/p/14586631.html 又写了两种投影技术的思路与比较：https://www.cnblogs.com/luyi07/p/14823838.html 在这个随笔中 阅读全文

摘要： 1. 回顾：角动量投影 之前写过一篇随笔，整理了角动量投影技术的推导：https://www.cnblogs.com/luyi07/p/14586631.html 投影算符定义为 \[ \hat{P}^J_{LK} | \Phi \rangle = \frac{ 2 J + 1 }{ 8 \pi^2 阅读全文

摘要： 1. mermaid：示意图 vscode 中安装 mermaid 插件以后，可以很容易地画示意图，例如： ```mermaid graph LR; A -->B; ''' 上面最后一行应为```，显示不了，所以用'''示意。效果如下： 还可以做出更多细致的效果，参考：https://blog.cs 阅读全文

摘要： 1. 矩阵：numpy.mat 1.1 设置 0 元素：numpy.zeros 例子： import numpy as np print(np.zeros((3,3))) 输出 [[0. 0. 0.] [0. 0. 0.] [0. 0. 0.]] 1.2 设置矩阵 例子： p = np.mat(np 阅读全文

摘要： 1. 动机：没有公网IP的痛 A：大家好，我是 A 机器，我爱 C 机器，但是 C 没有公网 IP，所以我看不见它，暗恋使我心焦+痛苦，谁替我捎给它秋天的菠菜？ C：我是 C 机器，我矜持、羞涩，但我也希望有人，给我带来秋天的菠菜。 B：兄弟们！我是 B 机器，我有公网 IP！我就是天上的鹊桥，是情 阅读全文