ContactSplit Class Reference

Split-based preconditioner for contact problems. More...

#include <ContactSplit.h>

Inheritance diagram for ContactSplit:

Public Member Functions
	ContactSplit (const InputParameters &params)
virtual void	setup (const std::string &prefix="-") override

Protected Attributes
std::vector< std::string >	_contact_master
std::vector< std::string >	_contact_slave
std::vector< int >	_contact_displaced
std::vector< std::string >	_uncontact_master
std::vector< std::string >	_uncontact_slave
std::vector< int >	_uncontact_displaced
bool	_include_all_contact_nodes

Detailed Description

Split-based preconditioner for contact problems.

Definition at line 24 of file ContactSplit.h.

Constructor & Destructor Documentation

ContactSplit()

ContactSplit::ContactSplit

(

const InputParameters &

params

)

Definition at line 43 of file ContactSplit.C.

  : Split(params),
    _contact_master(getParam<std::vector<std::string>>("contact_master")),
    _contact_slave(getParam<std::vector<std::string>>("contact_slave")),
    _contact_displaced(getParam<std::vector<int>>("contact_displaced")),
    _uncontact_master(getParam<std::vector<std::string>>("uncontact_master")),
    _uncontact_slave(getParam<std::vector<std::string>>("uncontact_slave")),
    _uncontact_displaced(getParam<std::vector<int>>("uncontact_displaced")),
    _include_all_contact_nodes(getParam<bool>("include_all_contact_nodes"))
{
  if (_contact_master.size() != _contact_slave.size())
  {
    std::ostringstream err;
    err << "Master and slave contact lists must have matching sizes: " << _contact_master.size()
        << " != " << _contact_slave.size();
    mooseError(err.str());
  }
  if (_contact_displaced.size() && _contact_master.size() != _contact_displaced.size())
  {
    std::ostringstream err;
    err << "Master and displaced contact lists must have matching sizes: " << _contact_master.size()
        << " != " << _contact_displaced.size();
    mooseError(err.str());
  }
  if (!_contact_displaced.size())
    _contact_displaced.resize(_contact_master.size());

  if (_uncontact_master.size() != _uncontact_slave.size())
  {
    std::ostringstream err;
    err << "Master and slave uncontact lists must have matching sizes: " << _uncontact_master.size()
        << " != " << _uncontact_slave.size();
    mooseError(err.str());
  }
  if (_uncontact_displaced.size() && _uncontact_master.size() != _uncontact_displaced.size())
  {
    std::ostringstream err;
    err << "Master and displaced uncontact lists must have matching sizes: "
        << _uncontact_master.size() << " != " << _uncontact_displaced.size();
    mooseError(err.str());
  }
  if (!_uncontact_displaced.size())
    _uncontact_displaced.resize(_uncontact_master.size());
}

Member Function Documentation

setup()

void ContactSplit::setup ( const std::string &  prefix = "-" )

overridevirtual

Definition at line 89 of file ContactSplit.C.

{
  // A reference to the PetscOptions
  Moose::PetscSupport::PetscOptions & po = _fe_problem.getPetscOptions();
  // prefix
  std::string dmprefix = prefix + "dm_moose_", opt, val;

  // contacts options
  if (_contact_master.size())
  {
    opt = dmprefix + "ncontacts";
    {
      std::ostringstream oval;
      oval << _contact_master.size();
      val = oval.str();
    }
    // push back PETSc options
    if (val == "")
      po.flags.push_back(opt);
    else
    {
      po.inames.push_back(opt);
      po.values.push_back(val);
    }
    for (unsigned int j = 0; j < _contact_master.size(); ++j)
    {
      std::ostringstream oopt;
      oopt << dmprefix << "contact_" << j;
      opt = oopt.str();
      val = _contact_master[j] + "," + _contact_slave[j];
      // push back PETSc options
      if (val == "")
        po.flags.push_back(opt);
      else
      {
        po.inames.push_back(opt);
        po.values.push_back(val);
      }
      if (_contact_displaced[j])
      {
        opt = opt + "_displaced";
        val = "yes";
        // push back PETSc options
        if (val == "")
          po.flags.push_back(opt);
        else
        {
          po.inames.push_back(opt);
          po.values.push_back(val);
        }
      }
    }
  }
  // uncontacts options
  if (_uncontact_master.size())
  {
    opt = dmprefix + "nuncontacts";
    {
      std::ostringstream oval;
      oval << _uncontact_master.size();
      val = oval.str();
    }
    // push back PETSc options
    if (val == "")
      po.flags.push_back(opt);
    else
    {
      po.inames.push_back(opt);
      po.values.push_back(val);
    }
    for (unsigned int j = 0; j < _uncontact_master.size(); ++j)
    {
      std::ostringstream oopt;
      oopt << dmprefix << "uncontact_" << j;
      opt = oopt.str();
      val = _uncontact_master[j] + "," + _uncontact_slave[j];
      // push back PETSc options
      if (val == "")
        po.flags.push_back(opt);
      else
      {
        po.inames.push_back(opt);
        po.values.push_back(val);
      }
      if (_uncontact_displaced[j])
      {
        opt = opt + "_displaced";
        val = "yes";
        // push back PETSc options
        if (val == "")
          po.flags.push_back(opt);
        else
        {
          po.inames.push_back(opt);
          po.values.push_back(val);
        }
      }
    }
  }

  // Whether to include all nodes on the contact surfaces
  // into the contact subsolver
  opt = dmprefix + "includeAllContactNodes";
  if (_include_all_contact_nodes)
    val = "yes";
  else
    val = "no";
  po.inames.push_back(opt);
  po.values.push_back(val);
  Split::setup(prefix);
}

Member Data Documentation

_contact_displaced

std::vector<int> ContactSplit::_contact_displaced

protected

Definition at line 34 of file ContactSplit.h.

Referenced by ContactSplit(), and setup().

_contact_master

std::vector<std::string> ContactSplit::_contact_master

protected

Definition at line 32 of file ContactSplit.h.

Referenced by ContactSplit(), and setup().

_contact_slave

std::vector<std::string> ContactSplit::_contact_slave

protected

Definition at line 33 of file ContactSplit.h.

Referenced by ContactSplit(), and setup().

_include_all_contact_nodes

bool ContactSplit::_include_all_contact_nodes

protected

Definition at line 38 of file ContactSplit.h.

Referenced by setup().

_uncontact_displaced

std::vector<int> ContactSplit::_uncontact_displaced

protected

Definition at line 37 of file ContactSplit.h.

Referenced by ContactSplit(), and setup().

_uncontact_master

std::vector<std::string> ContactSplit::_uncontact_master

protected

Definition at line 35 of file ContactSplit.h.

Referenced by ContactSplit(), and setup().

_uncontact_slave

std::vector<std::string> ContactSplit::_uncontact_slave

protected

Definition at line 36 of file ContactSplit.h.

Referenced by ContactSplit(), and setup().

---

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

• contact/include/ContactSplit.h
• contact/src/ContactSplit.C