adjoint_initial.C

Go to the documentation of this file.

// The libMesh Finite Element Library.
// Copyright (C) 2002-2025 Benjamin S. Kirk, John W. Peterson, Roy H. Stogner

// This library is free software; you can redistribute it and/or
// modify it under the terms of the GNU Lesser General Public
// License as published by the Free Software Foundation; either
// version 2.1 of the License, or (at your option) any later version.

// This library 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
// Lesser General Public License for more details.

// You should have received a copy of the GNU Lesser General Public
// License along with this library; if not, write to the Free Software
// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA



#include "adjoint_initial.h"

using namespace libMesh;

void adjoint_read_initial_parameters()
{
}

void adjoint_finish_initialization()
{
}



// Initial conditions
Number adjoint_initial_value(const Point & p,
                             const Parameters &,
                             const std::string &,
                             const std::string &)
{
  Real x = p(0), y = p(1);

  return (sin(M_PI * x) * sin(M_PI * y));
}



Gradient adjoint_initial_grad(const Point & p,
                              const Parameters &,
                              const std::string &,
                              const std::string &)
{
  Real x = p(0), y = p(1);

  return Gradient(Number(M_PI*cos(M_PI * x) * sin(M_PI * y)),
                  Number(M_PI*sin(M_PI * x) * cos(M_PI * y)));
}