From c07c3637da0eb2e6d394dd90fc0e84848431c728 Mon Sep 17 00:00:00 2001 From: Jaeryun Yim Date: Thu, 29 Sep 2016 01:22:39 +0900 Subject: [PATCH] Remove the parts for a graphical issue already resolved. --- source/fe/fe_p1nc.cc | 73 -------------------------------------------- 1 file changed, 73 deletions(-) diff --git a/source/fe/fe_p1nc.cc b/source/fe/fe_p1nc.cc index 6aa07ade45..310b4a5297 100644 --- a/source/fe/fe_p1nc.cc +++ b/source/fe/fe_p1nc.cc @@ -218,23 +218,6 @@ FE_P1NC::fill_fe_values (const Triangulation<2,2>::cell_iterator &cell } } } - - // When this function is called for the graphic out, - // MappingRelatedData does not work properly. - // In this case, the quadrature points on the real cell is computed in manual sense, using 'mapping' and 'quadrature'. - // This is a temporary solution. It needs to be fixed fundamentally. - if (n_q_points==0) - { - for (unsigned int i=0; idofs_per_cell; ++k) - { - Point<2> realquadrature ; - - realquadrature = mapping.transform_unit_to_real_cell(cell, quadrature.point(i)) ; - values[k] = coeffs[k][0]*realquadrature(0) + coeffs[k][1]*realquadrature(1) + coeffs[k][2] ; - output_data.shape_values[k][i] = values[k]; - } - } } @@ -281,34 +264,6 @@ FE_P1NC::fill_fe_face_values (const Triangulation<2,2>::cell_iterator } } - - // When this function is called for computation of facial jump residual, - // MappingRelatedData does not work properly. - // In this case, the quadrature points on the real cell is computed in manual sense, using 'mapping' and 'quadrature'. - // This is a temporary solution. It needs to be fixed fundamentally. - if (n_q_points==0) - { - Quadrature<2> cellquadrature = QProjector<2>::project_to_face(quadrature, face_no) ; - for (unsigned int i=0; idofs_per_cell; ++k) - { - if (flags & update_values) - { - Point<2> realquadrature ; - - realquadrature = mapping.transform_unit_to_real_cell(cell, cellquadrature.point(i)) ; - values[k] = coeffs[k][0]*realquadrature(0) + coeffs[k][1]*realquadrature(1) + coeffs[k][2] ; - output_data.shape_values[k][i] = values[k]; - } - - if (flags & update_gradients) - { - grads[k][0] = coeffs[k][0] ; - grads[k][1] = coeffs[k][1] ; - output_data.shape_gradients[k][i] = grads[k]; - } - } - } } @@ -355,34 +310,6 @@ FE_P1NC::fill_fe_subface_values (const Triangulation<2,2>::cell_iterator } } } - - // When this function is called for computation of facial jump residual, - // MappingRelatedData does not work properly. - // In this case, the quadrature points on the real cell is computed in manual sense, using 'mapping' and 'quadrature'. - // This is a temporary solution. It needs to be fixed fundamentally. - if (n_q_points==0) - { - Quadrature<2> cellquadrature = QProjector<2>::project_to_subface(quadrature, face_no, sub_no) ; - for (unsigned int i=0; idofs_per_cell; ++k) - { - if (flags & update_values) - { - Point<2> realquadrature ; - - realquadrature = mapping.transform_unit_to_real_cell(cell, cellquadrature.point(i)) ; - values[k] = coeffs[k][0]*realquadrature(0) + coeffs[k][1]*realquadrature(1) + coeffs[k][2] ; - output_data.shape_values[k][i] = values[k]; - } - - if (flags & update_gradients) - { - grads[k][0] = coeffs[k][0] ; - grads[k][1] = coeffs[k][1] ; - output_data.shape_gradients[k][i] = grads[k]; - } - } - } } -- 2.39.5