Télécharger la plaquette Cast3M


fr en

Télécharger mamu36.eso

Retour à la liste

Numérotation des lignes : 

mamu36

  1. C MAMU36    SOURCE    PV        20/02/22    21:15:13     10531          
  2.       subroutine mamu36(long,val11,lpl1,val1,lpl,pt)
  3.       implicit real*8 (a-h,o-z)
  4.       real*8 val11(*),val1(*)
  5.       real*8 pt(36)
  6.            id12=lpl1
  7.            id13=lpl1+ lpl1+1
  8.            id14=lpl1+ lpl1+1+ lpl1+2
  9.            id15=lpl1+ lpl1+1+ lpl1+2+ lpl1+3
  10.            id16=lpl1+ lpl1+1+ lpl1+2+ lpl1+3 +lpl1+4
  11.            id2=lpl
  12.            id3=lpl+ lpl+1
  13.            pt1=pt(1)
  14.            pt2=pt(2)
  15.            pt3=pt(3)
  16.            pt4=pt(4)
  17.            pt5=pt(5)
  18.            pt6=pt(6)
  19.            pt7=pt(7)
  20.            pt8=pt(8)
  21.            pt9=pt(9)
  22.            pt10=pt(10)
  23.            pt11=pt(11)
  24.            pt12=pt(12)
  25.            pt13=pt(13)
  26.            pt14=pt(14)
  27.            pt15=pt(15)
  28.            pt16=pt(16)
  29.            pt17=pt(17)
  30.            pt18=pt(18)
  31.            do 110 i=1,long
  32.               xval1=val1(i)
  33.               xval2=val1(i+id2)
  34.               xval3=val1(i+id3)
  35.               xval11=val11(i)
  36.               xval12=val11(i+id12)
  37.               xval13=val11(i+id13)
  38.               xval14=val11(i+id14)
  39.               xval15=val11(i+id15)
  40.               xval16=val11(i+id16)
  41.               pt1=pt1+xval1*xval11
  42.               pt2=pt2+xval1*xval12
  43.               pt3=pt3+xval1*xval13
  44.               pt4=pt4+xval1*xval14
  45.               pt5=pt5+xval1*xval15
  46.               pt6=pt6+xval1*xval16
  47.               pt7=pt7+xval2*xval11
  48.               pt8=pt8+xval2*xval12
  49.               pt9=pt9+xval2*xval13
  50.               pt10=pt10+xval2*xval14
  51.               pt11=pt11+xval2*xval15
  52.               pt12=pt12+xval2*xval16
  53.               pt13=pt13+xval3*xval11
  54.               pt14=pt14+xval3*xval12
  55.               pt15=pt15+xval3*xval13
  56.               pt16=pt16+xval3*xval14
  57.               pt17=pt17+xval3*xval15
  58.               pt18=pt18+xval3*xval16
  59.  110       continue
  60.            pt(1)=pt1
  61.            pt(2)=pt2
  62.            pt(3)=pt3
  63.            pt(4)=pt4
  64.            pt(5)=pt5
  65.            pt(6)=pt6
  66.            pt(7)=pt7
  67.            pt(8)=pt8
  68.            pt(9)=pt9
  69.            pt(10)=pt10
  70.            pt(11)=pt11
  71.            pt(12)=pt12
  72.            pt(13)=pt13
  73.            pt(14)=pt14
  74.            pt(15)=pt15
  75.            pt(16)=pt16
  76.            pt(17)=pt17
  77.            pt(18)=pt18
  78.  
  79.            end
  80.  
  81.  
  82.  
  83.  
  84.  
  85.  
  86.  
  87.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales