QPMS
/
qpms
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

yet some "real" tests 

Former-commit-id: f09d2f26f3f36cf2838deaadccd69cce7e8f1906
This commit is contained in:
Marek Nečada 2018-09-19 07:17:04 +03:00
parent 6a58a89628
commit 5426d5064c
3 changed files with 12 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

0
qpms/apps/hexlattice_ewald.c → qpms/apps/hexlattice_wrongewald.c
View File

0
qpms/apps/hexlattice_ewald_altin.c → qpms/apps/hexlattice_wrongewald_altin.c
View File

12
tests/ewaldshift.c
View File

@ -10,6 +10,14 @@
#include <stdio.h> #include <stdio.h>
#include <float.h> #include <float.h>
#include <gsl/gsl_sf_legendre.h> #include <gsl/gsl_sf_legendre.h>
#include <gsl/gsl_const_mksa.h>
static const double eV = GSL_CONST_MKSA_ELECTRON_CHARGE;
static const double c0 = GSL_CONST_MKSA_SPEED_OF_LIGHT;
static const double hbar = GSL_CONST_MKSA_PLANCKS_CONSTANT_HBAR;
typedef struct ewaldtest_triang_params { typedef struct ewaldtest_triang_params {
qpms_l_t lMax; qpms_l_t lMax;
point2d beta; point2d beta;
@ -39,6 +47,10 @@ typedef struct ewaldtest_triang_results {
ewaldtest_triang_params paramslist[] = { ewaldtest_triang_params paramslist[] = {
// lMax, beta, shift, k, a, eta, maxR, maxK, csphase, orientation // lMax, beta, shift, k, a, eta, maxR, maxK, csphase, orientation
{3, {2.7/576e-9, 1/576e-9}, {576e-9/2,0}, 4.1*M_PI/576e-9, 576e-9, 1/576e-9, 40*576e-9, 40/576e-9, 1., TRIANGULAR_VERTICAL},
{3, {2.7/576e-9, 1/576e-9}, {576e-9/2,0}, 4.1*M_PI/576e-9, 576e-9, 2/576e-9, 40*576e-9, 40/576e-9, 1., TRIANGULAR_VERTICAL},
{3, {2.7/576e-9, 1/576e-9}, {576e-9/2,0}, 4.1*M_PI/576e-9, 576e-9, 3/576e-9, 40*576e-9, 40/576e-9, 1., TRIANGULAR_VERTICAL},
{3, {2.7/576e-9, 1/576e-9}, {576e-9/2,0}, 4.1*M_PI/576e-9, 576e-9, 4/576e-9, 40*576e-9, 40/576e-9, 1., TRIANGULAR_VERTICAL},
/* /*
{ 2, {2.7, 1}, {0.5,0.1325}, 2.3, 0.97, 0.5, 20, 160, 1., TRIANGULAR_VERTICAL}, { 2, {2.7, 1}, {0.5,0.1325}, 2.3, 0.97, 0.5, 20, 160, 1., TRIANGULAR_VERTICAL},
{ 2, {2.7, 1}, {0.5,0.1325}, 2.3, 0.97, 1.5, 20, 160, 1., TRIANGULAR_VERTICAL}, { 2, {2.7, 1}, {0.5,0.1325}, 2.3, 0.97, 1.5, 20, 160, 1., TRIANGULAR_VERTICAL},