]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Update.
authorwolf <wolf@0785d39b-7218-0410-832d-ea1e28bc413d>
Fri, 6 Nov 1998 15:43:40 +0000 (15:43 +0000)
committerwolf <wolf@0785d39b-7218-0410-832d-ea1e28bc413d>
Fri, 6 Nov 1998 15:43:40 +0000 (15:43 +0000)
git-svn-id: https://svn.dealii.org/trunk@647 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/deal.II/source/fe/scripts/1d/lagrange
deal.II/deal.II/source/fe/scripts/1d/lagrange-cubic
deal.II/deal.II/source/fe/scripts/1d/lagrange-quartic

index c9a255a7113a45d06c9f8b96f37ed672a1f9d408..b57a028a09575127f035925ac450547f6cf6de31 100644 (file)
@@ -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;
index 662ff18e07de565af8cf97f10556f18d01edb52c..a33cbd26f2e4f5c51bfb3af905c56981cb2df32e 100644 (file)
@@ -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);
index 7b878c03d4a0a4ea22db86217cd900dd9ac13081..bf82c87237d068a5a88501d2ff8afffa2a2ffdae 100644 (file)
@@ -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);

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.