Télécharger ntap13.eso

Retour à la liste

Numérotation des lignes :

  1. C NTAP13 SOURCE CHAT 05/01/13 02:02:00 5004
  2. SUBROUTINE NTAP13 (MP,MCP,MCQ,M,N,MVDU,MVDL,KK,II)
  3. IMPLICIT INTEGER(I-N)
  4. IMPLICIT REAL*8(A-H,O-Z)
  5. *
  6. * MISE A JOUR DE LA MATRICE DE PROJECTION
  7. *
  8. -INC TMXMAT
  9. -INC SMLREEL
  10. -INC CCOPTIO
  11. MXMAT=MCP
  12. MXMA1=MCQ
  13. MXMA2=MVDU
  14. MXMA3=MVDL
  15. JG=M
  16. SEGINI MLREEL
  17. DO 1 J=1,M
  18. IF(XMAT(J,II).EQ.0.D0) THEN
  19. PROG(J)=MXMA1.XMAT(J,II)*MXMA2.XMAT(II,KK)*MXMA2.XMAT(II,KK-1)
  20. ELSE
  21. PROG(J)=XMAT(J,II)*MXMA3.XMAT(II,KK)*MXMA3.XMAT(II,KK-1)
  22. ENDIF
  23. 1 CONTINUE
  24. SEGINI MLREE1
  25. MXMAT=MP
  26. CALL MATVE1( XMAT,PROG,M,M,MLREE1.PROG,2)
  27. LDIM1=M
  28. LDIM2=M
  29. SEGINI MXMAT
  30. XX=0.D0
  31. DO 2 I=1,M
  32. XX = MLREE1.PROG(I)*MLREE1.PROG(I) + XX
  33. 2 CONTINUE
  34. CALL XITY(MLREE1.PROG,MLREE1.PROG,M,XMAT)
  35. XX=1./XX
  36. CALL REEMAT(XX,XMAT,M,M)
  37. MXMA1=MP
  38. CALL ADDIMA(MXMA1.XMAT,XMAT,M,M,2)
  39. MP=MXMA1
  40. SEGSUP MLREEL,MLREE1
  41. RETURN
  42. END
  43.  
  44.  

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