Real LevelSetAdvectionSUPG::computeQpJacobian() { computeQpVelocity(); Real tau = _current_elem->hmin() / (2 * _velocity.norm()); return (tau * _velocity * _grad_test[_i][_qp]) * (_velocity * _grad_phi[_j][_qp]); }
Real LevelSetAdvection::computeQpJacobian() { computeQpVelocity(); return _test[_i][_qp] * (_velocity * _grad_phi[_j][_qp]); }
Real LevelSetAdvection::computeQpResidual() { computeQpVelocity(); return _test[_i][_qp] * (_velocity * _grad_u[_qp]); }