/scilab-master-1333395999/modules/scicos_blocks/src/fortran/logblk.f

# · FORTRAN Legacy · 53 lines · 23 code · 5 blank · 25 comment · 4 complexity · 6d051a94b612edd413b260897a1193b6 MD5 · raw file 

    

  1. c     Scicos
    2. 

  2. c     
    3. 

  3. c     Copyright (C) INRIA - METALAU Project <scicos@inria.fr>
    4. 

  4. c     
    5. 

  5. c     This program is free software; you can redistribute it and/or modify
    6. 

  6. c     it under the terms of the GNU General Public License as published by
    7. 

  7. c     the Free Software Foundation; either version 2 of the License, or
    8. 

  8. c     (at your option) any later version.
    9. 

  9. c     
    10. 

  10. c     This program is distributed in the hope that it will be useful,
    11. 

  11. c     but WITHOUT ANY WARRANTY; without even the implied warranty of
    12. 

  12. c     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
    13. 

  13. c     GNU General Public License for more details.
    14. 

  14. c     
    15. 

  15. c     You should have received a copy of the GNU General Public License
    16. 

  16. c     along with this program; if not, write to the Free Software
    17. 

  17. c     Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
    18. 

  18. c     
    19. 

  19. c     See the file ./license.txt
    20. 

  20. c     
    21. 

  21. 

  22.       subroutine logblk(flag,nevprt,t,xd,x,nx,z,nz,tvec,ntvec,
    23. 

  23.      &     rpar,nrpar,ipar,nipar,u,nu,y,ny)
    24. 

  24. c     Copyright INRIA
    25. 

  25. 

  26. 

  27. c     Scicos block simulator
    28. 

  28. c     y=log(u)/log(rpar(1))
    29. 

  29. c
    30. 

  30.       double precision t,xd(*),x(*),z(*),tvec(*),rpar(*),u(*),y(*)
    31. 

  31.       integer flag,nevprt,nx,nz,ntvec,nrpar,ipar(*)
    32. 

  32.       integer nipar,nu,ny
    33. 

  33. 

  34. c
    35. 

  35.       if(flag.eq.1) then
    36. 

  36.          do 15 i=1,nu
    37. 

  37.             if(u(i).gt.0.0d0) then
    38. 

  38.                y(i)=log(u(i))/log(rpar(1))
    39. 

  39.             else
    40. 

  40.                flag=-2
    41. 

  41.                return
    42. 

  42.             endif
    43. 

  43.  15      continue
    44. 

  44.       endif
    45. 

  45.       if(flag.eq.6) then
    46. 

  46.          do 20 i=1,nu
    47. 

  47.             if(u(i).gt.0.0d0) then
    48. 

  48.                y(i)=log(u(i))/log(rpar(1))
    49. 

  49.             endif
    50. 

  50.  20      continue
    51. 

  51.       endif
    52. 

  52. 

  53.       end
    54.