19 #include "libmesh/dense_vector.h" 20 #include "metaphysicl/raw_type.h" 30 #ifdef MOOSE_KOKKOS_ENABLED 60 #ifdef MOOSE_KOKKOS_ENABLED 82 #ifdef MOOSE_KOKKOS_ENABLED 101 #ifdef MOOSE_KOKKOS_ENABLED 118 void useVectorTag(
const TagName & tag_name, VectorTagsKey);
120 void useMatrixTag(
const TagName & tag_name, MatrixTagsKey);
263 template <
typename Res
iduals,
typename Indices>
265 const Residuals & residuals,
266 const Indices & dof_indices,
267 Real scaling_factor);
272 template <
typename T,
typename Indices>
274 const DenseVector<T> & residuals,
275 const Indices & dof_indices,
276 Real scaling_factor);
282 template <
typename Res
iduals,
typename Indices>
284 const Residuals & residuals,
285 const Indices & dof_indices,
286 Real scaling_factor);
291 template <
typename Res
iduals,
typename Indices>
293 const Residuals & residuals,
294 const Indices & dof_indices,
295 Real scaling_factor);
318 template <
typename Res
iduals,
typename Indices>
320 const Residuals & residuals,
321 const Indices & dof_indices,
322 Real scaling_factor);
329 template <
typename Res
iduals,
typename Indices>
331 const Residuals & residuals,
332 const Indices & dof_indices,
333 Real scaling_factor);
340 template <
typename Res
iduals,
typename Indices>
342 const Residuals & residuals,
343 const Indices & dof_indices,
344 Real scaling_factor);
353 Real scaling_factor);
360 const std::vector<dof_id_type> & row_indices,
361 const std::vector<dof_id_type> & column_indices,
362 Real scaling_factor);
367 template <
typename T>
378 template <
typename SetRes
idualFunctor>
402 const std::set<TagID> & vector_tags,
403 const std::set<TagID> & absolute_value_vector_tags);
454 #define usingTaggingInterfaceMembers \ 455 using TaggingInterface::_subproblem; \ 456 using TaggingInterface::accumulateTaggedLocalResidual; \ 457 using TaggingInterface::accumulateTaggedLocalMatrix; \ 458 using TaggingInterface::prepareVectorTag; \ 459 using TaggingInterface::prepareMatrixTag; \ 460 using TaggingInterface::prepareVectorTagNeighbor; \ 461 using TaggingInterface::_local_re; \ 462 using TaggingInterface::prepareVectorTagLower; \ 463 using TaggingInterface::prepareMatrixTagNeighbor; \ 464 using TaggingInterface::prepareMatrixTagLower; \ 465 using TaggingInterface::_local_ke 467 template <
typename Res
iduals,
typename Indices>
470 const Residuals & residuals,
471 const Indices & dof_indices,
472 const Real scaling_factor)
490 template <
typename T,
typename Indices>
493 const DenseVector<T> & residuals,
494 const Indices & dof_indices,
495 const Real scaling_factor)
497 addResiduals(assembly, residuals.get_values(), dof_indices, scaling_factor);
500 template <
typename Res
iduals,
typename Indices>
503 const Residuals & residuals,
504 const Indices & dof_indices,
505 const Real scaling_factor)
523 template <
typename Res
iduals,
typename Indices>
526 const Residuals & residuals,
527 const Indices & dof_indices,
530 addResiduals(assembly, residuals, dof_indices, scaling_factor);
531 addJacobian(assembly, residuals, dof_indices, scaling_factor);
534 template <
typename Res
iduals,
typename Indices>
537 const Residuals & residuals,
538 const Indices & dof_indices,
545 template <
typename Res
iduals,
typename Indices>
548 const Residuals & residuals,
549 const Indices & dof_indices,
556 template <
typename Res
iduals,
typename Indices>
559 const Residuals & residuals,
560 const Indices & dof_indices,
572 const Real scaling_factor)
581 const std::vector<dof_id_type> & row_indices,
582 const std::vector<dof_id_type> & column_indices,
583 const Real scaling_factor)
590 template <
typename T>
607 template <
typename SetRes
idualFunctor>
613 set_residual_functor(sys.
getVector(tag_id));
Nonlinear eigenvalue system to be solved.
MetaPhysicL::DualNumber< V, D, asd > abs(const MetaPhysicL::DualNumber< V, D, asd > &a)
void accumulateTaggedNonlocalMatrix()
Nonlocal Jacobian blocks will be appended by adding the current nonlocal kernel Jacobian.
std::set< TagID > _ref_abs_vector_tags
A set of either size 1 or 0.
SubProblem & _subproblem
SubProblem that contains tag info.
void addResidualsAndJacobian(Assembly &assembly, const Residuals &residuals, const Indices &dof_indices, Real scaling_factor)
Add the provided incoming residuals and derivatives for the Jacobian, corresponding to the provided d...
const InputParameters & _tag_params
Parameters from moose object.
void accumulateTaggedLocalResidual()
Local residual blocks will be appended by adding the current local kernel residual.
bool hasVector(const std::string &tag_name) const
Check if the named vector exists in the system.
bool hasVectorTags() const
Keeps track of stuff related to assembling.
const std::vector< dof_id_type > & dof_indices
void assignTaggedLocalMatrix()
Local Jacobian blocks will assigned as the current local kernel Jacobian.
void addResiduals(Assembly &assembly, const Residuals &residuals, const Indices &dof_indices, Real scaling_factor)
Add the provided incoming residuals corresponding to the provided dof indices.
void assignTaggedLocalResidual()
Local residual blocks will assigned as the current local kernel residual.
std::vector< DenseVector< Number > * > _absre_blocks
Residual blocks for absolute value residual tags.
The base class for Kokkos residual objects.
const Real scaling_factor
ResidualTagType
Enumerate whether a (residual) vector tag is to be of a non-reference or reference tag type...
std::vector< DenseVector< Number > * > _re_blocks
Residual blocks Vectors For each Tag.
const DenseVector< ADReal > & residuals
std::set< TagID > _ref_vector_tags
A set of either size 1 or 0.
Base class for a system (of equations)
DenseMatrix< Number > _local_ke
Holds local Jacobian entries as they are accumulated by this Kernel.
std::set< TagID > _abs_vector_tags
The absolute value residual tag ids.
void cacheJacobian(GlobalDataKey)
Takes the values that are currently in _sub_Kee and appends them to the cached values.
Nonlinear system to be solved.
void prepareMatrixTagNeighbor(Assembly &assembly, unsigned int ivar, unsigned int jvar, Moose::DGJacobianType type)
Prepare data for computing element jacobian according to the active tags for DG and interface kernels...
TaggingInterface(const MooseObject *moose_object)
void addJacobian(Assembly &assembly, const Residuals &residuals, const Indices &dof_indices, Real scaling_factor)
Add the provided residual derivatives into the Jacobian for the provided dof indices.
InputParameters validParams()
std::set< TagID > _non_ref_abs_vector_tags
A set to hold absolute value vector tags excluding the reference residual tag.
void prepareMatrixTagNonlocal(Assembly &assembly, unsigned int ivar, unsigned int jvar)
Prepare data for computing nonlocal element jacobian according to the active tags.
Real value(unsigned n, unsigned alpha, unsigned beta, Real x)
Every object that can be built by the factory should be derived from this class.
std::set< TagID > _matrix_tags
The matrices this Kernel will contribute to.
void cacheJacobianWithoutConstraints(const Residuals &residuals, const Indices &row_indices, Real scaling_factor, LocalDataKey, const std::set< TagID > &matrix_tags)
Process the derivatives() data of a vector of ADReals.
void addResidualsAndJacobianWithoutConstraints(Assembly &assembly, const Residuals &residuals, const Indices &dof_indices, Real scaling_factor)
Add the provided incoming residuals and derivatives for the Jacobian, corresponding to the provided d...
std::set< TagID > _non_ref_vector_tags
A set to hold vector tags excluding the reference residual tag.
std::set< TagID > _vector_tags
The vector tag ids this Kernel will contribute to.
void accumulateTaggedLocalMatrix()
Local Jacobian blocks will be appended by adding the current local kernel Jacobian.
Utility structure for packaging up all of the residual object's information needed to add into the sy...
A storage container for MooseObjects that inherit from SetupInterface.
void cacheResiduals(const Residuals &residuals, const Indices &row_indices, Real scaling_factor, LocalDataKey, const std::set< TagID > &vector_tags)
Process the supplied residual values.
void useMatrixTag(const TagName &tag_name, MatrixTagsKey)
void insertNodalValue(libMesh::NumericVector< libMesh::Number > &residual, const OutputData &v)
Write a nodal value to the passed-in solution vector.
void prepareVectorTagNeighbor(Assembly &assembly, unsigned int ivar)
Prepare data for computing element residual the according to active tags for DG and interface kernels...
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
This is the common base class for objects that give residual contributions.
Generic class for solving transient nonlinear problems.
This is the common base class for objects that give contributions to a linear system.
void prepareMatrixTagLower(Assembly &assembly, unsigned int ivar, unsigned int jvar, Moose::ConstraintJacobianType type)
Prepare data for computing the jacobian according to the active tags for mortar.
DenseVector< Number > _local_re
Holds local residual entries as they are accumulated by this Kernel.
const MooseObject & _moose_object
Moose objct this tag works on.
const std::set< TagID > & getMatrixTags(MatrixTagsKey) const
void addResidualsWithoutConstraints(Assembly &assembly, const Residuals &residuals, const Indices &dof_indices, Real scaling_factor)
Add the provided incoming residuals corresponding to the provided dof indices.
virtual ~TaggingInterface()
virtual void set(const numeric_index_type i, const Number value)=0
void cacheJacobianBlock(DenseMatrix< Number > &jac_block, const std::vector< dof_id_type > &idof_indices, const std::vector< dof_id_type > &jdof_indices, Real scaling_factor, LocalDataKey, TagID tag)
Cache a local Jacobian block with the provided rows (idof_indices) and columns (jdof_indices) for eve...
std::vector< DenseMatrix< Number > * > _ke_blocks
Kernel blocks Vectors For each Tag.
void prepareVectorTagLower(Assembly &assembly, unsigned int ivar)
Prepare data for computing the residual according to active tags for mortar constraints.
void prepareVectorTagInternal(Assembly &assembly, unsigned int ivar, const std::set< TagID > &vector_tags, const std::set< TagID > &absolute_value_vector_tags)
Prepare data for computing element residual according to the specified tags Residual blocks for diffe...
void prepareVectorTag(Assembly &assembly, unsigned int ivar)
Prepare data for computing element residual according to active tags.
void addJacobianWithoutConstraints(Assembly &assembly, const Residuals &residuals, const Indices &dof_indices, Real scaling_factor)
Add the provided residual derivatives into the Jacobian for the provided dof indices.
void prepareMatrixTag(Assembly &assembly, unsigned int ivar, unsigned int jvar)
Prepare data for computing element jacobian according to the active tags.
void addJacobianElement(Assembly &assembly, Real value, dof_id_type row_index, dof_id_type column_index, Real scaling_factor)
Add into a single Jacobian element.
DenseMatrix< Number > _nonlocal_ke
Holds nonlocal Jacobian entries as they are accumulated by this Kernel.
static InputParameters validParams()
void setResidual(SystemBase &sys, const T &residual, MooseVariableFE< T > &var)
Set residual using the variables' insertion API.
virtual NumericVector< Number > & getVector(const std::string &name)
Get a raw NumericVector by name.
auto index_range(const T &sizable)
std::vector< Real > _absolute_residuals
A container to hold absolute values of residuals passed into addResiduals.
Key structure for APIs adding/caching local element residuals/Jacobians.
void useVectorTag(const TagName &tag_name, VectorTagsKey)
void cacheResidualsWithoutConstraints(const Residuals &residuals, const Indices &row_indices, Real scaling_factor, LocalDataKey, const std::set< TagID > &vector_tags)
Process the supplied residual values.
const std::set< TagID > & getVectorTags(VectorTagsKey) const