|
68 | 68 | "cell_type": "code", |
69 | 69 | "execution_count": null, |
70 | 70 | "id": "2", |
71 | | - "metadata": { |
72 | | - "vscode": { |
73 | | - "languageId": "python" |
74 | | - } |
75 | | - }, |
| 71 | + "metadata": {}, |
76 | 72 | "outputs": [], |
77 | 73 | "source": [ |
78 | 74 | "from mpi4py import MPI\n", |
|
105 | 101 | "cell_type": "code", |
106 | 102 | "execution_count": null, |
107 | 103 | "id": "4", |
108 | | - "metadata": { |
109 | | - "vscode": { |
110 | | - "languageId": "python" |
111 | | - } |
112 | | - }, |
| 104 | + "metadata": {}, |
113 | 105 | "outputs": [], |
114 | 106 | "source": [ |
115 | 107 | "from dolfinx.fem import functionspace\n", |
|
120 | 112 | "cell_type": "code", |
121 | 113 | "execution_count": null, |
122 | 114 | "id": "5", |
123 | | - "metadata": { |
124 | | - "vscode": { |
125 | | - "languageId": "python" |
126 | | - } |
127 | | - }, |
| 115 | + "metadata": {}, |
128 | 116 | "outputs": [], |
129 | 117 | "source": [ |
130 | 118 | "from dolfinx import fem\n", |
|
146 | 134 | "cell_type": "code", |
147 | 135 | "execution_count": null, |
148 | 136 | "id": "7", |
149 | | - "metadata": { |
150 | | - "vscode": { |
151 | | - "languageId": "python" |
152 | | - } |
153 | | - }, |
| 137 | + "metadata": {}, |
154 | 138 | "outputs": [], |
155 | 139 | "source": [ |
156 | 140 | "import numpy\n", |
|
180 | 164 | "cell_type": "code", |
181 | 165 | "execution_count": null, |
182 | 166 | "id": "9", |
183 | | - "metadata": { |
184 | | - "vscode": { |
185 | | - "languageId": "python" |
186 | | - } |
187 | | - }, |
| 167 | + "metadata": {}, |
188 | 168 | "outputs": [], |
189 | 169 | "source": [ |
190 | 170 | "boundary_dofs = fem.locate_dofs_topological(V, fdim, boundary_facets)\n", |
|
208 | 188 | "cell_type": "code", |
209 | 189 | "execution_count": null, |
210 | 190 | "id": "11", |
211 | | - "metadata": { |
212 | | - "vscode": { |
213 | | - "languageId": "python" |
214 | | - } |
215 | | - }, |
| 191 | + "metadata": {}, |
216 | 192 | "outputs": [], |
217 | 193 | "source": [ |
218 | 194 | "import ufl\n", |
|
233 | 209 | "cell_type": "code", |
234 | 210 | "execution_count": null, |
235 | 211 | "id": "13", |
236 | | - "metadata": { |
237 | | - "vscode": { |
238 | | - "languageId": "python" |
239 | | - } |
240 | | - }, |
| 212 | + "metadata": {}, |
241 | 213 | "outputs": [], |
242 | 214 | "source": [ |
243 | 215 | "from dolfinx import default_scalar_type\n", |
|
261 | 233 | "cell_type": "code", |
262 | 234 | "execution_count": null, |
263 | 235 | "id": "15", |
264 | | - "metadata": { |
265 | | - "vscode": { |
266 | | - "languageId": "python" |
267 | | - } |
268 | | - }, |
| 236 | + "metadata": {}, |
269 | 237 | "outputs": [], |
270 | 238 | "source": [ |
271 | 239 | "a = ufl.dot(ufl.grad(u), ufl.grad(v)) * ufl.dx\n", |
|
306 | 274 | "cell_type": "code", |
307 | 275 | "execution_count": null, |
308 | 276 | "id": "17", |
309 | | - "metadata": { |
310 | | - "vscode": { |
311 | | - "languageId": "python" |
312 | | - } |
313 | | - }, |
| 277 | + "metadata": {}, |
314 | 278 | "outputs": [], |
315 | 279 | "source": [ |
316 | 280 | "from dolfinx.fem.petsc import LinearProblem\n", |
|
332 | 296 | "cell_type": "code", |
333 | 297 | "execution_count": null, |
334 | 298 | "id": "19", |
335 | | - "metadata": { |
336 | | - "vscode": { |
337 | | - "languageId": "python" |
338 | | - } |
339 | | - }, |
| 299 | + "metadata": {}, |
340 | 300 | "outputs": [], |
341 | 301 | "source": [ |
342 | 302 | "V2 = fem.functionspace(domain, (\"Lagrange\", 2))\n", |
|
357 | 317 | "cell_type": "code", |
358 | 318 | "execution_count": null, |
359 | 319 | "id": "21", |
360 | | - "metadata": { |
361 | | - "vscode": { |
362 | | - "languageId": "python" |
363 | | - } |
364 | | - }, |
| 320 | + "metadata": {}, |
365 | 321 | "outputs": [], |
366 | 322 | "source": [ |
367 | 323 | "L2_error = fem.form(ufl.inner(uh - uex, uh - uex) * ufl.dx)\n", |
|
386 | 342 | "cell_type": "code", |
387 | 343 | "execution_count": null, |
388 | 344 | "id": "23", |
389 | | - "metadata": { |
390 | | - "vscode": { |
391 | | - "languageId": "python" |
392 | | - } |
393 | | - }, |
| 345 | + "metadata": {}, |
394 | 346 | "outputs": [], |
395 | 347 | "source": [ |
396 | 348 | "error_max = numpy.max(numpy.abs(uD.x.array-uh.x.array))\n", |
|
429 | 381 | "cell_type": "code", |
430 | 382 | "execution_count": null, |
431 | 383 | "id": "26", |
432 | | - "metadata": { |
433 | | - "vscode": { |
434 | | - "languageId": "python" |
435 | | - } |
436 | | - }, |
| 384 | + "metadata": {}, |
437 | 385 | "outputs": [], |
438 | 386 | "source": [ |
439 | 387 | "from dolfinx import plot\n", |
440 | | - "pyvista.start_xvfb()\n", |
441 | 388 | "domain.topology.create_connectivity(tdim, tdim)\n", |
442 | 389 | "topology, cell_types, geometry = plot.vtk_mesh(domain, tdim)\n", |
443 | 390 | "grid = pyvista.UnstructuredGrid(topology, cell_types, geometry)" |
|
464 | 411 | "cell_type": "code", |
465 | 412 | "execution_count": null, |
466 | 413 | "id": "29", |
467 | | - "metadata": { |
468 | | - "vscode": { |
469 | | - "languageId": "python" |
470 | | - } |
471 | | - }, |
| 414 | + "metadata": {}, |
472 | 415 | "outputs": [], |
473 | 416 | "source": [ |
474 | 417 | "plotter = pyvista.Plotter()\n", |
|
493 | 436 | "cell_type": "code", |
494 | 437 | "execution_count": null, |
495 | 438 | "id": "31", |
496 | | - "metadata": { |
497 | | - "vscode": { |
498 | | - "languageId": "python" |
499 | | - } |
500 | | - }, |
| 439 | + "metadata": {}, |
501 | 440 | "outputs": [], |
502 | 441 | "source": [ |
503 | 442 | "u_topology, u_cell_types, u_geometry = plot.vtk_mesh(V)" |
|
515 | 454 | "cell_type": "code", |
516 | 455 | "execution_count": null, |
517 | 456 | "id": "33", |
518 | | - "metadata": { |
519 | | - "vscode": { |
520 | | - "languageId": "python" |
521 | | - } |
522 | | - }, |
| 457 | + "metadata": {}, |
523 | 458 | "outputs": [], |
524 | 459 | "source": [ |
525 | 460 | "u_grid = pyvista.UnstructuredGrid(u_topology, u_cell_types, u_geometry)\n", |
|
544 | 479 | "cell_type": "code", |
545 | 480 | "execution_count": null, |
546 | 481 | "id": "35", |
547 | | - "metadata": { |
548 | | - "vscode": { |
549 | | - "languageId": "python" |
550 | | - } |
551 | | - }, |
| 482 | + "metadata": {}, |
552 | 483 | "outputs": [], |
553 | 484 | "source": [ |
554 | 485 | "warped = u_grid.warp_by_scalar()\n", |
|
571 | 502 | "cell_type": "code", |
572 | 503 | "execution_count": null, |
573 | 504 | "id": "37", |
574 | | - "metadata": { |
575 | | - "vscode": { |
576 | | - "languageId": "python" |
577 | | - } |
578 | | - }, |
| 505 | + "metadata": {}, |
579 | 506 | "outputs": [], |
580 | 507 | "source": [ |
581 | 508 | "from dolfinx import io\n", |
|
0 commit comments