use pointers in eval_gradient()

knodlseder
1 parent 3af07206
Showing 13 changed files with 15 additions and 15 deletions
include/GModel.hpp
include/GModelSpatial.hpp
include/GModelSpatialPtsrc.hpp
include/GModelSpectral.hpp
include/GModelSpectralPlaw.hpp
include/GModelTemporal.hpp
include/GModelTemporalConst.hpp
include/GModels.hpp
src/model/GModel.cpp
src/model/GModelSpatialPtsrc.cpp
src/model/GModelSpectralPlaw.cpp
src/model/GModelTemporalConst.cpp
src/model/GModels.cpp
@@ -56,7 +56,7 @@ public:
     void        name(const std::string& name) { m_name=name; return; }
     int         npars(void) const { return m_npars; }
     GModelPar*  par(int index) const;
-    void        eval_gradients(GSkyDir& dir, const double& energy, const double& time);
+    void        eval_gradients(GSkyDir* dir, double* energy, double* time);
  
 protected:
     // Protected methods
@@ -45,7 +45,7 @@ public:
     // Virtual methods
     virtual int        npars(void) const = 0;
     virtual GModelPar* par(int index) const = 0;
-    virtual void       eval_gradients(GSkyDir& dir) = 0;
+    virtual void       eval_gradients(GSkyDir* dir) = 0;
  
 protected:
     // Protected methods
@@ -49,7 +49,7 @@ public:
     // Methods
     int        npars(void) const { return m_npars; }
     GModelPar* par(int index) const;
-    void       eval_gradients(GSkyDir& dir);
+    void       eval_gradients(GSkyDir* dir);
     double     ra(void) const { return m_ra.real_value(); }
     double     dec(void) const { return m_dec.real_value(); }
  
@@ -45,7 +45,7 @@ public:
     // Virtual methods
     virtual int        npars(void) const = 0;
     virtual GModelPar* par(int index) const = 0;
-    virtual void       eval_gradients(const double& energy) = 0;
+    virtual void       eval_gradients(double* energy) = 0;
  
 protected:
     // Protected methods
@@ -48,7 +48,7 @@ public:
     // Methods
     int        npars(void) const { return m_npars; }
     GModelPar* par(int index) const;
-    void       eval_gradients(const double& energy);
+    void       eval_gradients(double* energy);
     void       autoscale(void);
     double     norm(void) const { return m_norm.real_value(); }
     double     index(void) const { return m_index.real_value(); }
@@ -45,7 +45,7 @@ public:
     // Virtual methods
     virtual int        npars(void) const = 0;
     virtual GModelPar* par(int index) const = 0;
-    virtual void       eval_gradients(const double& time) = 0;
+    virtual void       eval_gradients(double* time) = 0;
  
 protected:
     // Protected methods
@@ -49,7 +49,7 @@ public:
     // Methods
     int        npars(void) const { return m_npars; }
     GModelPar* par(int index) const;
-    void       eval_gradients(const double& time);
+    void       eval_gradients(double* time);
  
 protected:
     // Protected methods
@@ -48,7 +48,7 @@ public:
  
     // Methods
     void add(const GModel& model);
-    void eval_gradients(GSkyDir& dir, const double& energy, const double& time);
+    void eval_gradients(GSkyDir* dir, double* energy, double* time);
  
 protected:
     // Protected methods
@@ -167,7 +167,7 @@ GModelPar* GModel::par(int index) const
 /***********************************************************************//**
  * @brief Evaluate function gradients
  ***************************************************************************/
-void GModel::eval_gradients(GSkyDir& dir, const double& energy, const double& time)
+void GModel::eval_gradients(GSkyDir* dir, double* energy, double* time)
 {
     // Evaluate gradients in all components
     if (m_spatial  != NULL) m_spatial->eval_gradients(dir);
@@ -161,7 +161,7 @@ GModelPar* GModelSpatialPtsrc::par(int index) const
 /***********************************************************************//**
  * @brief Evaluate function gradients
  ***************************************************************************/
-void GModelSpatialPtsrc::eval_gradients(GSkyDir& dir)
+void GModelSpatialPtsrc::eval_gradients(GSkyDir* dir)
 {
     // Return
     return;
@@ -166,7 +166,7 @@ GModelPar* GModelSpectralPlaw::par(int index) const
  *
  * @param[in] energy Energy at which the gradients are to be computed.
  ***************************************************************************/
-void GModelSpectralPlaw::eval_gradients(const double& energy)
+void GModelSpectralPlaw::eval_gradients(double* energy)
 {
     // Set gradient to 0 is all parameters are fixed ...
     if ( !m_norm.isfree() && !m_index.isfree() && !m_pivot.isfree() ) {
@@ -179,7 +179,7 @@ void GModelSpectralPlaw::eval_gradients(const double&amp; energy)
     else {
  
         // Precompute some values
-        double e = energy/m_pivot.value();
+        double e = *energy/m_pivot.value();
         double p = pow(e, m_index.value());
  
         // Compute normalization gradient
@@ -190,7 +190,7 @@ void GModelSpectralPlaw::eval_gradients(const double&amp; energy)
  
         // Compute index gradient
         if (m_index.isfree()) {
-            double g = m_norm.value() * p * log(energy);
+            double g = m_norm.value() * p * log(*energy);
             m_index.gradient(g);
         }
         else
@@ -142,7 +142,7 @@ GModelPar* GModelTemporalConst::par(int index) const
 /***********************************************************************//**
  * @brief Evaluate function gradients
  ***************************************************************************/
-void GModelTemporalConst::eval_gradients(const double& time)
+void GModelTemporalConst::eval_gradients(double* time)
 {
     // Return
     return;
@@ -165,7 +165,7 @@ void GModels::add(const GModel&amp; model)
 /***********************************************************************//**
  * @brief Evaluate function gradients
  ***************************************************************************/
-void GModels::eval_gradients(GSkyDir& dir, const double& energy, const double& time)
+void GModels::eval_gradients(GSkyDir* dir, double* energy, double* time)
 {
     // Evaluate gradients for all models
     for (int i = 0; i < m_elements; ++i)