Esope source files

   1 : C CHNSNG    SOURCE    CHAT      05/01/12    22:00:42     5004           
   2 : *  impose par penalisation une relation dans am
   3 : *
   4 :       subroutine chnsng(am,nbnn,i,j,k)
   5 :       IMPLICIT INTEGER(I-N)
   6 :       IMPLICIT REAL*8 (A-H,O-Z)
   7 :       dimension am(nbnn,nbnn)
   8 : *      
   9 :       xpenl=am(i,i)*1d-3
  10 :       am(i,i)=am(i,i)-2*xpenl
  11 :       am(i,j)=am(i,j)+xpenl
  12 :       am(i,k)=am(i,k)+xpenl
  13 : *
  14 :       end
  15 :  
  16 :