Sources Esope | Cast3M

chnsng
C CHNSNG    SOURCE    PV        20/11/14    23:36:47     10781          
*  impose par penalisation une relation dans am
*
      subroutine chnsng(am,nbnn,i,j,k)
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)
      dimension am(nbnn,nbnn)
*
      xpenl=am(i,i)*1d-3
      am(i,i)=am(i,i)-2*xpenl
      am(i,j)=am(i,j)+xpenl
      am(i,k)=am(i,k)+xpenl
      am(j,i)=am(j,i)+xpenl
      am(k,i)=am(k,i)+xpenl
      am(j,j)=am(j,j)-xpenl/2
      am(k,k)=am(k,k)-xpenl/2
      am(j,k)=am(j,k)-xpenl/2
      am(k,j)=am(k,j)-xpenl/2
 
 
*
      end