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 /Trivac/src/VALPL.f | |
Initial commit from Polytechnique Montreal
Diffstat (limited to 'Trivac/src/VALPL.f')
| -rwxr-xr-x | Trivac/src/VALPL.f | 35 |
1 files changed, 35 insertions, 0 deletions
diff --git a/Trivac/src/VALPL.f b/Trivac/src/VALPL.f new file mode 100755 index 0000000..b66cb78 --- /dev/null +++ b/Trivac/src/VALPL.f @@ -0,0 +1,35 @@ +*DECK VALPL + FUNCTION VALPL(L,U) +* +*----------------------------------------------------------------------- +* +*Purpose: +* Return the Legendre function coefficients for the nodal collocation +* method. +* +*Copyright: +* Copyright (C) 2002 Ecole Polytechnique de Montreal +* This library is free software; you can redistribute it and/or +* modify it under the terms of the GNU Lesser General Public +* License as published by the Free Software Foundation; either +* version 2.1 of the License, or (at your option) any later version +* +*Author(s): A. Hebert +* +*Parameters: input +* L order of the Legendre polynomial. +* U indemendent variable. +* +*---------------------------------------------------------------------- +* +*---- +* SUBROUTINE ARGUMENTS +*---- + IF (L.EQ.0) P=1.0 + IF (L.EQ.1) P=2.0*U + IF (L.EQ.2) P=(6.0*U*U-0.5) + IF (L.EQ.3) P=(20.0*U*U-3.0)*U + IF (L.EQ.4) P=(U*U*(70.0*U*U-15.0)+0.375) + VALPL=SQRT(REAL(2*L+1))*P + RETURN + END |