ElementHDivSemiError.C

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://mooseframework.inl.gov
//*
//* All rights reserved, see COPYRIGHT for full restrictions
//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
//*
//* Licensed under LGPL 2.1, please see LICENSE for details
//* https://www.gnu.org/licenses/lgpl-2.1.html

#include "ElementHDivSemiError.h"
#include "Function.h"

registerMooseObject("MooseApp", ElementHDivSemiError);

InputParameters
ElementHDivSemiError::validParams()
{
  InputParameters params = ElementVectorL2Error::validParams();
  params.addClassDescription("Returns the H(div)-seminorm of the difference between a pair of "
                             "computed and analytical vector-valued solutions.");
  params.addRequiredParam<FunctionName>("function",
                                        "The vector analytical solution to compare against");
  params.addRequiredCoupledVar("variable", "The vector FE solution");
  params.suppressParameter<FunctionName>("function_x");
  params.suppressParameter<FunctionName>("function_y");
  params.suppressParameter<FunctionName>("function_z");
  params.suppressParameter<std::vector<VariableName>>("var_x");
  params.suppressParameter<std::vector<VariableName>>("var_y");
  params.suppressParameter<std::vector<VariableName>>("var_z");
  return params;
}

ElementHDivSemiError::ElementHDivSemiError(const InputParameters & parameters)
  : ElementVectorL2Error(parameters), _u_var(*getVectorVar("variable", 0)), _div_u(_u_var.divSln())
{
}

Real
ElementHDivSemiError::computeQpIntegral()
{
  Real diff = _div_u[_qp] - _func.div(_t, _q_point[_qp]);
  return diff * diff;
}