qpms/test_vswf_pitau.c

   1 #include "vswf.h"
   2 #include "indexing.h"
   3 #include <stdio.h>
   4 #include <gsl/gsl_math.h>
   5 const double dtheta = 0.01 * M_PI;
   6 const qpms_l_t lMax = 3;
   7 #ifdef CS
   8 #define CS_OFFSET QPMS_NORMALISATION_T_CSBIT
   9 #else
  10 #define CS_OFFSET 0
  11 #endif
  12 int main() {
  13         qpms_y_t nelem = qpms_lMax2nelem(lMax);
  14         for (double theta = 0.; theta <= M_PI; theta += dtheta) {
  15                 printf("%.5e %.5e ", theta, cos(theta));
  16                 for(qpms_normalisation_t norm = CS_OFFSET + 1; norm <= CS_OFFSET +  3; ++norm) {//fujka :D
  17                         qpms_pitau_t pt = qpms_pitau_get(theta, lMax, norm);
  18                         for (qpms_y_t y = 0; y < nelem; ++y)
  19                                 printf("%.5e %.5e %.5e ", pt.leg[y], pt.pi[y], pt.tau[y]);
  20                         qpms_pitau_free(pt);
  21                 }
  22                 printf("\n");
  23         }
  24 }
  25
  26
  27
  28
  29