diff --git a/test/blockstructured/poisson/CMakeLists.txt b/test/blockstructured/poisson/CMakeLists.txt index 7ce13eff78142fb61a2ea43c5c7cb45c388a0bc6..1b48ca5f343efeffe528a5066d3bc4073898b8a5 100644 --- a/test/blockstructured/poisson/CMakeLists.txt +++ b/test/blockstructured/poisson/CMakeLists.txt @@ -1,7 +1,7 @@ add_subdirectory(3d) -dune_symlink_to_source_files(FILES square.msh) +dune_symlink_to_source_files(FILES square.msh cube.msh) dune_add_formcompiler_system_test(UFLFILE poisson.ufl BASENAME blockstructured_poisson diff --git a/test/blockstructured/poisson/cube.msh b/test/blockstructured/poisson/cube.msh new file mode 100644 index 0000000000000000000000000000000000000000..0eac3436f979911323a1307eb10e1a0686f5b4a1 --- /dev/null +++ b/test/blockstructured/poisson/cube.msh @@ -0,0 +1,615 @@ +$MeshFormat +2.0 0 8 +$EndMeshFormat +$Nodes +366 +1 0 0.5 0 +2 0.244661 0.502233 0 +3 0.26544 0.715177 0 +4 0.122821 0.650966 0 +5 0.536558 0.28602 0 +6 0.512143 0.143732 0 +7 0.618254 0.14114 0 +8 0.681191 0.251832 0 +9 0.713611 0.720833 0 +10 0.669166 0.870357 0 +11 0.522464 0.854908 0 +12 0.540201 0.701589 0 +13 0.153432 0.365595 0 +14 0.261079 0.274632 0 +15 0.36626 0.375488 0 +16 0 0.333333 0 +17 0.5 0 0 +18 0.666667 0 0 +19 0.823879 0.54292 0 +20 1 0.5 0 +21 1 0.666667 0 +22 0.863154 0.676934 0 +23 0.5 1 0 +24 0.333333 1 0 +25 0.376237 0.855576 0 +26 0.365396 0.150972 0 +27 0.333333 0 0 +28 0.761616 0.400009 0 +29 1 0.333333 0 +30 0.666667 1 0 +31 0.546663 0.49302 0 +32 0.394587 0.70129 0 +33 0.393225 0.576529 0 +34 0.885313 0.770579 0 +35 0.821844 0.824177 0 +36 0.114959 0.238751 0 +37 0.177635 0.180744 0 +38 0.767996 0.888065 0 +39 0.119952 0.807903 0 +40 0.0837918 0.7004 0 +41 0.236315 0.116804 0 +42 0.702554 0.587171 0 +43 0.427914 0.282699 0 +44 0.273484 0.867522 0 +45 0.166667 1 0 +46 0.115282 0.116456 0 +47 0.166667 0 0 +48 0.885092 0.885433 0 +49 0.833333 1 0 +50 0 0.166667 0 +51 1 0.833333 0 +52 0 0.833333 0 +53 0 0.666667 0 +54 0.881426 0.118655 0 +55 1 0.166667 0 +56 0.806242 0.176795 0 +57 0.833333 0 0 +58 0 0 0 +59 1 1 0 +60 1 0 0 +61 0 1 0 +62 0.833333 1 0.2 +63 1 1 0.2 +64 0.833333 1 0.4 +65 1 1 0.4 +66 0.833333 1 0.6 +67 1 1 0.6 +68 0.833333 1 0.8 +69 1 1 0.8 +70 0.833333 1 1 +71 1 1 1 +72 0.666667 1 0.2 +73 0.666667 1 0.4 +74 0.666667 1 0.6 +75 0.666667 1 0.8 +76 0.666667 1 1 +77 0.5 1 0.2 +78 0.5 1 0.4 +79 0.5 1 0.6 +80 0.5 1 0.8 +81 0.5 1 1 +82 0.333333 1 0.2 +83 0.333333 1 0.4 +84 0.333333 1 0.6 +85 0.333333 1 0.8 +86 0.333333 1 1 +87 0.166667 1 0.2 +88 0.166667 1 0.4 +89 0.166667 1 0.6 +90 0.166667 1 0.8 +91 0.166667 1 1 +92 0 1 0.2 +93 0 1 0.4 +94 0 1 0.6 +95 0 1 0.8 +96 0 1 1 +97 0 0.833333 0.2 +98 0 0.833333 0.4 +99 0 0.833333 0.6 +100 0 0.833333 0.8 +101 0 0.833333 1 +102 0 0.666667 0.2 +103 0 0.666667 0.4 +104 0 0.666667 0.6 +105 0 0.666667 0.8 +106 0 0.666667 1 +107 0 0.5 0.2 +108 0 0.5 0.4 +109 0 0.5 0.6 +110 0 0.5 0.8 +111 0 0.5 1 +112 0 0.333333 0.2 +113 0 0.333333 0.4 +114 0 0.333333 0.6 +115 0 0.333333 0.8 +116 0 0.333333 1 +117 0 0.166667 0.2 +118 0 0.166667 0.4 +119 0 0.166667 0.6 +120 0 0.166667 0.8 +121 0 0.166667 1 +122 0 0 0.2 +123 0 0 0.4 +124 0 0 0.6 +125 0 0 0.8 +126 0 0 1 +127 0.166667 0 0.2 +128 0.166667 0 0.4 +129 0.166667 0 0.6 +130 0.166667 0 0.8 +131 0.166667 0 1 +132 0.333333 0 0.2 +133 0.333333 0 0.4 +134 0.333333 0 0.6 +135 0.333333 0 0.8 +136 0.333333 0 1 +137 0.5 0 0.2 +138 0.5 0 0.4 +139 0.5 0 0.6 +140 0.5 0 0.8 +141 0.5 0 1 +142 0.666667 0 0.2 +143 0.666667 0 0.4 +144 0.666667 0 0.6 +145 0.666667 0 0.8 +146 0.666667 0 1 +147 0.833333 0 0.2 +148 0.833333 0 0.4 +149 0.833333 0 0.6 +150 0.833333 0 0.8 +151 0.833333 0 1 +152 1 0 0.2 +153 1 0 0.4 +154 1 0 0.6 +155 1 0 0.8 +156 1 0 1 +157 1 0.166667 0.2 +158 1 0.166667 0.4 +159 1 0.166667 0.6 +160 1 0.166667 0.8 +161 1 0.166667 1 +162 1 0.333333 0.2 +163 1 0.333333 0.4 +164 1 0.333333 0.6 +165 1 0.333333 0.8 +166 1 0.333333 1 +167 1 0.5 0.2 +168 1 0.5 0.4 +169 1 0.5 0.6 +170 1 0.5 0.8 +171 1 0.5 1 +172 1 0.666667 0.2 +173 1 0.666667 0.4 +174 1 0.666667 0.6 +175 1 0.666667 0.8 +176 1 0.666667 1 +177 1 0.833333 0.2 +178 1 0.833333 0.4 +179 1 0.833333 0.6 +180 1 0.833333 0.8 +181 1 0.833333 1 +182 0.244661 0.502233 1 +183 0.26544 0.715177 1 +184 0.122821 0.650966 1 +185 0.536558 0.28602 1 +186 0.512143 0.143732 1 +187 0.618254 0.14114 1 +188 0.681191 0.251832 1 +189 0.713611 0.720833 1 +190 0.669166 0.870357 1 +191 0.522464 0.854908 1 +192 0.540201 0.701589 1 +193 0.153432 0.365595 1 +194 0.261079 0.274632 1 +195 0.36626 0.375488 1 +196 0.823879 0.54292 1 +197 0.863154 0.676934 1 +198 0.376237 0.855576 1 +199 0.365396 0.150972 1 +200 0.761616 0.400009 1 +201 0.546663 0.49302 1 +202 0.394587 0.70129 1 +203 0.393225 0.576529 1 +204 0.885313 0.770579 1 +205 0.821844 0.824177 1 +206 0.114959 0.238751 1 +207 0.177635 0.180744 1 +208 0.767996 0.888065 1 +209 0.119952 0.807903 1 +210 0.0837918 0.7004 1 +211 0.236315 0.116804 1 +212 0.702554 0.587171 1 +213 0.427914 0.282699 1 +214 0.273484 0.867522 1 +215 0.115282 0.116456 1 +216 0.885092 0.885433 1 +217 0.881426 0.118655 1 +218 0.806242 0.176795 1 +219 0.26544 0.715177 0.2 +220 0.122821 0.650966 0.2 +221 0.244661 0.502233 0.2 +222 0.26544 0.715177 0.4 +223 0.122821 0.650966 0.4 +224 0.244661 0.502233 0.4 +225 0.26544 0.715177 0.6 +226 0.122821 0.650966 0.6 +227 0.244661 0.502233 0.6 +228 0.26544 0.715177 0.8 +229 0.122821 0.650966 0.8 +230 0.244661 0.502233 0.8 +231 0.618254 0.14114 0.2 +232 0.681191 0.251832 0.2 +233 0.536558 0.28602 0.2 +234 0.512143 0.143732 0.2 +235 0.618254 0.14114 0.4 +236 0.681191 0.251832 0.4 +237 0.536558 0.28602 0.4 +238 0.512143 0.143732 0.4 +239 0.618254 0.14114 0.6 +240 0.681191 0.251832 0.6 +241 0.536558 0.28602 0.6 +242 0.512143 0.143732 0.6 +243 0.618254 0.14114 0.8 +244 0.681191 0.251832 0.8 +245 0.536558 0.28602 0.8 +246 0.512143 0.143732 0.8 +247 0.522464 0.854908 0.2 +248 0.540201 0.701589 0.2 +249 0.713611 0.720833 0.2 +250 0.669166 0.870357 0.2 +251 0.522464 0.854908 0.4 +252 0.540201 0.701589 0.4 +253 0.713611 0.720833 0.4 +254 0.669166 0.870357 0.4 +255 0.522464 0.854908 0.6 +256 0.540201 0.701589 0.6 +257 0.713611 0.720833 0.6 +258 0.669166 0.870357 0.6 +259 0.522464 0.854908 0.8 +260 0.540201 0.701589 0.8 +261 0.713611 0.720833 0.8 +262 0.669166 0.870357 0.8 +263 0.261079 0.274632 0.2 +264 0.36626 0.375488 0.2 +265 0.153432 0.365595 0.2 +266 0.261079 0.274632 0.4 +267 0.36626 0.375488 0.4 +268 0.153432 0.365595 0.4 +269 0.261079 0.274632 0.6 +270 0.36626 0.375488 0.6 +271 0.153432 0.365595 0.6 +272 0.261079 0.274632 0.8 +273 0.36626 0.375488 0.8 +274 0.153432 0.365595 0.8 +275 0.863154 0.676934 0.2 +276 0.823879 0.54292 0.2 +277 0.863154 0.676934 0.4 +278 0.823879 0.54292 0.4 +279 0.863154 0.676934 0.6 +280 0.823879 0.54292 0.6 +281 0.863154 0.676934 0.8 +282 0.823879 0.54292 0.8 +283 0.376237 0.855576 0.2 +284 0.376237 0.855576 0.4 +285 0.376237 0.855576 0.6 +286 0.376237 0.855576 0.8 +287 0.365396 0.150972 0.2 +288 0.365396 0.150972 0.4 +289 0.365396 0.150972 0.6 +290 0.365396 0.150972 0.8 +291 0.761616 0.400009 0.2 +292 0.761616 0.400009 0.4 +293 0.761616 0.400009 0.6 +294 0.761616 0.400009 0.8 +295 0.394587 0.70129 0.2 +296 0.393225 0.576529 0.2 +297 0.546663 0.49302 0.2 +298 0.394587 0.70129 0.4 +299 0.393225 0.576529 0.4 +300 0.546663 0.49302 0.4 +301 0.394587 0.70129 0.6 +302 0.393225 0.576529 0.6 +303 0.546663 0.49302 0.6 +304 0.394587 0.70129 0.8 +305 0.393225 0.576529 0.8 +306 0.546663 0.49302 0.8 +307 0.821844 0.824177 0.2 +308 0.885313 0.770579 0.2 +309 0.821844 0.824177 0.4 +310 0.885313 0.770579 0.4 +311 0.821844 0.824177 0.6 +312 0.885313 0.770579 0.6 +313 0.821844 0.824177 0.8 +314 0.885313 0.770579 0.8 +315 0.177635 0.180744 0.2 +316 0.114959 0.238751 0.2 +317 0.177635 0.180744 0.4 +318 0.114959 0.238751 0.4 +319 0.177635 0.180744 0.6 +320 0.114959 0.238751 0.6 +321 0.177635 0.180744 0.8 +322 0.114959 0.238751 0.8 +323 0.767996 0.888065 0.2 +324 0.767996 0.888065 0.4 +325 0.767996 0.888065 0.6 +326 0.767996 0.888065 0.8 +327 0.119952 0.807903 0.2 +328 0.0837918 0.7004 0.2 +329 0.119952 0.807903 0.4 +330 0.0837918 0.7004 0.4 +331 0.119952 0.807903 0.6 +332 0.0837918 0.7004 0.6 +333 0.119952 0.807903 0.8 +334 0.0837918 0.7004 0.8 +335 0.236315 0.116804 0.2 +336 0.236315 0.116804 0.4 +337 0.236315 0.116804 0.6 +338 0.236315 0.116804 0.8 +339 0.702554 0.587171 0.2 +340 0.702554 0.587171 0.4 +341 0.702554 0.587171 0.6 +342 0.702554 0.587171 0.8 +343 0.427914 0.282699 0.2 +344 0.427914 0.282699 0.4 +345 0.427914 0.282699 0.6 +346 0.427914 0.282699 0.8 +347 0.273484 0.867522 0.2 +348 0.273484 0.867522 0.4 +349 0.273484 0.867522 0.6 +350 0.273484 0.867522 0.8 +351 0.115282 0.116456 0.2 +352 0.115282 0.116456 0.4 +353 0.115282 0.116456 0.6 +354 0.115282 0.116456 0.8 +355 0.885092 0.885433 0.2 +356 0.885092 0.885433 0.4 +357 0.885092 0.885433 0.6 +358 0.885092 0.885433 0.8 +359 0.806242 0.176795 0.2 +360 0.881426 0.118655 0.2 +361 0.806242 0.176795 0.4 +362 0.881426 0.118655 0.4 +363 0.806242 0.176795 0.6 +364 0.881426 0.118655 0.6 +365 0.806242 0.176795 0.8 +366 0.881426 0.118655 0.8 +$EndNodes +$Elements +240 +1 5 0 1 2 3 4 107 221 219 220 +2 5 0 107 221 219 220 108 224 222 223 +3 5 0 108 224 222 223 109 227 225 226 +4 5 0 109 227 225 226 110 230 228 229 +5 5 0 111 184 183 182 110 229 228 230 +6 5 0 5 6 7 8 233 234 231 232 +7 5 0 233 234 231 232 237 238 235 236 +8 5 0 237 238 235 236 241 242 239 240 +9 5 0 241 242 239 240 245 246 243 244 +10 5 0 185 188 187 186 245 244 243 246 +11 5 0 9 10 11 12 249 250 247 248 +12 5 0 249 250 247 248 253 254 251 252 +13 5 0 253 254 251 252 257 258 255 256 +14 5 0 257 258 255 256 261 262 259 260 +15 5 0 189 192 191 190 261 260 259 262 +16 5 0 2 13 14 15 221 265 263 264 +17 5 0 221 265 263 264 224 268 266 267 +18 5 0 224 268 266 267 227 271 269 270 +19 5 0 227 271 269 270 230 274 272 273 +20 5 0 182 195 194 193 230 273 272 274 +21 5 0 1 16 13 2 107 112 265 221 +22 5 0 107 112 265 221 108 113 268 224 +23 5 0 108 113 268 224 109 114 271 227 +24 5 0 109 114 271 227 110 115 274 230 +25 5 0 111 182 193 116 110 230 274 115 +26 5 0 17 18 7 6 137 142 231 234 +27 5 0 137 142 231 234 138 143 235 238 +28 5 0 138 143 235 238 139 144 239 242 +29 5 0 139 144 239 242 140 145 243 246 +30 5 0 141 186 187 146 140 246 243 145 +31 5 0 20 21 22 19 167 172 275 276 +32 5 0 167 172 275 276 168 173 277 278 +33 5 0 168 173 277 278 169 174 279 280 +34 5 0 169 174 279 280 170 175 281 282 +35 5 0 171 196 197 176 170 282 281 175 +36 5 0 24 25 11 23 82 283 247 77 +37 5 0 82 283 247 77 83 284 251 78 +38 5 0 83 284 251 78 84 285 255 79 +39 5 0 84 285 255 79 85 286 259 80 +40 5 0 86 81 191 198 85 80 259 286 +41 5 0 17 6 26 27 137 234 287 132 +42 5 0 137 234 287 132 138 238 288 133 +43 5 0 138 238 288 133 139 242 289 134 +44 5 0 139 242 289 134 140 246 290 135 +45 5 0 141 136 199 186 140 135 290 246 +46 5 0 20 19 28 29 167 276 291 162 +47 5 0 167 276 291 162 168 278 292 163 +48 5 0 168 278 292 163 169 280 293 164 +49 5 0 169 280 293 164 170 282 294 165 +50 5 0 171 166 200 196 170 165 294 282 +51 5 0 30 23 11 10 72 77 247 250 +52 5 0 72 77 247 250 73 78 251 254 +53 5 0 73 78 251 254 74 79 255 258 +54 5 0 74 79 255 258 75 80 259 262 +55 5 0 76 190 191 81 75 262 259 80 +56 5 0 33 31 12 32 296 297 248 295 +57 5 0 296 297 248 295 299 300 252 298 +58 5 0 299 300 252 298 302 303 256 301 +59 5 0 302 303 256 301 305 306 260 304 +60 5 0 203 202 192 201 305 304 260 306 +61 5 0 9 22 34 35 249 275 308 307 +62 5 0 249 275 308 307 253 277 310 309 +63 5 0 253 277 310 309 257 279 312 311 +64 5 0 257 279 312 311 261 281 314 313 +65 5 0 189 205 204 197 261 313 314 281 +66 5 0 13 36 37 14 265 316 315 263 +67 5 0 265 316 315 263 268 318 317 266 +68 5 0 268 318 317 266 271 320 319 269 +69 5 0 271 320 319 269 274 322 321 272 +70 5 0 193 194 207 206 274 272 321 322 +71 5 0 9 35 38 10 249 307 323 250 +72 5 0 249 307 323 250 253 309 324 254 +73 5 0 253 309 324 254 257 311 325 258 +74 5 0 257 311 325 258 261 313 326 262 +75 5 0 189 190 208 205 261 262 326 313 +76 5 0 40 4 3 39 328 220 219 327 +77 5 0 328 220 219 327 330 223 222 329 +78 5 0 330 223 222 329 332 226 225 331 +79 5 0 332 226 225 331 334 229 228 333 +80 5 0 210 209 183 184 334 333 228 229 +81 5 0 14 37 41 26 263 315 335 287 +82 5 0 263 315 335 287 266 317 336 288 +83 5 0 266 317 336 288 269 319 337 289 +84 5 0 269 319 337 289 272 321 338 290 +85 5 0 194 199 211 207 272 290 338 321 +86 5 0 42 19 22 9 339 276 275 249 +87 5 0 339 276 275 249 340 278 277 253 +88 5 0 340 278 277 253 341 280 279 257 +89 5 0 341 280 279 257 342 282 281 261 +90 5 0 212 189 197 196 342 261 281 282 +91 5 0 33 2 15 31 296 221 264 297 +92 5 0 296 221 264 297 299 224 267 300 +93 5 0 299 224 267 300 302 227 270 303 +94 5 0 302 227 270 303 305 230 273 306 +95 5 0 203 201 195 182 305 306 273 230 +96 5 0 31 15 43 5 297 264 343 233 +97 5 0 297 264 343 233 300 267 344 237 +98 5 0 300 267 344 237 303 270 345 241 +99 5 0 303 270 345 241 306 273 346 245 +100 5 0 201 185 213 195 306 245 346 273 +101 5 0 42 31 28 19 339 297 291 276 +102 5 0 339 297 291 276 340 300 292 278 +103 5 0 340 300 292 278 341 303 293 280 +104 5 0 341 303 293 280 342 306 294 282 +105 5 0 212 196 200 201 342 282 294 306 +106 5 0 39 3 44 45 327 219 347 87 +107 5 0 327 219 347 87 329 222 348 88 +108 5 0 329 222 348 88 331 225 349 89 +109 5 0 331 225 349 89 333 228 350 90 +110 5 0 209 91 214 183 333 90 350 228 +111 5 0 46 47 41 37 351 127 335 315 +112 5 0 351 127 335 315 352 128 336 317 +113 5 0 352 128 336 317 353 129 337 319 +114 5 0 353 129 337 319 354 130 338 321 +115 5 0 215 207 211 131 354 321 338 130 +116 5 0 48 49 38 35 355 62 323 307 +117 5 0 355 62 323 307 356 64 324 309 +118 5 0 356 64 324 309 357 66 325 311 +119 5 0 357 66 325 311 358 68 326 313 +120 5 0 216 205 208 70 358 313 326 68 +121 5 0 50 46 37 36 117 351 315 316 +122 5 0 117 351 315 316 118 352 317 318 +123 5 0 118 352 317 318 119 353 319 320 +124 5 0 119 353 319 320 120 354 321 322 +125 5 0 121 206 207 215 120 322 321 354 +126 5 0 48 35 34 51 355 307 308 177 +127 5 0 355 307 308 177 356 309 310 178 +128 5 0 356 309 310 178 357 311 312 179 +129 5 0 357 311 312 179 358 313 314 180 +130 5 0 216 181 204 205 358 180 314 313 +131 5 0 32 25 44 3 295 283 347 219 +132 5 0 295 283 347 219 298 284 348 222 +133 5 0 298 284 348 222 301 285 349 225 +134 5 0 301 285 349 225 304 286 350 228 +135 5 0 202 183 214 198 304 228 350 286 +136 5 0 32 12 11 25 295 248 247 283 +137 5 0 295 248 247 283 298 252 251 284 +138 5 0 298 252 251 284 301 256 255 285 +139 5 0 301 256 255 285 304 260 259 286 +140 5 0 202 198 191 192 304 286 259 260 +141 5 0 31 5 8 28 297 233 232 291 +142 5 0 297 233 232 291 300 237 236 292 +143 5 0 300 237 236 292 303 241 240 293 +144 5 0 303 241 240 293 306 245 244 294 +145 5 0 201 200 188 185 306 294 244 245 +146 5 0 53 40 39 52 102 328 327 97 +147 5 0 102 328 327 97 103 330 329 98 +148 5 0 103 330 329 98 104 332 331 99 +149 5 0 104 332 331 99 105 334 333 100 +150 5 0 106 101 209 210 105 100 333 334 +151 5 0 55 29 56 54 157 162 359 360 +152 5 0 157 162 359 360 158 163 361 362 +153 5 0 158 163 361 362 159 164 363 364 +154 5 0 159 164 363 364 160 165 365 366 +155 5 0 161 217 218 166 160 366 365 165 +156 5 0 54 56 18 57 360 359 142 147 +157 5 0 360 359 142 147 362 361 143 148 +158 5 0 362 361 143 148 364 363 144 149 +159 5 0 364 363 144 149 366 365 145 150 +160 5 0 217 151 146 218 366 150 145 365 +161 5 0 50 58 47 46 117 122 127 351 +162 5 0 117 122 127 351 118 123 128 352 +163 5 0 118 123 128 352 119 124 129 353 +164 5 0 119 124 129 353 120 125 130 354 +165 5 0 121 215 131 126 120 354 130 125 +166 5 0 48 51 59 49 355 177 63 62 +167 5 0 355 177 63 62 356 178 65 64 +168 5 0 356 178 65 64 357 179 67 66 +169 5 0 357 179 67 66 358 180 69 68 +170 5 0 216 70 71 181 358 68 69 180 +171 5 0 15 14 26 43 264 263 287 343 +172 5 0 264 263 287 343 267 266 288 344 +173 5 0 267 266 288 344 270 269 289 345 +174 5 0 270 269 289 345 273 272 290 346 +175 5 0 195 213 199 194 273 346 290 272 +176 5 0 42 9 12 31 339 249 248 297 +177 5 0 339 249 248 297 340 253 252 300 +178 5 0 340 253 252 300 341 257 256 303 +179 5 0 341 257 256 303 342 261 260 306 +180 5 0 212 201 192 189 342 306 260 261 +181 5 0 27 26 41 47 132 287 335 127 +182 5 0 132 287 335 127 133 288 336 128 +183 5 0 133 288 336 128 134 289 337 129 +184 5 0 134 289 337 129 135 290 338 130 +185 5 0 136 131 211 199 135 130 338 290 +186 5 0 29 28 8 56 162 291 232 359 +187 5 0 162 291 232 359 163 292 236 361 +188 5 0 163 292 236 361 164 293 240 363 +189 5 0 164 293 240 363 165 294 244 365 +190 5 0 166 218 188 200 165 365 244 294 +191 5 0 30 10 38 49 72 250 323 62 +192 5 0 72 250 323 62 73 254 324 64 +193 5 0 73 254 324 64 74 258 325 66 +194 5 0 74 258 325 66 75 262 326 68 +195 5 0 76 70 208 190 75 68 326 262 +196 5 0 53 1 4 40 102 107 220 328 +197 5 0 102 107 220 328 103 108 223 330 +198 5 0 103 108 223 330 104 109 226 332 +199 5 0 104 109 226 332 105 110 229 334 +200 5 0 106 210 184 111 105 334 229 110 +201 5 0 16 50 36 13 112 117 316 265 +202 5 0 112 117 316 265 113 118 318 268 +203 5 0 113 118 318 268 114 119 320 271 +204 5 0 114 119 320 271 115 120 322 274 +205 5 0 116 193 206 121 115 274 322 120 +206 5 0 24 45 44 25 82 87 347 283 +207 5 0 82 87 347 283 83 88 348 284 +208 5 0 83 88 348 284 84 89 349 285 +209 5 0 84 89 349 285 85 90 350 286 +210 5 0 86 198 214 91 85 286 350 90 +211 5 0 21 51 34 22 172 177 308 275 +212 5 0 172 177 308 275 173 178 310 277 +213 5 0 173 178 310 277 174 179 312 279 +214 5 0 174 179 312 279 175 180 314 281 +215 5 0 176 197 204 181 175 281 314 180 +216 5 0 56 8 7 18 359 232 231 142 +217 5 0 359 232 231 142 361 236 235 143 +218 5 0 361 236 235 143 363 240 239 144 +219 5 0 363 240 239 144 365 244 243 145 +220 5 0 218 146 187 188 365 145 243 244 +221 5 0 33 32 3 2 296 295 219 221 +222 5 0 296 295 219 221 299 298 222 224 +223 5 0 299 298 222 224 302 301 225 227 +224 5 0 302 301 225 227 305 304 228 230 +225 5 0 203 182 183 202 305 230 228 304 +226 5 0 5 43 26 6 233 343 287 234 +227 5 0 233 343 287 234 237 344 288 238 +228 5 0 237 344 288 238 241 345 289 242 +229 5 0 241 345 289 242 245 346 290 246 +230 5 0 185 186 199 213 245 246 290 346 +231 5 0 55 54 57 60 157 360 147 152 +232 5 0 157 360 147 152 158 362 148 153 +233 5 0 158 362 148 153 159 364 149 154 +234 5 0 159 364 149 154 160 366 150 155 +235 5 0 161 156 151 217 160 155 150 366 +236 5 0 52 39 45 61 97 327 87 92 +237 5 0 97 327 87 92 98 329 88 93 +238 5 0 98 329 88 93 99 331 89 94 +239 5 0 99 331 89 94 100 333 90 95 +240 5 0 101 96 91 209 100 95 90 333 +$EndElements diff --git a/test/blockstructured/poisson/poisson_tensor.mini b/test/blockstructured/poisson/poisson_tensor.mini index 9eb55a9a1fdafcfc2e862d51b08365700e552b28..c2922469f5d631d8647edbee40f04133b5908b8d 100644 --- a/test/blockstructured/poisson/poisson_tensor.mini +++ b/test/blockstructured/poisson/poisson_tensor.mini @@ -1,8 +1,16 @@ __name = blockstructured_poisson_tensor_{__exec_suffix} -__exec_suffix = nonvec, vec | expand vectorized +__exec_suffix = {grid_suffix}_{vec_suffix}_{dim_suffix} -cells = 10 10 -extension = 1. 1. +dim = 2, 3 | expand dimension + +grid_suffix = structured, unstructured | expand unstructured +vec_suffix = nonvec, vec | expand vectorized +dim_suffix = 2d, 3d | expand dimension + +cells = 10, 5 | expand dimension | repeat {dim} +extension = 1. | repeat {dim} + +gmshFile = square.msh, cube.msh | expand dimension [wrapper.vtkcompare] name = {__name} @@ -11,10 +19,15 @@ extension = vtu [formcompiler] compare_l2errorsquared = 1e-7 +grid_unstructured = 0, 1 | expand unstructured [formcompiler.r] matrix_free = 1 vectorization_blockstructured = 0, 1 | expand vectorized generate_jacobians = 0 blockstructured = 1 -number_of_blocks = 8 +number_of_blocks = 8, 4 | expand dimension + +[formcompiler.ufl_variants] +cell = quadrilateral, hexahedron | expand dimension +dim = 2, 3 | expand dimension \ No newline at end of file diff --git a/test/blockstructured/poisson/poisson_tensor.ufl b/test/blockstructured/poisson/poisson_tensor.ufl index 239e14fca9550031cd53994d9528cb0c41e424a6..568e697ee377fa93b3a31468d4af3009daf7366b 100644 --- a/test/blockstructured/poisson/poisson_tensor.ufl +++ b/test/blockstructured/poisson/poisson_tensor.ufl @@ -1,12 +1,10 @@ -cell = quadrilateral - x = SpatialCoordinate(cell) -I = Identity(2) -A = as_matrix([[x[i]*x[j] + I[i,j] for j in range(2)] for i in range(2)]) -g = x[0]**2 + x[1]**2 -c = 8. -f = -4. +I = Identity(dim) +A = as_matrix([[x[i]*x[j] + I[i,j] for j in range(dim)] for i in range(dim)]) +g = sum([c**2 for c in x]) +c = 6. + (dim - 1) * 2 +f = - 2. * dim V = FiniteElement("CG", cell, 1) u = TrialFunction(V) diff --git a/test/blockstructured/poisson/square.geo b/test/blockstructured/poisson/square.geo new file mode 100644 index 0000000000000000000000000000000000000000..8a630d3011a5c0ff219a141d546086af864f3655 --- /dev/null +++ b/test/blockstructured/poisson/square.geo @@ -0,0 +1,28 @@ +// Gmsh project created on Tue Oct 23 11:21:07 2018 +SetFactory("OpenCASCADE"); +//+ +Point(1) = {0, 0, 0, 0.23}; +//+ +Point(2) = {0, 1, 0, 0.23}; +//+ +Point(3) = {1, 1, 0, 0.23}; +//+ +Point(4) = {1, 0, 0, 0.23}; +//+ +Line(1) = {1, 4}; +//+ +Line(2) = {4, 3}; +//+ +Line(3) = {3, 2}; +//+ +Line(4) = {2, 1}; +//+ +Line Loop(1) = {3, 4, 1, 2}; +//+ +Plane Surface(1) = {1}; +//+ +Recombine Surface {1}; +//+ +Extrude {0, 0, 1} { + Surface{1}; Layers{5}; Recombine; +}