share/hompack/fortran/rho.f

   1       SUBROUTINE RHO(A,LAMBDA,X,V,PAR,IPAR)
   2       DOUBLE PRECISION A(*),LAMBDA,X(*),V(*),PAR(*)
   3       INTEGER IPAR(*)
   4 C
   5 C PAR(1:*) AND IPAR(1:*) ARE ARRAYS FOR (OPTIONAL) USER PARAMETERS,
   6 C    WHICH ARE SIMPLY PASSED THROUGH TO THE USER WRITTEN SUBROUTINES
   7 C    RHO, RHOJAC.
   8 C
   9 C EVALUATE  RHO(A,LAMBDA,X)  AND RETURN IN THE VECTOR  V .
  10 C
  11 C THE FOLLOWING CODE IS SPECIFICALLY FOR THE POLYNOMIAL SYSTEM DRIVER
  12 C  POLSYS , AND SHOULD BE USED VERBATUM WITH  POLSYS .  IF THE USER IS
  13 C CALLING  FIXP??  OR   STEP??  DIRECTLY, HE MUST SUPPLY APPROPRIATE
  14 C REPLACEMENT CODE HERE.
  15       INTEGER J,N
  16 C FORCE PREDICTED POINT TO HAVE  LAMBDA .GE. 0  .
  17       IF (LAMBDA .LT. 0.0) LAMBDA=0.0
  18       CALL HFUNP(A,LAMBDA,X,PAR,IPAR)
  19       N=IPAR(1)
  20       DO 10 J=1,2*N
  21         V(J)=PAR(IPAR(3 + (4-1)) + (J-1))
  22    10 CONTINUE
  23 C
  24       RETURN
  25       END