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

Timestamp:

Feb 23, 2001, 12:26:48 PM (25 years ago)

Author:

lemeur

Message:

surcharge d'operateurs =, +=, *= etc...

File:

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

trunk/SophyaLib/SkyMap/spherethetaphi.cc

@@ -55 +55 @@
 //    m is the number of slices in theta on an hemisphere (the polar cap
 //    forms the first slice).
-//    pet is a dummy parameter at the moment. 
 //--
 {
@@ -129 +128 @@
   Theta_.CloneOrShare(a.Theta_);
   pixels_.CloneOrShare(a.pixels_);
+  if (mInfo_) {delete mInfo_; mInfo_ = NULL;}
+  if (a.mInfo_) mInfo_ = new DVList(*(a.mInfo_));
+}
+template<class T>
+void  SphereThetaPhi<T>::Share(const  SphereThetaPhi<T>& a)
+{
+
+  NTheta_= a.NTheta_;
+  NPix_  = a.NPix_;
+  Omega_ = a.Omega_;
+  NPhi_.Share(a.NPhi_);
+  TNphi_.Share(a.TNphi_);
+  Theta_.Share(a.Theta_);
+  pixels_.Share(a.pixels_);
+  if (mInfo_) {delete mInfo_; mInfo_ = NULL;}
+  if (a.mInfo_) mInfo_ = new DVList(*(a.mInfo_));
 }
 
@@ -710 +725 @@
   os << endl;
 }
+
+//   ...... Operations de calcul  ......
+
+
+//! Fill a SphereThetaPhi with a constant value \b a
+template <class T>
+SphereThetaPhi<T>& SphereThetaPhi<T>::SetT(T a) 
+{
+  if (NbPixels() < 1) 
+    throw RangeCheckError("SphereThetaPhi<T>::SetT(T )  - SphereThetaPhi not dimensionned ! ");
+  pixels_ = a;
+  return (*this);
+}
+
+/*! Add a constant value \b x to a SphereThetaPhi */
+template <class T>
+SphereThetaPhi<T>& SphereThetaPhi<T>::Add(T a)
+ {
+  if (NbPixels()< 1) 
+    throw RangeCheckError("SphereThetaPhi<T>::Add(T )  - SphereThetaPhi not dimensionned ! ");
+  pixels_ += a; 
+  return (*this);
+}
+
+/*! Substract a constant value \b a to a SphereThetaPhi */
+template <class T>
+SphereThetaPhi<T>& SphereThetaPhi<T>::Sub(T a) 
+{
+  if (NbPixels()< 1) 
+    throw RangeCheckError("SphereThetaPhi<T>::Sub(T )  - SphereThetaPhi not dimensionned ! ");
+  pixels_ -= a; 
+  return (*this);
+} 
+
+/*! multiply a SphereThetaPhi by a constant value \b a */
+template <class T>
+SphereThetaPhi<T>& SphereThetaPhi<T>::Mul(T a) 
+{
+  if (NbPixels()< 1) 
+    throw RangeCheckError("SphereThetaPhi<T>::Mul(T )  - SphereThetaPhi not dimensionned ! ");
+  pixels_ *= a; 
+  return (*this);
+}
+
+/*! divide a SphereThetaPhi by a constant value \b a */
+template <class T>
+SphereThetaPhi<T>& SphereThetaPhi<T>::Div(T a) 
+{
+  if (NbPixels()< 1) 
+    throw RangeCheckError("SphereThetaPhi<T>::Div(T )  - SphereThetaPhi not dimensionned ! ");
+  pixels_ /= a; 
+  return (*this);
+} 
+
+//  >>>> Operations avec 2nd membre de type SphereThetaPhi
+//! Add two SphereThetaPhi
+
+template <class T>
+SphereThetaPhi<T>& SphereThetaPhi<T>::AddElt(const SphereThetaPhi<T>& a)
+{
+  if (NbPixels()!= a.NbPixels())
+    {
+    throw(SzMismatchError("SphereThetaPhi<T>::AddElt(const SphereThetaPhi<T>&) SizeMismatch")) ;
+    }
+  pixels_ += a.pixels_;
+  return (*this);
+}
+
+//! Substract two SphereThetaPhi
+template <class T>
+SphereThetaPhi<T>& SphereThetaPhi<T>::SubElt(const SphereThetaPhi<T>& a)
+{
+  if (NbPixels()!= a.NbPixels())
+    {
+    throw(SzMismatchError("SphereThetaPhi<T>::SubElt(const SphereThetaPhi<T>&) SizeMismatch")) ;
+    }
+  pixels_ -= a.pixels_;
+  return (*this);
+}
+
+//! Multiply two SphereThetaPhi (elements by elements)
+template <class T>
+SphereThetaPhi<T>& SphereThetaPhi<T>::MulElt(const SphereThetaPhi<T>& a)
+{
+  if (NbPixels()!= a.NbPixels())
+    {
+    throw(SzMismatchError("SphereThetaPhi<T>::SubElt(const SphereThetaPhi<T>&) SizeMismatch")) ;
+    }
+  pixels_ *= a.pixels_;
+  return (*this);
+}
+