Changeset 3510 in Sophya for trunk/SophyaLib/Samba/sphericaltransformserver.cc

Timestamp:

Aug 8, 2008, 3:11:45 PM (17 years ago)

Author:

ansari

Message:

Nettoyage alm.h .cc , utilisation RandomGenerator par Alm<T> et SphericalTransformServer , Reza 08/08/2008

File:

• trunk/SophyaLib/Samba/sphericaltransformserver.cc (modified) (4 diffs)

trunk/SophyaLib/Samba/sphericaltransformserver.cc

@@ -6 +6 @@
 #include "sphericaltransformserver.h"
 #include "tvector.h"
-#include "srandgen.h"
+#include "stsrand.h"
 #include "nbmath.h"
 #include "timing.h"
@@ -152 +152 @@
 
  */
+
+//!  Default constructor - Creates a non thread-safe RandomGenerator to be used by GenerateFromCl 
+template<class T>
+SphericalTransformServer<T>::SphericalTransformServer()
+: rg_(1, false) 
+{
+  fftIntfPtr_=new FFTPackServer(true); // preserveinput = true
+  fftIntfPtr_->setNormalize(false);
+}
+
+//!  Constructor with the specification of a RandomGenerator object to be used by GenerateFromCl 
+template<class T>
+SphericalTransformServer<T>::SphericalTransformServer(RandomGenerator const & rg)
+: rg_(rg) 
+{
+  fftIntfPtr_=new FFTPackServer(true); // preserveinput = true
+  fftIntfPtr_->setNormalize(false);
+}
+
+template<class T>
+SphericalTransformServer<T>::~SphericalTransformServer()
+{
+  if (fftIntfPtr_!=NULL) delete fftIntfPtr_; 
+}
+
+/*!
+ 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.
+ \warning The FFTServerInterface object should NOT overwrite the input arrays
+*/
+template<class T>
+void SphericalTransformServer<T>::SetFFTServer(FFTServerInterface* srv) 
+{
+  if (fftIntfPtr_!=NULL) delete fftIntfPtr_;
+  fftIntfPtr_=srv;
+  fftIntfPtr_->setNormalize(false);
+}
+
 
  /*! \fn void SOPHYA::SphericalTransformServer::GenerateFromAlm( SphericalMap<T>& map, int_4 pixelSizeIndex, const Alm<T>& alm) const
@@ -1275 +1313 @@
   //  Alm<T> a2lme = almFromCl(Cle, fwhm);
   // Alm<T> a2lmb = almFromCl(Clb, fwhm);
-  Alm<T> a2lme(Cle, fwhm);
-  Alm<T> a2lmb(Clb, fwhm);
+  Alm<T> a2lme(Cle, fwhm, rg_);
+  Alm<T> a2lmb(Clb, fwhm, rg_);
 
   GenerateFromAlm(sphq,sphu,pixelSizeIndex,a2lme,a2lmb); 
@@ -1293 +1331 @@
 {
 
-  Alm<T> alm(Cl, fwhm);
+  Alm<T> alm(Cl, fwhm, rg_);
   GenerateFromAlm(sph,pixelSizeIndex, alm ); 
 }