Add new examples

Signed-off-by: Tom Freudenberg <[email protected]>
TomF98 · Oct 25, 2023 · 0c02911 · 0c02911
1 parent 76cbd65
commit 0c02911
Show file tree

Hide file tree

Showing 13 changed files with 6,028 additions and 10 deletions.
diff --git a/examples/backup/BackUp_1.ipynb b/examples/backup/BackUp_1.ipynb
diff --git a/examples/backup/BackUp_2.ipynb b/examples/backup/BackUp_2.ipynb
diff --git a/examples/backup/BackUp_Sol_1.ipynb b/examples/backup/BackUp_Sol_1.ipynb
diff --git a/examples/backup/BackUp_Sol_2.ipynb b/examples/backup/BackUp_Sol_2.ipynb
diff --git a/examples/new_examples/DrillingData/coordinates.npy b/examples/new_examples/DrillingData/coordinates.npy
diff --git a/examples/new_examples/DrillingData/model_drilling.pt b/examples/new_examples/DrillingData/model_drilling.pt
diff --git a/examples/new_examples/DrillingData/solution.npy b/examples/new_examples/DrillingData/solution.npy
diff --git a/examples/new_examples/DrillingData/time_points.npy b/examples/new_examples/DrillingData/time_points.npy
diff --git a/examples/new_examples/Exercise_2.ipynb b/examples/new_examples/Exercise_2.ipynb
@@ -69,7 +69,7 @@
         "import torchphysics as tp\n",
         "X = tp.spaces.R2('x')\n",
         "U = tp.spaces.R1('u')\n",
-        "T = tp.spaces. # <- Add the time variable \"t\" of dimension 1 (e.g. R1)"
+        "T = tp.spaces. # TODO: Add the time variable \"t\" of dimension 1 (e.g. R1)"
       ]
     },
     {
@@ -91,8 +91,8 @@
       "outputs": [],
       "source": [
         "omega = tp.domains.Parallelogram(X, [0,0], [1,0], [0,1])\n",
-        "time_interval = tp.domains.Interval(T, ) # <-add the bounds of the Interval (0, 2)\n",
-        "product_domain = time_interval # <- products are define with: * omega"
+        "time_interval = tp.domains.Interval(T, ) # TODO: Add the bounds of the Interval (0, 2)\n",
+        "product_domain = time_interval # TODO: Create the product domain. Products are define with: * omega"
       ]
     },
     {
@@ -113,13 +113,13 @@
       },
       "outputs": [],
       "source": [
-        "# Add the product domain of time and space:\n",
+        "# TODO: Add the product domain of time and space:\n",
         "inner_sampler = tp.samplers.RandomUniformSampler(, n_points=25000) \n",
         "\n",
         "# The boundary sampler is done already.\n",
         "bound_sampler = tp.samplers.RandomUniformSampler(time_interval*omega.boundary, n_points=10000)\n",
         "\n",
-        "# Currently only the left interval side {0} is passed in the initial sampler, create the product domain with omega:\n",
+        "# Currently only the left interval side {0} is passed in the initial sampler. TODO: Create the product domain with omega:\n",
         "initial_sampler = tp.samplers.RandomUniformSampler(time_interval.boundary_left, n_points=5000)"
       ]
     },
@@ -143,6 +143,7 @@
       },
       "outputs": [],
       "source": [
+        "# TODO: Add spaces\n",
         "model = tp.models.FCN(input_space=, output_space=, hidden=(30,30,30))"
       ]
     },
@@ -167,7 +168,7 @@
       "outputs": [],
       "source": [
         "def pde_residual(u, x, t):\n",
-        "    # in the differential operators you have to pass in the correct variables \n",
+        "    # TODO: Pass in the correct variables \n",
         "    # for the derivative computation as the second argument:\n",
         "    return tp.utils.grad(u, ) - 0.1*tp.utils.laplacian(u, ) - 1.0\n",
         "\n",
@@ -212,11 +213,11 @@
       "metadata": {},
       "outputs": [],
       "source": [
-        "# Implement the residual for the initial condition:\n",
+        "# TODO: Implement the residual for the initial condition:\n",
         "def initial_residual(u):\n",
         "    return \n",
         "\n",
-        "initial_cond = tp.conditions.PINNCondition() # add the model, sampler and residual"
+        "initial_cond = tp.conditions.PINNCondition() # TODO: Add the model, sampler and residual"
       ]
     },
     {
@@ -238,7 +239,7 @@
       "outputs": [],
       "source": [
         "optim = tp.OptimizerSetting(torch.optim.Adam, lr=0.005)\n",
-        "solver = tp.solver.Solver([.,.], optimizer_setting=optim)"
+        "solver = tp.solver.Solver([.,.], optimizer_setting=optim) # TODO: Collect all conditions"
       ]
     },
     {
@@ -282,7 +283,7 @@
       "metadata": {},
       "outputs": [],
       "source": [
-        "plot_sampler = tp.samplers.PlotSampler(plot_domain=omega, n_points=2000, data_for_other_variables={\"t\": 0.0})\n",
+        "plot_sampler = tp.samplers.PlotSampler(plot_domain=omega, n_points=2000, data_for_other_variables={\"t\": 0.1})\n",
         "fig = tp.utils.plot(model, lambda u : u, plot_sampler)\n",
         "\n",
         "\n",