Sources Esope | Cast3M

jslc3d
C JSLC3D    SOURCE    PV090527  23/01/27    21:15:45     11574          
      subroutine jslc3d(NBRINC,NDIMG,NINC,K,G,FRAC,JSE,ERR1,affiche)
 
c     Jacibienne d sigma / depse en espace generalisee
 
c     organisation de la Jacobienne
c     colonne : [eie1, eie2, eie3, ere1, ere2, ere3, eoe1, eoe2, eoe3, eme1, eme2, eme3]
c                 1     2     3     4     5     6     7     8     9     10    11    12
c     ligne = transpose de colonne en remplacant e par s 
c     et onajoute s 13-15 = s 10-12 et s 16-18= f* s 1-3 +(1-f)=s 10-12
      implicit real*8 (a-h,o-z)
      implicit integer (i-n)
 
      integer NBRINC,NINC,ERR1            
      real*8 JSE(NDIMG,NDIMG)
      real*8 K(0:NBRINC),G(0:NBRINC),FRAC(0:NBRINC)
      logical affiche
 
      if (NINC.ne.1) then
        print*,'Pb dans JSLC3D car NINC ne 1'
        err1=1
        return
      end if
      f=FRAC(1)
 
      t2 = K(1) + 0.4D1 / 0.3D1 * G(1)
      t4 = K(1) - 0.2D1 / 0.3D1 * G(1)
      t6 = K(0) + 0.4D1 / 0.3D1 * G(0)
      t8 = K(0) - 0.2D1 / 0.3D1 * G(0)
      t11 = K(0) / 0.2D1 - G(0) / 0.3D1
      JSE(1,1) = t2
      JSE(1,2) = t4
      JSE(1,3) = t4
      JSE(1,4) = 0.0D0
      JSE(1,5) = 0.0D0
      JSE(1,6) = 0.0D0
      JSE(1,7) = 0.0D0
      JSE(1,8) = 0.0D0
      JSE(1,9) = 0.0D0
      JSE(1,10) = 0.0D0
      JSE(1,11) = 0.0D0
      JSE(1,12) = 0.0D0
      JSE(2,1) = t4
      JSE(2,2) = t2
      JSE(2,3) = t4
      JSE(2,4) = 0.0D0
      JSE(2,5) = 0.0D0
      JSE(2,6) = 0.0D0
      JSE(2,7) = 0.0D0
      JSE(2,8) = 0.0D0
      JSE(2,9) = 0.0D0
      JSE(2,10) = 0.0D0
      JSE(2,11) = 0.0D0
      JSE(2,12) = 0.0D0
      JSE(3,1) = t4
      JSE(3,2) = t4
      JSE(3,3) = t2
      JSE(3,4) = 0.0D0
      JSE(3,5) = 0.0D0
      JSE(3,6) = 0.0D0
      JSE(3,7) = 0.0D0
      JSE(3,8) = 0.0D0
      JSE(3,9) = 0.0D0
      JSE(3,10) = 0.0D0
      JSE(3,11) = 0.0D0
      JSE(3,12) = 0.0D0
      JSE(4,1) = 0.0D0
      JSE(4,2) = 0.0D0
      JSE(4,3) = 0.0D0
      JSE(4,4) = t6
      JSE(4,5) = 0.0D0
      JSE(4,6) = 0.0D0
      JSE(4,7) = 0.0D0
      JSE(4,8) = t8
      JSE(4,9) = t8
      JSE(4,10) = 0.0D0
      JSE(4,11) = 0.0D0
      JSE(4,12) = 0.0D0
      JSE(5,1) = 0.0D0
      JSE(5,2) = 0.0D0
      JSE(5,3) = 0.0D0
      JSE(5,4) = 0.0D0
      JSE(5,5) = t6
      JSE(5,6) = 0.0D0
      JSE(5,7) = t8
      JSE(5,8) = 0.0D0
      JSE(5,9) = t8
      JSE(5,10) = 0.0D0
      JSE(5,11) = 0.0D0
      JSE(5,12) = 0.0D0
      JSE(6,1) = 0.0D0
      JSE(6,2) = 0.0D0
      JSE(6,3) = 0.0D0
      JSE(6,4) = 0.0D0
      JSE(6,5) = 0.0D0
      JSE(6,6) = t6
      JSE(6,7) = t8
      JSE(6,8) = t8
      JSE(6,9) = 0.0D0
      JSE(6,10) = 0.0D0
      JSE(6,11) = 0.0D0
      JSE(6,12) = 0.0D0
      JSE(7,1) = 0.0D0
      JSE(7,2) = 0.0D0
      JSE(7,3) = 0.0D0
      JSE(7,4) = 0.0D0
      JSE(7,5) = t11
      JSE(7,6) = t11
      JSE(7,7) = t6
      JSE(7,8) = t11
      JSE(7,9) = t11
      JSE(7,10) = 0.0D0
      JSE(7,11) = 0.0D0
      JSE(7,12) = 0.0D0
      JSE(8,1) = 0.0D0
      JSE(8,2) = 0.0D0
      JSE(8,3) = 0.0D0
      JSE(8,4) = t11
      JSE(8,5) = 0.0D0
      JSE(8,6) = t11
      JSE(8,7) = t11
      JSE(8,8) = t6
      JSE(8,9) = t11
      JSE(8,10) = 0.0D0
      JSE(8,11) = 0.0D0
      JSE(8,12) = 0.0D0
      JSE(9,1) = 0.0D0
      JSE(9,2) = 0.0D0
      JSE(9,3) = 0.0D0
      JSE(9,4) = t11
      JSE(9,5) = t11
      JSE(9,6) = 0.0D0
      JSE(9,7) = t11
      JSE(9,8) = t11
      JSE(9,9) = t6
      JSE(9,10) = 0.0D0
      JSE(9,11) = 0.0D0
      JSE(9,12) = 0.0D0
      JSE(10,1) = 0.0D0
      JSE(10,2) = 0.0D0
      JSE(10,3) = 0.0D0
      JSE(10,4) = 0.0D0
      JSE(10,5) = 0.0D0
      JSE(10,6) = 0.0D0
      JSE(10,7) = 0.0D0
      JSE(10,8) = 0.0D0
      JSE(10,9) = 0.0D0
      JSE(10,10) = t6
      JSE(10,11) = t8
      JSE(10,12) = t8
      JSE(11,1) = 0.0D0
      JSE(11,2) = 0.0D0
      JSE(11,3) = 0.0D0
      JSE(11,4) = 0.0D0
      JSE(11,5) = 0.0D0
      JSE(11,6) = 0.0D0
      JSE(11,7) = 0.0D0
      JSE(11,8) = 0.0D0
      JSE(11,9) = 0.0D0
      JSE(11,10) = t8
      JSE(11,11) = t6
      JSE(11,12) = t8
      JSE(12,1) = 0.0D0
      JSE(12,2) = 0.0D0
      JSE(12,3) = 0.0D0
      JSE(12,4) = 0.0D0
      JSE(12,5) = 0.0D0
      JSE(12,6) = 0.0D0
      JSE(12,7) = 0.0D0
      JSE(12,8) = 0.0D0
      JSE(12,9) = 0.0D0
      JSE(12,10) = t8
      JSE(12,11) = t8
      JSE(12,12) = t6
 
      if(affiche) then
        print*,'Dans JSE1inc3d'
        do i=1,12
            write(*,'(12(e10.3,1X))')
     #      (JSE(i,j),j=1,12)
        end do     
      end if
 
      return
      end