19 #include "libmesh/dense_vector.h" 20 #include "metaphysicl/raw_type.h" 97 void useVectorTag(
const TagName & tag_name, VectorTagsKey);
99 void useMatrixTag(
const TagName & tag_name, MatrixTagsKey);
242 template <
typename Res
iduals,
typename Indices>
244 const Residuals & residuals,
245 const Indices & dof_indices,
246 Real scaling_factor);
251 template <
typename T,
typename Indices>
253 const DenseVector<T> & residuals,
254 const Indices & dof_indices,
255 Real scaling_factor);
261 template <
typename Res
iduals,
typename Indices>
263 const Residuals & residuals,
264 const Indices & dof_indices,
265 Real scaling_factor);
270 template <
typename Res
iduals,
typename Indices>
272 const Residuals & residuals,
273 const Indices & dof_indices,
274 Real scaling_factor);
297 template <
typename Res
iduals,
typename Indices>
299 const Residuals & residuals,
300 const Indices & dof_indices,
301 Real scaling_factor);
308 template <
typename Res
iduals,
typename Indices>
310 const Residuals & residuals,
311 const Indices & dof_indices,
312 Real scaling_factor);
319 template <
typename Res
iduals,
typename Indices>
321 const Residuals & residuals,
322 const Indices & dof_indices,
323 Real scaling_factor);
332 Real scaling_factor);
339 const std::vector<dof_id_type> & row_indices,
340 const std::vector<dof_id_type> & column_indices,
341 Real scaling_factor);
346 template <
typename T>
357 template <
typename SetRes
idualFunctor>
381 const std::set<TagID> & vector_tags,
382 const std::set<TagID> & absolute_value_vector_tags);
433 #define usingTaggingInterfaceMembers \ 434 using TaggingInterface::_subproblem; \ 435 using TaggingInterface::accumulateTaggedLocalResidual; \ 436 using TaggingInterface::accumulateTaggedLocalMatrix; \ 437 using TaggingInterface::prepareVectorTag; \ 438 using TaggingInterface::prepareMatrixTag; \ 439 using TaggingInterface::prepareVectorTagNeighbor; \ 440 using TaggingInterface::_local_re; \ 441 using TaggingInterface::prepareVectorTagLower; \ 442 using TaggingInterface::prepareMatrixTagNeighbor; \ 443 using TaggingInterface::prepareMatrixTagLower; \ 444 using TaggingInterface::_local_ke 446 template <
typename Res
iduals,
typename Indices>
449 const Residuals & residuals,
450 const Indices & dof_indices,
451 const Real scaling_factor)
469 template <
typename T,
typename Indices>
472 const DenseVector<T> & residuals,
473 const Indices & dof_indices,
474 const Real scaling_factor)
476 addResiduals(assembly, residuals.get_values(), dof_indices, scaling_factor);
479 template <
typename Res
iduals,
typename Indices>
482 const Residuals & residuals,
483 const Indices & dof_indices,
484 const Real scaling_factor)
502 template <
typename Res
iduals,
typename Indices>
505 const Residuals & residuals,
506 const Indices & dof_indices,
509 addResiduals(assembly, residuals, dof_indices, scaling_factor);
510 addJacobian(assembly, residuals, dof_indices, scaling_factor);
513 template <
typename Res
iduals,
typename Indices>
516 const Residuals & residuals,
517 const Indices & dof_indices,
524 template <
typename Res
iduals,
typename Indices>
527 const Residuals & residuals,
528 const Indices & dof_indices,
535 template <
typename Res
iduals,
typename Indices>
538 const Residuals & residuals,
539 const Indices & dof_indices,
551 const Real scaling_factor)
560 const std::vector<dof_id_type> & row_indices,
561 const std::vector<dof_id_type> & column_indices,
562 const Real scaling_factor)
569 template <
typename T>
586 template <
typename SetRes
idualFunctor>
592 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.
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
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