From: guido Date: Thu, 30 May 2002 18:53:38 +0000 (+0000) Subject: access to internals of FEValues eliminated X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=415bd2c718bc4df73a2fffc64863d8290821d66f;p=dealii-svn.git access to internals of FEValues eliminated git-svn-id: https://svn.dealii.org/trunk@5947 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/examples/step-12/step-12.cc b/deal.II/examples/step-12/step-12.cc index 399e4fbc4e..07dfedcbdc 100644 --- a/deal.II/examples/step-12/step-12.cc +++ b/deal.II/examples/step-12/step-12.cc @@ -273,8 +273,6 @@ void DGTransportEquation::assemble_cell_term( // First we ask ``fe_v'' for the // shape gradients, shape values and // quadrature weights, - const std::vector > > &grad_v = fe_v.get_shape_grads (); - const FullMatrix &v = fe_v.get_shape_values (); const std::vector &JxW = fe_v.get_JxW_values (); // Then the flow field beta and the @@ -294,11 +292,11 @@ void DGTransportEquation::assemble_cell_term( for (unsigned int i=0; i::assemble_boundary_term( // function, we ask the ``FEValues'' // object for the shape values and // the quadrature weights - const FullMatrix &v = fe_v.get_shape_values (); const std::vector &JxW = fe_v.get_JxW_values (); // but here also for the normals. const std::vector > &normals = fe_v.get_normal_vectors (); @@ -356,8 +353,8 @@ void DGTransportEquation::assemble_boundary_term( for (unsigned int i=0; i::assemble_boundary_term( for (unsigned int i=0; i::assemble_face_term1( // objects for the shape values, // the quadrature weights and the // normals - const FullMatrix &v = fe_v.get_shape_values (); - const FullMatrix &v_neighbor = fe_v_neighbor.get_shape_values (); const std::vector &JxW = fe_v.get_JxW_values (); const std::vector > &normals = fe_v.get_normal_vectors (); @@ -444,8 +439,8 @@ void DGTransportEquation::assemble_face_term1( for (unsigned int i=0; i::assemble_face_term1( for (unsigned int i=0; i::assemble_face_term2( FullMatrix &un_vn_matrix) const { // the first few lines are the same - const FullMatrix &v = fe_v.get_shape_values (); - const FullMatrix &v_neighbor = fe_v_neighbor.get_shape_values (); const std::vector &JxW = fe_v.get_JxW_values (); const std::vector > &normals = fe_v.get_normal_vectors (); @@ -506,8 +499,8 @@ void DGTransportEquation::assemble_face_term2( for (unsigned int i=0; i::assemble_face_term2( for (unsigned int k=0; k::assemble_face_term2( for (unsigned int i=0; i::assemble_face_term2( for (unsigned int k=0; k