From: wolf Date: Fri, 6 Nov 1998 15:43:40 +0000 (+0000) Subject: Update. X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=c060898ffa261aab5f609146757e932915f4fb73;p=dealii-svn.git Update. git-svn-id: https://svn.dealii.org/trunk@647 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/deal.II/source/fe/scripts/1d/lagrange b/deal.II/deal.II/source/fe/scripts/1d/lagrange index c9a255a711..b57a028a09 100644 --- a/deal.II/deal.II/source/fe/scripts/1d/lagrange +++ b/deal.II/deal.II/source/fe/scripts/1d/lagrange @@ -10,6 +10,7 @@ phi_polynom := array(0..n_functions-1); grad_phi_polynom := array(0..n_functions-1); + grad_grad_phi_polynom := array(0..n_functions-1); local_mass_matrix := array(0..n_functions-1, 0..n_functions-1); for i from 0 to n_functions-1 do @@ -28,6 +29,7 @@ phi_polynom[i] := interp (shifted_support_points, values, xi); grad_phi_polynom[i] := diff(phi_polynom[i], xi); + grad_grad_phi_polynom[i] := diff(grad_phi_polynom[i], xi); od; phi:= proc(i,x) subs(xi=x, phi_polynom[i]); end; diff --git a/deal.II/deal.II/source/fe/scripts/1d/lagrange-cubic b/deal.II/deal.II/source/fe/scripts/1d/lagrange-cubic index 662ff18e07..a33cbd26f2 100644 --- a/deal.II/deal.II/source/fe/scripts/1d/lagrange-cubic +++ b/deal.II/deal.II/source/fe/scripts/1d/lagrange-cubic @@ -15,6 +15,7 @@ readlib(C); C(phi_polynom, filename=cubic1d_shape_value); C(grad_phi_polynom, filename=cubic1d_shape_grad); + C(grad_grad_phi_polynom, filename=cubic1d_shape_grad_grad); C(prolongation, filename=cubic1d_prolongation); C(restriction, filename=cubic1d_restriction); C(local_mass_matrix, optimized, filename=cubic1d_massmatrix); diff --git a/deal.II/deal.II/source/fe/scripts/1d/lagrange-quartic b/deal.II/deal.II/source/fe/scripts/1d/lagrange-quartic index 7b878c03d4..bf82c87237 100644 --- a/deal.II/deal.II/source/fe/scripts/1d/lagrange-quartic +++ b/deal.II/deal.II/source/fe/scripts/1d/lagrange-quartic @@ -14,8 +14,9 @@ # write data to files readlib(C); - C(phi_polynom, filename=cubic1d_shape_value); - C(grad_phi_polynom, filename=cubic1d_shape_grad); - C(prolongation, filename=cubic1d_prolongation); - C(restriction, filename=cubic1d_restriction); - C(local_mass_matrix, optimized, filename=cubic1d_massmatrix); + C(phi_polynom, filename=quartic1d_shape_value); + C(grad_phi_polynom, filename=quartic1d_shape_grad); + C(grad_grad_phi_polynom, filename=quartic1d_shape_grad_grad); + C(prolongation, filename=quartic1d_prolongation); + C(restriction, filename=quartic1d_restriction); + C(local_mass_matrix, optimized, filename=quartic1d_massmatrix);