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

Fix invalid pointer in qpms_apply_tmatrix. 

Former-commit-id: ba9400c3e1a39d472cfdebf7e61ab175c5c8fb6e
This commit is contained in:
Marek Nečada 2019-11-14 13:36:45 +02:00
parent 1821a0d8f5
commit c9a5661519
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/scatsystem.h
View File

@ -230,7 +230,7 @@ complex double *qpms_scatsys_irrep_unpack_vector(complex double *target_full,
/// Global translation matrix. /// Global translation matrix.
/** /**
* The diagonal (particle self-) block are filled with zeros. * The diagonal (particle self-) block are filled with zeros (even for regular Bessel waves).
* This may change in the future. * This may change in the future.
*/ */
complex double *qpms_scatsys_build_translation_matrix_full( complex double *qpms_scatsys_build_translation_matrix_full(

2
qpms/tmatrices.c
View File

@ -616,7 +616,7 @@ complex double *qpms_apply_tmatrix(
QPMS_CRASHING_CALLOC(f, n, sizeof(complex double)); QPMS_CRASHING_CALLOC(f, n, sizeof(complex double));
const complex double one = 1; const complex double one = 1;
const complex double zero = 0; const complex double zero = 0;
cblas_zgemv(CblasRowMajor, CblasNoTrans, n, n, &one, T->m, n, a, 1, &zero, NULL, 1); cblas_zgemv(CblasRowMajor, CblasNoTrans, n, n, &one, T->m, n, a, 1, &zero, f, 1);
return f; return f;
} }