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

Fix hexlattice ewald kin option 

Former-commit-id: 306c547053b53c594029142095ff3649380e7ab0
This commit is contained in:
Marek Nečada 2018-09-19 22:02:20 +03:00
parent e5412da897
commit dba2b9d104
1 changed files with 3 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
qpms/apps/hexlattice_ewald.c
View File

@ -3,6 +3,7 @@
// c99 -o ew -Wall -I ../.. -O2 -ggdb -DLATTICESUMS32 hexlattice_ewald.c ../translations.c ../ewald.c ../ewaldsf.c ../gaunt.c ../lattices2d.c -lgsl -lm -lblas // c99 -o ew -Wall -I ../.. -O2 -ggdb -DLATTICESUMS32 hexlattice_ewald.c ../translations.c ../ewald.c ../ewaldsf.c ../gaunt.c ../lattices2d.c -lgsl -lm -lblas
#include <stdio.h> #include <stdio.h>
#include <math.h> #include <math.h>
#include <string.h>
#include <qpms/translations.h> #include <qpms/translations.h>
#include <qpms/lattices.h> #include <qpms/lattices.h>
#include <gsl/gsl_const_mksa.h> #include <gsl/gsl_const_mksa.h>
@ -68,7 +69,7 @@ int main (int argc, char **argv) {
size_t kcount = 0; size_t kcount = 0;
size_t klist_capacity = MAXKCOUNT; size_t klist_capacity = MAXKCOUNT;
cart2_t *klist = malloc(sizeof(cart2_t) * klist_capacity); cart2_t *klist = malloc(sizeof(cart2_t) * klist_capacity);
while (fscanf(f, "%lf %lf", &(klist[kcount].x), &(klist[kcount].y)) == 2) { while (scanf("%lf %lf", &(klist[kcount].x), &(klist[kcount].y)) == 2) {
++kcount; ++kcount;
if(kcount >= klist_capacity) { if(kcount >= klist_capacity) {
klist_capacity *= 2; klist_capacity *= 2;
@ -204,7 +205,7 @@ int main (int argc, char **argv) {
if(err) fclose(err); if(err) fclose(err);
#ifdef KSTDIN #ifdef KSTDIN
free(klist) free(klist);
#endif #endif
triangular_lattice_gen_free(Klg); triangular_lattice_gen_free(Klg);