Home
Publications
Presentations
Support
Download
Examples
Developers
Class Docs
Namespaces
Classes
Files
Mailing Lists
Gallery
Wiki

fem_system.h

Go to the documentation of this file.

// The libMesh Finite Element Library.
// Copyright (C) 2002-2012 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



#ifndef LIBMESH_FEM_SYSTEM_H
#define LIBMESH_FEM_SYSTEM_H

// Local Includes
#include "libmesh/diff_system.h"
#include "libmesh/fem_physics.h"

// C++ includes
#include <cstddef>

namespace libMesh
{

// Forward Declarations
class DiffContext;
class FEMContext;


// ------------------------------------------------------------
// FEMSystem class definition

class FEMSystem : public DifferentiableSystem,
                  public FEMPhysics
{
public:

  FEMSystem (EquationSystems& es,
                 const std::string& name,
                 const unsigned int number);

  virtual ~FEMSystem ();

  typedef FEMSystem sys_type;

  typedef DifferentiableSystem Parent;

  virtual void clear ();

  virtual void assembly (bool get_residual, bool get_jacobian);

  virtual void solve ();

  void mesh_position_get();

  void mesh_position_set();

  virtual AutoPtr<DiffContext> build_context();

  /*
   * Prepares the result of a build_context() call for use.
   *
   * Most FEMSystem-based problems will need to reimplement this in order to
   * call FE::get_*() as their particular physics requires.
   */
  virtual void init_context(DiffContext &);

  virtual void postprocess ();

  virtual void assemble_qoi
    (const QoISet& indices = QoISet());

  virtual void assemble_qoi_derivative
    (const QoISet& indices = QoISet());

  bool fe_reinit_during_postprocess;

  Real numerical_jacobian_h;

  Real verify_analytic_jacobians;

  typedef bool (TimeSolver::*TimeSolverResPtr)(bool, DiffContext&);

  void numerical_jacobian (TimeSolverResPtr res,
                           FEMContext &context);

  void numerical_elem_jacobian (FEMContext &context);

  void numerical_side_jacobian (FEMContext &context);

protected:
  virtual void init_data ();
};


} // namespace libMesh


#endif // LIBMESH_FEM_SYSTEM_H

---

Site Created By: libMesh Developers
Last modified: February 05 2013 19:54:47 UTC

Hosted By: