weights.c

Go to the documentation of this file.

 
#include "s2kit/weights.h"
 
#include <math.h>
 
void GenerateWeightsForDLT(const int bw, double* weights) {
    double coeff = M_PI / (4. * bw);
 
    for (int i = 0; i < 2 * bw; ++i) {
        double sum = 0.;
        double k = 2. * i + 1.;
 
        for (int j = 0; j < bw; ++j)
            sum += 1. / (2. * j + 1.) * sin(k * (2. * j + 1.) * coeff);
 
        sum *= 2. * sin(k * coeff) / bw;
 
        weights[i] = sum;
        weights[i + 2 * bw] = sum * sin(k * coeff);
    }
}