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

Correct orbit particle index type 

Former-commit-id: bf7c240d62ea661e22e320cb26be7ee06a058832
This commit is contained in:
Marek Nečada 2019-08-03 12:51:27 +03:00
parent 8baf3a259c
commit 3f25161e21
1 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
qpms/scatsystem.c
View File

@ -461,9 +461,9 @@ complex double *qpms_orbit_action_matrix(complex double *target,
memset(target, 0, n*n*N*N*sizeof(complex double)); memset(target, 0, n*n*N*N*sizeof(complex double));
complex double tmp[n][n]; // this is the 'per-particle' action complex double tmp[n][n]; // this is the 'per-particle' action
qpms_irot3_uvswfi_dense(tmp[0], bspec, sym->rep3d[g]); qpms_irot3_uvswfi_dense(tmp[0], bspec, sym->rep3d[g]);
for(qpms_gmi_t Col = 0; Col < ot->size; ++Col) { for(qpms_ss_orbit_pi_t Col = 0; Col < ot->size; ++Col) {
// Row is the 'destination' of the symmetry operation, Col is the 'source' // Row is the 'destination' of the symmetry operation, Col is the 'source'
const qpms_gmi_t Row = ot->action[sym->order * Col + g]; const qpms_ss_orbit_pi_t Row = ot->action[sym->order * Col + g];
for(size_t row = 0; row < bspec->n; ++row) for(size_t row = 0; row < bspec->n; ++row)
for(size_t col = 0; col < bspec->n; ++col) for(size_t col = 0; col < bspec->n; ++col)
target[n*n*N*Row + n*Col + n*N*row + col] = conj(tmp[row][col]); //CHECKCONJ target[n*n*N*Row + n*Col + n*N*row + col] = conj(tmp[row][col]); //CHECKCONJ