OptFEM3DP1 Toolbox  V1.0
Matlab/Octave Optimized P1-Lagrange Finite Element Method in 3D
 All Files Functions Variables Pages
valid_FEMmatrices.m
Go to the documentation of this file.
1 function Test=valid_FEMmatrices()
2 % automatic generation with sage
3 % Copyright:
4 % See \ref license
5 i=1;
6 Test(i).u=@(x,y,z) x + y + z;
7 Test(i).cu='x + y + z';
8 Test(i).v=@(x,y,z) x - y - z;
9 Test(i).cv='x - y - z';
10 Test(i).w=@(x,y,z) 2.*x + y - 3.*z;
11 Test(i).cw='2*x + y - 3*z';
12 Test(i).Mass=-0.8333333333333333;
13 Test(i).MassW=0.5000000000000000;
14 Test(i).Stiff=-1.000000000000000;
15 Test(i).degree=1;
16 i=2;
17 Test(i).u=@(x,y,z) 3.*x + 2.*y - z - 1;
18 Test(i).cu='3*x + 2*y - z - 1';
19 Test(i).v=@(x,y,z) 2.*x - 2.*y + 2.*z + 1;
20 Test(i).cv='2*x - 2*y + 2*z + 1';
21 Test(i).w=@(x,y,z) 2.*x - 5.*y - 3.*z;
22 Test(i).cw='2*x - 5*y - 3*z';
23 Test(i).Mass=2.000000000000000;
24 Test(i).MassW=-5.500000000000000;
25 Test(i).Stiff=0.0000000000000000;
26 Test(i).degree=1;
27 i=3;
28 Test(i).u=@(x,y,z) 3.*x.^2 - x.*y + 2.*y.^2 + y.*z - z.^2 - 3;
29 Test(i).cu='3*x^2 - x*y + 2*y^2 + y*z - z^2 - 3';
30 Test(i).v=@(x,y,z) 2.*x.^2 + x.*y - 3.*y.^2 - x.*z - y;
31 Test(i).cv='2*x^2 + x*y - 3*y^2 - x*z - y';
32 Test(i).w=@(x,y,z) x.*y + y.*z;
33 Test(i).cw='x*y + y*z';
34 Test(i).Mass=1.236111111111111;
35 Test(i).MassW=0.7134259259259259;
36 Test(i).Stiff=-1.916666666666667;
37 Test(i).degree=2;
38 i=4;
39 Test(i).u=@(x,y,z) 5.*x.^3 - z.^3 - x.*y + 2.*y.^2 + y.*z - 3;
40 Test(i).cu='5*x^3 - z^3 - x*y + 2*y^2 + y*z - 3';
41 Test(i).v=@(x,y,z) 2.*x.^2.*y + x.*y.*z - 3.*y.^2 - x.*z - y.*z;
42 Test(i).cv='2*x^2*y + x*y*z - 3*y^2 - x*z - y*z';
43 Test(i).w=@(x,y,z) x.*y.*z + x.*y - 2;
44 Test(i).cw='x*y*z + x*y - 2';
45 Test(i).Mass=1.220833333333333;
46 Test(i).MassW=-2.220205026455026;
47 Test(i).Stiff=-1.416666666666667;
48 Test(i).degree=3;
49 i=5;
50 Test(i).u=@(x,y,z) -x.*y.^2.*z + 2.*y.^2.*z.^2 - x.^3 - z.^3 + y.*z - 3;
51 Test(i).cu='-x*y^2*z + 2*y^2*z^2 - x^3 - z^3 + y*z - 3';
52 Test(i).v=@(x,y,z) x.^2.*y.*z + x.*y.*z.^2 - y.*z.^2 - 3.*y.^2 - x.*z;
53 Test(i).cv='x^2*y*z + x*y*z^2 - y*z^2 - 3*y^2 - x*z';
54 Test(i).w=@(x,y,z) x.*y.*z.^2 + x.*y.*z - 2;
55 Test(i).cw='x*y*z^2 + x*y*z - 2';
56 Test(i).Mass=3.646990740740741;
57 Test(i).MassW=-6.175414021164021;
58 Test(i).Stiff=-3.338888888888889;
59 Test(i).degree=4;
60 i=6;
61 Test(i).u=@(x,y,z) 4.*x.^5 - x.*y.^2.*z + 2.*y.^2.*z.^2 - z.^3 + y.*z - 3;
62 Test(i).cu='4*x^5 - x*y^2*z + 2*y^2*z^2 - z^3 + y*z - 3';
63 Test(i).v=@(x,y,z) x.^2.*y.*z.^2 + x.^2.*y.*z - 2.*x.*y.*z.^2 + y.*z.^2 - 3.*y.^2 - x.*z;
64 Test(i).cv='x^2*y*z^2 + x^2*y*z - 2*x*y*z^2 + y*z^2 - 3*y^2 - x*z';
65 Test(i).w=@(x,y,z) -2.*x.^2.*y.^3 + x.*y.*z.^2 + x.*y.*z;
66 Test(i).cw='-2*x^2*y^3 + x*y*z^2 + x*y*z';
67 Test(i).Mass=2.157953042328042;
68 Test(i).MassW=0.07763668430335097;
69 Test(i).Stiff=-3.102777777777778;
70 Test(i).degree=5;
71 i=7;
72 Test(i).u=@(x,y,z) cos(x - 2./3.*y + 4./3.*z - 2./3);
73 Test(i).cu='cos(x - 2/3*y + 4/3*z - 2/3)';
74 Test(i).v=@(x,y,z) x.*y + x.*z - z.^2;
75 Test(i).cv='x*y + x*z - z^2';
76 Test(i).w=@(x,y,z) cos(1./2.*x - 3./4.*y + 1./4.*z);
77 Test(i).cw='cos(1/2*x - 3/4*y + 1/4*z)';
78 Test(i).Mass=0.1437236424227304;
79 Test(i).MassW=0.1407206851537448;
80 Test(i).Stiff=0.1639024563527328;
81 Test(i).degree=-1;