FEFaceValues<dim> fe_values (fe, q,
update_values |
update_gradients |
- update_second_derivatives |
+ update_hessians |
update_quadrature_points |
update_jacobians);
fe_values.reinit (dof_handler.begin_active(), 0);
:
fe_values (fe, quadrature,
update_values |
- update_second_derivatives |
+ update_hessians |
update_quadrature_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_quadrature_points | update_JxW_values);
FEValues<dim> neighbor_fe_values (fe, quadrature,
:
fe_values (fe, quadrature,
update_values |
- update_second_derivatives |
+ update_hessians |
update_quadrature_points |
update_JxW_values),
right_hand_side (&right_hand_side)
hp::QCollection<dim> quadrature(q);
hp::FEValues<dim> fe_values (fe, quadrature,
update_values | update_gradients |
- update_second_derivatives |
+ update_hessians |
update_quadrature_points | update_JxW_values);
hp::FEValues<dim> neighbor_fe_values (fe, quadrature,
QTrapez<dim> q;
// QIterated<dim> q(q_trapez, div);
FEValues<dim> fe(mapping, finel, q, UpdateFlags(update_gradients
- | update_second_derivatives));
+ | update_hessians));
fe.reinit(c);
unsigned int k=0;
QTrapez<dim> q;
// QIterated<dim> q(q_trapez, div);
FEValues<dim> fe(mapping, finel, q, UpdateFlags(update_gradients
- | update_second_derivatives));
+ | update_hessians));
fe.reinit(c);
unsigned int k=0;
QGauss<dim-1> q(1);
// QIterated<dim> q(q_trapez, div);
FEFaceValues<dim> fe(mapping, finel, q, UpdateFlags(update_gradients
- | update_second_derivatives));
+ | update_hessians));
for (unsigned int face=0; face<GeometryInfo<dim>::faces_per_cell; ++face)
{
fe.reinit(c, face);
QIterated<dim> q(q_trapez, div);
FEValues<dim> fe(mapping, finel, q, UpdateFlags(update_values
| update_gradients
- | update_second_derivatives));
+ | update_hessians));
sprintf(fname, "Cell%dd-%s", dim, name);
// cerr << "\n" << fname << "\n";
const char *name,
UpdateFlags uflags = UpdateFlags(update_values
| update_gradients
- | update_second_derivatives))
+ | update_hessians))
{
Triangulation<dim> tr;
DoFHandler<dim> dof(tr);
dof.distribute_dofs(fe);
const UpdateFlags update_all = (update_values | update_gradients |
- update_second_derivatives);
+ update_hessians);
// first check this for FEValues
// objects
UpdateFlags(update_values |
update_JxW_values |
update_gradients |
- update_second_derivatives));
+ update_hessians));
typename DoFHandler<dim>::cell_iterator c = dof.begin();
fe_values.reinit(c);
const QGauss<1> quad (fe_degree+1);
typename MatrixFree<dim,number>::AdditionalData data;
data.tasks_parallel_scheme = MatrixFree<dim,number>::AdditionalData::none;
- data.mapping_update_flags = update_gradients | update_second_derivatives;
+ data.mapping_update_flags = update_gradients | update_hessians;
mf_data.reinit (dof, constraints, quad, data);
}
MappingQ<dim> mapping (2);
typename MatrixFree<dim,number>::AdditionalData data;
data.tasks_parallel_scheme = MatrixFree<dim,number>::AdditionalData::none;
- data.mapping_update_flags = update_gradients | update_second_derivatives;
+ data.mapping_update_flags = update_gradients | update_hessians;
mf_data.reinit (mapping, dof, constraints, quad, data);
MatrixFreeTest<dim,fe_degree,fe_degree+1,number> mf (mf_data, mapping);
mf.test_functions (solution);
const QGauss<1> quad (fe_degree+1);
typename MatrixFree<dim,number>::AdditionalData data;
data.tasks_parallel_scheme = MatrixFree<dim,number>::AdditionalData::none;
- data.mapping_update_flags = update_gradients | update_second_derivatives;
+ data.mapping_update_flags = update_gradients | update_hessians;
mf_data.reinit (mapping, dof, constraints, quad, data);
}
MappingQ<dim> mapping (4);
typename MatrixFree<dim,number>::AdditionalData data;
data.tasks_parallel_scheme = MatrixFree<dim,number>::AdditionalData::none;
- data.mapping_update_flags = update_gradients | update_second_derivatives;
+ data.mapping_update_flags = update_gradients | update_hessians;
mf_data.reinit (mapping, dof, constraints, quad, data);
MatrixFreeTest<dim,fe_degree,fe_degree+1,number> mf (mf_data, mapping);
mf.test_functions (solution);
MappingQ<dim> mapping (2);
typename MatrixFree<dim,number>::AdditionalData data;
data.tasks_parallel_scheme = MatrixFree<dim,number>::AdditionalData::none;
- data.mapping_update_flags = update_gradients | update_second_derivatives;
+ data.mapping_update_flags = update_gradients | update_hessians;
mf_data.reinit (mapping, dof, constraints, quad, data);
MatrixFreeTest<dim,fe_degree,n_q_points_1d,number> mf (mf_data, mapping);
mf.test_functions (solution);
const QGauss<1> quad (fe_degree+1);
typename MatrixFree<dim>::AdditionalData data;
data.tasks_parallel_scheme = MatrixFree<dim>::AdditionalData::none;
- data.mapping_update_flags = update_gradients | update_second_derivatives;
+ data.mapping_update_flags = update_gradients | update_hessians;
mf_data.reinit (dof, constraints, quad, data);
}
MatrixFree<dim,number>::AdditionalData::partition_partition;
}
data.tasks_block_size = 7;
- data.mapping_update_flags |= update_second_derivatives;
+ data.mapping_update_flags |= update_hessians;
mf_data.reinit (dof, constraints, quad, data);
}
FEValues<dim> fe_values (dof.get_fe(), quadrature_formula,
update_values | update_gradients |
- update_JxW_values | update_second_derivatives);
+ update_JxW_values | update_hessians);
const unsigned int dofs_per_cell = dof.get_fe().dofs_per_cell;
const unsigned int n_q_points = quadrature_formula.size();
const QGauss<1> quad (fe_degree+1);
typename MatrixFree<dim,number>::AdditionalData data;
data.tasks_parallel_scheme = MatrixFree<dim,number>::AdditionalData::none;
- data.mapping_update_flags = update_gradients | update_second_derivatives;
+ data.mapping_update_flags = update_gradients | update_hessians;
data.initialize_indices = false;
mf_data.reinit (dof, constraints, quad, data);
}
const QGauss<1> quad (fe_degree+1);
typename MatrixFree<dim,double>::AdditionalData data;
data.tasks_parallel_scheme = MatrixFree<dim,double>::AdditionalData::none;
- data.mapping_update_flags = update_gradients | update_second_derivatives;
+ data.mapping_update_flags = update_gradients | update_hessians;
if (i==1)
data.initialize_indices = false;
mf_data.reinit (mapping, dof, constraints, quad, data);
FEValues<dim> fe2 (mapping, element, quadrature,
update_values | update_gradients
- | update_second_derivatives
+ | update_hessians
| update_quadrature_points | update_JxW_values);
double l2 = 0.;