Uploaded image for project: 'Data Management'
  1. Data Management
  2. DM-14510

Implement line search

    Details

    • Templates:
    • Story Points:
      8
    • Sprint:
      AP S18-6, AP F18-1
    • Team:
      Alert Production

      Description

      Pierre Astier suggests that we should implement a line search in jointcal, to deal with the non-linearities in the photometric model (and because we may eventually have those in astrometry, as the astrometry model gets more complex). He suggested using Brent's method, which is available in GSL. It might be easier to just take the code directly from GSL, instead of writing an interface layer, given the way the nature of the jointcal models.

      At the top of the outlier rejection loop in FitterBase.minimize(), this would look something like:

      Eigen::VectorXd delta = chol.solve(grad);
      scale = lineSearch(delta);
      offsetParams(delta*scale));
      

        Attachments

          Issue Links

            Activity

              People

              • Assignee:
                Parejkoj John Parejko
                Reporter:
                Parejkoj John Parejko
                Reviewers:
                Pim Schellart [X] (Inactive)
                Watchers:
                Dominique Boutigny, Jim Bosch, John Parejko, John Swinbank, Pierre Astier, Pim Schellart [X] (Inactive), Tim Jenness
              • Votes:
                0 Vote for this issue
                Watchers:
                7 Start watching this issue

                Dates

                • Created:
                  Updated:
                  Resolved:

                  Summary Panel