130
130
const Material * material,
133
134
// call base class to get electron_temperature terms
134
135
WeakEquationElectronTemperature::N_integrand(fields, grad_fields, material, flux);
137
139
material->electric_field(fields, grad_fields, _E_);
138
140
jouleHeating = _J_[0].mult_by_element(_E_[0]);
139
141
for (DENS_MAT_VEC::size_type i=1; i < _J_.size(); i++)
140
142
jouleHeating += _J_[i].mult_by_element(_E_[i]);