ElementH1SemiError.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 "ElementH1SemiError.h"
#include "Function.h"
 
registerMooseObject("MooseApp", ElementH1SemiError);
 
InputParameters
ElementH1SemiError::validParams()
{
  InputParameters params = ElementIntegralVariablePostprocessor::validParams();
  params.addRequiredParam<FunctionName>("function", "The analytic solution to compare against");
  params.addClassDescription("Returns the gradient difference norm part of the H1 error");
  return params;
}
 
ElementH1SemiError::ElementH1SemiError(const InputParameters & parameters)
  : ElementIntegralVariablePostprocessor(parameters), _func(getFunction("function"))
{
}
 
Real
ElementH1SemiError::getValue() const
{
  return std::sqrt(ElementIntegralVariablePostprocessor::getValue());
}
 
Real
ElementH1SemiError::computeQpIntegral()
{
  RealGradient diff = _grad_u[_qp] - _func.gradient(_t, _q_point[_qp]);
  return diff * diff;
}