scnc
/
phonopy
mirror of https://github.com/phonopy/phonopy.git
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

Add special case of sinh function that is used for counting phonon occupation

This commit is contained in:
Atsushi Togo 2014-04-29 16:00:16 +09:00
parent 692a1b7707
commit 12ca7b4a4a
1 changed files with 6 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
c/anharmonic/phonoc_utils.c
View File

@ -290,7 +290,12 @@ double bose_einstein(const double x, const double t)
double gaussian(const double x, const double sigma)
{
return INVSQRT2PI / sigma * exp(-x * x / 2 / sigma / sigma);
}
}
double inv_sinh_occupation(const double x, const double t)
{
return 1.0 / sinh(x * THZTOEVPARKB / 2 / t);
}
static int collect_undone_grid_points(int *undone,
char *phonon_done,