Télécharger kprior.eso

Retour à la liste

Numérotation des lignes :

  1. C KPRIOR SOURCE CHAT 05/01/13 01:04:48 5004
  2. SUBROUTINE KPRIOR(NES,NSA,NSB,A,B,U,V,IVU)
  3. IMPLICIT INTEGER(I-N)
  4. IMPLICIT REAL*8 (A-H,O-Z)
  5. C
  6. DIMENSION A(NES,NSA),B(NES,NSB),U(1),V(1)
  7. C
  8. EPS = 1E-5
  9. C
  10. LS1 = 0
  11. DO 1 IS=1,NSB
  12. F=0.
  13. DO 11 IES = 1,NES
  14. F = F+B(IES,IS)*U(IES)
  15. 11 CONTINUE
  16. F = F + U(NES+1)
  17. LS2=KSIG(F,EPS)
  18. IF (IS.GE.2.AND.LS1.NE.LS2.AND.LS2.NE.0.AND.LS1.NE.0) GOTO 4
  19. IF(LS2.NE.0) LS1=LS2
  20. 1 CONTINUE
  21. C
  22. IF(LS1.EQ.1) THEN
  23. LS1 = 0
  24. DO 2 IS=1,NSA
  25. F=0.
  26. DO 21 IES = 1,NES
  27. F = F+A(IES,IS)*V(IES)
  28. 21 CONTINUE
  29. F = F + V(NES+1)
  30. LS2=KSIG(F,EPS)
  31. IF (IS.GE.2.AND.LS1.NE.LS2.AND.LS2.NE.0.AND.LS1.NE.0) THEN
  32. IVU = 2
  33. RETURN
  34. ENDIF
  35. IF(LS2.NE.0) LS1=LS2
  36. 2 CONTINUE
  37. C
  38. IF(LS1.EQ.1) THEN
  39. IVU=1
  40. RETURN
  41. ELSE
  42. IVU = 0
  43. RETURN
  44. ENDIF
  45. C
  46. ELSE
  47. IVU = 0
  48. RETURN
  49. ENDIF
  50. C
  51. 4 CONTINUE
  52. F=0.
  53. DO 41 IES = 1,NES
  54. F = F+A(IES,1)*V(IES)
  55. 41 CONTINUE
  56. F = F + V(NES+1)
  57. LS2 = KSIG(F,EPS)
  58. IF (LS2.EQ.-1) THEN
  59. IVU = 0
  60. RETURN
  61. ELSE
  62. IVU = 3
  63. RETURN
  64. ENDIF
  65. END 00000070
  66.  
  67.  

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