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

---
 Donjon/src/DETSPLI2.f | 31 +++++++++++++++++++++++++++++++
 1 file changed, 31 insertions(+)
 create mode 100644 Donjon/src/DETSPLI2.f

(limited to 'Donjon/src/DETSPLI2.f')

diff --git a/Donjon/src/DETSPLI2.f b/Donjon/src/DETSPLI2.f
new file mode 100644
index 0000000..1d37d25
--- /dev/null
+++ b/Donjon/src/DETSPLI2.f
@@ -0,0 +1,31 @@
+*DECK DETSPLI2
+      SUBROUTINE DETSPLI2(XA,YA,Y2A,N,X,Y)
+*
+*Parameters: 
+* XA       
+* YA       
+* Y2A      
+* N        
+* X         
+* Y    
+*     
+      DIMENSION XA(N),YA(N),Y2A(N)
+      KLO=1
+      KHI=N
+1     IF (KHI-KLO.GT.1) THEN
+        K=(KHI+KLO)/2
+        IF(XA(K).GT.X)THEN
+          KHI=K
+        ELSE
+          KLO=K
+        ENDIF
+      GOTO 1
+      ENDIF
+      H=XA(KHI)-XA(KLO)
+      IF (H.EQ.0.) CALL XABORT('DETSPLI2: BAD XA INPUT.')
+      A=(XA(KHI)-X)/H
+      B=(X-XA(KLO))/H
+      Y=A*YA(KLO)+B*YA(KHI)+
+     *      ((A**3-A)*Y2A(KLO)+(B**3-B)*Y2A(KHI))*(H**2)/6.
+      RETURN
+      END
-- 
cgit v1.2.3