pages v0.2.3

Universite-Gustave-Eiffel · Apr 23, 2024 · 6663d4a · 6663d4a
1 parent 24e28a6
commit 6663d4a
Show file tree

Hide file tree

Showing 143 changed files with 2,703 additions and 9,276 deletions.
diff --git a/docs/.buildinfo b/docs/.buildinfo
@@ -1,4 +1,4 @@
 # Sphinx build info version 1
 # This file hashes the configuration used when building these files. When it is not found, a full rebuild will be done.
-config: 19a2012b399ad64330f3f1b4a80396e2
+config: 8404aef1a2ecbed613652486840678e2
 tags: 645f666f9bcd5a90fca523b33c5a78b7
diff --git a/docs/.doctrees/api/elastodynamicsx.plot.doctree b/docs/.doctrees/api/elastodynamicsx.plot.doctree
diff --git a/docs/.doctrees/demos/_ln_demo/eigs_3D_AluminumCube.doctree b/docs/.doctrees/demos/_ln_demo/eigs_3D_AluminumCube.doctree
diff --git a/docs/.doctrees/demos/_ln_demo/eigs_3D_ElasticBeam.doctree b/docs/.doctrees/demos/_ln_demo/eigs_3D_ElasticBeam.doctree
diff --git a/docs/.doctrees/demos/_ln_demo/freq_2D-PSV_FullSpace.doctree b/docs/.doctrees/demos/_ln_demo/freq_2D-PSV_FullSpace.doctree
diff --git a/docs/.doctrees/demos/_ln_demo/freq_2D-SH_FullSpace.doctree b/docs/.doctrees/demos/_ln_demo/freq_2D-SH_FullSpace.doctree
diff --git a/docs/.doctrees/demos/_ln_demo/tdsdyn_3D_ElasticBeam.doctree b/docs/.doctrees/demos/_ln_demo/tdsdyn_3D_ElasticBeam.doctree
diff --git a/docs/.doctrees/demos/_ln_demo/weq_2D-PSV_FullSpace.doctree b/docs/.doctrees/demos/_ln_demo/weq_2D-PSV_FullSpace.doctree
diff --git a/docs/.doctrees/demos/_ln_demo/weq_2D-PSV_HalfSpace_Lamb_KomatitschVilotte_BSSA1998.doctree b/docs/.doctrees/demos/_ln_demo/weq_2D-PSV_HalfSpace_Lamb_KomatitschVilotte_BSSA1998.doctree
diff --git a/docs/.doctrees/demos/_ln_demo/weq_2D-SH_FullSpace.doctree b/docs/.doctrees/demos/_ln_demo/weq_2D-SH_FullSpace.doctree
diff --git a/docs/.doctrees/demos/_ln_demo/weqnl_1D-PSV_Murnaghan_Pwave.doctree b/docs/.doctrees/demos/_ln_demo/weqnl_1D-PSV_Murnaghan_Pwave.doctree
diff --git a/docs/.doctrees/environment.pickle b/docs/.doctrees/environment.pickle
diff --git a/docs/.doctrees/gettingstarted/installation.doctree b/docs/.doctrees/gettingstarted/installation.doctree
diff --git a/docs/.doctrees/gettingstarted/overview.doctree b/docs/.doctrees/gettingstarted/overview.doctree
diff --git a/docs/.doctrees/nbsphinx/demos/_ln_demo/eigs_3D_AluminumCube.ipynb b/docs/.doctrees/nbsphinx/demos/_ln_demo/eigs_3D_AluminumCube.ipynb
@@ -2,7 +2,7 @@
  "cells": [
   {
    "cell_type": "markdown",
-   "id": "05e1f859",
+   "id": "5a86fc7c",
    "metadata": {},
    "source": [
     "# Resonances of an aluminum cube\n",
@@ -15,13 +15,13 @@
   {
    "cell_type": "code",
    "execution_count": 1,
-   "id": "b912b17a",
+   "id": "6f9110bd",
    "metadata": {
     "execution": {
-     "iopub.execute_input": "2024-03-29T23:22:34.318202Z",
-     "iopub.status.busy": "2024-03-29T23:22:34.317988Z",
-     "iopub.status.idle": "2024-03-29T23:22:35.177352Z",
-     "shell.execute_reply": "2024-03-29T23:22:35.176736Z"
+     "iopub.execute_input": "2024-04-23T10:14:19.924815Z",
+     "iopub.status.busy": "2024-04-23T10:14:19.924106Z",
+     "iopub.status.idle": "2024-04-23T10:14:20.751559Z",
+     "shell.execute_reply": "2024-04-23T10:14:20.749779Z"
     }
    },
    "outputs": [],
@@ -30,15 +30,14 @@
     "\n",
     "from dolfinx import mesh, fem, default_scalar_type\n",
     "from mpi4py import MPI\n",
-    "from petsc4py import PETSc\n",
     "\n",
     "from elastodynamicsx.pde import material, PDE\n",
     "from elastodynamicsx.solvers import EigenmodesSolver"
    ]
   },
   {
    "cell_type": "markdown",
-   "id": "dc112e91",
+   "id": "eb755d1d",
    "metadata": {},
    "source": [
     "### FE domain"
@@ -47,13 +46,13 @@
   {
    "cell_type": "code",
    "execution_count": 2,
-   "id": "76ee3541",
+   "id": "07f01203",
    "metadata": {
     "execution": {
-     "iopub.execute_input": "2024-03-29T23:22:35.179592Z",
-     "iopub.status.busy": "2024-03-29T23:22:35.179398Z",
-     "iopub.status.idle": "2024-03-29T23:22:35.190592Z",
-     "shell.execute_reply": "2024-03-29T23:22:35.190034Z"
+     "iopub.execute_input": "2024-04-23T10:14:20.757856Z",
+     "iopub.status.busy": "2024-04-23T10:14:20.757314Z",
+     "iopub.status.idle": "2024-04-23T10:14:20.777691Z",
+     "shell.execute_reply": "2024-04-23T10:14:20.776946Z"
     }
    },
    "outputs": [],
@@ -68,7 +67,7 @@
   },
   {
    "cell_type": "markdown",
-   "id": "5cd792a1",
+   "id": "4ecf6459",
    "metadata": {},
    "source": [
     "### Define the material law\n",
@@ -85,13 +84,13 @@
   {
    "cell_type": "code",
    "execution_count": 3,
-   "id": "2c58771d",
+   "id": "6e80bbd6",
    "metadata": {
     "execution": {
-     "iopub.execute_input": "2024-03-29T23:22:35.192818Z",
-     "iopub.status.busy": "2024-03-29T23:22:35.192620Z",
-     "iopub.status.idle": "2024-03-29T23:22:35.196589Z",
-     "shell.execute_reply": "2024-03-29T23:22:35.196236Z"
+     "iopub.execute_input": "2024-04-23T10:14:20.780235Z",
+     "iopub.status.busy": "2024-04-23T10:14:20.780035Z",
+     "iopub.status.idle": "2024-04-23T10:14:20.784583Z",
+     "shell.execute_reply": "2024-04-23T10:14:20.784113Z"
     }
    },
    "outputs": [],
@@ -108,7 +107,7 @@
   },
   {
    "cell_type": "markdown",
-   "id": "c6990aa7",
+   "id": "3d6faf79",
    "metadata": {},
    "source": [
     "### Assemble the PDE"
@@ -117,13 +116,13 @@
   {
    "cell_type": "code",
    "execution_count": 4,
-   "id": "e75f2d28",
+   "id": "6af3cf34",
    "metadata": {
     "execution": {
-     "iopub.execute_input": "2024-03-29T23:22:35.198529Z",
-     "iopub.status.busy": "2024-03-29T23:22:35.198228Z",
-     "iopub.status.idle": "2024-03-29T23:22:35.200508Z",
-     "shell.execute_reply": "2024-03-29T23:22:35.200166Z"
+     "iopub.execute_input": "2024-04-23T10:14:20.786580Z",
+     "iopub.status.busy": "2024-04-23T10:14:20.786367Z",
+     "iopub.status.idle": "2024-04-23T10:14:20.788623Z",
+     "shell.execute_reply": "2024-04-23T10:14:20.788275Z"
     }
    },
    "outputs": [],
@@ -133,7 +132,7 @@
   },
   {
    "cell_type": "markdown",
-   "id": "8fb5ce2e",
+   "id": "c3d511c4",
    "metadata": {},
    "source": [
     "### Solve"
@@ -142,13 +141,13 @@
   {
    "cell_type": "code",
    "execution_count": 5,
-   "id": "b13578fa",
+   "id": "3048ea69",
    "metadata": {
     "execution": {
-     "iopub.execute_input": "2024-03-29T23:22:35.202267Z",
-     "iopub.status.busy": "2024-03-29T23:22:35.202100Z",
-     "iopub.status.idle": "2024-03-29T23:22:35.275852Z",
-     "shell.execute_reply": "2024-03-29T23:22:35.275358Z"
+     "iopub.execute_input": "2024-04-23T10:14:20.790639Z",
+     "iopub.status.busy": "2024-04-23T10:14:20.790295Z",
+     "iopub.status.idle": "2024-04-23T10:14:20.864133Z",
+     "shell.execute_reply": "2024-04-23T10:14:20.863455Z"
     }
    },
    "outputs": [],
@@ -163,13 +162,13 @@
   {
    "cell_type": "code",
    "execution_count": 6,
-   "id": "3362721b",
+   "id": "bba2151a",
    "metadata": {
     "execution": {
-     "iopub.execute_input": "2024-03-29T23:22:35.277849Z",
-     "iopub.status.busy": "2024-03-29T23:22:35.277688Z",
-     "iopub.status.idle": "2024-03-29T23:22:36.291776Z",
-     "shell.execute_reply": "2024-03-29T23:22:36.291234Z"
+     "iopub.execute_input": "2024-04-23T10:14:20.866178Z",
+     "iopub.status.busy": "2024-04-23T10:14:20.866022Z",
+     "iopub.status.idle": "2024-04-23T10:14:21.867734Z",
+     "shell.execute_reply": "2024-04-23T10:14:21.867167Z"
     }
    },
    "outputs": [],
@@ -182,13 +181,13 @@
   {
    "cell_type": "code",
    "execution_count": 7,
-   "id": "5ff05e17",
+   "id": "6178177a",
    "metadata": {
     "execution": {
-     "iopub.execute_input": "2024-03-29T23:22:36.294292Z",
-     "iopub.status.busy": "2024-03-29T23:22:36.293963Z",
-     "iopub.status.idle": "2024-03-29T23:22:37.400534Z",
-     "shell.execute_reply": "2024-03-29T23:22:37.400133Z"
+     "iopub.execute_input": "2024-04-23T10:14:21.869567Z",
+     "iopub.status.busy": "2024-04-23T10:14:21.869415Z",
+     "iopub.status.idle": "2024-04-23T10:14:22.935297Z",
+     "shell.execute_reply": "2024-04-23T10:14:22.934882Z"
     },
     "lines_to_next_cell": 2
    },
@@ -211,12 +210,12 @@
     "eigenfreqs = eps.getEigenfrequencies()\n",
     "# eigenmodes = eps.getEigenmodes()\n",
     "\n",
-    "eps.plot(V, slice(6,6+9), wireframe=True, factor=30)  # Avoid the first 6 rigid body modes"
+    "eps.plot(V, slice(6, 6+9), wireframe=True, factor=30)  # Avoid the first 6 rigid body modes"
    ]
   },
   {
    "cell_type": "markdown",
-   "id": "afb89b60",
+   "id": "778492eb",
    "metadata": {},
    "source": [
     "### Compare with literature values\n",
@@ -230,13 +229,13 @@
   {
    "cell_type": "code",
    "execution_count": 8,
-   "id": "ad17287f",
+   "id": "4ba16f69",
    "metadata": {
     "execution": {
-     "iopub.execute_input": "2024-03-29T23:22:37.406883Z",
-     "iopub.status.busy": "2024-03-29T23:22:37.406717Z",
-     "iopub.status.idle": "2024-03-29T23:22:37.412365Z",
-     "shell.execute_reply": "2024-03-29T23:22:37.411992Z"
+     "iopub.execute_input": "2024-04-23T10:14:22.941887Z",
+     "iopub.status.busy": "2024-04-23T10:14:22.941671Z",
+     "iopub.status.idle": "2024-04-23T10:14:22.946718Z",
+     "shell.execute_reply": "2024-04-23T10:14:22.946232Z"
     }
    },
    "outputs": [
@@ -246,19 +245,19 @@
      "text": [
       "Eigenfrequencies: comparison with litterature values\n",
       "  FE   \tOgi et al, calc.\t Ogi et al, exp. \t(kHz)\n",
-      "116.414\t     116.32\t\t     116.716\n",
+      "116.414\t     116.320\t\t     116.716\n",
       "143.605\t     143.186\t\t     143.783\n",
-      "158.517\t     158.44\t\t     158.081\n",
-      "166.27\t     166.113\t\t     166.5\n",
+      "158.517\t     158.440\t\t     158.081\n",
+      "166.270\t     166.113\t\t     166.500\n",
       "169.457\t     169.338\t\t     169.523\n",
-      "178.607\t     178.36\t\t     177.846\n",
-      "184.935\t     184.57\t\t     183.875\n",
+      "178.607\t     178.360\t\t     177.846\n",
+      "184.935\t     184.570\t\t     183.875\n",
       "185.384\t     185.078\t\t     186.047\n",
       "190.293\t     190.206\t\t     190.341\n",
       "198.003\t     197.692\t\t     197.386\n",
       "202.098\t     201.462\t\t     201.133\n",
       "207.531\t     207.096\t\t     207.386\n",
-      "211.323\t     211.0\t\t     209.836\n",
+      "211.323\t     211.000\t\t     209.836\n",
       "216.115\t     215.613\t\t     214.753\n"
      ]
     }
@@ -271,17 +270,18 @@
     "                               310.109, 316.197, 317.392, 326.462, 329.034, 332.441, 333.364, 336.65,\n",
     "                               337.359, 338.276])\n",
     "\n",
-    "freqs_OgiEtAl_calc = np.array([116.32 , 143.186, 158.44 , 166.113, 169.338, 178.36 , 184.57 , 185.078, \\\n",
-    "                               190.206, 197.692, 201.462, 207.096, 211    , 215.613, 223.219, 230.804, \\\n",
-    "                               233.329, 234.758, 250.777, 251.038, 252.303, 256.849, 258.064, 258.874, \\\n",
-    "                               259.203, 267.746, 276.736, 279.144, 282.773, 293.016, 304.593, 305.316, \\\n",
-    "                               309.591, 315.775, 317.931, 326.556, 329.369, 332.732, 332.271, 336.218, \\\n",
+    "freqs_OgiEtAl_calc = np.array([116.32 , 143.186, 158.44 , 166.113, 169.338, 178.36 , 184.57 , 185.078,\n",
+    "                               190.206, 197.692, 201.462, 207.096, 211    , 215.613, 223.219, 230.804,\n",
+    "                               233.329, 234.758, 250.777, 251.038, 252.303, 256.849, 258.064, 258.874,\n",
+    "                               259.203, 267.746, 276.736, 279.144, 282.773, 293.016, 304.593, 305.316,\n",
+    "                               309.591, 315.775, 317.931, 326.556, 329.369, 332.732, 332.271, 336.218,\n",
     "                               337.511, 337.71])\n",
     "\n",
     "print('Eigenfrequencies: comparison with litterature values')\n",
     "print('  FE   \\tOgi et al, calc.\\t Ogi et al, exp. \\t(kHz)')\n",
-    "for fFE, fOgi_calc, fOgi_exp in zip(eigenfreqs[6:]*1e3, freqs_OgiEtAl_calc, freqs_OgiEtAl_exp):  # *1e3 to convert MHz into kHz\n",
-    "    print(str(round(fFE, 3)) +\"\\t     \"+ str(round(fOgi_calc, 3)) +\"\\t\\t     \"+ str(round(fOgi_exp, 3)))"
+    "for fFE, fOgi_calc, fOgi_exp in zip(eigenfreqs[6:] * 1e3, freqs_OgiEtAl_calc, freqs_OgiEtAl_exp):\n",
+    "    # *1e3 to convert MHz into kHz\n",
+    "    print(f\"{fFE:.3f}\\t     {fOgi_calc:.3f}\\t\\t     {fOgi_exp:.3f}\")"
    ]
   }
  ],