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
RevLine 
[729]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
[746]11namespace SOPHYA {
[729]12
[866]13
[729]14template <class T>
15class SphericalTransformServer
16{
17
18 public:
19
[1218]20 SphericalTransformServer()
21 {
22 fftIntfPtr_=new FFTPackServer;
23 fftIntfPtr_->setNormalize(false);
24 };
25 ~SphericalTransformServer(){ if (fftIntfPtr_!=NULL) delete fftIntfPtr_;};
26
27 /*!
[729]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*/
[1218]30 void SetFFTServer(FFTServerInterface* srv=NULL)
31 {
32 if (fftIntfPtr_!=NULL) delete fftIntfPtr_;
33 fftIntfPtr_=srv;
34 fftIntfPtr_->setNormalize(false);
35 }
[729]36
[1218]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,
[729]41 const TVector<T>& Cl, const r_8 fwhm) const;
[1218]42 void GenerateFromCl(SphericalMap<T>& sphq, SphericalMap<T>& sphu,
[729]43 int_4 pixelSizeIndex,
44 const TVector<T>& Cle, const TVector<T>& Clb,
45 const r_8 fwhm) const;
46
47
[1683]48void DecomposeToAlm(SphericalMap<T>& map, Alm<T>& alm, int_4 nlmax, r_8 cos_theta_cut, int iterationOrder = 0) const;
[729]49
[1683]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;
[729]57
[1683]58TVector<T> DecomposeToCl(SphericalMap<T>& sph,
59 int_4 nlmax, r_8 cos_theta_cut, int iterationOrder = 0) const;
[729]60
61
[1218]62 private:
[1683]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,
[729]75 int_4 nph, r_8 phi0) const;
[1683]76TVector<T> RfourierSynthesisFromB(const Bm<complex<T> >& b_m,
[729]77 int_4 nph, r_8 phi0) const;
[1683]78TVector< complex<T> > CFromFourierAnalysis(int_4 mmax,
[729]79 const TVector<complex<T> > datain,
80 r_8 phi0) const;
[1683]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;
[729]87
[1218]88 void almFromWX(int_4 nlmax, int_4 nmmax, r_8 phi0,
[729]89 r_8 domega, r_8 theta,
90 const TVector<T>& dataq, const TVector<T>& datau,
91 Alm<T>& alme, Alm<T>& almb) const;
[1218]92 void almFromPM(int_4 nph, int_4 nlmax, int_4 nmmax,
[729]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
[1218]97 void mapFromWX(int_4 nlmax, int_4 nmmax,
[729]98 SphericalMap<T>& mapq, SphericalMap<T>& mapu,
99 const Alm<T>& alme, const Alm<T>& almb) const;
100
101
[1218]102 void mapFromPM(int_4 nlmax, int_4 nmmax,
[729]103 SphericalMap<T>& mapq, SphericalMap<T>& mapu,
104 const Alm<T>& alme, const Alm<T>& almb) const;
105
106
107
108 FFTServerInterface* fftIntfPtr_;
[1683]109
[729]110};
[746]111} // Fin du namespace
[729]112
113
114#endif
Note: See TracBrowser for help on using the repository browser.