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

Fix fatal typo in translations.c 

Former-commit-id: 54d8ff3b8cd8f443584bd567e8ba068539e5656f
This commit is contained in:
Marek Nečada 2018-10-01 17:03:56 +03:00
parent 533c3f70bb
commit 8507b09332
2 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
qpms/tmatrices.py
View File

@ -229,7 +229,7 @@ def apply_ndmatrix_right(tensor, matrix, axes):
Multiplies a tensor with a 2N-dimensional matrix, conserving the axis order.
"""
N = len(axes)
matrix = np.tensordot(tensor, matrix, axes = (range(N), axes))
matrix = np.tensordot(tensor, matrix, axes = (axes, range(N)))
matrix = np.moveaxis(matrix, [-N+axn for axn in range(N)], axes)
return matrix

2
qpms/translations.c
View File

@ -1241,7 +1241,7 @@ int qpms_trans_calculator_get_AB_arrays_e32_both_points_and_shift(const qpms_tra
for(qpms_l_t q = 0; q <= qmax; ++q) {
const qpms_l_t p_ = n + nu - 2*q + 1;
const qpms_y_t y_sc = qpms_mn2y_sc(mu_m, p_);
const complex double multiplier = c->A_multipliers[i][q-BQ_OFFSET];
const complex double multiplier = c->B_multipliers[i][q-BQ_OFFSET];
complex double sigma = sigmas_total[y_sc];
ckahanadd(&Bsum, &Bsumc, multiplier * sigma);
if (Berr) kahanadd(&Bsumerr, &Bsumerrc, multiplier * serr_total[y_sc]);