UELThread.h

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://mooseframework.inl.gov
//*
//* All rights reserved, see COPYRIGHT for full restrictions
//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
//*
//* Licensed under LGPL 2.1, please see LICENSE for details
//* https://www.gnu.org/licenses/lgpl-2.1.html
 
#pragma once
 
#include "ThreadedElementLoop.h"
#include "AbaqusUserElement.h"
#include "libmesh/elem_range.h"
#include "AuxiliarySystem.h"
 
class UELThread : public ThreadedElementLoop<ConstElemRange>
{
public:
  UELThread(FEProblemBase & fe_problem, AbaqusUserElement & uel_uo);
 
  // Splitting Constructor
  UELThread(UELThread & x, Threads::split split);
 
  void join(const UELThread &) {}
 
protected:
  virtual void subdomainChanged() override final;
  virtual void onElement(const Elem * elem) override final;
 
private:
  NonlinearSystemBase & _sys;
 
  AuxiliarySystem * _aux_sys;
 
  std::vector<dof_id_type> _var_dof_indices;
  std::vector<dof_id_type> _all_dof_indices;
  std::vector<Real> _all_dof_values;
  std::vector<Real> _all_dof_increments;
  std::vector<Real> _all_udot_dof_values;
  std::vector<Real> _all_udotdot_dof_values;
 
  // vectors for aux variables
  std::vector<dof_id_type> _aux_var_dof_indices;
  std::vector<dof_id_type> _all_aux_var_dof_indices;
  std::vector<Real> _all_aux_var_dof_increments;
  std::vector<Real> _all_aux_var_dof_values;
  std::vector<Real> _aux_var_values_to_uel;
 
  DenseVector<Real> _local_re;
  DenseMatrix<Real> _local_ke;
 
  const std::vector<const MooseVariableFieldBase *> & _variables;
 
  const std::vector<const MooseVariableFieldBase *> & _aux_variables;
 
  AbaqusUserElement & _uel_uo;
 
  const AbaqusUserElement::uel_t & _uel;
 
  std::array<int, 5> _lflags;
  std::vector<Real> _coords;
 
  std::vector<Real> _statev_copy;
};