Télécharger kmp3.eso

Retour à la liste

Numérotation des lignes :

  1. C KMP3 SOURCE CHAT 05/01/13 00:58:04 5004
  2. SUBROUTINE KMP3(AM,B,GX,IPADL,NUM,NP,NEL,NPT)
  3. C
  4. CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  5. C C
  6. C PARTIE NUMERIQUE INTENSE DE CAMTP2, CAS TRIDIMENSIONNEL C
  7. C CAMTP C
  8. C C
  9. C POUR OPTIMISER LES ACCES MEMOIRE ON BOUCLE SUR LES C
  10. C ELEMENTS A L'INTERIEUR DE LA BOUCLE SUR LES SOMMETS. C
  11. C C
  12. C ON AMELIORE ENCORE LA RESIDENCE DES DONNEES EN TRAITANT C
  13. C PAR PAQUETS D'ELEMENTS. LES PAQUETS NE DOIVENT PAS ETRE C
  14. C TROP GROS. C
  15. C C
  16. C C'EST LA TECHNIQUE UTILISEE DANS XCV**. C
  17. C C
  18. CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  19. C
  20. IMPLICIT INTEGER(I-N)
  21. IMPLICIT REAL*8 (A-H,O-Z)
  22. C
  23. DIMENSION AM(NEL,NP,3),GX(NPT,3),B(1)
  24. DIMENSION NUM(NP,1),IPADL(1)
  25. C
  26. LPACK=48
  27. NNN=MOD(NEL,LPACK)
  28. IF(NNN.EQ.0) NPACK=NEL/LPACK
  29. IF(NNN.NE.0) NPACK=1+(NEL-NNN)/LPACK
  30.  
  31. DO 8 IPACK=1,NPACK
  32. KDEB=(IPACK-1)*LPACK+1
  33. KFIN=KDEB+LPACK-1
  34. IF(KFIN.GT.NEL) KFIN=NEL
  35. DO 7 I=1,NP
  36. DO 6 K=KDEB,KFIN
  37. IU=IPADL(NUM(I,K))
  38. GX(IU,1)=GX(IU,1)+AM(K,I,1)*B(K)
  39. GX(IU,2)=GX(IU,2)+AM(K,I,2)*B(K)
  40. GX(IU,3)=GX(IU,3)+AM(K,I,3)*B(K)
  41. 6 CONTINUE
  42. 7 CONTINUE
  43.  
  44. 8 CONTINUE
  45. C
  46. RETURN
  47. END
  48.  
  49.  
  50.  

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