FORTRAN(cutコマンドを作るその6)

-bオプションに対応したcutコマンドを作成しました。

takk@deb9:~$ cat -n cut6.f
     1        PROGRAM MAIN
     2          CHARACTER*100 LINE
     3          INTEGER DISP_LOC(100),N
     4          CHARACTER*100 ARGS
     5
     6          N = IARGC()
     7          DO 300 I=1,N
     8            CALL GETARG(I,ARGS)
     9
    10            IF(ARGS(1:2).EQ.'-b') THEN
    11              CALL B_OPT(ARGS(3:),DISP_LOC)
    12            ELSE
    13              OPEN(UNIT=10,FILE=ARGS(1:INDEX(ARGS,' ')-1),
    14       1      STATUS='OLD',
    15       2      ACCESS='SEQUENTIAL')
    16     10       READ(10,*,END=100) LINE
    17              CALL CUT(LINE,DISP_LOC)
    18              GOTO 10
    19    100       CLOSE(10,STATUS='KEEP')
    20            END IF
    21    300 CONTINUE
    22            STOP
    23        END
    24
    25        SUBROUTINE B_OPT(B_OPT_STRING,DISP_LOC)
    26          CHARACTER*100 B_OPT_STRING
    27          CHARACTER*100 SSTRING1(10),SSTRING2(10)
    28          INTEGER I,SSIZE1,SSIZE2,K,L,M
    29          INTEGER DISP_LOC(100)
    30
    31          DO 10000 I=1,100
    32            DISP_LOC(I) = 0
    33  10000   CONTINUE
    34
    35          CALL SPLIT(B_OPT_STRING,',', SSTRING1,SSIZE1)
    36          DO 10200 I=1,SSIZE1
    37            CALL SPLIT(SSTRING1(I),'-', SSTRING2,SSIZE2)
    38            IF(SSIZE2.EQ.1) THEN
    39              READ(SSTRING1(I),*) L
    40              DISP_LOC(L) = 1
    41            ELSE
    42              READ(SSTRING2(1),*) L
    43              READ(SSTRING2(2),*) M
    44              DO 10100 K=L,M
    45                DISP_LOC(K) = 1
    46  10100       CONTINUE
    47            END IF
    48  10200   CONTINUE
    49        END
    50
    51        SUBROUTINE SPLIT(SRC,SEP,DEST,SPLIT_SIZE)
    52          CHARACTER*100 SRC
    53          CHARACTER*1 SEP
    54          CHARACTER*100 DEST(10)
    55          INTEGER IND,S1,S2,RSLT,SPLIT_SIZE
    56          IND = 0
    57          S2=-1
    58          I = 1
    59  20100   RSLT = INDEX(SRC(IND+1:),SEP)
    60          IF(RSLT.NE.0) THEN
    61            S1 = IND+1
    62            IND = IND + RSLT
    63            S2 = IND-1
    64            DEST(I) = SRC(S1:S2)
    65            I = I + 1
    66            GOTO 20100
    67          END IF
    68          DEST(I) = SRC(S2+2:)
    69          SPLIT_SIZE = I
    70          RETURN
    71        END
    72
    73        SUBROUTINE CUT(LINE,DISP_LOC)
    74          INTEGER DISP_LOC(100)
    75          CHARACTER*100 LINE
    76          DO 30100 I=1,100
    77            IF(DISP_LOC(I).EQ.1) THEN
    78              WRITE(*,31000) LINE(I:I)
    79            END IF
    80            IF(LINE(I:I).EQ.' ') THEN
    81              WRITE(*,*) ''
    82              GOTO 30200
    83            END IF
    84  30100   CONTINUE
    85  30200   RETURN
    86  31000   FORMAT(A$)
    87        END
takk@deb9:~$

73行目のCUT関数が、列抽出する関数です。100列分配列を確認して、1がはいっていたら、部分文字列の抽出をします。空白が出現したら、文字列の終端とみなして、改行するためにWRITE(*,*) ”を実行しています。

動くでしょうか。GDBを使ってみます。

takk@deb9:~$ gfortran -g -fdollar-ok cut6.f
takk@deb9:~$ gdb a.out

パラメータを設定し、パラメータ数の取得関数前まで実行します。

(gdb) set args -b2-3,5 x00
(gdb) l
1             PROGRAM MAIN
2               CHARACTER*100 LINE
3               INTEGER DISP_LOC(100),N
4               CHARACTER*100 ARGS
5
6               N = IARGC()
7               DO 300 I=1,N
8                 CALL GETARG(I,ARGS)
9
10                IF(ARGS(1:2).EQ.'-b') THEN
(gdb) b 6
Breakpoint 1 at 0x147b: file cut6.f, line 6.
(gdb) run
Starting program: /home/takk/a.out -b2-3,5 x00

Breakpoint 1, MAIN__ () at cut6.f:6
6               N = IARGC()
(gdb)

パラメータ数は2ですので、合ってます。

(gdb) n
7               DO 300 I=1,N
(gdb) p N
$1 = 2
(gdb)

1を設定する配列の位置も、-bで指定した通りとなっています。

(gdb) n
8                 CALL GETARG(I,ARGS)
(gdb)
10                IF(ARGS(1:2).EQ.'-b') THEN
(gdb)
11                  CALL B_OPT(ARGS(3:),DISP_LOC)
(gdb)
7               DO 300 I=1,N
(gdb) p DISP_LOC
$2 = (0, 1, 1, 0, 1, 0, 0, ~
(gdb)

ファイルから読み込んだ文字列もOKです。

(gdb) n
8                 CALL GETARG(I,ARGS)
(gdb)
10                IF(ARGS(1:2).EQ.'-b') THEN
(gdb)
15           2      ACCESS='SEQUENTIAL')
(gdb)
16         10       READ(10,*,END=100) LINE
(gdb)
17                  CALL CUT(LINE,DISP_LOC)
(gdb) p line
$3 = '10010', ' ' <repeats 95 times>
(gdb)
(gdb) n
000
18                  GOTO 10
(gdb)
16         10       READ(10,*,END=100) LINE
(gdb)
17                  CALL CUT(LINE,DISP_LOC)
(gdb) s
cut (line=..., disp_loc=..., _line=100) at cut6.f:76
76              DO 30100 I=1,100
(gdb) n
77                IF(DISP_LOC(I).EQ.1) THEN
(gdb)
80                IF(LINE(I:I).EQ.' ') THEN
(gdb)
76              DO 30100 I=1,100
(gdb)
77                IF(DISP_LOC(I).EQ.1) THEN
(gdb)
78                  WRITE(*,31000) LINE(I:I)
(gdb)
080               IF(LINE(I:I).EQ.' ') THEN
(gdb)

実行結果を見てしまった方が早そうです。

takk@deb9:~$ ./a.out -b1,2,3,4 x00
1001
1011
1021
1031
1041
1051
1061
1071
1081
1091
takk@deb9:~$ ./a.out -b1-4 x00
1001
1011
1021
1031
1041
1051
1061
1071
1081
1091
takk@deb9:~$ ./a.out -b2-3,5 x00
000
010
020
030
040
050
060
070
080
090
takk@deb9:~$

上手く行きました。

Leave a Reply

Your email address will not be published. Required fields are marked *

CAPTCHA