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
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;
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);
# 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);