Extract substring of Fortran string array

后端 未结 3 1172
感动是毒
感动是毒 2020-12-19 18:33

How does one extract a substring of a Fortran string array? For example

program testcharindex
    implicit none
    character(len=10), dimension(5) :: s
             


        
相关标签:
3条回答
  • 2020-12-19 18:55

    You have a couple of problems here. One of which is easily addressed (and has been in other questions: you can find these for more detail).

    The line1

    n = s(1:i-1)
    

    about which the compiler is complaining is an attempt to reference a section of the array s, not an array of substrings of elements of array s. To access the substrings of the array you will need

    n = s(:)(1:i-1)
    

    However, this is related to your second problem. As the compiler complains for accessing the array section, the i must be a scalar. This is also true for the case of accessing substrings of an array. The above line will still not work.

    Essentially, if you wish to access substrings of an array, each substring has to have exactly the same structure. That is, in s(:)(i:j) both i and j must be scalar integer expressions. This is motived by the desire to have every element of the returned array being the same length.

    You will, then, need to use a loop.


    1 As High Performance Mark once commented, there's also a problem with the assignment itself. I considered simply the expression on the right-hand side. Even corrected for a valid array substring, the expression is still a character array, which cannot be assigned to the integer scalar n as desired.

    If you want the literal answer about selecting substrings then read as above. If you simply care about "converting part of a character array to an integer array" then another answer covers things well.

    0 讨论(0)
  • 2020-12-19 19:07

    In case loop is to be avoided and there is no other (simple) methods, it might be useful to define an elemental substring function and apply it to an array of strings. For example,

    module str_mod
        implicit none
    contains
        elemental function substr( s, a, b ) result( res )
            character(*), intent(in) :: s
            integer,      intent(in) :: a, b
            character(len(s)) :: res
    
            res = s( a : b )
        endfunction
    endmodule
    
    program main
        use str_mod
        implicit none
        character(10) :: s( 5 )
        integer, allocatable :: ind(:)
        character(len(s)), allocatable :: comp(:)
    
        s = [ '1_E ', '2_S ', '3_E ', '14_E', '25_S' ]
        ! s = [ character(len(s)) :: '1_E', '2_S', '3_E', '14_E', '25_S' ]
    
        print *, "test(scalar) : ", substr( s(1), 1, 2 )
        print *, "test(array ) : ", substr( s,    1, 2 )
    
        ind = index( s, '_' )
        comp = substr( s, 1, ind-1 )
    
        print *
        print *, "string (all)    : ", s
        print *, "string before _ : ", comp
        print *, "string after _  : ", substr( s, ind+1, len(s) )
    endprogram
    

    which gives (with gfortran-7.3)

     test(scalar) : 1_        
     test(array ) : 1_        2_        3_        14        25        
    
     string (all)    : 1_E       2_S       3_E       14_E      25_S      
     string before _ : 1         2         3         14        25        
     string after _  : E         S         E         E         S     
    
    0 讨论(0)
  • 2020-12-19 19:10

    @francescalus has already explained the error, here's my contribution towards the question OP seems really to be tackling, ie how to read the integers from a string array such as

    s = (/ '1_E ', '2_S ', '3_E ', '14_E', '25_S' /)
    

    OP wants to do this without loops, and @roygvib points us towards using an elemental function. Here's my version of such a function to read the integer from a string. This ignores any leading spaces, so should cope with strings such as 12_e. It then stops scanning at the first non-digit character (so reads 12 from a string such as 12_3).

    ELEMENTAL INTEGER FUNCTION read_int(str)
      CHARACTER(*), INTENT(in) :: str
      CHARACTER(:), ALLOCATABLE :: instr
    
      instr = adjustl(str)
      instr = instr(1:VERIFY(instr,'0123456789')-1)
      ! if the string doesn't have a leading digit instr will be empty, return a guard value
      IF(instr=='') instr = '-999'
      READ(instr,*) read_int
    END FUNCTION read_int
    

    I trust that this is clear enough. OP could then write

    n = read_int(s)
    
    0 讨论(0)
提交回复
热议问题