lisa.pizzo/LisaPizzoExercises
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

.

Browse source
This commit is contained in:
Lisa Pizzo 2026-01-10 15:28:10 +01:00
commit 3355c00bcf
5 changed files with 879 additions and 0 deletions
Download patch file Download diff file Expand all files Collapse all files

42
Sheet7/E14/jacob_template/square.m Normal file
View file

@ -0,0 +1,42 @@
% Square:
% flatpak run org.octave.Octave <filename>
% or
% octave --no-window-system --no-gui -qf <filename>
clear all
clc
% %% L-shape
% g=[2 0 2 0 0 1 0; % #vertices,v_1x, v_2x, v_1y, v_2y, subdomain_left, subdomain_right
% 2 2 2 0 1 1 0;
% 2 2 1 1 0.5 1 0;
% 2 1 1 0.5 2 1 0;
% 2 1 0 2 2 1 0;
% 2 0 0 2 0 1 0]';
%% square
g=[2 0 1 0 0 1 0; % #vertices,v_1x, v_2x, v_1y, v_2y, subdomain_left, subdomain_right
2 1 1 0 1 1 0;
2 1 0 1 1 1 0;
2 0 0 1 0 1 0]';
%[p,e,t] = initmesh(g,'hmax',0.01);
[p,e,t] = initmesh(g,'hmax',0.5);
pdemesh(p,e,t)
%% GH
% output from <https://de.mathworks.com/help/pde/ug/initmesh.html initmesh>
%
% coordinates p: [2][nnode]
% connectivity t: [4][nelem] with t(4,:) are the subdomain numbers
% edges e: [7][nedges] boundary edges
% e([1,2],:) - start/end vertex of edge
% e([3,4],:) - start/end values
% e(5,:) - segment number
% e([6,7],:) - left/right subdomain
ascii_write_mesh( p, t, e, mfilename);
% tmp=t(1:3,:)

726
Sheet7/E14/jacob_template/square_bb_4.txt Normal file
View file

@ -0,0 +1,726 @@
209
2
384
3
0 0
1 0
1 1
0 1
2 0
2 1
1 2
0 2
2 2
1.5 1.5
0.5 1.5
1.5 0.5
0.5 0.5
1 0.5
1.5 1
0.5 1
1 1.5
0.5 0
1 0.25
0.75 1
0 0.5
1.5 0
2 0.5
1.75 1
1 1.25
0.5 2
0 1.5
2 1.5
1.5 2
1 0.75
1.25 1
0.25 1
1 1.75
1.75 1.25
1.25 1.75
1.75 1.75
0.25 1.25
0.75 1.75
0.25 1.75
1.25 0.25
1.75 0.25
1.75 0.75
0.25 0.25
0.75 0.25
0.25 0.75
1.25 0.5
0.75 0.5
1.5 1.25
1.5 0.75
1.25 0.75
0.5 1.25
0.5 0.75
0.75 0.75
1.25 1.5
0.75 1.5
1.25 1.25
0.75 1.25
0.25 0
1 0.125
0.875 1
0 0.75
1.25 0
2 0.25
1.875 1
1 1.125
0.75 2
0 1.75
2 1.25
1.75 2
1 0.625
1.375 1
0.375 1
1 1.625
0.75 0
1 0.375
0.625 1
0 0.25
1.75 0
2 0.75
1.625 1
1 1.375
0.25 2
0 1.25
2 1.75
1.25 2
1 0.875
1.125 1
0.125 1
1 1.875
0.875 1.125
0.875 0.875
1.125 1.125
1.125 0.875
1.875 1.125
1.625 1.375
1.75 1.125
1.875 1.375
1.125 1.875
1.375 1.625
1.375 1.875
1.125 1.75
1.875 1.875
1.625 1.625
1.875 1.625
1.625 1.875
1.75 1.5
1.5 1.75
0.125 1.125
0.375 1.375
0.125 1.375
0.25 1.125
0.875 1.875
0.625 1.625
0.625 1.875
0.875 1.75
0.125 1.875
0.375 1.625
0.375 1.875
0.125 1.625
0.25 1.5
0.5 1.75
1.125 0.125
1.375 0.375
1.125 0.25
1.375 0.125
1.875 0.125
1.625 0.375
1.625 0.125
1.875 0.375
1.5 0.25
1.875 0.875
1.625 0.625
1.875 0.625
1.75 0.875
1.75 0.5
0.125 0.125
0.375 0.375
0.375 0.125
0.125 0.375
0.875 0.125
0.625 0.375
0.625 0.125
0.875 0.25
0.5 0.25
0.125 0.875
0.375 0.625
0.125 0.625
0.25 0.875
0.25 0.5
1.375 0.5
1.125 0.5
1.125 0.375
1.25 0.375
0.625 0.5
0.875 0.5
0.875 0.375
0.75 0.375
1.5 1.375
1.5 1.125
1.625 1.125
1.625 1.25
1.5 0.625
1.5 0.875
1.625 0.875
1.625 0.75
1.375 0.625
1.125 0.625
1.375 0.875
1.125 0.75
1.25 0.875
1.25 0.625
1.375 0.75
0.5 1.375
0.5 1.125
0.375 1.125
0.375 1.25
0.5 0.625
0.5 0.875
0.375 0.875
0.375 0.75
0.625 0.625
0.875 0.625
0.625 0.875
0.75 0.875
0.875 0.75
0.75 0.625
0.625 0.75
1.375 1.5
1.125 1.5
1.125 1.625
1.25 1.625
1.375 1.375
0.625 1.5
0.875 1.5
0.875 1.625
0.75 1.625
0.625 1.375
1.375 1.125
1.125 1.375
1.125 1.25
1.25 1.125
1.375 1.25
1.25 1.375
0.625 1.125
0.875 1.375
0.75 1.125
0.875 1.25
0.625 1.25
0.75 1.375
4 61 145
1 58 136
13 141 154
17 205 81
16 183 76
8 67 116
10 99 188
2 62 122
9 69 102
6 68 94
7 66 112
17 199 189
5 63 126
12 132 162
11 109 173
6 94 64
2 122 59
7 112 89
15 168 163
14 167 70
4 145 88
14 182 155
15 198 71
16 204 174
21 77 139
18 74 142
44 140 143
57 204 206
53 182 185
27 83 110
35 98 101
22 78 128
29 85 100
28 84 104
26 82 118
56 198 202
23 79 133
42 131 134
37 108 111
34 95 161
40 123 153
38 113 196
50 167 171
50 168 170
45 146 180
53 183 187
56 199 200
57 205 209
45 149 146
43 144 137
47 156 155
25 90 65
20 91 60
39 120 117
54 190 189
40 130 123
36 107 103
34 106 95
38 121 113
54 192 188
41 135 127
49 164 163
51 175 174
24 160 80
19 152 75
33 195 73
49 166 162
30 93 86
32 179 72
47 181 154
31 92 87
51 197 173
43 149 139
44 144 142
19 156 143
20 90 206
30 91 185
37 120 110
33 190 101
41 130 128
35 107 100
36 106 104
39 121 118
48 192 202
42 135 133
24 164 134
32 175 111
48 160 161
46 152 153
55 195 196
46 166 171
31 93 170
52 179 180
52 181 187
25 92 200
55 197 209
21 147 61
18 138 58
44 157 141
57 207 205
53 184 183
27 119 67
35 191 99
22 125 62
29 105 69
28 97 68
26 114 66
56 203 199
23 129 63
42 165 132
37 176 109
34 96 94
40 124 122
38 115 112
50 172 168
50 169 167
45 148 145
53 186 182
56 201 198
57 208 204
1 136 77
2 140 74
2 59 140
16 76 204
14 70 182
4 108 83
7 89 98
5 126 78
7 98 85
9 102 84
8 116 82
15 159 198
6 131 79
6 64 131
4 88 108
10 158 95
12 150 123
11 193 113
14 151 167
15 71 168
13 177 146
16 178 183
17 81 199
17 194 205
43 137 149
44 141 144
19 75 156
20 60 90
30 86 91
37 109 120
33 73 190
41 127 130
35 99 107
36 103 106
39 117 121
48 158 192
42 132 135
24 80 164
32 72 175
48 159 160
46 151 152
55 194 195
46 150 166
31 87 93
52 178 179
52 177 181
25 65 92
55 193 197
45 147 149
43 138 144
47 157 156
25 207 90
20 184 91
39 119 120
54 191 190
40 125 130
36 105 107
34 97 106
38 114 121
54 203 192
41 129 135
49 165 164
51 176 175
24 96 160
19 124 152
33 115 195
49 172 166
30 169 93
32 148 179
47 186 181
31 201 92
51 208 197
45 145 147
43 136 138
47 154 157
25 81 207
20 76 184
39 116 119
54 188 191
40 122 125
36 102 105
34 94 97
38 112 114
54 189 203
41 126 129
49 162 165
51 173 176
24 64 96
19 59 124
33 89 115
49 163 172
30 70 169
32 88 148
47 155 186
31 71 201
51 174 208
43 139 136
44 142 140
19 143 59
20 206 76
30 185 70
37 110 108
33 101 89
41 128 126
35 100 98
36 104 102
39 118 116
48 202 159
42 133 131
24 134 64
32 111 88
48 161 158
46 153 150
55 196 193
46 171 151
31 170 71
52 180 177
52 187 178
25 200 81
55 209 194
13 146 137
13 137 141
14 155 75
3 65 60
3 60 86
11 117 109
17 189 73
12 123 127
10 103 99
10 95 103
11 113 117
10 188 158
12 127 132
15 163 80
16 174 72
15 80 159
14 75 151
17 73 194
12 162 150
3 86 87
16 72 178
13 154 177
3 87 65
11 173 193
21 139 147
18 142 138
44 143 157
57 206 207
53 185 184
27 110 119
35 101 191
22 128 125
29 100 105
28 104 97
26 118 114
56 202 203
23 133 129
42 134 165
37 111 176
34 161 96
40 153 124
38 196 115
50 171 172
50 170 169
45 180 148
53 187 186
56 200 201
57 209 208
147 145 61
138 136 58
157 154 141
207 81 205
184 76 183
119 116 67
191 188 99
125 122 62
105 102 69
97 94 68
114 112 66
203 189 199
129 126 63
165 162 132
176 173 109
96 64 94
124 59 122
115 89 112
172 163 168
169 70 167
148 88 145
186 155 182
201 71 198
208 174 204
136 139 77
140 142 74
59 143 140
76 206 204
70 185 182
108 110 83
89 101 98
126 128 78
98 100 85
102 104 84
116 118 82
159 202 198
131 133 79
64 134 131
88 111 108
158 161 95
150 153 123
193 196 113
151 171 167
71 170 168
177 180 146
178 187 183
81 200 199
194 209 205
137 146 149
141 137 144
75 155 156
60 65 90
86 60 91
109 117 120
73 189 190
127 123 130
99 103 107
103 95 106
117 113 121
158 188 192
132 127 135
80 163 164
72 174 175
159 80 160
151 75 152
194 73 195
150 162 166
87 86 93
178 72 179
177 154 181
65 87 92
193 173 197
147 139 149
138 142 144
157 143 156
207 206 90
184 185 91
119 110 120
191 101 190
125 128 130
105 100 107
97 104 106
114 118 121
203 202 192
129 133 135
165 134 164
176 111 175
96 161 160
124 153 152
115 196 195
172 171 166
169 170 93
148 180 179
186 187 181
201 200 92
208 209 197
64
1 58
2 59
3 60
4 61
2 62
5 63
6 64
3 65
7 66
8 67
6 68
9 69
14 70
15 71
16 72
17 73
18 74
19 75
20 76
21 77
22 78
23 79
24 80
25 81
26 82
27 83
28 84
29 85
30 86
31 87
32 88
33 89
58 18
59 19
60 20
61 21
62 22
63 23
64 24
65 25
66 26
67 27
68 28
69 29
70 30
71 31
72 32
73 33
74 2
75 14
76 16
77 1
78 5
79 6
80 15
81 17
82 8
83 4
84 9
85 7
86 3
87 3
88 4
89 7
1 0
1 2
1 3
1 0
2 0
2 0
2 4
3 4
3 0
3 0
4 0
4 0
1 2
2 4
1 3
3 4
1 0
1 2
1 3
1 0
2 0
2 0
2 4
3 4
3 0
3 0
4 0
4 0
1 2
2 4
1 3
3 4
1 0
1 2
1 3
1 0
2 0
2 0
2 4
3 4
3 0
3 0
4 0
4 0
1 2
2 4
1 3
3 4
1 0
1 2
1 3
1 0
2 0
2 0
2 4
3 4
3 0
3 0
4 0
4 0
1 2
2 4
1 3
3 4

BIN
Sheet7/E14/jacob_template/square_tiny.pdf Normal file
View file

Binary file not shown.

95
Sheet7/E14/jacob_template/square_tiny.txt Normal file
View file

@ -0,0 +1,95 @@
13
2
16
3
0
0
1
0
1
1
0
1
0.5
0
1
0.5
0.5
1
0
0.5
0.4999999999999999
0.4999999999999999
0.3333333333333333
0.6666666666666666
0.6666666666666666
0.6666666666666666
0.6666666666666666
0.3333333333333333
0.3333333333333333
0.3333333333333333
8
1
13
5
2
12
6
3
11
7
4
10
1
5
13
10
8
13
2
6
12
3
7
11
4
8
10
12
9
13
10
9
11
7
10
11
11
9
12
6
11
12
9
10
13
5
12
13
8
1
5
5
2
2
6
6
3
3
7
7
4
4
8
8
1

16
Sheet7/E14/jacob_template/userset.cpp Normal file
View file

@ -0,0 +1,16 @@
#include "userset.h"
#include <cmath>
double FunctF(double const x , double const y)
{
// return std::sin(3.14159*1*x)*std::sin(3.14159*1*y);
// return 16.0*1024. ;
// return (double)1.0 ;
return x * x * std::sin(2.5 * 3.14159 * y);
}
double FunctU(const double /* x */, double const /* y */)
{
return 1.0 ;
}