gradient_problem_evaluator.h 4.0 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103
  1. // Ceres Solver - A fast non-linear least squares minimizer
  2. // Copyright 2023 Google Inc. All rights reserved.
  3. // http://ceres-solver.org/
  4. //
  5. // Redistribution and use in source and binary forms, with or without
  6. // modification, are permitted provided that the following conditions are met:
  7. //
  8. // * Redistributions of source code must retain the above copyright notice,
  9. // this list of conditions and the following disclaimer.
  10. // * Redistributions in binary form must reproduce the above copyright notice,
  11. // this list of conditions and the following disclaimer in the documentation
  12. // and/or other materials provided with the distribution.
  13. // * Neither the name of Google Inc. nor the names of its contributors may be
  14. // used to endorse or promote products derived from this software without
  15. // specific prior written permission.
  16. //
  17. // THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
  18. // AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
  19. // IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
  20. // ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
  21. // LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
  22. // CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
  23. // SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
  24. // INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
  25. // CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
  26. // ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
  27. // POSSIBILITY OF SUCH DAMAGE.
  28. //
  29. // Author: sameeragarwal@google.com (Sameer Agarwal)
  30. #ifndef CERES_INTERNAL_GRADIENT_PROBLEM_EVALUATOR_H_
  31. #define CERES_INTERNAL_GRADIENT_PROBLEM_EVALUATOR_H_
  32. #include <map>
  33. #include <memory>
  34. #include <string>
  35. #include "ceres/evaluator.h"
  36. #include "ceres/execution_summary.h"
  37. #include "ceres/gradient_problem.h"
  38. #include "ceres/internal/disable_warnings.h"
  39. #include "ceres/internal/export.h"
  40. #include "ceres/sparse_matrix.h"
  41. #include "ceres/wall_time.h"
  42. namespace ceres::internal {
  43. class CERES_NO_EXPORT GradientProblemEvaluator final : public Evaluator {
  44. public:
  45. explicit GradientProblemEvaluator(const GradientProblem& problem)
  46. : problem_(problem) {}
  47. std::unique_ptr<SparseMatrix> CreateJacobian() const final { return nullptr; }
  48. bool Evaluate(const EvaluateOptions& /*evaluate_options*/,
  49. const double* state,
  50. double* cost,
  51. double* /*residuals*/,
  52. double* gradient,
  53. SparseMatrix* jacobian) final {
  54. CHECK(jacobian == nullptr);
  55. ScopedExecutionTimer total_timer("Evaluator::Total", &execution_summary_);
  56. // The reason we use Residual and Jacobian here even when we are
  57. // only computing the cost and gradient has to do with the fact
  58. // that the line search minimizer code is used by both the
  59. // GradientProblemSolver and the main CeresSolver coder where the
  60. // Evaluator evaluates the Jacobian, and these magic strings need
  61. // to be consistent across the code base for the time accounting
  62. // to work.
  63. ScopedExecutionTimer call_type_timer(
  64. gradient == nullptr ? "Evaluator::Residual" : "Evaluator::Jacobian",
  65. &execution_summary_);
  66. return problem_.Evaluate(state, cost, gradient);
  67. }
  68. bool Plus(const double* state,
  69. const double* delta,
  70. double* state_plus_delta) const final {
  71. return problem_.Plus(state, delta, state_plus_delta);
  72. }
  73. int NumParameters() const final { return problem_.NumParameters(); }
  74. int NumEffectiveParameters() const final {
  75. return problem_.NumTangentParameters();
  76. }
  77. int NumResiduals() const final { return 1; }
  78. std::map<std::string, internal::CallStatistics> Statistics() const final {
  79. return execution_summary_.statistics();
  80. }
  81. private:
  82. const GradientProblem& problem_;
  83. ::ceres::internal::ExecutionSummary execution_summary_;
  84. };
  85. } // namespace ceres::internal
  86. #include "ceres/internal/reenable_warnings.h"
  87. #endif // CERES_INTERNAL_GRADIENT_PROBLEM_EVALUATOR_H_