diff --git a/Sheet7/E14/jacob_template/square.m b/Sheet7/E14/jacob_template/square.m new file mode 100644 index 0000000..f793583 --- /dev/null +++ b/Sheet7/E14/jacob_template/square.m @@ -0,0 +1,42 @@ +% Square: +% flatpak run org.octave.Octave +% or +% octave --no-window-system --no-gui -qf + +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 +% +% 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,:) + diff --git a/Sheet7/E14/jacob_template/square_bb_4.txt b/Sheet7/E14/jacob_template/square_bb_4.txt new file mode 100644 index 0000000..b3136b8 --- /dev/null +++ b/Sheet7/E14/jacob_template/square_bb_4.txt @@ -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 diff --git a/Sheet7/E14/jacob_template/square_tiny.pdf b/Sheet7/E14/jacob_template/square_tiny.pdf new file mode 100644 index 0000000..079189b Binary files /dev/null and b/Sheet7/E14/jacob_template/square_tiny.pdf differ diff --git a/Sheet7/E14/jacob_template/square_tiny.txt b/Sheet7/E14/jacob_template/square_tiny.txt new file mode 100644 index 0000000..97c91b1 --- /dev/null +++ b/Sheet7/E14/jacob_template/square_tiny.txt @@ -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 diff --git a/Sheet7/E14/jacob_template/userset.cpp b/Sheet7/E14/jacob_template/userset.cpp new file mode 100644 index 0000000..e11dc39 --- /dev/null +++ b/Sheet7/E14/jacob_template/userset.cpp @@ -0,0 +1,16 @@ +#include "userset.h" +#include + + +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 ; +}