From 1e9709e86af027f8e8b74237b2e12934200545bd Mon Sep 17 00:00:00 2001 From: =?utf8?q?Micha=C5=82=20Wichrowski?= Date: Fri, 2 May 2025 18:40:31 +0200 Subject: [PATCH] resolve comments --- include/deal.II/lac/full_matrix.h | 8 ++++---- include/deal.II/numerics/tensor_product_matrix_creator.h | 5 ++--- 2 files changed, 6 insertions(+), 7 deletions(-) diff --git a/include/deal.II/lac/full_matrix.h b/include/deal.II/lac/full_matrix.h index 2f9ae851be..6838eaec1c 100644 --- a/include/deal.II/lac/full_matrix.h +++ b/include/deal.II/lac/full_matrix.h @@ -984,10 +984,10 @@ public: * m x n matrix A and a p x q matrix B is an (m*p) x (n*q) matrix defined as: * * ``` - * A ⊗ B = | a11*B a12*B ... a1n*B | - * | a21*B a22*B ... a2n*B | - * | ... ... ... ... | - * | am1*B am2*B ... amn*B | + * A ⊗ B = | a11*B a12*B ... a1n*B | + * | a21*B a22*B ... a2n*B | + * | ... ... ... ... | + * | am1*B am2*B ... amn*B | * ``` * * where aij are the elements of the matrix A. diff --git a/include/deal.II/numerics/tensor_product_matrix_creator.h b/include/deal.II/numerics/tensor_product_matrix_creator.h index 60b672eb59..ab5159cdd4 100644 --- a/include/deal.II/numerics/tensor_product_matrix_creator.h +++ b/include/deal.II/numerics/tensor_product_matrix_creator.h @@ -556,7 +556,6 @@ namespace TensorProductMatrixCreator const unsigned int degree = fe.degree; const unsigned int n_dofs_per_cell = fe.dofs_per_cell; - const Number &JxW = h; QGauss<1> quadrature(degree + 1); FullMatrix cell_matrix(n_dofs_per_cell, n_dofs_per_cell); @@ -573,7 +572,7 @@ namespace TensorProductMatrixCreator cell_matrix(i - shift, j - shift) += (fe.shape_value(numbering[i], quadrature.point(q)) * fe.shape_value(numbering[j], quadrature.point(q))) * - JxW * quadrature.weight(q); + (h * quadrature.weight(q)); return cell_matrix; } @@ -620,7 +619,7 @@ namespace TensorProductMatrixCreator cell_matrix(i - shift, j - shift) += (fe.shape_grad(numbering[i], quadrature.point(q)) / h * fe.shape_grad(numbering[j], quadrature.point(q))) / - h * JxW * quadrature.weight(q); + h * (h * quadrature.weight(q)); return cell_matrix; } -- 2.39.5