SparsityBasedContactConstraint Class Reference

#include <SparsityBasedContactConstraint.h>

Inheritance diagram for SparsityBasedContactConstraint:

Public Member Functions
	SparsityBasedContactConstraint (const InputParameters &parameters)
virtual	~SparsityBasedContactConstraint ()
virtual Real	computeQpSlaveValue ()
virtual Real	computeQpResidual (Moose::ConstraintType)
virtual Real	computeQpJacobian (Moose::ConstraintJacobianType)
virtual void	getConnectedDofIndices ()
	Gets the indices for all dofs conected to the constraint Get indices for all the slave Jacobian columns, not only those based on the mesh connectivity. More...

Detailed Description

Definition at line 21 of file SparsityBasedContactConstraint.h.

Constructor & Destructor Documentation

SparsityBasedContactConstraint()

SparsityBasedContactConstraint::SparsityBasedContactConstraint ( const InputParameters &  parameters )

inline

Definition at line 24 of file SparsityBasedContactConstraint.h.

    : NodeFaceConstraint(parameters){};

~SparsityBasedContactConstraint()

virtual SparsityBasedContactConstraint::~SparsityBasedContactConstraint ( )

inlinevirtual

Definition at line 26 of file SparsityBasedContactConstraint.h.

{}

Member Function Documentation

computeQpJacobian()

virtual Real SparsityBasedContactConstraint::computeQpJacobian ( Moose::ConstraintJacobianType  )

inlinevirtual

Reimplemented in GluedContactConstraint.

Definition at line 40 of file SparsityBasedContactConstraint.h.

  {
    mooseError("Unimplemented pure virtual method.  SparsityBasedContactConstraint should only be "
               "used as a base class");
    return 0;
  }

computeQpResidual()

virtual Real SparsityBasedContactConstraint::computeQpResidual ( Moose::ConstraintType  )

inlinevirtual

Reimplemented in GluedContactConstraint.

Definition at line 34 of file SparsityBasedContactConstraint.h.

  {
    mooseError("Unimplemented pure virtual method.  SparsityBasedContactConstraint should only be "
               "used as a base class");
    return 0;
  }

computeQpSlaveValue()

virtual Real SparsityBasedContactConstraint::computeQpSlaveValue ( )

inlinevirtual

Reimplemented in GluedContactConstraint.

Definition at line 28 of file SparsityBasedContactConstraint.h.

  {
    mooseError("Unimplemented pure virtual method.  SparsityBasedContactConstraint should only be "
               "used as a base class");
    return 0;
  }

getConnectedDofIndices()

void SparsityBasedContactConstraint::getConnectedDofIndices ( )

virtual

Gets the indices for all dofs conected to the constraint Get indices for all the slave Jacobian columns, not only those based on the mesh connectivity.

Definition at line 29 of file SparsityBasedContactConstraint.C.

{
#if defined(LIBMESH_HAVE_PETSC) && !PETSC_VERSION_LESS_THAN(3, 3, 0)
  _connected_dof_indices.clear();
 
  // An ugly hack: have to extract the row and look at it's sparsity structure, since otherwise I
  // won't get the dofs connected to this row by virtue of intervariable coupling
  // This is easier than sifting through all coupled variables, selecting those active on the
  // current subdomain, dealing with the scalar variables, etc.
  // Also, importantly, this will miss coupling to variables that might have introduced by prior
  // constraints similar to this one!
  PetscMatrix<Number> * petsc_jacobian = dynamic_cast<PetscMatrix<Number> *>(_jacobian);
  mooseAssert(petsc_jacobian, "Expected a PETSc matrix");
  Mat jac = petsc_jacobian->mat();
  PetscErrorCode ierr;
  PetscInt ncols;
  const PetscInt * cols;
  ierr = MatGetRow(jac, static_cast<PetscInt>(_var.nodalDofIndex()), &ncols, &cols, PETSC_NULL);
  CHKERRABORT(_communicator.get(), ierr);
  bool debug = false;
  if (debug)
  {
    libMesh::out << "_connected_dof_indices: adding " << ncols << " dofs from Jacobian row["
                 << _var.nodalDofIndex() << "] = [";
  }
  for (PetscInt i = 0; i < ncols; ++i)
  {
    if (debug)
    {
      libMesh::out << cols[i] << " ";
    }
    _connected_dof_indices.push_back(cols[i]);
  }
  if (debug)
  {
    libMesh::out << "]\n";
  }
  ierr = MatRestoreRow(jac, static_cast<PetscInt>(_var.nodalDofIndex()), &ncols, &cols, PETSC_NULL);
  CHKERRABORT(_communicator.get(), ierr);
#else
  NodeFaceConstraint::getConnectedDofIndices();
#endif
}

---

The documentation for this class was generated from the following files:

• contact/include/constraints/SparsityBasedContactConstraint.h
• contact/src/constraints/SparsityBasedContactConstraint.C