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

Fix k-argument related bugs etc. 

Former-commit-id: 72c955f31bcd1bfd9cd714c5b19d038f9c7ec6e3
This commit is contained in:
Marek Nečada 2020-01-20 17:30:36 +02:00
parent 71852aa017
commit a2a51d0de6
2 changed files with 3 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
qpms/qpms_c.pyx
View File

@ -622,7 +622,7 @@ cdef class _ScatteringSystemAtOmega:
cdef qpms_scatsys_at_omega_t *rawpointer(self): cdef qpms_scatsys_at_omega_t *rawpointer(self):
return self.ssw return self.ssw
def scatter_solver(self, double k, iri=None): def scatter_solver(self, iri=None):
self.check() self.check()
return ScatteringMatrix(self, iri) return ScatteringMatrix(self, iri)
@ -637,7 +637,7 @@ cdef class _ScatteringSystemAtOmega:
def modeproblem_matrix_full(self): def modeproblem_matrix_full(self):
self.check() self.check()
cdef size_t flen = self.s[0].fecv_size cdef size_t flen = self.ss_pyref.s[0].fecv_size
cdef np.ndarray[np.complex_t, ndim=2] target = np.empty( cdef np.ndarray[np.complex_t, ndim=2] target = np.empty(
(flen,flen),dtype=complex, order='C') (flen,flen),dtype=complex, order='C')
cdef cdouble[:,::1] target_view = target cdef cdouble[:,::1] target_view = target

2
qpms/scatsystem.c
View File

@ -1675,7 +1675,7 @@ complex double *qpms_scatsys_build_translation_matrix_e_irrep_packed(
pthread_mutex_t opistartR_mutex; pthread_mutex_t opistartR_mutex;
QPMS_ENSURE_SUCCESS(pthread_mutex_init(&opistartR_mutex, NULL)); QPMS_ENSURE_SUCCESS(pthread_mutex_init(&opistartR_mutex, NULL));
const struct qpms_scatsys_build_translation_matrix_e_irrep_packed_parallelR_thread_arg const struct qpms_scatsys_build_translation_matrix_e_irrep_packed_parallelR_thread_arg
arg = {ss, &opistartR, &opistartR_mutex, iri, target_packed, J}; arg = {ss, &opistartR, &opistartR_mutex, iri, target_packed, k, J};
// FIXME THIS IS NOT PORTABLE: // FIXME THIS IS NOT PORTABLE:
long nthreads; long nthreads;