Esope source files

   1 : C CHOLE3I   SOURCE    PV        13/10/31    21:15:02     7850           
   2 :       subroutine chole3i(ithr)
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL*8 (A-H,O-Z)
   5 : -INC SMMATRI
   6 : -INC CCHOLE
   7 :       SEGMENT IMASQ(LMASQ)
   8 :       SEGMENT ITMASQ(NBLIG)
   9 :       segment immx(io)
  10 : *     write (6,*) ' ith dans chole3i ',ith
  11 :       ic=-1
  12 :       do 100 jbib=il1,il2,nbthr
  13 :       ic=mod(ic+1,nbthr)
  14 :       jbi=jbib+mod(ithr+ic-1,nbthr)
  15 :       if (jbi.gt.il2) return
  16 :       LIGN=ILIGN(JBI)
  17 :       im=imm
  18 :       if (im.gt.ider) then
  19 : *      write (6,*) ' il1 il2 ider jbib im ',il1,il2,ider,jbib,im
  20 :        return
  21 :       endif
  22 :       im=immm(1)
  23 :       do ih=2,immm(/1)
  24 :        im=min(im,immm(ih))
  25 :       enddo
  26 :       imasq=itmasq(jbi)
  27 :       do 10 ip=max(im,iper),ider
  28 :       lig1=ilign(ip)
  29 :       ippr=lig1.iprel
  30 :       iddr=lig1.iderl
  31 : 
  31 :       CALL CHOLE3(IPREL,IDERL,IPPVV(1),IPPR,IDDR,IVPO(1),
  32 :      > LIG1.IPPVV(1),VAL(1),LIG1.VAL(1),LIG1.IVPO(1),
  33 :      > imasq(1),precc)
  34 : *  en multithread il peut y avoir n'importe quoi dans oov(1) du
  35 : *  aux acces simultanés et ca crache gemat. donc :
  36 :        oov(1)=0
  37 :   10  continue
  38 :  100  continue
  39 :       end
  40 : 
  41 : 
  42 : 
  43 : 
  44 : 
  45 : 
  46 : 
  47 :  
  48 :  
  49 :  
  50 :  
  51 :  
  52 :  
  53 :  
  54 :