From 5927840bb4b810bacbdef77c0c191691d67eb6fd Mon Sep 17 00:00:00 2001
From: Sylwia Majchrowska <s.majchrowska@neurosys.com>
Date: Sat, 31 Jul 2021 17:25:54 +0200
Subject: [PATCH] add possibility to use Debye Fit in gprMax initial script by
 #python

---
 user_libs/DebyeFit/__init__.py              | 11 +++++++++
 user_models/HavriliakNegami_material_box.in | 26 +++++++++++++++++++++
 2 files changed, 37 insertions(+)
 create mode 100644 user_libs/DebyeFit/__init__.py
 create mode 100644 user_models/HavriliakNegami_material_box.in

diff --git a/user_libs/DebyeFit/__init__.py b/user_libs/DebyeFit/__init__.py
new file mode 100644
index 00000000..1b70d5ce
--- /dev/null
+++ b/user_libs/DebyeFit/__init__.py
@@ -0,0 +1,11 @@
+import os, sys; sys.path.append(os.path.dirname(os.path.realpath(__file__)))
+from .Debye_Fit import (HavriliakNegami,
+                       Jonscher,
+                       Crim,
+                       Rawdata)
+
+__all__ = [
+    'HavriliakNegami', 'Jonscher', 'Crim',
+    'Rawdata'
+    ]
+
diff --git a/user_models/HavriliakNegami_material_box.in b/user_models/HavriliakNegami_material_box.in
new file mode 100644
index 00000000..2f38af48
--- /dev/null
+++ b/user_models/HavriliakNegami_material_box.in
@@ -0,0 +1,26 @@
+#domain: 0.7 0.5 0.004
+#dx_dy_dz: 0.004 0.004 0.004
+#time_window: 3000
+
+#python:
+from user_libs.DebyeFit import HavriliakNegami
+setup = HavriliakNegami(f_min=1e4, f_max=1e11,
+                        alpha=0.3, beta=1,
+                        e_inf=3.4, de=2.7, tau_0=.8e-10,
+                        sigma=0.45e-3, mu=1, mu_sigma=0,
+                        material_name="dry_sand",
+                        number_of_debye_poles=5, f_n=100,
+                        plot=False, save=False,
+                        optimizer_options={'swarmsize':30,
+                                           'maxiter':100,
+                                           'omega':0.5,
+                                           'phip':1.4,
+                                           'phig':1.4,
+                                           'minstep':1e-8,
+                                           'minfun':1e-8,
+                                           'seed': 111,
+                                           'pflag': False})
+setup.run()
+#end_python:
+
+#box: 0 0 0 0.07 0.4 0.004 dry_sand