← Back to branch summary

~ubuntu-branches/ubuntu/vivid/grass/vivid-proposed

« back to all changes in this revision

Viewing changes to imagery/i.atcorr/Gauss.h

  • Committer: Package Import Robot
  • Author(s): Bas Couwenberg
  • Date: 2015-02-20 23:12:08 UTC
  • mfrom: (8.2.6 experimental)
  • Revision ID: package-import@ubuntu.com-20150220231208-1u6qvqm84v430b10
Tags: 7.0.0-1~exp1
* New upstream release.
* Update python-ctypes-ternary.patch to use if/else instead of and/or.
* Drop check4dev patch, rely on upstream check.
* Add build dependency on libpq-dev to grass-dev for libpq-fe.h.
* Drop patches applied upstream, refresh remaining patches.
* Update symlinks for images switched from jpg to png.

Show diffs side-by-side

added added

removed removed

Lines of Context:
imagery/i.atcorr/Gauss.h
1
 
#ifndef MY_GAUSS_H
2
 
#define MY_GAUSS_H
3
 
 
4
 
const long int mu       = 25;
5
 
const long int np       = 49;
6
 
 
7
 
/* index standard interval [-1,1] */
8
 
#define STDI(X) ((X)+mu)
9
 
 
10
 
struct Gauss
11
 
{
12
 
private:
13
 
        static float angmu[10];
14
 
        static float angphi[13];
15
 
 
16
 
public:
17
 
        /* [a,b] = [0,2*Pi] */
18
 
        float rp[np];                   /* gaussian angles */
19
 
        float gp[np];                   /* gaussian weights */
20
 
 
21
 
        // [a,b] = [-1,1]
22
 
        float rm[2*mu+1];               /* shifted gaussian angles */
23
 
        float gb[2*mu+1];               /* shifted gaussian weights */
24
 
                                        /* with the ends zeroed as well as the center */
25
 
                                        /* [0 ? ? ? ? 0 ? ? ? ? 0] */
26
 
 
27
 
    /*  preliminary computations for gauss integration */
28
 
        void init();
29
 
 
30
 
        /*      Compute for a given n, the gaussian quadrature (the n gaussian angles and the
31
 
        their respective weights). The gaussian quadrature is used in numerical integration involving the
32
 
        cosine of emergent or incident direction zenith angle. */
33
 
        static void gauss (float a, float b, float *x, float *w, long int n);
34
 
};
35
 
 
36
 
#endif /* MY_GAUSS_H */