From 7dfcc480ba1e19bd3232349fc733caef94034292 Mon Sep 17 00:00:00 2001 From: stainer_t Date: Mon, 8 Sep 2025 13:48:49 +0200 Subject: Initial commit from Polytechnique Montreal --- Donjon/src/DETPAR.f | 31 +++++++++++++++++++++++++++++++ 1 file changed, 31 insertions(+) create mode 100644 Donjon/src/DETPAR.f (limited to 'Donjon/src/DETPAR.f') diff --git a/Donjon/src/DETPAR.f b/Donjon/src/DETPAR.f new file mode 100644 index 0000000..d4feb2a --- /dev/null +++ b/Donjon/src/DETPAR.f @@ -0,0 +1,31 @@ +*DECK DETPAR + SUBROUTINE DETPAR(X1,X2,X3,Y1,Y2,Y3,AS,BS,CS) +* +*---------------------------------------------------------------------- +*Purpose: +* Routine calculating the parabolic coefficients needed for +* a parabolic interpolation Y = AS*X*X + BS*X + CS +* +*Author(s): +* M. Beaudet +* +*Parameters: +* X1 +* X2 +* X3 +* Y1 +* Y2 +* Y3 +* AS +* BS +* CS +* +*---------------------------------------------------------------------- +* + ANUM = Y1*(X2-X3)+Y3*(X1-X2)+Y2*(X3-X1) + ADEN = (X1-X2)*(X1-X3)*(X2-X3) + AS = ANUM/ADEN + BS = (Y2-Y3-AS*(X2*X2-X3*X3))/(X2-X3) + CS = Y1-BS*X1-AS*X1*X1 + RETURN + END -- cgit v1.2.3