Télécharger versym.eso

Retour à la liste

Numérotation des lignes :

  1. C VERSYM SOURCE CB215821 19/03/18 21:16:32 10160
  2. *
  3. * verification symetrie de la matrice
  4. *
  5. subroutine versym(re,ire1,ire2,ire3)
  6. implicit real*8 (a-h,o-z)
  7. -INC CCOPTIO
  8. -INC CCREEL
  9. real*8 re(ire1,ire2,ire3)
  10. * facteur de blocage pour optimiser le cache
  11. PARAMETER (NINC=32)
  12. *
  13. if (ire1.ne.ire2) call erreur(1044)
  14. xzref=(xpetit/xzprec)
  15. xzprla=xzprec*5.5D0
  16. mdim=ire1
  17. do iel=1,ire3
  18. do i=1,mdim,ninc
  19. do jr=1,mdim
  20. do ir=i,min(jr-1,i+ninc-1)
  21. re1=re(ir,jr,iel)
  22. re2=re(jr,ir,iel)
  23. xn=abs(re1-re2)
  24. xd=abs((re1+re2) / 2.D0)+xzref
  25. if (xn.gt.xd*xzprla) then
  26. reaerr(1)=re1
  27. reaerr(2)=re2
  28. reaerr(3)=xn
  29. call erreur(1044)
  30. return
  31. endif
  32. enddo
  33. enddo
  34. enddo
  35. enddo
  36. end
  37.  
  38.  

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