From f7c6cc271c7f5109ac7f910920362cec41414e81 Mon Sep 17 00:00:00 2001
From: kittel_c <christoph.kittel@psi.ch>
Date: Wed, 12 Oct 2022 15:46:33 +0200
Subject: [PATCH] changed folder structure, added and renamed files

---
 .gitignore                                    |  5 ++++
 phases/{parallel.py => full_polarization.py}  |  0
 phases/models/__init__.py                     |  0
 .../anti-parallel_model.py}                   |  0
 phases/models/parallel_model.py               | 27 +++++++++++++++++++
 5 files changed, 32 insertions(+)
 rename phases/{parallel.py => full_polarization.py} (100%)
 create mode 100644 phases/models/__init__.py
 rename phases/{model.py => models/anti-parallel_model.py} (100%)
 create mode 100644 phases/models/parallel_model.py

diff --git a/.gitignore b/.gitignore
index 73358ad..031224e 100644
--- a/.gitignore
+++ b/.gitignore
@@ -1,3 +1,8 @@
+# editor related files
+*.swp
+*~
+*.bak
+
 # Byte-compiled / optimized / DLL files
 __pycache__/
 *.py[cod]
diff --git a/phases/parallel.py b/phases/full_polarization.py
similarity index 100%
rename from phases/parallel.py
rename to phases/full_polarization.py
diff --git a/phases/models/__init__.py b/phases/models/__init__.py
new file mode 100644
index 0000000..e69de29
diff --git a/phases/model.py b/phases/models/anti-parallel_model.py
similarity index 100%
rename from phases/model.py
rename to phases/models/anti-parallel_model.py
diff --git a/phases/models/parallel_model.py b/phases/models/parallel_model.py
new file mode 100644
index 0000000..5f3e969
--- /dev/null
+++ b/phases/models/parallel_model.py
@@ -0,0 +1,27 @@
+import numpy as np
+from scipy.optimize import fsolve
+
+
+def parallel2gap(K, phi, undudict):
+    gLH = K2gap(K, undudict['K-value_LH'])
+    if phi >= 0.0:
+        gLV = K2gap(K, undudict['K-value_LV+'])
+        gC = K2gap(K, undudict['K-value_C+'])
+        dgLV = gLV - gLH
+        dgC = gC - gLH - dgLV/2
+    else:
+        gLV = K2gap(K, undudict['K-value_LV-'])
+        gC = K2gap(K, undudict['K-value_C-'])
+        dgLV = gLV - gLH
+        dgC = gC - gLH - dgLV/2
+    return gLH + dgLV * np.sin(0.5 * phi)**2 + dgC * np.sin(phi)**2
+
+
+def K2gap(Kval, fitparam):
+    g2K_func = np.poly1d(fitparam[::-1])
+    tau_init = 1.0
+    k_log = float(np.log(Kval))
+    return float(fsolve(k_log - g2K_func, tau_init))
+
+
+