From 7dfcc480ba1e19bd3232349fc733caef94034292 Mon Sep 17 00:00:00 2001
From: stainer_t <thomas.stainer@oecd-nea.org>
Date: Mon, 8 Sep 2025 13:48:49 +0200
Subject: Initial commit from Polytechnique Montreal

---
 PyGan/data/Equivalence_SPH_APEX.py | 29 +++++++++++++++++++++++++++++
 1 file changed, 29 insertions(+)
 create mode 100644 PyGan/data/Equivalence_SPH_APEX.py

(limited to 'PyGan/data/Equivalence_SPH_APEX.py')

diff --git a/PyGan/data/Equivalence_SPH_APEX.py b/PyGan/data/Equivalence_SPH_APEX.py
new file mode 100644
index 0000000..0ddcf03
--- /dev/null
+++ b/PyGan/data/Equivalence_SPH_APEX.py
@@ -0,0 +1,29 @@
+#
+# Equivalence_SPH_APEX: incorporation of SPH factors in an APEX file
+#
+import lifo
+import lcm
+import cle2000
+from assertS import *
+import numpy as np
+
+# recover the macro-geometry
+ipLifo=lifo.new()
+ipLifo.pushEmpty("Geom", "LCM")
+ipLifo.push('UOX_5x5_TG6_sym8_multiDom.h5')
+ipLifo.push(int) # ncals
+GeomProc = cle2000.new('Umzone5x5',ipLifo,1)
+GeomProc.exec()
+ipGeom = ipLifo.node("Geom")
+ncals = ipLifo.node(2)
+print("test Equivalence_SPH_APEX: number of calculations=",ncals)
+ipGeom.lib()
+
+# perform transport-diffusion SPH equivalence
+SPH_RT_Proc = cle2000.new('SPH_RT_Proc',ipLifo,1)
+SPH_RT_Proc.exec()
+
+# perform transport-SN SPH equivalence
+SPH_RT_Proc = cle2000.new('SPH_SN_Proc',ipLifo,1)
+SPH_RT_Proc.exec()
+print("test Equivalence_SPH_APEX completed")
-- 
cgit v1.2.3