void
ComputeResidualThread::subdomainChanged()
{
_fe_problem.subdomainSetup(_subdomain, _tid);
_sys._kernels[_tid].updateActiveKernels(_subdomain);
if (_sys._doing_dg)
_sys._dg_kernels[_tid].updateActiveDGKernels(_fe_problem.time(), _fe_problem.dt());
std::set<MooseVariable *> needed_moose_vars;
const std::vector<KernelBase *> & kernels = _sys._kernels[_tid].active();
for (std::vector<KernelBase *>::const_iterator it = kernels.begin(); it != kernels.end(); ++it)
{
const std::set<MooseVariable *> & mv_deps = (*it)->getMooseVariableDependencies();
needed_moose_vars.insert(mv_deps.begin(), mv_deps.end());
}
// Boundary Condition Dependencies
const std::set<unsigned int> & subdomain_boundary_ids = _mesh.getSubdomainBoundaryIds(_subdomain);
for (std::set<unsigned int>::const_iterator id_it = subdomain_boundary_ids.begin();
id_it != subdomain_boundary_ids.end();
++id_it)
{
std::vector<IntegratedBC *> bcs = _sys._bcs[_tid].activeIntegrated(*id_it);
if (bcs.size() > 0)
{
for (std::vector<IntegratedBC *>::iterator it = bcs.begin(); it != bcs.end(); ++it)
{
IntegratedBC * bc = (*it);
if (bc->shouldApply())
{
const std::set<MooseVariable *> & mv_deps = bc->getMooseVariableDependencies();
needed_moose_vars.insert(mv_deps.begin(), mv_deps.end());
}
}
}
}
// DG Kernel dependencies
{
std::vector<DGKernel *> dgks = _sys._dg_kernels[_tid].active();
for (std::vector<DGKernel *>::iterator it = dgks.begin(); it != dgks.end(); ++it)
{
const std::set<MooseVariable *> & mv_deps = (*it)->getMooseVariableDependencies();
needed_moose_vars.insert(mv_deps.begin(), mv_deps.end());
}
}
_fe_problem.setActiveElementalMooseVariables(needed_moose_vars, _tid);
_fe_problem.prepareMaterials(_subdomain, _tid);
}
