source: Sophya/trunk/SophyaLib/Samba/sphericaltransformserver.h@ 1683

Last change on this file since 1683 was 1683, checked in by lemeur, 24 years ago

methode iterative pour analyse harmonique
File size: 3.3 KB
Line 
1#ifndef SPHERICALTRANFORMSERVER_SEEN
2#define SPHERICALTRANFORMSERVER_SEEN
3
4#include "sphericalmap.h"
5#include "fftservintf.h"
6#include "fftpserver.h"
7#include "alm.h"
8#include "lambdaBuilder.h"
9
10
11namespace SOPHYA {
12
13
14template <class T>
15class SphericalTransformServer
16{
17
18 public:
19
20 SphericalTransformServer()
21 {
22 fftIntfPtr_=new FFTPackServer;
23 fftIntfPtr_->setNormalize(false);
24 };
25 ~SphericalTransformServer(){ if (fftIntfPtr_!=NULL) delete fftIntfPtr_;};
26
27 /*!
28 Set a fft server. The constructor sets a default fft server (fft-pack). So it is not necessary to call this method for a standard use.
29*/
30 void SetFFTServer(FFTServerInterface* srv=NULL)
31 {
32 if (fftIntfPtr_!=NULL) delete fftIntfPtr_;
33 fftIntfPtr_=srv;
34 fftIntfPtr_->setNormalize(false);
35 }
36
37 void GenerateFromAlm( SphericalMap<T>& map, int_4 pixelSizeIndex, const Alm<T>& alm) const;
38 void GenerateFromAlm(SphericalMap<T>& mapq, SphericalMap<T>& mapu, int_4 pixelSizeIndex, const Alm<T>& alme, const Alm<T>& almb) const;
39
40 void GenerateFromCl(SphericalMap<T>& sph, int_4 pixelSizeIndex,
41 const TVector<T>& Cl, const r_8 fwhm) const;
42 void GenerateFromCl(SphericalMap<T>& sphq, SphericalMap<T>& sphu,
43 int_4 pixelSizeIndex,
44 const TVector<T>& Cle, const TVector<T>& Clb,
45 const r_8 fwhm) const;
46
47
48void DecomposeToAlm(SphericalMap<T>& map, Alm<T>& alm, int_4 nlmax, r_8 cos_theta_cut, int iterationOrder = 0) const;
49
50void DecomposeToAlm(SphericalMap<T>& mapq,
51 SphericalMap<T>& mapu,
52 Alm<T>& alme,
53 Alm<T>& almb,
54 int_4 nlmax,
55 r_8 cos_theta_cut,
56 int iterationOrder = 0) const;
57
58TVector<T> DecomposeToCl(SphericalMap<T>& sph,
59 int_4 nlmax, r_8 cos_theta_cut, int iterationOrder = 0) const;
60
61
62 private:
63
64void carteVersAlm(const SphericalMap<T>& map, int_4 nlmax, r_8 cos_theta_cut, Alm<T>& alm) const;
65
66void carteVersAlm(const SphericalMap<T>& mapq,
67 const SphericalMap<T>& mapu,
68 Alm<T>& alme,
69 Alm<T>& almb,
70 int_4 nlmax,
71 r_8 cos_theta_cut) const;
72
73
74TVector< complex<T> > fourierSynthesisFromB(const Bm<complex<T> >& b_m,
75 int_4 nph, r_8 phi0) const;
76TVector<T> RfourierSynthesisFromB(const Bm<complex<T> >& b_m,
77 int_4 nph, r_8 phi0) const;
78TVector< complex<T> > CFromFourierAnalysis(int_4 mmax,
79 const TVector<complex<T> > datain,
80 r_8 phi0) const;
81 // TVector< complex<T> > CFromFourierAnalysis(int_4 mmax,
82 // const TVector<T> datain,
83 // r_8 phi0) const;
84 void CFromFourierAnalysis(int_4 mmax, const TVector<T> datain,
85 TVector< complex<T> >& dataout,
86 r_8 phi0) const;
87
88 void almFromWX(int_4 nlmax, int_4 nmmax, r_8 phi0,
89 r_8 domega, r_8 theta,
90 const TVector<T>& dataq, const TVector<T>& datau,
91 Alm<T>& alme, Alm<T>& almb) const;
92 void almFromPM(int_4 nph, int_4 nlmax, int_4 nmmax,
93 r_8 phi0, r_8 domega, r_8 theta,
94 const TVector<T>& dataq, const TVector<T>& datau,
95 Alm<T>& alme, Alm<T>& almb) const;
96
97 void mapFromWX(int_4 nlmax, int_4 nmmax,
98 SphericalMap<T>& mapq, SphericalMap<T>& mapu,
99 const Alm<T>& alme, const Alm<T>& almb) const;
100
101
102 void mapFromPM(int_4 nlmax, int_4 nmmax,
103 SphericalMap<T>& mapq, SphericalMap<T>& mapu,
104 const Alm<T>& alme, const Alm<T>& almb) const;
105
106
107
108 FFTServerInterface* fftIntfPtr_;
109
110};
111} // Fin du namespace
112
113
114#endif
Note: See TracBrowser for help on using the repository browser.