18 #include "libmesh/threads.h" 19 #include "libmesh/quadrature.h" 31 mooseWarning(
"NonlocalKernel is a computationally expensive experimental capability used only " 32 "for integral terms.");
54 for (
unsigned int i = 0; i < rows; i++)
58 var->sys().solution().add_vector(diag, var->dofIndices());
75 const auto phi_size = jvar.dofIndices().size();
78 for (
_j = 0;
_j < phi_size;
98 unsigned int n_total_dofs = var_alldofindices.size();
101 for (
_k = 0;
_k < n_total_dofs;
105 auto it = local_dofindices.find(var_alldofindices[
_k]);
106 if (it == local_dofindices.end())
134 const std::vector<dof_id_type> & jv_alldofindices = jv.
allDofIndices();
135 unsigned int n_total_dofs = jv_alldofindices.size();
138 for (
_k = 0;
_k < n_total_dofs;
142 auto it = local_dofindices.find(jv_alldofindices[
_k]);
143 if (it == local_dofindices.end())
MooseMesh & _mesh
Reference to this Kernel's mesh object.
void accumulateTaggedNonlocalMatrix()
Nonlocal Jacobian blocks will be appended by adding the current nonlocal kernel Jacobian.
static InputParameters validParams()
MooseVariable & _var
This is a regular kernel so we cast to a regular MooseVariable.
virtual void getUserObjectJacobian(unsigned int, dof_id_type)
Optimization option for getting jocobinas from userobject once per dof.
const MooseArray< Real > & _JxW
The current quadrature point weight value.
virtual void precalculateOffDiagJacobian(unsigned int)
unsigned int number() const
Get variable number coming from libMesh.
virtual Real computeQpOffDiagJacobian(unsigned int)
For coupling standard variables.
This class provides an interface for common operations on field variables of both FE and FV types wit...
virtual void computeJacobian() override
computeJacobian and computeQpOffDiagJacobian methods are almost same as Kernel except for few additio...
THREAD_ID _tid
The thread ID for this kernel.
bool _has_diag_save_in
The aux variables to save the diagonal Jacobian contributions to.
DenseMatrix< Number > _local_ke
Holds local Jacobian entries as they are accumulated by this Kernel.
const MooseVariableFieldBase & getVariable(unsigned int jvar_num) const
Retrieve the variable object from our system associated with jvar_num.
void mooseWarning(Args &&... args) const
Emits a warning prefixed with object name and type.
NonlocalKernel(const InputParameters ¶meters)
SystemBase & _sys
Reference to the EquationSystem object.
void errorIfDistributedMesh(std::string name) const
Generate a unified error message if the underlying libMesh mesh is a DistributedMesh.
void prepareMatrixTagNonlocal(Assembly &assembly, unsigned int ivar, unsigned int jvar)
Prepare data for computing nonlocal element jacobian according to the active tags.
unsigned int size() const
The number of elements that can currently be stored in the array.
const VariableTestValue & _test
the current test function
const std::vector< dof_id_type > & dofIndices() const final
Get local DoF indices.
static InputParameters validParams()
virtual void precalculateJacobian()
std::vector< MooseVariableFEBase * > _diag_save_in
const QBase *const & _qrule
active quadrature rule
virtual bool globalDoFEnabled(MooseVariableFEBase &, dof_id_type)
optimization option for executing nonlocal jacobian calculation only for nonzero elements ...
virtual Real computeQpJacobian()
Compute this Kernel's contribution to the Jacobian at the current quadrature point.
virtual const std::vector< dof_id_type > & dofIndices() const
Get local DoF indices.
unsigned int _i
current index for the test function
void accumulateTaggedLocalMatrix()
Local Jacobian blocks will be appended by adding the current local kernel Jacobian.
const MooseArray< Real > & _coord
The scaling factor to convert from cartesian to another coordinate system (e.g rz, spherical, etc.)
Assembly & _assembly
Reference to this Kernel's assembly object.
virtual void computeOffDiagJacobian(unsigned int jvar) override
Computes d-residual / d-jvar... storing the result in Ke.
virtual Real computeQpNonlocalOffDiagJacobian(unsigned int, dof_id_type)
unsigned int _j
current index for the shape function
virtual void computeNonlocalOffDiagJacobian(unsigned int jvar) override
Computes Jacobian entries corresponding to nonlocal dofs of the jvar.
const std::vector< dof_id_type > & allDofIndices() const
Get all global dofindices for the variable.
MooseVariableFieldBase & getVariable(THREAD_ID tid, const std::string &var_name) const
Gets a reference to a variable of with specified name.
virtual void computeNonlocalJacobian() override
computeNonlocalJacobian and computeNonlocalOffDiagJacobian methods are introduced for providing the j...
void prepareMatrixTag(Assembly &assembly, unsigned int ivar, unsigned int jvar)
Prepare data for computing element jacobian according to the active tags.
DenseMatrix< Number > _nonlocal_ke
Holds nonlocal Jacobian entries as they are accumulated by this Kernel.
virtual Real computeQpNonlocalJacobian(dof_id_type)
Compute this Kernel's contribution to the Jacobian corresponding to nolocal dof at the current quadra...
const VariablePhiValue & _phi
the current shape functions
unsigned int _qp
The current quadrature point index.