]>
Commit | Line | Data |
---|---|---|
ef94df36 | 1 | *CMZ : 23/08/93 13.30.11 by Jonathan Butterworth |
2 | *-- Author : | |
3 | SUBROUTINE DELBND (N,NACTV,IACTV,NFREE,IFREE,IBDEL,IBTRUE,GACTV,NA | |
4 | 1DIM,AHESS,GFREE,GNORM) | |
5 | INTEGER N, NACTV, NFREE, IBDEL, IBTRUE, NADIM | |
6 | INTEGER IACTV(N), IFREE(N) | |
7 | DOUBLE PRECISION GNORM | |
8 | DOUBLE PRECISION AHESS(NADIM, N), GFREE(N), GACTV(N) | |
9 | INTEGER IPOS, J | |
10 | DOUBLE PRECISION RVAL | |
11 | SAVE | |
12 | CALL DELETE(NACTV,IACTV,IBTRUE,IPOS) | |
13 | NACTV=NACTV-1 | |
14 | IBTRUE=ABS(IBTRUE) | |
15 | IFREE(NFREE+1)=IBTRUE | |
16 | IF(NFREE.EQ.0) GOTO 20 | |
17 | DO 10 J=1,NFREE | |
18 | AHESS(NFREE+1,J)=0.0D+0 | |
19 | 10 CONTINUE | |
20 | 20 NFREE=NFREE+1 | |
21 | AHESS(NFREE,NFREE)=1.0D+0 | |
22 | GFREE(NFREE)=GACTV(IBDEL) | |
23 | RVAL=GNORM*GNORM+GFREE(NFREE)**2 | |
24 | GNORM=SQRT(RVAL) | |
25 | RETURN | |
26 | END |