* * $Id: delbnd.F,v 1.1.1.1 1996/04/01 15:03:27 mclareni Exp $ * * $Log: delbnd.F,v $ * Revision 1.1.1.1 1996/04/01 15:03:27 mclareni * Mathlib gen * * #include "gen/pilot.h" SUBROUTINE DELBND (N,NACTV,IACTV,NFREE,IFREE,IBDEL,IBTRUE,GACTV,NA 1DIM,AHESS,GFREE,GNORM) INTEGER N, NACTV, NFREE, IBDEL, IBTRUE, NADIM INTEGER IACTV(N), IFREE(N) DOUBLE PRECISION GNORM DOUBLE PRECISION AHESS(NADIM, N), GFREE(N), GACTV(N) INTEGER IPOS, J DOUBLE PRECISION RVAL CALL DELETE(NACTV,IACTV,IBTRUE,IPOS) NACTV=NACTV-1 IBTRUE=ABS(IBTRUE) IFREE(NFREE+1)=IBTRUE IF(NFREE.EQ.0) GOTO 20 DO 10 J=1,NFREE AHESS(NFREE+1,J)=0.0D+0 10 CONTINUE 20 NFREE=NFREE+1 AHESS(NFREE,NFREE)=1.0D+0 GFREE(NFREE)=GACTV(IBDEL) RVAL=GNORM*GNORM+GFREE(NFREE)**2 GNORM=SQRT(RVAL) RETURN END