This shows you the differences between two versions of the page.
Both sides previous revision Previous revision | Last revision Both sides next revision | ||
codes [2014/08/08 15:07] bogner |
codes [2014/08/08 15:08] bogner |
||
---|---|---|---|
Line 6: | Line 6: | ||
* {{:coulomboscrelmeV1.f90.zip|}} Some f90 subroutines (NOT a standalone code) from a larger code that calculates relative two-body matrix element $\langle nl|V|n'l\rangle$ of the Coulomb potential in the HO basis. The main purpose of this code is to illustrate the algorithm to calculate generalized Laguerre polynomials $L^{l+1/2}_n$ that enter the definition of the HO wf's $R_{nl}$. {{makefile.gz|Here}} is the Makefile that was missing. | * {{:coulomboscrelmeV1.f90.zip|}} Some f90 subroutines (NOT a standalone code) from a larger code that calculates relative two-body matrix element $\langle nl|V|n'l\rangle$ of the Coulomb potential in the HO basis. The main purpose of this code is to illustrate the algorithm to calculate generalized Laguerre polynomials $L^{l+1/2}_n$ that enter the definition of the HO wf's $R_{nl}$. {{makefile.gz|Here}} is the Makefile that was missing. | ||
+ | |||
+ | UPDATE: One of the files (renorm-modules.f90) in the TBME package had a subtle but important bug where matrix elements in the spin-triplet channel, which should vanish due to the structure of the Minnesota potential, were non-zero. Here is a fixed version of the file in question {{:renorm-modules.f90.gz|}}. | ||
* {{:eigen_basis.zip|}} C++ code adapted from Dick Furnstahl's Computational Physics course at Ohio State to calculate the Coulomb hamiltonian on a HO basis for $l=0$. This calculates the HO wf's a bit differently than Fortran code above. | * {{:eigen_basis.zip|}} C++ code adapted from Dick Furnstahl's Computational Physics course at Ohio State to calculate the Coulomb hamiltonian on a HO basis for $l=0$. This calculates the HO wf's a bit differently than Fortran code above. |