Updated documentation

Author: robbietuk

src/include/stir/recon_buildblock/PoissonLogLikelihoodWithLinearModelForMean.h

@@ -121,6 +121,10 @@ public  GeneralisedObjectiveFunction<TargetT>
   //PoissonLogLikelihoodWithLinearModelForMean(); 
 
   //! Computes the gradient of the data fit term
+  /*!
+      This function is implemented in terms of \c actual_compute_sub_gradient_without_penalty()
+      by setting do_subtraction = true
+   */
   virtual void 
     compute_sub_gradient_without_penalty(TargetT& gradient, 
                                          const TargetT &current_estimate, 
@@ -130,12 +134,8 @@ public  GeneralisedObjectiveFunction<TargetT>
   /*! 
     This function is used for instance by OSMAPOSL.
 
-    This computes
-    \f[ {\partial L \over \partial \lambda_v} + P_v =
-      \sum_b P_{bv} {y_b \over Y_b}
-      \f]
-    (see the class general documentation).
-    The sum will however be restricted to a subset.
+    This function is implemented in terms of \c actual_compute_sub_gradient_without_penalty()
+    by setting do_subtraction = false
    */
   virtual void 
     compute_sub_gradient_without_penalty_plus_sensitivity(TargetT& gradient, 
@@ -253,7 +253,20 @@ public  GeneralisedObjectiveFunction<TargetT>
   */
   void compute_sensitivities();
 
+  //! computes the objective function subset gradient without the penalty
+  /*!
+    If do_subtraction is false, this computes
+    \f[ {\partial L \over \partial \lambda_v} + P_v =
+      \sum_b P_{bv} {y_b \over Y_b}
+      \f]
+    (see the class general documentation).
+    The sum will however be restricted to a subset.
 
+    However, if do_subtraction is true, this function will instead compute
+    \f[ {\partial L \over \partial \lambda_v} =
+      \sum_b P_{bv} ({y_b \over Y_b} - 1)
+    \f]
+  */
   virtual void
   actual_compute_sub_gradient_without_penalty(TargetT& gradient,
                                               const TargetT &current_estimate,