final rev1.0

pescap · Aug 2, 2023 · c369ec3 · c369ec3
1 parent d99e1b8
commit c369ec3
Show file tree

Hide file tree

Showing 18 changed files with 110 additions and 47 deletions.
diff --git a/examples/fichera/mc.py b/examples/fichera/mc.py
@@ -127,7 +127,7 @@ def set_random_seed(seed):
 
 config = bemppUQ.config.set_case("MC")
 
-precision = 10
+precision = 30
 
 h = 2.0 * np.pi / (precision * config["k_int"])
 
@@ -142,7 +142,7 @@ def set_random_seed(seed):
 
 Ut, _ = evaluate_far_field(transmission_operators_t, config)
 
-np.savetxt("results/out" + str(seed) + ".txt", Ut.view(float))
+np.savetxt("results1/out" + str(seed) + ".txt", Ut.view(float))
 
 t0 = time.time() - t0
 print("Execution time: ", t0)
diff --git a/examples/fichera/plot_mc.ipynb b/examples/fichera/plot_mc.ipynb
@@ -9,6 +9,15 @@
     "#for run in {1..20}; do python3 mc.py&python3 mc.py&python3 mc.py&python3 mc.py&python3 mc.py&python3 mc.py&python3 mc.py&python3 mc.py&python3 mc.py&python3 mc.py; done"
    ]
   },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "for run in {1..20}; do python3 mc.py&python3 mc.py&python3 mc.py&python3 mc.py&python3 mc.py; done"
+   ]
+  },
   {
    "cell_type": "code",
    "execution_count": 22,

diff --git a/examples/fichera/plots/fichera0.pdf b/examples/fichera/plots/fichera0.pdf
diff --git a/examples/fichera/plots/fichera1.pdf b/examples/fichera/plots/fichera1.pdf
diff --git a/examples/fichera/plots/fichera2.pdf b/examples/fichera/plots/fichera2.pdf
diff --git a/examples/fichera/plots_mc.ipynb b/examples/fichera/plots_mc.ipynb
diff --git a/examples/fichera/results/full/1_2.txt.npy b/examples/fichera/results/full/1_2.txt.npy
diff --git a/examples/fichera/results/full/6_10.txt.npy b/examples/fichera/results/full/6_10.txt.npy
diff --git a/examples/fichera/solver.py b/examples/fichera/solver.py
@@ -12,7 +12,8 @@
 )
 
 bempp.api.global_parameters.assembly.potential_operator_assembly_type = "dense"
-bempp.api.global_parameters.hmat.eps = 1e-4
+bempp.api.global_parameters.hmat.eps = 1e-6
+print(bempp.api.global_parameters.hmat.eps)
 
 ta = time.time()
 
@@ -94,16 +95,22 @@ def get_base_grid(h):
 parser.add_argument("--l0", default=0, type=int)
 parser.add_argument("--nrank", default=1, type=int)
 
+parser.add_argument("--prec", default=10, type=int)
+
 args = parser.parse_args()
 l0 = args.l0
 nrank = args.nrank
+prec = args.prec
 
 print(l0, "l0")
 
 
 precision_list = [2, 5, 10]
 precision = precision_list[l0]
 
+if prec != 10:
+    precision = prec
+
 h = 2.0 * np.pi / (precision * config["k_int"])
 
 if l0 == 2:

diff --git a/examples/kite/plot_mc.ipynb b/examples/kite/plot_mc.ipynb
diff --git a/examples/kite/plots/kite.pdf b/examples/kite/plots/kite.pdf
diff --git a/examples/tp/.ipynb_checkpoints/plots-transmission_problem-checkpoint.ipynb b/examples/tp/.ipynb_checkpoints/plots-transmission_problem-checkpoint.ipynb
@@ -2,7 +2,7 @@
  "cells": [
   {
    "cell_type": "code",
-   "execution_count": null,
+   "execution_count": 1,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -13,7 +13,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 1,
+   "execution_count": 2,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -96,7 +96,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 24,
+   "execution_count": 3,
    "metadata": {},
    "outputs": [
     {
@@ -139,7 +139,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 19,
+   "execution_count": 4,
    "metadata": {},
    "outputs": [
     {
@@ -225,7 +225,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 25,
+   "execution_count": 5,
    "metadata": {},
    "outputs": [
     {
@@ -275,7 +275,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 26,
+   "execution_count": 6,
    "metadata": {},
    "outputs": [
     {
@@ -336,7 +336,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 27,
+   "execution_count": 7,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -352,7 +352,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 33,
+   "execution_count": 8,
    "metadata": {},
    "outputs": [
     {
@@ -394,7 +394,28 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 36,
+   "execution_count": 9,
+   "metadata": {},
+   "outputs": [
+    {
+     "data": {
+      "text/plain": [
+       "array([ 1.        ,  6.78149363,  9.10193553, 35.3964499 , 33.05938767,\n",
+       "       27.65405125, 17.87109752, 14.33922143])"
+      ]
+     },
+     "execution_count": 9,
+     "metadata": {},
+     "output_type": "execute_result"
+    }
+   ],
+   "source": [
+    "1/compression_hat"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 10,
    "metadata": {},
    "outputs": [
     {
@@ -404,7 +425,7 @@
        "       70.99225214, 45.23645742, 37.85853709])"
       ]
      },
-     "execution_count": 36,
+     "execution_count": 10,
      "metadata": {},
      "output_type": "execute_result"
     }
@@ -415,7 +436,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 30,
+   "execution_count": 11,
    "metadata": {},
    "outputs": [
     {

diff --git a/examples/tp/plots-transmission_problem.ipynb b/examples/tp/plots-transmission_problem.ipynb
@@ -2,7 +2,7 @@
  "cells": [
   {
    "cell_type": "code",
-   "execution_count": null,
+   "execution_count": 1,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -13,7 +13,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 1,
+   "execution_count": 2,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -96,7 +96,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 24,
+   "execution_count": 3,
    "metadata": {},
    "outputs": [
     {
@@ -139,7 +139,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 19,
+   "execution_count": 4,
    "metadata": {},
    "outputs": [
     {
@@ -225,7 +225,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 25,
+   "execution_count": 5,
    "metadata": {},
    "outputs": [
     {
@@ -275,7 +275,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 26,
+   "execution_count": 6,
    "metadata": {},
    "outputs": [
     {
@@ -336,7 +336,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 27,
+   "execution_count": 7,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -352,7 +352,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 33,
+   "execution_count": 8,
    "metadata": {},
    "outputs": [
     {
@@ -394,7 +394,28 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 36,
+   "execution_count": 9,
+   "metadata": {},
+   "outputs": [
+    {
+     "data": {
+      "text/plain": [
+       "array([ 1.        ,  6.78149363,  9.10193553, 35.3964499 , 33.05938767,\n",
+       "       27.65405125, 17.87109752, 14.33922143])"
+      ]
+     },
+     "execution_count": 9,
+     "metadata": {},
+     "output_type": "execute_result"
+    }
+   ],
+   "source": [
+    "1/compression_hat"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 10,
    "metadata": {},
    "outputs": [
     {
@@ -404,7 +425,7 @@
        "       70.99225214, 45.23645742, 37.85853709])"
       ]
      },
-     "execution_count": 36,
+     "execution_count": 10,
      "metadata": {},
      "output_type": "execute_result"
     }
@@ -415,7 +436,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 30,
+   "execution_count": 11,
    "metadata": {},
    "outputs": [
     {

diff --git a/examples/tp/plots/FirstSecond.pdf b/examples/tp/plots/FirstSecond.pdf
diff --git a/examples/tp/plots/SecondDof.pdf b/examples/tp/plots/SecondDof.pdf
diff --git a/examples/tp/plots/SecondEfficiency.pdf b/examples/tp/plots/SecondEfficiency.pdf
diff --git a/examples/tp/plots/Secondh.pdf b/examples/tp/plots/Secondh.pdf
diff --git a/examples/tp/plots/combVsSym0.pdf b/examples/tp/plots/combVsSym0.pdf