| fe4da5cc |
1 | * |
| 2 | * $Id$ |
| 3 | * |
| 4 | * $Log$ |
| 5 | * Revision 1.1.1.1 1996/04/01 15:03:25 mclareni |
| 6 | * Mathlib gen |
| 7 | * |
| 8 | * |
| 9 | #include "gen/pilot.h" |
| 10 | SUBROUTINE TSTEXT (N,XEXTR,XLOW,XUP,EDGFAC,FRACT,NCUT,NCDIM,ICUT,D |
| 11 | 1ELPLS,DELNEG) |
| 12 | INTEGER N, NCDIM, NCUT |
| 13 | INTEGER ICUT(NCDIM) |
| 14 | DOUBLE PRECISION EDGFAC, FRACT |
| 15 | DOUBLE PRECISION XEXTR(N), XLOW(N), XUP(N), DELPLS(N), DELNEG(N) |
| 16 | INTEGER I |
| 17 | DOUBLE PRECISION DIFX |
| 18 | NCUT=0 |
| 19 | DO 10 I=1,N |
| 20 | DIFX=XUP(I)-XEXTR(I) |
| 21 | DELPLS(I)=DIFX |
| 22 | IF(DIFX.LT.EDGFAC*(XUP(I)-XLOW(I))) GOTO 10 |
| 23 | DELPLS(I)=FRACT*DIFX |
| 24 | NCUT=NCUT+1 |
| 25 | ICUT(NCUT)=I |
| 26 | 10 CONTINUE |
| 27 | DO 20 I=1,N |
| 28 | DIFX=XEXTR(I)-XLOW(I) |
| 29 | DELNEG(I)=DIFX |
| 30 | IF(DIFX.LT.EDGFAC*(XUP(I)-XLOW(I))) GOTO 20 |
| 31 | DELNEG(I)=FRACT*DIFX |
| 32 | NCUT=NCUT+1 |
| 33 | ICUT(NCUT)=-I |
| 34 | 20 CONTINUE |
| 35 | RETURN |
| 36 | END |
git://git.uio.no / u / mrichter / AliRoot.git / blame