Télécharger la plaquette Cast3M


fr en

Télécharger mamu45.eso

Retour à la liste

Numérotation des lignes : 

mamu45

  1. C MAMU45    SOURCE    PV        20/02/22    21:15:17     10531          
  2.       subroutine mamu45(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.            id2=lpl
  11.            id3=lpl+ lpl+1
  12.            id4=lpl+ lpl+1 +lpl+2
  13.            pt1=pt(1)
  14.            pt2=pt(2)
  15.            pt3=pt(3)
  16.            pt4=pt(4)
  17.            pt5=pt(5)
  18.            pt7=pt(6)
  19.            pt8=pt(7)
  20.            pt9=pt(8)
  21.            pt10=pt(9)
  22.            pt11=pt(10)
  23.            pt13=pt(11)
  24.            pt14=pt(12)
  25.            pt15=pt(13)
  26.            pt16=pt(14)
  27.            pt17=pt(15)
  28.            pt19=pt(16)
  29.            pt20=pt(17)
  30.            pt21=pt(18)
  31.            pt22=pt(19)
  32.            pt23=pt(20)
  33.            do 110 i=1,long
  34.               xval1=val1(i)
  35.               xval2=val1(i+id2)
  36.               xval3=val1(i+id3)
  37.               xval4=val1(i+id4)
  38.               xval11=val11(i)
  39.               xval12=val11(i+id12)
  40.               xval13=val11(i+id13)
  41.               xval14=val11(i+id14)
  42.               xval15=val11(i+id15)
  43.               pt1=pt1+xval1*xval11
  44.               pt2=pt2+xval1*xval12
  45.               pt3=pt3+xval1*xval13
  46.               pt4=pt4+xval1*xval14
  47.               pt5=pt5+xval1*xval15
  48.               pt7=pt7+xval2*xval11
  49.               pt8=pt8+xval2*xval12
  50.               pt9=pt9+xval2*xval13
  51.               pt10=pt10+xval2*xval14
  52.               pt11=pt11+xval2*xval15
  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.               pt19=pt19+xval4*xval11
  59.               pt20=pt20+xval4*xval12
  60.               pt21=pt21+xval4*xval13
  61.               pt22=pt22+xval4*xval14
  62.               pt23=pt23+xval4*xval15
  63.  110       continue
  64.            pt(1)=pt1
  65.            pt(2)=pt2
  66.            pt(3)=pt3
  67.            pt(4)=pt4
  68.            pt(5)=pt5
  69.            pt(6)=pt7
  70.            pt(7)=pt8
  71.            pt(8)=pt9
  72.            pt(9)=pt10
  73.            pt(10)=pt11
  74.            pt(11)=pt13
  75.            pt(12)=pt14
  76.            pt(13)=pt15
  77.            pt(14)=pt16
  78.            pt(15)=pt17
  79.            pt(16)=pt19
  80.            pt(17)=pt20
  81.            pt(18)=pt21
  82.            pt(19)=pt22
  83.            pt(20)=pt23
  84.  
  85.            end
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  

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