luyi07的笔记

摘要：(未经允许,不得转载) 参考文献：J. M. Thijssen, 《Computational Physics》 1. Formalism 为了简化，设一维谐振子的哈密顿量为 $$ \hat{H} = - \frac{d^2}{dx^2} + x^2 $$ 引入拟设波函数 $$ \psi(x, \a 阅读全文

摘要：这是上次一个小文献笔记(https://www.cnblogs.com/luyi07/p/15442971.html)里一个定理的实践。 1. 实数反对称矩阵 \(M\) 所有矩阵元为实数，并且有反对称性 \(M^\top = - M\)。 2. 反对称阵的正则形式 如果反对称矩阵 \(M\) 是块 阅读全文

摘要：随便记在这吧 LAPkb3g.sh #!/bin/bash # ./LAPkb3g.sh 1 1 example/kb3g.jsp example/kb3g.pn y 10 0.001 0.001 projection_pf0202_p.dat projection_pf0202_n.dat nuc 阅读全文

摘要：我们需要做很多很多小矩阵相乘（维数只有几十），但是次数很多，所以用哪个矩阵库的函数对我们很重要。这里写一个很小的测试代码，测试lapack(包含着朴素的blas)，cblas，还有手写函数，对比它们做小矩阵相乘的效率。 对于给定的维数，这三种办法，每种都做１０００次方阵相乘$AB = C$，每次相乘 阅读全文

摘要：Ising模型哈密顿量为 \[ \hat{H} = - J \sum_{i,j相邻} s_i s_j - B \sum_i s_i \] 若取 \(J=1, B=0\)，做 4 x 4 一共 16 个 spin 的 Ising 模型模拟，可以得到很平缓的相变曲线，如下图 曲线最左侧有些点低于低于1. 阅读全文