// $Id$
// Version: $Name$
//
-// Copyright (C) 2004, 2005, 2006 by the deal.II authors
+// Copyright (C) 2004, 2005, 2006, 2007 by the deal.II authors
//
// This file is subject to QPL and may not be distributed
// without copyright and license information. Please refer
* Mapping object used in
* conjunction with this finite
* element. <li> if
- * <tt>update_second_derivatives</tt>
+ * <tt>update_hessians</tt>
* is set, the result will
* contain
- * <tt>update_second_derivatives</tt>
+ * <tt>update_hessians</tt>
* and
* <tt>update_covariant_transformation</tt>.
* The rationale is the same as
// $Id$
// Version: $Name$
//
-// Copyright (C) 2006 by the deal.II authors
+// Copyright (C) 2006, 2007 by the deal.II authors
//
// This file is subject to QPL and may not be distributed
// without copyright and license information. Please refer
if (flags & update_gradients)
out |= update_gradients | update_covariant_transformation;
- if (flags & update_second_derivatives)
- out |= update_second_derivatives | update_covariant_transformation;
+ if (flags & update_hessians)
+ out |= update_hessians | update_covariant_transformation;
return out;
}
// finite differencing is required,
// then initialize some objects for
// that
- if (flags & update_second_derivatives)
+ if (flags & update_hessians)
data->initialize_2nd (this, mapping, quadrature);
// next already fill those fields
}
const typename QProjector<dim>::DataSetDescriptor dsd;
- if (flags & update_second_derivatives)
+ if (flags & update_hessians)
this->compute_2nd (mapping, cell, dsd.cell(),
mapping_data, fe_data, data);
}
}
}
- if (flags & update_second_derivatives)
+ if (flags & update_hessians)
this->compute_2nd (mapping, cell, offset, mapping_data, fe_data, data);
}
}
}
- if (flags & update_second_derivatives)
+ if (flags & update_hessians)
this->compute_2nd (mapping, cell, offset, mapping_data, fe_data, data);
}
// $Id$
// Version: $Name$
//
-// Copyright (C) 1999, 2000, 2001, 2002, 2003, 2004, 2005, 2006 by the deal.II authors
+// Copyright (C) 1999, 2000, 2001, 2002, 2003, 2004, 2005, 2006, 2007 by the deal.II authors
//
// This file is subject to QPL and may not be distributed
// without copyright and license information. Please refer
* second derivatives shall
* be computed.
*/
- bool compute_second_derivatives;
+ bool compute_hessians;
/**
* Gives write-access to the
| update_JxW_values;
//TODO: Set update flags appropriately and figure out, how the second
// derivatives for the RT elements can be computed correctly.
- if (flags & update_second_derivatives)
- out |= update_second_derivatives | update_contravariant_transformation;
+ if (flags & update_hessians)
+ out |= update_hessians | update_contravariant_transformation;
return out;
}
out |= update_values | update_covariant_transformation;
if (flags & update_gradients)
out |= update_gradients | update_covariant_transformation;
- if (flags & update_second_derivatives)
- out |= update_second_derivatives | update_covariant_transformation;
+ if (flags & update_hessians)
+ out |= update_hessians | update_covariant_transformation;
return out;
}
// finite differencing is required,
// then initialize some objects for
// that
- if (flags & update_second_derivatives)
+ if (flags & update_hessians)
data->initialize_2nd (this, mapping, quadrature);
// next already fill those fields
};
}
- if (flags & update_second_derivatives)
+ if (flags & update_hessians)
this->compute_2nd (mapping, cell,
QProjector<dim>::DataSetDescriptor::cell(),
mapping_data, fe_data, data);
};
}
- if (flags & update_second_derivatives)
+ if (flags & update_hessians)
this->compute_2nd (mapping, cell, offset, mapping_data, fe_data, data);
}
};
}
- if (flags & update_second_derivatives)
+ if (flags & update_hessians)
this->compute_2nd (mapping, cell, offset, mapping_data, fe_data, data);
}
// finite differencing is required,
// then initialize some objects for
// that
- if (flags & update_second_derivatives)
+ if (flags & update_hessians)
{
// grad_grads.resize (this->dofs_per_cell);
data->initialize_2nd (this, mapping, quadrature);
}
}
- if (flags & update_second_derivatives)
+ if (flags & update_hessians)
this->compute_2nd (mapping, cell,
typename QProjector<dim>::DataSetDescriptor().cell(),
mapping_data, fe_data, data);
}
}
- if (flags & update_second_derivatives)
+ if (flags & update_hessians)
this->compute_2nd (mapping, cell, offset, mapping_data, fe_data, data);
}
}
}
- if (flags & update_second_derivatives)
+ if (flags & update_hessians)
this->compute_2nd (mapping, cell, offset, mapping_data, fe_data, data);
}
out |= update_values | update_covariant_transformation;
if (flags & update_gradients)
out |= update_gradients | update_covariant_transformation;
- if (flags & update_second_derivatives)
- out |= update_second_derivatives | update_covariant_transformation;
+ if (flags & update_hessians)
+ out |= update_hessians | update_covariant_transformation;
return out;
}
| update_JxW_values;
//TODO: Set update flags appropriately and figure out, how the second
// derivatives for the RT elements can be computed correctly.
- if (flags & update_second_derivatives)
- out |= update_second_derivatives | update_contravariant_transformation;
+ if (flags & update_hessians)
+ out |= update_hessians | update_contravariant_transformation;
return out;
}
// therefore have to set the
// respective flag since the base
// elements don't have them
- if (flags & update_second_derivatives)
- out |= update_second_derivatives | update_covariant_transformation;
+ if (flags & update_hessians)
+ out |= update_hessians | update_covariant_transformation;
return out;
}
// finite differencing are required,
// then initialize some objects for
// that
- data->compute_second_derivatives = flags & update_second_derivatives;
- if (data->compute_second_derivatives)
+ data->compute_hessians = flags & update_hessians;
+ if (data->compute_hessians)
{
// delete
- // update_second_derivatives
+ // update_hessians
// from flags list
- sub_flags = UpdateFlags (sub_flags ^ update_second_derivatives);
+ sub_flags = UpdateFlags (sub_flags ^ update_hessians);
data->initialize_2nd (this, mapping, quadrature);
}
// make sure that *we* compute
// second derivatives, base
// elements should not do it
- Assert (!(base_fe_data->update_each & update_second_derivatives),
+ Assert (!(base_fe_data->update_each & update_hessians),
ExcInternalError());
- Assert (!(base_fe_data->update_once & update_second_derivatives),
+ Assert (!(base_fe_data->update_once & update_hessians),
ExcInternalError());
// The FEValuesData @p{data}
// should not
// have computed
// them!
- Assert (!(base_flags & update_second_derivatives),
+ Assert (!(base_flags & update_hessians),
ExcInternalError());
};
};
}
}
- if (fe_data.compute_second_derivatives)
+ if (fe_data.compute_hessians)
{
unsigned int offset = 0;
if (face_no != invalid_face_number)
if (update_flags & update_gradients)
this->dof_data[dataset]->get_function_gradients (fe_patch_values,
data.patch_gradients);
- if (update_flags & update_second_derivatives)
+ if (update_flags & update_hessians)
this->dof_data[dataset]->get_function_second_derivatives (fe_patch_values,
data.patch_second_derivatives);
postprocessor->compute_derived_quantities(data.postprocessed_values[dataset],
if (update_flags & update_gradients)
this->dof_data[dataset]->get_function_gradients (fe_patch_values,
data.patch_gradients_system);
- if (update_flags & update_second_derivatives)
+ if (update_flags & update_hessians)
this->dof_data[dataset]->get_function_second_derivatives (fe_patch_values,
data.patch_second_derivatives_system);
postprocessor->compute_derived_quantities(data.postprocessed_values[dataset],
// $Id$
// Version: $Name$
//
-// Copyright (C) 2000, 2001, 2002, 2003, 2004, 2005, 2006 by the deal.II authors
+// Copyright (C) 2000, 2001, 2002, 2003, 2004, 2005, 2006, 2007 by the deal.II authors
//
// This file is subject to QPL and may not be distributed
// without copyright and license information. Please refer
if (update_flags & update_gradients)
this->dof_data[dataset]->get_function_gradients (fe_patch_values,
data.patch_gradients);
- if (update_flags & update_second_derivatives)
+ if (update_flags & update_hessians)
this->dof_data[dataset]->get_function_second_derivatives (fe_patch_values,
data.patch_second_derivatives);
postprocessor->compute_derived_quantities(data.postprocessed_values[dataset],
if (update_flags & update_gradients)
this->dof_data[dataset]->get_function_gradients (fe_patch_values,
data.patch_gradients_system);
- if (update_flags & update_second_derivatives)
+ if (update_flags & update_hessians)
this->dof_data[dataset]->get_function_second_derivatives (fe_patch_values,
data.patch_second_derivatives_system);
postprocessor->compute_derived_quantities(data.postprocessed_values[dataset],
if (update_flags & update_gradients)
this->dof_data[dataset]->get_function_gradients (fe_patch_values,
data.patch_gradients);
- if (update_flags & update_second_derivatives)
+ if (update_flags & update_hessians)
this->dof_data[dataset]->get_function_second_derivatives (fe_patch_values,
data.patch_second_derivatives);
postprocessor->compute_derived_quantities(data.postprocessed_values[dataset],
if (update_flags & update_gradients)
this->dof_data[dataset]->get_function_gradients (fe_patch_values,
data.patch_gradients_system);
- if (update_flags & update_second_derivatives)
+ if (update_flags & update_hessians)
this->dof_data[dataset]->get_function_second_derivatives (fe_patch_values,
data.patch_second_derivatives_system);
postprocessor->compute_derived_quantities(data.postprocessed_values[dataset],
const UpdateFlags DerivativeApproximation::SecondDerivative<dim>::update_flags = update_gradients;
template <int dim>
-const UpdateFlags DerivativeApproximation::ThirdDerivative<dim>::update_flags = update_second_derivatives;
+const UpdateFlags DerivativeApproximation::ThirdDerivative<dim>::update_flags = update_hessians;
:
fe_values (fe, quadrature,
update_values |
- update_second_derivatives |
+ update_hessians |
update_q_points |
update_JxW_values),
right_hand_side (&right_hand_side)
QTrapez<dim> quadrature;
FEValues<dim> fe_values (fe, quadrature,
update_values | update_gradients |
- update_second_derivatives |
+ update_hessians |
update_q_points | update_JxW_values);
// The error indicator formula presented in