From: Jonathan Robey Date: Thu, 1 Sep 2016 00:48:58 +0000 (-0700) Subject: Change update_second_derivatives to update_hessians X-Git-Tag: v8.5.0-rc1~717^2~2 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=194392f096277023d4a18f485a5fc871e604b2c5;p=dealii.git Change update_second_derivatives to update_hessians --- diff --git a/tests/bits/rt_1.cc b/tests/bits/rt_1.cc index e4fc58f0c4..42ec1686e3 100644 --- a/tests/bits/rt_1.cc +++ b/tests/bits/rt_1.cc @@ -49,7 +49,7 @@ void test (const unsigned int degree, FEFaceValues 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); diff --git a/tests/bits/step-14.cc b/tests/bits/step-14.cc index 9193f5588e..c24272fc0f 100644 --- a/tests/bits/step-14.cc +++ b/tests/bits/step-14.cc @@ -1474,7 +1474,7 @@ namespace LaplaceSolver : fe_values (fe, quadrature, update_values | - update_second_derivatives | + update_hessians | update_quadrature_points | update_JxW_values), right_hand_side (&right_hand_side) diff --git a/tests/bits/step-15.cc b/tests/bits/step-15.cc index fc6dbfba8f..5af0ed76eb 100644 --- a/tests/bits/step-15.cc +++ b/tests/bits/step-15.cc @@ -397,7 +397,7 @@ void MinimizationProblem<1>::refine_grid () QTrapez quadrature; FEValues fe_values (fe, quadrature, update_values | update_gradients | - update_second_derivatives | + update_hessians | update_quadrature_points | update_JxW_values); FEValues neighbor_fe_values (fe, quadrature, diff --git a/tests/fail/hp-step-14.cc b/tests/fail/hp-step-14.cc index 00e71c3fc0..c895147643 100644 --- a/tests/fail/hp-step-14.cc +++ b/tests/fail/hp-step-14.cc @@ -1472,7 +1472,7 @@ namespace LaplaceSolver : fe_values (fe, quadrature, update_values | - update_second_derivatives | + update_hessians | update_quadrature_points | update_JxW_values), right_hand_side (&right_hand_side) diff --git a/tests/fail/hp-step-15.cc b/tests/fail/hp-step-15.cc index c6ce14c34d..804b5590af 100644 --- a/tests/fail/hp-step-15.cc +++ b/tests/fail/hp-step-15.cc @@ -399,7 +399,7 @@ void MinimizationProblem<1>::refine_grid () hp::QCollection quadrature(q); hp::FEValues fe_values (fe, quadrature, update_values | update_gradients | - update_second_derivatives | + update_hessians | update_quadrature_points | update_JxW_values); hp::FEValues neighbor_fe_values (fe, quadrature, diff --git a/tests/fe/derivatives.cc b/tests/fe/derivatives.cc index e3ab65ee09..f63d6a75ff 100644 --- a/tests/fe/derivatives.cc +++ b/tests/fe/derivatives.cc @@ -51,7 +51,7 @@ plot_derivatives(Mapping &mapping, QTrapez q; // QIterated q(q_trapez, div); FEValues fe(mapping, finel, q, UpdateFlags(update_gradients - | update_second_derivatives)); + | update_hessians)); fe.reinit(c); unsigned int k=0; diff --git a/tests/fe/derivatives_bernstein.cc b/tests/fe/derivatives_bernstein.cc index 323d01dba7..caa7e4fc31 100644 --- a/tests/fe/derivatives_bernstein.cc +++ b/tests/fe/derivatives_bernstein.cc @@ -51,7 +51,7 @@ plot_derivatives(Mapping &mapping, QTrapez q; // QIterated q(q_trapez, div); FEValues fe(mapping, finel, q, UpdateFlags(update_gradients - | update_second_derivatives)); + | update_hessians)); fe.reinit(c); unsigned int k=0; diff --git a/tests/fe/derivatives_face.cc b/tests/fe/derivatives_face.cc index 180d981c3e..bee7d926d0 100644 --- a/tests/fe/derivatives_face.cc +++ b/tests/fe/derivatives_face.cc @@ -49,7 +49,7 @@ plot_derivatives(Mapping &mapping, QGauss q(1); // QIterated q(q_trapez, div); FEFaceValues fe(mapping, finel, q, UpdateFlags(update_gradients - | update_second_derivatives)); + | update_hessians)); for (unsigned int face=0; face::faces_per_cell; ++face) { fe.reinit(c, face); diff --git a/tests/fe/shapes.h b/tests/fe/shapes.h index 5b80b6c898..c81b99584c 100644 --- a/tests/fe/shapes.h +++ b/tests/fe/shapes.h @@ -44,7 +44,7 @@ plot_shape_functions(Mapping &mapping, QIterated q(q_trapez, div); FEValues 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"; @@ -106,7 +106,7 @@ plot_face_shape_functions( const char *name, UpdateFlags uflags = UpdateFlags(update_values | update_gradients - | update_second_derivatives)) + | update_hessians)) { Triangulation tr; DoFHandler dof(tr); @@ -389,7 +389,7 @@ void test_compute_functions (const Mapping &mapping, dof.distribute_dofs(fe); const UpdateFlags update_all = (update_values | update_gradients | - update_second_derivatives); + update_hessians); // first check this for FEValues // objects diff --git a/tests/mappings/mapping_q1_eulerian.cc b/tests/mappings/mapping_q1_eulerian.cc index 3edfd4b17b..6e28f74da8 100644 --- a/tests/mappings/mapping_q1_eulerian.cc +++ b/tests/mappings/mapping_q1_eulerian.cc @@ -65,7 +65,7 @@ show_values(FiniteElement &fe, UpdateFlags(update_values | update_JxW_values | update_gradients | - update_second_derivatives)); + update_hessians)); typename DoFHandler::cell_iterator c = dof.begin(); fe_values.reinit(c); diff --git a/tests/matrix_free/get_functions_common.h b/tests/matrix_free/get_functions_common.h index 039dd4bf8d..ddc289fff1 100644 --- a/tests/matrix_free/get_functions_common.h +++ b/tests/matrix_free/get_functions_common.h @@ -251,7 +251,7 @@ void do_test (const DoFHandler &dof, const QGauss<1> quad (fe_degree+1); typename MatrixFree::AdditionalData data; data.tasks_parallel_scheme = MatrixFree::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); } diff --git a/tests/matrix_free/get_functions_gl.cc b/tests/matrix_free/get_functions_gl.cc index 43ceb7b845..64fc569e69 100644 --- a/tests/matrix_free/get_functions_gl.cc +++ b/tests/matrix_free/get_functions_gl.cc @@ -85,7 +85,7 @@ void test () MappingQ mapping (2); typename MatrixFree::AdditionalData data; data.tasks_parallel_scheme = MatrixFree::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 mf (mf_data, mapping); mf.test_functions (solution); diff --git a/tests/matrix_free/get_functions_mappingq.cc b/tests/matrix_free/get_functions_mappingq.cc index 9b7f97d792..2165aa45af 100644 --- a/tests/matrix_free/get_functions_mappingq.cc +++ b/tests/matrix_free/get_functions_mappingq.cc @@ -85,7 +85,7 @@ void test () const QGauss<1> quad (fe_degree+1); typename MatrixFree::AdditionalData data; data.tasks_parallel_scheme = MatrixFree::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); } diff --git a/tests/matrix_free/get_functions_q_hierarchical.cc b/tests/matrix_free/get_functions_q_hierarchical.cc index 3434d2fba8..9cac15c86f 100644 --- a/tests/matrix_free/get_functions_q_hierarchical.cc +++ b/tests/matrix_free/get_functions_q_hierarchical.cc @@ -88,7 +88,7 @@ void test () MappingQ mapping (4); typename MatrixFree::AdditionalData data; data.tasks_parallel_scheme = MatrixFree::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 mf (mf_data, mapping); mf.test_functions (solution); diff --git a/tests/matrix_free/get_functions_rect.cc b/tests/matrix_free/get_functions_rect.cc index e875f40f6f..f42029c259 100644 --- a/tests/matrix_free/get_functions_rect.cc +++ b/tests/matrix_free/get_functions_rect.cc @@ -41,7 +41,7 @@ void sub_test (const DoFHandler &dof, MappingQ mapping (2); typename MatrixFree::AdditionalData data; data.tasks_parallel_scheme = MatrixFree::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 mf (mf_data, mapping); mf.test_functions (solution); diff --git a/tests/matrix_free/get_functions_variants.cc b/tests/matrix_free/get_functions_variants.cc index 072c6940af..f6bbf76f3c 100644 --- a/tests/matrix_free/get_functions_variants.cc +++ b/tests/matrix_free/get_functions_variants.cc @@ -178,7 +178,7 @@ void test () const QGauss<1> quad (fe_degree+1); typename MatrixFree::AdditionalData data; data.tasks_parallel_scheme = MatrixFree::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); } diff --git a/tests/matrix_free/matrix_vector_18.cc b/tests/matrix_free/matrix_vector_18.cc index 87a09d32e8..7f70b6a420 100644 --- a/tests/matrix_free/matrix_vector_18.cc +++ b/tests/matrix_free/matrix_vector_18.cc @@ -120,7 +120,7 @@ void do_test (const DoFHandler &dof, MatrixFree::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); } @@ -155,7 +155,7 @@ void do_test (const DoFHandler &dof, FEValues 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(); diff --git a/tests/matrix_free/no_index_initialize.cc b/tests/matrix_free/no_index_initialize.cc index 462f2bac30..17297273ca 100644 --- a/tests/matrix_free/no_index_initialize.cc +++ b/tests/matrix_free/no_index_initialize.cc @@ -135,7 +135,7 @@ void do_test (const DoFHandler &dof, const QGauss<1> quad (fe_degree+1); typename MatrixFree::AdditionalData data; data.tasks_parallel_scheme = MatrixFree::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); } diff --git a/tests/matrix_free/update_mapping_only.cc b/tests/matrix_free/update_mapping_only.cc index 1472c1ce08..85ebeb3717 100644 --- a/tests/matrix_free/update_mapping_only.cc +++ b/tests/matrix_free/update_mapping_only.cc @@ -82,7 +82,7 @@ void test () const QGauss<1> quad (fe_degree+1); typename MatrixFree::AdditionalData data; data.tasks_parallel_scheme = MatrixFree::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); diff --git a/tests/numerics/derivatives.cc b/tests/numerics/derivatives.cc index 94ab3da865..a596f7b155 100644 --- a/tests/numerics/derivatives.cc +++ b/tests/numerics/derivatives.cc @@ -118,7 +118,7 @@ check (const unsigned int level, FEValues fe2 (mapping, element, quadrature, update_values | update_gradients - | update_second_derivatives + | update_hessians | update_quadrature_points | update_JxW_values); double l2 = 0.;