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

Error on (temporarily) unsupported normalisations. 

Former-commit-id: b7ed7054cd4f717a430f9ef546f1d9967748edd7
This commit is contained in:
Marek Nečada 2019-07-22 11:22:21 +03:00
parent d01ab445fe
commit c6e666f4e6
1 changed files with 19 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

19
qpms/translations.c
View File

@ -59,6 +59,24 @@ static inline void TROPS_ONLY_EIMF_IMPLEMENTED(qpms_normalisation_t norm) {
QPMS_NOT_IMPLEMENTED("Translation operators for real or inverse complex spherical harmonics based waves are not implemented."); QPMS_NOT_IMPLEMENTED("Translation operators for real or inverse complex spherical harmonics based waves are not implemented.");
} }
// Use if only the symmetric form [[A, B], [B, A]] (without additional factors) of translation operator is allowed.
static inline void TROPS_ONLY_AB_SYMMETRIC_NORMS_IMPLEMENTED(qpms_normalisation_t norm) {
switch (norm & QPMS_NORMALISATION_NORM_BITS) {
case QPMS_NORMALISATION_NORM_SPHARM:
case QPMS_NORMALISATION_NORM_POWER:
break; // OK
default:
QPMS_NOT_IMPLEMENTED("Only spherical harmonic and power normalisation supported.");
}
if (
( !(norm & QPMS_NORMALISATION_N_I) != !(norm & QPMS_NORMALISATION_M_I) )
||
( !(norm & QPMS_NORMALISATION_N_MINUS) != !(norm & QPMS_NORMALISATION_M_MINUS) )
)
QPMS_NOT_IMPLEMENTED("Only normalisations without a phase factors between M and N waves are supported.");
}
/* /*
* References: * References:
* [Xu_old] Yu-Lin Xu, Journal of Computational Physics 127, 285–298 (1996) * [Xu_old] Yu-Lin Xu, Journal of Computational Physics 127, 285–298 (1996)
@ -1094,6 +1112,7 @@ qpms_errno_t qpms_trans_calculator_get_trans_array(const qpms_trans_calculator *
const qpms_vswf_set_spec_t *srcspec, size_t srcstride, const qpms_vswf_set_spec_t *srcspec, size_t srcstride,
sph_t kdlj, bool r_ge_d, qpms_bessel_t J) sph_t kdlj, bool r_ge_d, qpms_bessel_t J)
{ {
TROPS_ONLY_AB_SYMMETRIC_NORMS_IMPLEMENTED(c->normalisation);
assert(c->normalisation == destspec->norm && c->normalisation == srcspec->norm); assert(c->normalisation == destspec->norm && c->normalisation == srcspec->norm);
assert(c->lMax >= destspec->lMax && c->lMax >= srcspec->lMax); assert(c->lMax >= destspec->lMax && c->lMax >= srcspec->lMax);
assert(destspec->lMax_L < 0 && srcspec->lMax_L < 0); assert(destspec->lMax_L < 0 && srcspec->lMax_L < 0);