Unverified Commit 6a973745 authored by rossial's avatar rossial Committed by GitHub
Browse files

Delete tnormi.for

parent ba81adc9
C ----------------------------------------------------------------------
C TNORMI generates a truncated normal random number
C through the inversion method
C Developed by A.Rossi, C.Planas and G.Fiorentini
C
C Copyright (C) 2010-2014 European Commission
C
C This file is part of Program DMM
C
C DMM is free software developed at the Joint Research Centre of the
C European Commission: you can redistribute it and/or modify it under
C the terms of the GNU General Public License as published by
C the Free Software Foundation, either version 3 of the License, or
C (at your option) any later version.
C
C DMM is distributed in the hope that it will be useful,
C but WITHOUT ANY WARRANTY; without even the implied warranty of
C MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
C GNU General Public License for more details.
C
C You should have received a copy of the GNU General Public License
C along with DMM. If not, see <http://www.gnu.org/licenses/>.
C ----------------------------------------------------------------------
DOUBLE PRECISION FUNCTION TNORMI(PHIL,PHIP)
C INPUT
DOUBLE PRECISION PHIL,PHIP,T
C EXTERNAL FUNCTIONS
DOUBLE PRECISION genunf,INVNORMCDF !PPND16
T = PHIL+genunf(0.d0,1.d0)*(PHIP-PHIL) ! Rescaling U(PHIL,PHIP)
C TNORMI = G01FAF('L',T,IFAIL) ! INVERSE of N(0,1)
C TNORMI = PPND16(T,IFAIL)
TNORMI = INVNORMCDF(T)
RETURN
END
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment