// This may look like C code, but it is really -*- C++ -*- #ifndef LOCALMAP_SEEN #define LOCALMAP_SEEN #include "pixelmap.h" #include "sphericalmap.h" #include "ndatablock.h" #include "anydataobj.h" #include "ppersist.h" // ***************** Class LocalMap ***************************** namespace SOPHYA { template class FIO_LocalMap; template class FITS_LocalMap; /* Class LocalMap */ template class LocalMap : public PixelMap { // friend declaration for classes which handle persistence and FITS IO friend class FIO_LocalMap; friend class FITS_LocalMap; public: LocalMap(); LocalMap(int_4 nx, int_4 ny); LocalMap(int_4 nx, int_4 ny, double angleX,double angleY, double theta0,double phi0,int_4 x0,int_4 y0,double angle=0.); LocalMap(int_4 nx, int_4 ny, double angleX,double angleY, double theta0,double phi0, double angle=0.); LocalMap(const LocalMap& lm, bool share); LocalMap(const LocalMap& lm); virtual ~LocalMap(); inline virtual bool IsTemp(void) const { return pixels_.IsTemp();} /*! Setting blockdata to temporary (see ndatablock documentation) */ inline virtual void SetTemp(bool temp=false) const {pixels_.SetTemp(temp);}; // ---------- Overloading of () to access pixel number k ---- inline T& operator()(int_4 k) {return(PixVal(k));} inline T const& operator()(int_4 k) const {return(PixVal(k));} inline T& operator()(int_4 ix, int_4 iy) {return PixVal(iy*nSzX_+ix);}; inline T const& operator()(int_4 ix, int_4 iy) const {return PixVal(iy*nSzX_+ix);}; // Acces a un pixel inline T& operator()(double theta,double phi) {return(PixValSph(theta,phi));}; inline T const& operator()(double theta,double phi) const {return(PixValSph(theta,phi));}; // ---------- Definition of PixelMap abstract methods ------- virtual int_4 NbPixels() const; // D.Y. int change en int_4 rationalisation Mac virtual T& PixVal(int_4 k); virtual T const& PixVal(int_4 k) const; //! if fg = true, will generate exceptions , when theta-phi out of range inline void SetThrowExceptionWhenOutofMapFlag(bool fg = false) { exc_outofmap_ = fg; } //! true --> exceptions wil be active when theta-phi out of range inline bool GetThrowExceptionWhenOutofMapFlag() { return exc_outofmap_; } virtual bool ContainsSph(double theta, double phi) const; //! return pixel number correspoinding to (theta,phi) virtual int_4 PixIndexSph(double theta,double phi) const; //! return (theta,phi) correspoinding to pixel number k virtual void PixThetaPhi(int_4 k,double& theta,double& phi) const; //! set all pixel values equal to v virtual T SetPixels(T v); // return the solid angle covered by pixel k (steradians) virtual double PixSolAngle(int_4 k) const; // ---------- Specific methods ------------------------------ void ReSize(int_4 nx, int_4 ny); void SetOrigin(double theta=90.,double phi=0.,double angle=0.); void SetOrigin(double theta,double phi,int_4 x0,int_4 y0,double angle=0.); void SetSize(double angleX,double angleY); /*! \fn TypeOfMap() return the string "LOCAL" */ inline virtual string TypeOfMap() const {return string("LOCAL");}; /*! Check to see if the local mapping is done */ inline bool LocalMap_isDone() const {return( originSet_ && angleSized_ && pixelSized_);}; void PixThetaPhi(int_4 ip,int_4 it, double& theta,double& phi) const; void ProjectionToSphere(SphericalMap&) const; /* There should be a more complex algorithm somewhere to combine *several* local maps to a full sphere. -> static method, or separate class */ /*! provides a integer characterizing the pixelization refinement (here : number of pixels) */ inline virtual int_4 SizeIndex() const {return(nPix_);} inline int_4 Size_x() const {return nSzX_;} inline int_4 Size_y() const {return nSzY_;} /* Je rajoute ces 2 fonctions inlines pour compatibilite d'interface avec TArray - Reza 30/8/2000 */ inline int_4 SizeX() const {return nSzX_;} inline int_4 SizeY() const {return nSzY_;} inline void Origin(double& theta,double& phi,int_4& x0,int_4& y0,double& angle) const {theta= thetaDegresC_; phi= phiDegresC_; x0= x0_; y0= y0_;angle= angleDegres_;} /*! total aperture in theta and phi, in degrees ( from SetSize() ) */ inline void Aperture(double& anglex,double& angley) const {anglex= angleDegresX_; angley= angleDegresY_;} /* Acces to the DataBlock */ inline NDataBlock& DataBlock() {return pixels_.DataBlock();} inline const NDataBlock& DataBlock() const {return pixels_.DataBlock();} /* Acces to the matrix */ //! access to matrix /*! \warning : a pixel is defined by the pair of a phi-like index (x-axis) and a theta-like index (y-axis). The phi-like index denotes the column number in the matrix ; the theta-like index denotes the row number. */ inline TMatrix& Matrix() {return pixels_;} inline const TMatrix& Matrix() const {return pixels_;} /* impression */ void print(ostream& os) const; // Operations diverses = , +=, ... virtual LocalMap& Set(const LocalMap& a); inline LocalMap& operator = (const LocalMap& a) {return Set(a);} // A += -= *= /= x (ajoute, soustrait, ... x a tous les elements) //! Fill LocalMap with all elements equal to \b x virtual LocalMap& SetT(T a); inline LocalMap& operator = (T a) {return SetT(a);} // Add \b x to all elements virtual LocalMap& Add(T a); inline LocalMap& operator += (T x) { return Add(x); } // Substract \b x to all elements virtual LocalMap& Sub(T a,bool fginv=false); inline LocalMap& operator -= (T x) { return Sub(x); } // Multiply all elements by \b x virtual LocalMap& Mul(T a); inline LocalMap& operator *= (T x) { return Mul(x); } // Divide all elements by \b x virtual LocalMap& Div(T a); inline LocalMap& operator /= (T x) { return Div(x); } // A += -= (ajoute, soustrait element par element les deux spheres ) //! Operator LocalMap += LocalMap virtual LocalMap& AddElt(const LocalMap& a); inline LocalMap& operator += (const LocalMap& a) { return AddElt(a); } virtual LocalMap& SubElt(const LocalMap& a); //! Operator LocalMap -= LocalMap inline LocalMap& operator -= (const LocalMap& a) { return SubElt(a); } // Multiplication, division element par element les deux LocalMap virtual LocalMap& MulElt(const LocalMap& a); inline LocalMap& operator *= (const LocalMap& a) { return MulElt(a); } virtual LocalMap& DivElt(const LocalMap& a); inline LocalMap& operator /= (const LocalMap& a) { return DivElt(a); } void CloneOrShare(const LocalMap& a); void Share(const LocalMap& a); LocalMap& CopyElt(const LocalMap& a); // ---------- Méthodes internes ----------------------------- private : void InitNul(); void SetCoorC(double theta0, double phi0); TMatrix CalculMatricePassage(); void Getij(int_4 k,int_4& i,int_4& j) const; void PixToSphereC(int_4 ip, int_4 it, double& XP, double& YP, double& ZP) const; void recopierVariablesSimples(const LocalMap& lm); void initializationError() const; // ---------- Variables internes ---------------------------- // variables suffisantes pour reconstruire l'objet bool pixelSized_; bool angleSized_; bool originSet_; int_4 nSzX_; int_4 nSzY_; double angleDegresX_; double angleDegresY_; double thetaDegresC_; double phiDegresC_; int_4 x0_; int_4 y0_; double angleDegres_; // NDataBlock pixels_; TMatrix pixels_; // variables derivees (redondantes, precalculees pour ameliorer // les performances) int_4 nPix_; double thetaC_; double phiC_; double csthC_; double snthC_; double csphC_; double snphC_; double XC_; double YC_; double ZC_; double angle_; double cosAngle_; double sinAngle_; double deltaPhi_; double deltaTheta_; bool exc_outofmap_; }; //////////////////////////////////////////////////////////////// // Surcharge d'operateurs A (+,-,*,/) (T) x /*! \ingroup SkyMap \fn operator+(const LocalMap&,T) \brief Operator LocalMap = LocalMap + constant */ template inline LocalMap operator + (const LocalMap& a, T b) {LocalMap result; result.CloneOrShare(a); result.SetTemp(true); result.Add(b); return result;} /*! \ingroup SkyMap \fn operator+(T,const LocalMap&) \brief Operator LocalMap = constant + LocalMap */ template inline LocalMap operator + (T b,const LocalMap& a) {LocalMap result; result.CloneOrShare(a); result.SetTemp(true); result.Add(b); return result;} /*! \ingroup SkyMap\fn operator-(const LocalMap&,T) \brief Operator LocalMap = LocalMap - constant */ template inline LocalMap operator - (const LocalMap& a, T b) {LocalMap result; result.CloneOrShare(a); result.SetTemp(true); result.Sub(b); return result;} /*! \ingroup \fn operator-(T,const LocalMap&) \brief Operator LocalMap = constant - LocalMap */ template inline LocalMap operator - (T b,const LocalMap& a) {LocalMap result; result.CloneOrShare(a); result.SetTemp(true); result.Sub(b,true); return result;} /*! \ingroup SkyMap \fn operator*(const LocalMap&,T) \brief Operator LocalMap = LocalMap * constant */ template inline LocalMap operator * (const LocalMap& a, T b) {LocalMap result; result.CloneOrShare(a); result.SetTemp(true); result.Mul(b); return result;} /*! \ingroup SkyMap \fn operator*(T,const LocalMap&) \brief Operator LocalMap = constant * LocalMap */ template inline LocalMap operator * (T b,const LocalMap& a) {LocalMap result; result.CloneOrShare(a); result.SetTemp(true); result.Mul(b); return result;} /*! \ingroup SkyMap \fn operator/(const LocalMap&,T) \brief Operator LocalMap = LocalMap / constant */ template inline LocalMap operator / (const LocalMap& a, T b) {LocalMap result; result.CloneOrShare(a); result.SetTemp(true); result.Div(b); return result;} /*! \ingroup SkyMap \fn operator/(T,const LocalMap&) \brief Operator LocalMap = constant / LocalMap */ template inline LocalMap operator / (T b, const LocalMap& a) {LocalMap result; result.CloneOrShare(a); result.SetTemp(true); result.Div(b, true); return result;} //////////////////////////////////////////////////////////////// // Surcharge d'operateurs C = A (+,-) B /*! \ingroup SkyMap \fn operator+(const LocalMap&,const LocalMap&) \brief Operator LocalMap = LocalMap + LocalMap */ template inline LocalMap operator + (const LocalMap& a,const LocalMap& b) { LocalMap result; result.SetTemp(true); if (b.IsTemp()) { result.Share(b); result.AddElt(a); } else { result.CloneOrShare(a); result.AddElt(b); } return result; } /*! \ingroup SkyMap \fn operator-(const LocalMap&,const LocalMap&) \brief Operator LocalMap = LocalMap - LocalMap */ template inline LocalMap operator - (const LocalMap& a,const LocalMap& b) { LocalMap result; result.SetTemp(true); if (b.IsTemp()) { result.Share(b); result.SubElt(a); } else { result.CloneOrShare(a); result.SubElt(b); } return result; } //////////////////////////////////////////////////////////////// // Surcharge d'operateurs C = A (*,/) B /*! \ingroup SkyMap \fn operator*(const LocalMap&,const LocalMap&) \brief Operator LocalMap = LocalMap * LocalMap (pixel by pixel multiply)*/ template inline LocalMap operator * (const LocalMap& a,const LocalMap& b) { LocalMap result; result.SetTemp(true); if (b.IsTemp()) { result.Share(b); result.MulElt(a); } else { result.CloneOrShare(a); result.MulElt(b); } return result; } /*! \ingroup SkyMap \fn operator/(const LocalMap&,const LocalMap&) \brief Operator LocalMap = LocalMap / LocalMap (pixel by pixel divide)*/ template inline LocalMap operator / (const LocalMap& a,const LocalMap& b) { LocalMap result; result.SetTemp(true); if (b.IsTemp()) { result.Share(b); result.DivElt(a); } else { result.CloneOrShare(a); result.DivElt(b); } return result; } } // Fin du namespace #endif