Changeset 980 in Sophya for trunk/SophyaLib/SkyMap/spherethetaphi.cc

Timestamp:

Apr 28, 2000, 4:12:44 PM (25 years ago)

Author:

ansari

Message:

modifs constructeurs copie et operateur =

File:

• trunk/SophyaLib/SkyMap/spherethetaphi.cc (modified) (6 diffs)

trunk/SophyaLib/SkyMap/spherethetaphi.cc

@@ -39 +39 @@
 template <class T>
 SphereThetaPhi<T>::SphereThetaPhi()
+  : NPhi_(), TNphi_(), Theta_(), pixels_() 
 
 //--
 {
   InitNul();
-  pixels_.Reset();
 }
 
@@ -67 +67 @@
     pixels_(s.pixels_ , share) 
 {
-  if(s.mInfo_) mInfo_= new DVList(*s.mInfo_);
 
   NTheta_= s.NTheta_;
   NPix_  = s.NPix_;
   Omega_ = s.Omega_;
+  if(s.mInfo_) mInfo_= new DVList(*s.mInfo_);
 }
 
@@ -78 +78 @@
   : NPhi_(s.NPhi_), TNphi_(s.TNphi_), Theta_(s.Theta_), pixels_(s.pixels_) 
 {
-  if(s.mInfo_) mInfo_= new DVList(*s.mInfo_);
 
   NTheta_= s.NTheta_;
   NPix_  = s.NPix_;
   Omega_ = s.Omega_;
+  if(s.mInfo_) mInfo_= new DVList(*s.mInfo_);
 }
 
@@ -122 +122 @@
 {
 
-    NPhi_.CloneOrShare(a.NPhi_);
-    TNphi_.CloneOrShare(a.TNphi_);
-    Theta_.CloneOrShare(a.Theta_);
-    pixels_.CloneOrShare(a.pixels_);
-}
-
+  NTheta_= a.NTheta_;
+  NPix_  = a.NPix_;
+  Omega_ = a.Omega_;
+  NPhi_.CloneOrShare(a.NPhi_);
+  TNphi_.CloneOrShare(a.TNphi_);
+  Theta_.CloneOrShare(a.Theta_);
+  pixels_.CloneOrShare(a.pixels_);
+}
+
+////////////////////////// methodes de copie/share
 template<class T>
 SphereThetaPhi<T>& SphereThetaPhi<T>::Set(const SphereThetaPhi<T>& a)
-   {
+{
   if (this != &a) 
     { 
-
-
-      NTheta_= a.NTheta_;
-      NPix_  = a.NPix_;
-      Omega_ = a.Omega_;
-      CloneOrShare(a);
+      
+      
+      if (a.NbPixels() < 1) 
+        throw RangeCheckError("SphereThetaPhi<T>::Set(a ) - Array a not allocated ! ");
+      if (NbPixels() < 1) CloneOrShare(a);
+      else CopyElt(a);
       if (mInfo_) delete mInfo_;
       mInfo_ = NULL;
@@ -144 +148 @@
     }
   return(*this);
-   }
-
-
+}
+
+template<class T>
+SphereThetaPhi<T>& SphereThetaPhi<T>::CopyElt(const SphereThetaPhi<T>& a)
+{
+  if (NbPixels() < 1) 
+    throw RangeCheckError("SphereThetaPhi<T>::CopyElt(const SphereThetaPhi<T>& )  - Not Allocated Array ! ");
+  if (NbPixels() != a.NbPixels()) 
+    throw(SzMismatchError("SphereThetaPhi<T>::CopyElt(const SphereThetaPhi<T>&) SizeMismatch")) ;
+
+  NTheta_= a.NTheta_;
+  NPix_  = a.NPix_;
+  Omega_ = a.Omega_;
+  int k;
+  for (k=0; k< NPix_; k++) pixels_(k) = a.pixels_(k);
+  for (k=0; k< a.NPhi_.Size(); k++) NPhi_(k) = a.NPhi_(k);
+  for (k=0; k< a.TNphi_.Size(); k++) TNphi_(k) = a.TNphi_(k);
+  for (k=0; k< a.Theta_.Size(); k++) Theta_(k) = a.Theta_(k);
+  return(*this);
+
+}
 
 /* --Methode-- */
@@ -524 +546 @@
   // Le nombre de pixels en phi de chacune des bandes en theta
   //  NPhi_ = new int_4[m];
-  NPhi_.ReSize(m);
+  // une taille de m suffit, mais je mets m+1 pour que les 3 tableaux aient 
+  // la meme taille pour une manipulation plus faciles par la librairie 
+  // cfitsio -- GLM (13-04-00)
+  NPhi_.ReSize(m+1);
   
   // Le nombre/Deuxpi total des pixels contenus dans les bandes de z superieur a une