Changeset 2198 in Sophya for trunk/SophyaLib/SkyMap/localmap.cc

Timestamp:

Sep 24, 2002, 4:05:46 PM (23 years ago)

Author:

lemeur

Message:

toilette d'ete

File:

• trunk/SophyaLib/SkyMap/localmap.cc (modified) (15 diffs)

trunk/SophyaLib/SkyMap/localmap.cc

@@ -65 +65 @@
 }
 
+//! partial constructor of the local map
+/*!
+  \param nx : number of pixels in x direction
+  \param ny : number of pixels in y direction
+
+must be followed by calling the methodsa setSize() and setOrigin() in order the object become usable
+
+ */
+
+template<class T>
+LocalMap<T>::LocalMap(int_4 nx, int_4 ny)  
+{
+  InitNul();
+  ReSize(nx, ny);
+  //  nSzX_ = nx;
+  //  nSzY_ = ny;
+  //  nPix_= nx*ny;
+  //  pixels_.ReSize(ny,nx);
+}
+
 //! full constructor of the local map
 /*!
@@ -82 +102 @@
 {
   InitNul();
-  nSzX_ = nx;
-  nSzY_ = ny;
-  nPix_= nx*ny;
-  pixels_.ReSize(ny,nx);
+  //  nSzX_ = nx;
+  // nSzY_ = ny;
+  // nPix_= nx*ny;
+  // pixels_.ReSize(ny,nx);
+  ReSize(nx, ny);
   SetSize(angleX, angleY);
   SetOrigin(theta0, phi0, x0, y0, angle);
@@ -107 +128 @@
 {
   InitNul();
+  //  nSzX_ = nx;
+  //  nSzY_ = ny;
+  // nPix_= nx*ny;
+  //  pixels_.ReSize(ny,nx);
+  ReSize(nx, ny);
+  SetSize(angleX, angleY);
+  SetOrigin(theta0, phi0, angle);
+}
+
+//! copy constructor 
+/*!
+  \param share : if true, share data. If false copy data
+ */
+template<class T>
+LocalMap<T>::LocalMap(const LocalMap<T>& lm, bool share)
+  : pixels_(lm.pixels_, share)
+{
+
+  if(lm.mInfo_) mInfo_= new DVList(*lm.mInfo_);
+  recopierVariablesSimples(lm);
+}
+
+//! copy constructor 
+/*!
+  \warning datas are \b SHARED with \b lm.
+  \sa TMatrix::TMatrix(const TMatrix<T>&)
+*/
+template<class T>
+LocalMap<T>::LocalMap(const LocalMap<T>& lm)
+  : pixels_(lm.pixels_)
+
+{
+
+  if(lm.mInfo_) mInfo_= new DVList(*lm.mInfo_);
+  recopierVariablesSimples(lm);
+}
+
+template<class T>
+LocalMap<T>::~LocalMap() {;}
+
+
+
+/*!  \fn void SOPHYA::LocalMap::ReSize(int_4 nx, int_4 ny) 
+
+  Resize storage area for pixels 
+*/
+template<class T>
+void LocalMap<T>::ReSize(int_4 nx, int_4 ny) 
+{
+    // angles par pixel,  en radians 
+  deltaPhi_   *=  nSzX_/nx;   
+  deltaTheta_ *=  nSzY_/ny;
+
   nSzX_ = nx;
   nSzY_ = ny;
   nPix_= nx*ny;
   pixels_.ReSize(ny,nx);
-  cout << " tableau pixels initialise " << endl;
-  SetSize(angleX, angleY);
-  cout << " fin du setsize " << endl;
-  SetOrigin(theta0, phi0, angle);
-  cout << " fin du set oorigine " << endl;
-}
-
-//! copy constructor 
-/*!
-  \param share : if true, share data. If false copy data
- */
-template<class T>
-LocalMap<T>::LocalMap(const LocalMap<T>& lm, bool share)
-  : pixels_(lm.pixels_, share)
-{
-
-  if(lm.mInfo_) mInfo_= new DVList(*lm.mInfo_);
-  recopierVariablesSimples(lm);
-}
-
-//! copy constructor 
-/*!
-  \warning datas are \b SHARED with \b lm.
-  \sa TMatrix::TMatrix(const TMatrix<T>&)
-*/
-template<class T>
-LocalMap<T>::LocalMap(const LocalMap<T>& lm)
-  : pixels_(lm.pixels_)
-
-{
-
-  if(lm.mInfo_) mInfo_= new DVList(*lm.mInfo_);
-  recopierVariablesSimples(lm);
-}
-
-template<class T>
-LocalMap<T>::~LocalMap() {;}
-
-
-
-/*!  \fn void SOPHYA::LocalMap::ReSize(int_4 nx, int_4 ny) 
-
-  Resize storage area for pixels 
-*/
-template<class T>
-void LocalMap<T>::ReSize(int_4 nx, int_4 ny) 
-{
-    // angles par pixel,  en radians 
-  deltaPhi_   *=  nSzX_/nx;   
-  deltaTheta_ *=  nSzY_/ny;
-
-  nSzX_ = nx;
-  nSzY_ = ny;
-  nPix_= nx*ny;
-  pixels_.ReSize(ny,nx);
+  pixelSized_ = true;
 }
 
@@ -285 +305 @@
   int_4 i,j;
 
+  if (! LocalMap_isDone()) initializationError();
   double csTheta = cos(theta);
   double snTheta = sin(theta);
@@ -296 +317 @@
   if ( denom == 0.) 
     {
-      throw PException("LocalMap::PixIndexSph : vector radius parallel to map plane ");
     }
   double lambda = 1./denom;
@@ -428 +448 @@
 void LocalMap<T>::PixToSphereC(int_4 ip, int_4 it, double& XP, double& YP, double& ZP) const
 {
+  if (! LocalMap_isDone()) initializationError();
   double X= double(ip-x0_)* deltaPhi_;
   double Y= double(it-y0_)*deltaTheta_; 
@@ -466 +487 @@
 void LocalMap<T>::SetOrigin(double theta0,double phi0,int_4 x0,int_4 y0,double angle)
 {
-  // if (originFlag_)
-  //   {
-  //     throw PException("LocalMap::SetOrigin : redefining origin is not allowed at the moment ");
-  //  }
 
   if ( theta0 < 0. || theta0 > 180. || phi0 < 0. || phi0 > 360. || angle < -90. || angle > 90.)
@@ -483 +500 @@
   cosAngle_= cos(angle_);
   sinAngle_= sin(angle_);
+
+   if (originSet_)
+     {
+       cout << " WARNING:  LocalMap<T>::SetOrigin() : reset origine of a local map which had already an origin " << endl;
+     }
+   originSet_ = true;
 } 
 
@@ -510 +533 @@
 {
   cout << " calcul matrice de passage " << endl;
+  if (! LocalMap_isDone()) initializationError();
   TMatrix<double> passage(3,3);
   double cos_th_axeZ;
@@ -578 +602 @@
 void LocalMap<T>::SetSize(double angleX,double angleY)
 {
+  if (!pixelSized_) initializationError();
+
   if ( angleX <= 0. || angleX > 180. || angleY <= 0. || angleY > 180.)
     {
@@ -589 +615 @@
   deltaPhi_   = angleX*Pi/(180.*nSzX_);
   deltaTheta_ = angleY*Pi/(180.*nSzY_);
-
+ 
+  if (angleSized_)
+     {
+       cout << " WARNING LocalMap::SetSize  ; angle resizing of an already sized local map " << endl;
+     }
+  angleSized_  = true;
 
 }
@@ -624 +655 @@
 //--
 { 
+
+  pixelSized_  = false;
+  angleSized_  = false;
+  originSet_ = false;
+
   nSzX_ = 0;
   nSzY_ = 0;
@@ -635 +671 @@
 
 
+
   nPix_ = 0;
 
@@ -699 +736 @@
 {
 
+  pixelSized_ = lm.pixelSized_;
+  angleSized_ = lm.angleSized_;
+  originSet_  = lm.originSet_ ;
+
+
   nSzX_ = lm.nSzX_;
   nSzY_ = lm.nSzY_;
@@ -721 +763 @@
   deltaTheta_ = lm.deltaTheta_;
 }
+
+template<class T>
+void LocalMap<T>::initializationError() const
+{
+  cout << endl;
+  if (!pixelSized_) cout << " LocalMap::initializationError() ; pixel sizing is missing " << endl;
+  if (!angleSized_ ) cout << " LocalMap::initializationError() ; angle sizing is missing " << endl;
+  if (!originSet_ ) cout << " LocalMap::initializationError() ; origin of map is not set " << endl;
+      throw PException("LocalMap::initializationError() : bad initialization of the local map ");
+
+
+}
+
 
 //   ...... Operations de calcul  ......