IP_SSAHardavTaoTikhonovProblem.hh

Go to the documentation of this file.

// Copyright (C) 2013, 2014, 2015, 2016, 2023  David Maxwell and Constantine Khroulev
//
// This file is part of PISM.
//
// PISM is free software; you can redistribute it and/or modify it under the
// terms of the GNU General Public License as published by the Free Software
// Foundation; either version 3 of the License, or (at your option) any later
// version.
//
// PISM is distributed in the hope that it will be useful, but WITHOUT ANY
// WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
// FOR A PARTICULAR PURPOSE.  See the GNU General Public License for more
// details.
//
// You should have received a copy of the GNU General Public License
// along with PISM; if not, write to the Free Software
// Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA  02110-1301  USA
 
 
#ifndef IP_SSAHARDAVTAOTIKHONOVPROBLEM_HH_HB8UWICX
#define IP_SSAHARDAVTAOTIKHONOVPROBLEM_HH_HB8UWICX
 
#include "pism/inverse/IPTaoTikhonovProblem.hh"
#include "pism/inverse/IP_SSAHardavForwardProblem.hh"
 
#include "pism/inverse/TaoUtil.hh"
#include "pism/inverse/functional/IPFunctional.hh"
 
namespace pism {
namespace inverse {
 
//! Defines an IPTaoTikhonovProblem for inversion of basal yeild stresses \f$\tau_c\f$ from %SSA velocities.
/*! The forward problem for the inversion is defined by an IP_SSAHardavForwardProblem.  The problem itself
  is solved with a TaoBasicSolver as described by the class-level documentation for IPTaoTikhonovProblem.  It
  is a reduced space method, inasmuch as we are performing unconstrained minimization on a Tikhonov functional
  that depends on a design variable (the value of \f$\tau_c\f$). It is compatible with any of the elementary
  TAO minimization algorithms, e.g. tao_cg, tao_lmvm.  If the minimization algorithm tao_blmvm is selected,
  the values of \f$\tau_c\f$ will be constrained by the config variables \a inverse.ssa.hardav_min
  and \a inverse.ssa.hardav_max.
*/
class IP_SSAHardavTaoTikhonovProblem: public IPTaoTikhonovProblem<IP_SSAHardavForwardProblem> {
public:
 
  IP_SSAHardavTaoTikhonovProblem(IP_SSAHardavForwardProblem &forward, 
                                  IP_SSAHardavTaoTikhonovProblem::DesignVec &d0, 
                                  IP_SSAHardavTaoTikhonovProblem::StateVec &u_obs, double eta, 
                                  IPFunctional<IP_SSAHardavTaoTikhonovProblem::DesignVec>&designFunctional, 
                                  IPFunctional<IP_SSAHardavTaoTikhonovProblem::StateVec>&stateFunctional) :
    IPTaoTikhonovProblem<IP_SSAHardavForwardProblem>(forward,d0,u_obs,eta,designFunctional,stateFunctional) {};
 
  virtual ~IP_SSAHardavTaoTikhonovProblem() {};
 
  virtual void connect(Tao tao);
 
  //! Callback to TAO to set bounds on \f$\tau_c\f$ for constrained minimization algorithms.
  virtual void getVariableBounds(Tao tao, Vec lo, Vec hi); 
 
};
 
} // end of namespace inverse
} // end of namespace pism
 
#endif /* end of include guard: IP_SSAHARDAVTIKHONOVPROBLEM_HH_HB8UWICX */