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

Tag symmetries.c with non-implemented macros etc. 

Former-commit-id: 0cf1dc127b43bdded6471501a666ce325e5528f5
This commit is contained in:
Marek Nečada 2019-07-14 09:33:30 +03:00
parent 4c5cd88598
commit 45b8233862
1 changed files with 25 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

32
qpms/symmetries.c
View File

@ -2,38 +2,46 @@
#include "tiny_inlines.h" #include "tiny_inlines.h"
#include "indexing.h" #include "indexing.h"
#include "wigner.h" #include "wigner.h"
#include "qpms_error.h"
// TODO at some point, maybe support also other norms. // TODO at some point, maybe support also other norms.
// (perhaps just use qpms_normalisation_t_factor() at the right places) // (perhaps just use qpms_normalisation_t_factor() at the right places)
static inline void check_norm_compat(const qpms_vswf_set_spec_t *s) static inline void check_norm_compat(const qpms_vswf_set_spec_t *s)
{ {
switch (qpms_normalisation_t_normonly(s->norm)) { switch (s->norm & QPMS_NORMALISATION_NORM_BITS) {
case QPMS_NORMALISATION_POWER: case QPMS_NORMALISATION_NORM_POWER:
break; break;
case QPMS_NORMALISATION_SPHARM: case QPMS_NORMALISATION_NORM_SPHARM:
break; break;
default: default:
abort(); // Only SPHARM and POWER norms are supported right now. QPMS_NOT_IMPLEMENTED("At the moment, only spherical harmonics of spherical harmonics or power normalisations implemented.");
} }
} }
static inline void ONLY_EIMF_IMPLEMENTED(const qpms_normalisation_t norm)
{
if (norm & QPMS_NORMALISATION_SPHARM_REAL)
QPMS_NOT_IMPLEMENTED("Support for real spherical harmonics not implemented yet.");
}
// Used in the functions below to ensure memory allocation and checks for bspec validity // Used in the functions below to ensure memory allocation and checks for bspec validity
static inline complex double *ensure_alloc(complex double *target, static inline complex double *ensure_alloc(complex double *target,
const qpms_vswf_set_spec_t *bspec) { const qpms_vswf_set_spec_t *bspec) {
check_norm_compat(bspec); check_norm_compat(bspec);
const size_t n = bspec->n; const size_t n = bspec->n;
if (target == NULL) if (target == NULL)
target = malloc(n * n * sizeof(complex double)); QPMS_CRASHING_MALLOC(target, n * n * sizeof(complex double));
if (target == NULL) abort();
return target; return target;
} }
complex double *qpms_zflip_uvswi_dense( complex double *qpms_zflip_uvswi_dense(
complex double *target, complex double *target,
const qpms_vswf_set_spec_t *bspec) const qpms_vswf_set_spec_t *bspec)
{ {
check_norm_compat(bspec);
target = ensure_alloc(target, bspec); target = ensure_alloc(target, bspec);
const size_t n = bspec->n; const size_t n = bspec->n;
@ -69,6 +77,8 @@ complex double *qpms_yflip_uvswi_dense(
complex double *target, complex double *target,
const qpms_vswf_set_spec_t *bspec) const qpms_vswf_set_spec_t *bspec)
{ {
check_norm_compat(bspec);
ONLY_EIMF_IMPLEMENTED(bspec->norm);
target = ensure_alloc(target, bspec); target = ensure_alloc(target, bspec);
const size_t n = bspec->n; const size_t n = bspec->n;
@ -104,6 +114,8 @@ complex double *qpms_xflip_uvswi_dense(
complex double *target, complex double *target,
const qpms_vswf_set_spec_t *bspec) const qpms_vswf_set_spec_t *bspec)
{ {
check_norm_compat(bspec);
ONLY_EIMF_IMPLEMENTED(bspec->norm);
target = ensure_alloc(target, bspec); target = ensure_alloc(target, bspec);
const size_t n = bspec->n; const size_t n = bspec->n;
@ -142,6 +154,9 @@ complex double *qpms_zrot_uvswi_dense(
double phi ///< Rotation angle double phi ///< Rotation angle
) )
{ {
QPMS_UNTESTED; // not sure about the C.-S. phase. Don't forget documenting it as well.
check_norm_compat(bspec);
ONLY_EIMF_IMPLEMENTED(bspec->norm);
target = ensure_alloc(target, bspec); target = ensure_alloc(target, bspec);
const size_t n = bspec->n; const size_t n = bspec->n;
@ -182,6 +197,9 @@ complex double *qpms_irot3_uvswfi_dense(
const qpms_vswf_set_spec_t *bspec, const qpms_vswf_set_spec_t *bspec,
const qpms_irot3_t t) const qpms_irot3_t t)
{ {
QPMS_UNTESTED; // not sure about the C.-S. phase. Don't forget documenting it as well.
check_norm_compat(bspec);
ONLY_EIMF_IMPLEMENTED(bspec->norm);
target = ensure_alloc(target, bspec); target = ensure_alloc(target, bspec);
const size_t n = bspec->n; const size_t n = bspec->n;