CHInterface::CHInterface(const std::string & name, InputParameters parameters) :
DerivativeMaterialInterface<JvarMapInterface<Kernel> >(name, parameters),
_kappa(getMaterialProperty<Real>("kappa_name")),
_M(getMaterialProperty<Real>("mob_name")),
_dMdc(getMaterialPropertyDerivative<Real>("mob_name", _var.name())),
_d2Mdc2(getMaterialPropertyDerivative<Real>("mob_name", _var.name(), _var.name())),
_second_u(second()),
_second_test(secondTest()),
_second_phi(secondPhi()),
_nvar(_coupled_moose_vars.size()),
_dMdarg(_nvar),
_d2Mdcdarg(_nvar),
_d2Mdargdarg(_nvar),
_coupled_grad_vars(_nvar)
{
// Iterate over all coupled variables
for (unsigned int i = 0; i < _nvar; ++i)
{
//Set material property values
_dMdarg[i] = &getMaterialPropertyDerivative<Real>("mob_name", _coupled_moose_vars[i]->name());
_d2Mdcdarg[i] = &getMaterialPropertyDerivative<Real>("mob_name", _var.name(), _coupled_moose_vars[i]->name());
_d2Mdargdarg[i].resize(_nvar);
for (unsigned int j = 0; j < _nvar; ++j)
_d2Mdargdarg[i][j] = &getMaterialPropertyDerivative<Real>("mob_name", _coupled_moose_vars[i]->name(), _coupled_moose_vars[j]->name());
//Set coupled variable gradients
_coupled_grad_vars[i] = &coupledGradient("args", i);
}
}
CHInterface::CHInterface(const std::string & name, InputParameters parameters) :
Kernel(name, parameters),
_kappa_name(getParam<std::string>("kappa_name")),
_mob_name(getParam<std::string>("mob_name")),
_Dmob_name(getParam<std::string>("Dmob_name")),
_grad_mob_name(getParam<std::string>("grad_mob_name")),
_kappa(getMaterialProperty<Real>(_kappa_name)),
_M(getMaterialProperty<Real>(_mob_name)),
_has_MJac(getParam<bool>("has_MJac")),
_DM(_has_MJac ? &getMaterialProperty<Real>(_Dmob_name) : NULL),
_grad_M(getMaterialProperty<RealGradient>(_grad_mob_name)),
_Dgrad_Mnp(_has_MJac ? &getMaterialProperty<RealGradient>("Dgrad_Mnp") : NULL),
_Dgrad_Mngp(_has_MJac ? &getMaterialProperty<Real>("Dgrad_Mngp") : NULL),
_second_u(second()),
_second_test(secondTest()),
_second_phi(secondPhi())
{
}
bool RC4Test::doTest()
{
bool ret1 = firstTest();
bool ret2 = secondTest();
return ret1 && ret2;
}