「10日くらいでできる!プログラミング言語自作入門」の続編#1-5

  • (by K, 2021.04.22)

(1) HL-15

  • 今回の最適化は、a = b + c - d; みたいな式に関するものです。HL-14aでは、こうなります。
    >codedump 1
    
    (len=0) 
    
    >a=b+c-d
    8b 05 58 0b 42 00 03 05 5c 0b 42 00 89 05 90 0a 42 00 8b 05 90 0a 42 00 2b 05 60 0b 42 00 89 05 94 0a 42 00 8b 05 94 0a 42 00 89 05 54 0b 42 00
    (len=48)
  • これを見やすくするとこうなっています。
    8b 05 58 0b 42 00; // EAX = b;
    03 05 5c 0b 42 00; // EAX += c;
    89 05 90 0a 42 00; // _t0 = EAX;
    8b 05 90 0a 42 00; // EAX = _t0;
    2b 05 60 0b 42 00; // EAX -= d;
    89 05 94 0a 42 00; // _t1 = EAX;
    8b 05 94 0a 42 00; // EAX = _t1;
    89 05 54 0b 42 00; // a = EAX;
  • はい、これはひどいです。_t0や_t1に書いたり読み込んだり、そんなの全然必要ないことです。掛け算や割り算など、異なる優先順位の演算がまざってもいいように問答無用で一時変数に代入するようにしていることが、この無駄の原因になっています。
  • これに対して、理想形はこうなります。
    8b 05 58 0b 42 00; // EAX = b;
    03 05 5c 0b 42 00; // EAX += c;
    2b 05 60 0b 42 00; // EAX -= d;
    89 05 54 0b 42 00; // a = EAX;
  • っていうか、普通こうなるべきですよね。

  • もう一つ気になっているのは、このコンパイル結果です。
    >if (a+b<c-d) { x=0; }
    
        8b 05 54 0b 42 00; // EAX = a;
        03 05 58 0b 42 00; // EAX += b;
        89 05 90 0a 42 00; // _t0 = EAX;
        8b 05 5c 0b 42 00; // EAX = c;
        2b 05 60 0b 42 00; // EAX -= d;
        89 05 94 0a 42 00; // _t1 = EAX;
        8b 05 90 0a 42 00; // EAX = _t0;
        3b 05 94 0a 42 00; // if (EAX < _t1) { EAX = 1: } (SETL)
        0f 9c c0;
        0f b6 c0;
        89 05 98 0a 42 00; // _t2 = EAX;
        8b 05 98 0a 42 00; // EAX = _t2;
        85 c0;             // if (EAX == 0) goto skip;
        0f 84 0c 00 00 00;
        8b 05 0c 0a 42 00; // EAX = 0;
        89 05 64 0b 42 00; // x = EAX;
    skip:
  • これがだめだなあと思うのは、_t2への無駄な読み書きと、0f_9c_c0;命令で比較結果を計算した後で、またその結果に対して85_c0;で比較していることです。理想はこうです。
        8b 05 54 0b 42 00; // EAX = a;
        03 05 58 0b 42 00; // EAX += b;
        89 05 90 0a 42 00; // _t0 = EAX;
        8b 05 5c 0b 42 00; // EAX = c;
        2b 05 60 0b 42 00; // EAX -= d;
        89 05 94 0a 42 00; // _t1 = EAX;
        8b 05 90 0a 42 00; // EAX = _t0;
        3b 05 94 0a 42 00; // if (EAX >= _t1) goto skip;
        0f 8d 0c 00 00 00;
        8b 05 0c 0a 42 00; // EAX = 0;
        89 05 64 0b 42 00; // x = EAX;
    skip:

  • ということで、HL-15はこの2つの問題を解決しようと思います。
  • [1]regVar()関数の後、putIcX86_sub()関数の前に以下を追加
    unsigned char *icq0, *icq1, *icqSet;
    
    void optimizerX86()
    {
        int i;
        if (icq0 != icq) {
            if (icq0[0] == 0x0f && 0x90 <= icq0[1] && icq0[1] <= 0x9f) {	// SETcc
                icqSet = icq0;
            }
            if (icq0[0] == 0x8b && icq1 != 0 && icq1[0] == 0x89 && icq0[1] == 0x05 && icq1[1] == 0x05 && get32(icq0 + 2) == get32(icq1 + 2)) {
                icq = icq0;	// 	8b命令は削除.
                i = (IntP) get32(icq1 + 2) - var;
                if (TcTmp0 <= i && i <= TcTmp9) {
                    icq = icq1;	// 89命令も削除.
                }
            }
            icq1 = icq0;
            icq0 = icq;
        }
        if (icqSet + 14 == icq && memcmp(&icqSet[2], "\xc0\x0f\xb6\xc0\x85\xc0\x0f", 7) == 0 && 0x84 <= icqSet[9] && icqSet[9] <= 0x85) {
            memcpy(&icqSet[2], &icqSet[10], 4);
            icqSet[1] -= 0x10;	// SETcc → Jcc.
            if (icqSet[9] == 0x84) {
                icqSet[1] ^= 1; // 条件反転.
            }
            icq0 = icq = icqSet + 6;
            icq1 = icqSet = 0;
            jmps[jp - 1] = icq - 4 - ic;
        }
    }
  • [2]putIcX86_sub()関数を一部改造
    void putIcX86_sub(String s, IntP a[])
    {
        int i, j, k;
        unsigned char addVal = 0, r = 0;
        for (i = 0; s[i] != 0; ) {
    !       if (s[i] == ' ' || s[i] == '\t' || s[i] == '_' || s[i] == ':') {
                i++;
    +       } else if (s[i] == ';') {
    +           i++;
    +           optimizerX86();
           } else if (getHex(s[i]) >= 0 && getHex(s[i + 1]) >= 0) { // 16進数2桁.
           (中略)
    }
  • [3]defLabel()関数に一行追加
    void defLabel(int i)
    {
        var[i] = icq - ic;
    +   icq1 = icqSet = 0;
    }
  • compile()関数を少し改造
    int compile(String s)
    {
        int pc, pc1, i, j;
        unsigned char *icq1, *icp;
        pc1 = lexer(s, tc);
        tc[pc1++] = TcSemi;	// 末尾に「;」を付け忘れることが多いので、付けてあげる.
        tc[pc1] = tc[pc1 + 1] = tc[pc1 + 2] = tc[pc1 + 3] = TcDot;	// エラー表示用のために末尾にピリオドを登録しておく.
    !   icq = icq0 = ic;
        jp = 0;
    +   icq1 = icqSet = 0;
        putIcX86("60; 83_ec_7c;", 0, 0, 0, 0); // PUSHAD(); SUB(ESP,124);
        regVarSaveLoad(RvLoad);
        dump0 = icq;
        (中略)
    }

  • 以上すべての改造を終えると、プログラムは942行になります。
  • さて、うまく動くでしょうか。
    >a=b+c-d
    8b 05 30 0b 42 00 03 05 34 0b 42 00 2b 05 38 0b 42 00 89 05 2c 0b 42 00
    (len=24)
  • うまくいっています。ちゃんと上記の理想通りになっています。
  • if文のほうはどうでしょうか。
    >if (a+b<c-d) { x=0; }
    8b 05 2c 0b 42 00 03 05 30 0b 42 00 89 05 a0 0a 42 00 8b 05 34 0b 42 00 2b 05 38 0b 42 00 89 05 a4 0a 42 00
    8b 05 a0 0a 42 00 3b 05 a4 0a 42 00 0f 9c c0 83 e0 01 85 c0 0f 84 0c 00 00 00 8b 05 1c 0a 42 00
    89 05 60 0b 42 00
    (len=74)
  • これもうまくいっています!

(2) プログラムの説明

  • 今回の改造の中心は、optimizerX86()関数です。
  • この関数は、putIcX86()でセミコロンが来たタイミングで呼ばれます。
  • まず、3つのグローバル変数があります。
    • icq0 : 直前のセミコロンの位置(=今の命令の開始位置)
    • icq1 : さらにその前のセミコロンの位置(=1つ前の命令の開始位置)
    • icqSet : SETcc命令を見つけたら、その先頭を記憶しておく
  • icqとicq0とicq1を使えば、連続した2命令をチェックすることができます。これを使って、直前に書いたところをすぐに読み込む場合があったら、その読み込み命令を消します。また、書き込み対象が一時変数なら、書き込み命令も消します。
  • またSETcc命令の後に最適化すべきパターンが来ていたら、前述の理想形になるように、命令列を加工しています。また条件分岐命令の位置が変わるので、 jmps[jp - 1] を上書きします。
  • 最適化対象の機械語フレーズの中にラベルがあると、どこからかそこに分岐して入ってくる可能性があるので、勝手な最適化はできません。ということで、defLabel()では最適化を妨げるためにicq1とicqSetを0クリアしています。

(3) 成果の比較

  • 上記みたいなのはすごく特別でレアケースなのか、それともよくあるメジャーケースなのか、たぶん簡単には判断できないと思うので、いくつかのプログラムでcodedumpしてバイト数を比較してみました。
  • サイズが小さくなったからいいとか悪いとかではなく、「サイズが小さくなった=それだけ適用された」ということだと考えてください。
HL-13aHL-14HL-14aHL-15
mandel.c10871007942768
maze.c2192219221921770
kcube.c4623462346233695
invader.c3260326032602750

次回に続く

こめんと欄


コメントお名前NameLink

トップ   編集 凍結 差分 バックアップ 添付 複製 名前変更 リロード   新規 一覧 単語検索 最終更新   ヘルプ   最終更新のRSS
Last-modified: 2021-04-26 (月) 22:51:33 (180d)