Télécharger la plaquette Cast3M


fr en

Télécharger mamu63.eso

Retour à la liste

Numérotation des lignes : 

mamu63

  1. C MAMU63    SOURCE    PV        20/02/22    21:15:23     10531          
  2.       subroutine mamu63(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.            id2=lpl
  9.            id3=lpl+ lpl+1
  10.            id4=lpl+ lpl+1 +lpl+2
  11.            id5=lpl+ lpl+1 +lpl+2 +lpl+3
  12.            id6=lpl+ lpl+1 +lpl+2 +lpl+3 +lpl+4
  13.            pt1=pt(1)
  14.            pt2=pt(2)
  15.            pt3=pt(3)
  16.            pt7=pt(4)
  17.            pt8=pt(5)
  18.            pt9=pt(6)
  19.            pt13=pt(7)
  20.            pt14=pt(8)
  21.            pt15=pt(9)
  22.            pt19=pt(10)
  23.            pt20=pt(11)
  24.            pt21=pt(12)
  25.            pt25=pt(13)
  26.            pt26=pt(14)
  27.            pt27=pt(15)
  28.            pt31=pt(16)
  29.            pt32=pt(17)
  30.            pt33=pt(18)
  31.            do 110 i=1,long
  32.               xval1=val1(i)
  33.               xval2=val1(i+id2)
  34.               xval3=val1(i+id3)
  35.               xval4=val1(i+id4)
  36.               xval5=val1(i+id5)
  37.               xval6=val1(i+id6)
  38.               xval11=val11(i)
  39.               xval12=val11(i+id12)
  40.               xval13=val11(i+id13)
  41.               pt1=pt1+xval1*xval11
  42.               pt2=pt2+xval1*xval12
  43.               pt3=pt3+xval1*xval13
  44.               pt7=pt7+xval2*xval11
  45.               pt8=pt8+xval2*xval12
  46.               pt9=pt9+xval2*xval13
  47.               pt13=pt13+xval3*xval11
  48.               pt14=pt14+xval3*xval12
  49.               pt15=pt15+xval3*xval13
  50.               pt19=pt19+xval4*xval11
  51.               pt20=pt20+xval4*xval12
  52.               pt21=pt21+xval4*xval13
  53.               pt25=pt25+xval5*xval11
  54.               pt26=pt26+xval5*xval12
  55.               pt27=pt27+xval5*xval13
  56.               pt31=pt31+xval6*xval11
  57.               pt32=pt32+xval6*xval12
  58.               pt33=pt33+xval6*xval13
  59.  110       continue
  60.            pt(1)=pt1
  61.            pt(2)=pt2
  62.            pt(3)=pt3
  63.            pt(4)=pt7
  64.            pt(5)=pt8
  65.            pt(6)=pt9
  66.            pt(7)=pt13
  67.            pt(8)=pt14
  68.            pt(9)=pt15
  69.            pt(10)=pt19
  70.            pt(11)=pt20
  71.            pt(12)=pt21
  72.            pt(13)=pt25
  73.            pt(14)=pt26
  74.            pt(15)=pt27
  75.            pt(16)=pt31
  76.            pt(17)=pt32
  77.            pt(18)=pt33
  78.  
  79.            end
  80.  
  81.  
  82.  
  83.  
  84.  
  85.  
  86.  
  87.  

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