declareEspace.edp 2.08 KB
Newer Older
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126



int ndofs=0;

mesh Th0=readmesh("th0.msh");
fespace vh0(Th0,P1);
vh0 XX0=x;
vh0 YY0=y;
//SaveVec(XX0[],"X");
//SaveVec(YY0[],"Y");
ndofs+=vh0.ndof;

mesh Th1=readmesh("th1.msh");
fespace vh1(Th1,P1);
vh1 XX1=x;
vh1 YY1=y;
//SaveVec(XX1[],"X");
//SaveVec(YY1[],"Y");
ndofs+=vh1.ndof;


mesh Th2=readmesh("th2.msh");
fespace vh2(Th2,P1);
vh2 XX2=x;
vh2 YY2=y;
//SaveVec(XX2[],"X");
//SaveVec(YY2[],"Y");
ndofs+=vh2.ndof;

mesh Th3=readmesh("th3.msh");
fespace vh3(Th3,P1);
vh3 XX3=x;
vh3 YY3=y;
//SaveVec(XX3[],"X");
//SaveVec(YY3[],"Y");
ndofs+=vh3.ndof;

mesh Th4=readmesh("th4.msh");
fespace vh4(Th4,P1);
vh4 XX4=x;
vh4 YY4=y;
//SaveVec(XX4[],"X");
//SaveVec(YY4[],"Y");
ndofs+=vh4.ndof;

mesh Th5=readmesh("th5.msh");
fespace vh5(Th5,P1);
vh5 XX5=x;
vh5 YY5=y;
//SaveVec(XX5[],"X");
//SaveVec(YY5[],"Y");
ndofs+=vh5.ndof;

mesh Th6=readmesh("th6.msh");
fespace vh6(Th6,P1);
vh6 XX6=x;
vh6 YY6=y;
//SaveVec(XX6[],"X");
//SaveVec(YY6[],"Y");
ndofs+=vh6.ndof;

mesh Th7=readmesh("th7.msh");
fespace vh7(Th7,P1);
vh7 XX7=x;
vh7 YY7=y;
//SaveVec(XX7[],"X");
//SaveVec(YY7[],"Y");
ndofs+=vh7.ndof;

mesh Th8=readmesh("th8.msh");
fespace vh8(Th8,P1);
vh8 XX8=x;
vh8 YY8=y;
//SaveVec(XX8[],"X");
//SaveVec(YY8[],"Y");
ndofs+=vh8.ndof;

mesh Th9=readmesh("th9.msh");
fespace vh9(Th9,P1);
vh9 XX9=x;
vh9 YY9=y;
//SaveVec(XX9[],"X");
//SaveVec(YY9[],"Y");
ndofs+=vh9.ndof;

mesh Th10=readmesh("th10.msh");
fespace vh10(Th10,P1);
vh10 XX10=x;
vh10 YY10=y;
//SaveVec(XX10[],"X");
//SaveVec(YY10[],"Y");
ndofs+=vh10.ndof;

mesh Th11=readmesh("th11.msh");
fespace vh11(Th11,P1);
vh11 XX11=x;
vh11 YY11=y;
//SaveVec(XX11[],"X");
//SaveVec(YY11[],"Y");
ndofs+=vh11.ndof;

mesh Th12=readmesh("th12.msh");
fespace vh12(Th12,P1);
vh12 XX12=x;
vh12 YY12=y;
//SaveVec(XX12[],"X");
//SaveVec(YY12[],"Y");
ndofs+=vh12.ndof;

mesh Th13=readmesh("th13.msh");
fespace vh13(Th13,P1);
vh13 XX13=x;
vh13 YY13=y;
//SaveVec(XX13[],"X");
//SaveVec(YY13[],"Y");
ndofs+=vh13.ndof;

mesh Th14=readmesh("th14.msh");
fespace vh14(Th14,P1);
vh14 XX14=x;
vh14 YY14=y;
//SaveVec(XX14[],"X");
//SaveVec(YY14[],"Y");
ndofs+=vh14.ndof;
cout<<"ndofs total="<<ndofs<<endl;;