ssht_dl.h

Go to the documentation of this file.

// SSHT package to perform spin spherical harmonic transforms
// Copyright (C) 2011  Jason McEwen
// See LICENSE.txt for license details
 
 
#ifndef SSHT_DL
#define SSHT_DL
 
 
#ifdef __cplusplus
extern "C"{
#endif
 
typedef enum {SSHT_DL_QUARTER = 0, 
              SSHT_DL_QUARTER_EXTENDED, 
              SSHT_DL_HALF, 
              SSHT_DL_FULL} ssht_dl_size_t;
 
typedef enum {SSHT_DL_RISBO = 0, 
              SSHT_DL_TRAPANI} ssht_dl_method_t;
 
 
double* ssht_dl_calloc(int L, ssht_dl_size_t dl_size);
int ssht_dl_get_offset(int L, ssht_dl_size_t dl_size);
int ssht_dl_get_stride(int L, ssht_dl_size_t dl_size);
 
void ssht_dl_beta_risbo_full_table(double *dl, double beta, int L, 
                                   ssht_dl_size_t dl_size,
                                   int el, double *sqrt_tbl);
void ssht_dl_beta_risbo_half_table(double *dl, double beta, int L, 
                                                                           ssht_dl_size_t dl_size,
                                                                           int el, double *sqrt_tbl,
                                                                           double *signs);
void ssht_dl_beta_risbo_eighth_table(double *dl, double beta, int L, 
                                     ssht_dl_size_t dl_size,
                                      int el, double *sqrt_tbl,
                                      double *signs);
void ssht_dl_beta_risbo_fill_eighth2quarter_table(double *dl, 
                                                  double *dl8,
                                                  int L,
                                                  ssht_dl_size_t dl_size,
                                                  ssht_dl_size_t dl8_size,
                                                  int el, 
                                                  double *signs);
 
void ssht_dl_beta_kostelec_full_table(double *dlm1p1, double *dl, 
                                 double beta, int L, 
                                 ssht_dl_size_t dl_size,
                                 int el, 
                                 double *sqrt_tbl, double *signs);
void ssht_dl_beta_kostelec_line_table(double *dlm1p1_line, double *dl_line, 
                                      double beta, int L, int mm, int el, 
                                      double *sqrt_tbl, double *signs);
void ssht_dl_beta_kostelec_halfline_table(double *dlm1p1_line, double *dl_line, 
                                          double beta, int L, int mm, int el, 
                                          double *sqrt_tbl, double *signs);
 
void ssht_dl_halfpi_trapani_eighth_table(double *dl, int L, 
                                         ssht_dl_size_t dl_size,
                                         int el, double *sqrt_tbl);
void ssht_dl_halfpi_trapani_quarter_table(double *dl, int L, 
                                          ssht_dl_size_t dl_size,
                                          int el, double *sqrt_tbl);
void ssht_dl_halfpi_trapani_fill_eighth2righthalf_table(double *dl, int L,
                                                        ssht_dl_size_t dl_size,
                                                        int el, double *signs);
void ssht_dl_halfpi_trapani_fill_eighth2quarter_table(double *dl, int L,
                                                     ssht_dl_size_t dl_size,
                                                     int el, double *signs);
 
#ifdef __cplusplus
}
#endif
 
#endif