From 26d6ae5f93f23028ce988f8cada3a64c150a6932 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Beno=C3=AEt=20Sierro?= <benoit.sierro@bluewin.ch>
Date: Fri, 17 Mar 2023 13:45:45 +0100
Subject: [PATCH] test vincetti

---
 testing/test_vincetti.py | 21 +++++++++++++++++++++
 1 file changed, 21 insertions(+)
 create mode 100644 testing/test_vincetti.py

diff --git a/testing/test_vincetti.py b/testing/test_vincetti.py
new file mode 100644
index 0000000..1a1bfa5
--- /dev/null
+++ b/testing/test_vincetti.py
@@ -0,0 +1,21 @@
+import matplotlib.pyplot as plt
+import numpy as np
+
+import scgenerator as sc
+
+wl = np.linspace(200e-9, 2e-6, 2048)
+w = sc.units.m(wl)
+wl0 = 800e-9
+gas = sc.materials.Gas("argon")
+ng2 = gas.sellmeier.n_gas_2(wl, pressure=1e5)
+
+n = sc.fiber.n_eff_vincetti(wl, wl0, ng2, 1e-6, 20e-6, 5e-6, 7)
+b2 = sc.fiber.beta2(w, n)
+
+bcap = sc.capillary_dispersion(
+    wl, sc.fiber.core_radius_from_capillaries(20e-6, 5e-6, 7), "argon", pressure=1e5
+)
+
+plt.plot(wl, b2)
+plt.plot(wl, bcap)
+plt.show()