OptFEM2DP1 Toolbox  V1.2b3
Matlab/Octave Optimized P1-Lagrange Finite Element Method in 2D
 All Files Functions Pages
valid_StiffElas.m
Go to the documentation of this file.
1 function Test=valid_StiffElas()
2 % Copyright:
3 % See \ref license
4 i=1;
5 Test(i).u={@(x,y) x - 2.*y, @(x,y) x + y};
6 Test(i).cu={'x - 2*y','x + y'};
7 Test(i).v={@(x,y) x + 2.*y, @(x,y) 2.*x - y};
8 Test(i).cv={'x + 2*y','2*x - y'};
9 Test(i).Mass=-5/12;
10 Test(i).Stiff=-8;
11 Test(i).lambda=1;
12 Test(i).mu=2;
13 Test(i).degree=2;
14 i=2;
15 Test(i).u={@(x,y) x.^2 - 2.*x.*y, @(x,y) y.^2 + x};
16 Test(i).cu={'x^2 - 2*x*y','y^2 + x'};
17 Test(i).v={@(x,y) x.^2 + 2.*y.^2, @(x,y) 2.*x.^2 - x.*y};
18 Test(i).cv={'x^2 + 2*y^2','2*x^2 - x*y'};
19 Test(i).Mass=37/360;
20 Test(i).Stiff=-4/3;
21 Test(i).lambda=1;
22 Test(i).mu=2;
23 Test(i).degree=2;
24 i=3;
25 Test(i).u={@(x,y) x.^2 - y.^2, @(x,y) x.^3 - y.^3};
26 Test(i).cu={'x^2 - y^2','x^3 - y^3'};
27 Test(i).v={@(x,y) 2.*y.^3 + x, @(x,y) 2.*x.^2 - x.*y};
28 Test(i).cv={'2*y^3 + x','2*x^2 - x*y'};
29 Test(i).Mass=1/12;
30 Test(i).Stiff=3;
31 Test(i).lambda=1;
32 Test(i).mu=1;
33 Test(i).degree=2;
34 i=4;
35 Test(i).u={@(x,y) x.^2 - y.^4, @(x,y) x.^4 - y.^3};
36 Test(i).cu={'x^2 - y^4','x^4 - y^3'};
37 Test(i).v={@(x,y) x.^2.*y.^2 + 2.*y.^3, @(x,y) 2.*x.^2 - x.*y};
38 Test(i).cv={'x^2*y^2 + 2*y^3','2*x^2 - x*y'};
39 Test(i).Mass=1/14;
40 Test(i).Stiff=5/9;
41 Test(i).lambda=3;
42 Test(i).mu=1;
43 Test(i).degree=2;
44 i=5;
45 Test(i).u={@(x,y) x.^3.*y.^2 - x.^2 - y.^4, @(x,y) x.^4 - 2.*y.^5 - y.^3};
46 Test(i).cu={'x^3*y^2 - x^2 - y^4','x^4 - 2*y^5 - y^3'};
47 Test(i).v={@(x,y) x.^2.*y.^3 + 2.*y.^3, @(x,y) -x.^4.*y + 2.*x.^2};
48 Test(i).cv={'x^2*y^3 + 2*y^3','-x^4*y + 2*x^2'};
49 Test(i).Mass=-1927/4200;
50 Test(i).Stiff=101/35;
51 Test(i).lambda=1;
52 Test(i).mu=3;
53 Test(i).degree=2;
54 i=6;
55 Test(i).u={@(x,y) x.^6 + x.^3.*y.^2 - x.^2 - y.^4, @(x,y) -2.*x.*y.^5 + x.^5 - y.^3};
56 Test(i).cu={'x^6 + x^3*y^2 - x^2 - y^4','-2*x*y^5 + x^5 - y^3'};
57 Test(i).v={@(x,y) x.^3.*y.^3 + 2.*y.^3, @(x,y) -x.^4.*y.^2 + 2.*x.^2};
58 Test(i).cv={'x^3*y^3 + 2*y^3','-x^4*y^2 + 2*x^2'};
59 Test(i).Mass=-367/1120;
60 Test(i).Stiff=1073/504;
61 Test(i).lambda=1;
62 Test(i).mu=2;
63 Test(i).degree=2;