diff options
| author | stainer_t <thomas.stainer@oecd-nea.org> | 2025-09-08 13:48:49 +0200 |
|---|---|---|
| committer | stainer_t <thomas.stainer@oecd-nea.org> | 2025-09-08 13:48:49 +0200 |
| commit | 7dfcc480ba1e19bd3232349fc733caef94034292 (patch) | |
| tree | 03ee104eb8846d5cc1a981d267687a729185d3f3 /Donjon/src/DETPAR.f | |
Initial commit from Polytechnique Montreal
Diffstat (limited to 'Donjon/src/DETPAR.f')
| -rw-r--r-- | Donjon/src/DETPAR.f | 31 |
1 files changed, 31 insertions, 0 deletions
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 |