summaryrefslogtreecommitdiff
path: root/Donjon/src/DETPAR2.f
diff options
context:
space:
mode:
authorstainer_t <thomas.stainer@oecd-nea.org>2025-09-08 13:48:49 +0200
committerstainer_t <thomas.stainer@oecd-nea.org>2025-09-08 13:48:49 +0200
commit7dfcc480ba1e19bd3232349fc733caef94034292 (patch)
tree03ee104eb8846d5cc1a981d267687a729185d3f3 /Donjon/src/DETPAR2.f
Initial commit from Polytechnique Montreal
Diffstat (limited to 'Donjon/src/DETPAR2.f')
-rw-r--r--Donjon/src/DETPAR2.f31
1 files changed, 31 insertions, 0 deletions
diff --git a/Donjon/src/DETPAR2.f b/Donjon/src/DETPAR2.f
new file mode 100644
index 0000000..94f0db9
--- /dev/null
+++ b/Donjon/src/DETPAR2.f
@@ -0,0 +1,31 @@
+*DECK DETPAR2
+ SUBROUTINE DETPAR2(V1,V2,V3,U1,U2,U3,AS,BS,CS)
+*
+*----------------------------------------------------------------------
+*Purpose: routine de HQSIMEX
+*
+*Author(s):
+* M. Beaudet
+*
+*Parameters:
+* V1
+* V2
+* V3
+* U1
+* U2
+* U3
+* AS
+* BS
+* CS
+*
+*----------------------------------------------------------------------
+*
+ CHARACTER*6 CLNAME
+ CLNAME = 'PAR'
+ ANUM = U1*(V2-V3)+U3*(V1-V2)+U2*(V3-V1)
+ ADEN = (V1-V2)*(V1-V3)*(V2-V3)
+ AS = ANUM/ADEN
+ BS = (U2-U3-AS*(V2*V2-V3*V3))/(V2-V3)
+ CS = U1-BS*V1-AS*V1*V1
+ RETURN
+ END