void MyMatrixMultiply( vFloat A[4], vFloat B[4], vFloat C[4] )

{

    vFloat A1 = vLoad( A );         // Aの行1

    vFloat A2 = vLoad( A + 1 );     // Aの行2

    vFloat A3 = vLoad( A + 2 );     // Aの行3

    vFloat A4 = vLoad( A + 3);      // Aの行4

    vFloat C1 = vZero();         // Cの行1、ゼロに初期化

    vFloat C2 = vZero();         // Cの行2、ゼロに初期化

    vFloat C3 = vZero();         // Cの行3、ゼロに初期化

    vFloat C4 = vZero();         // Cの行4、ゼロに初期化

    vFloat B1 = vLoad( B );         // Bの行1

    vFloat B2 = vLoad( B + 1);      // Bの行2

    vFloat B3 = vLoad( B + 2);      // Bの行3

    vFloat B4 = vLoad( B + 3);      // Bの行4

    // Bの第1行とAの第1列を乗算する（合計はしない）

    C1 = vMADD( vSplat( A1, 0 ), B1, C1 );

    C2 = vMADD( vSplat( A2, 0 ), B1, C2 );

    C3 = vMADD( vSplat( A3, 0 ), B1, C3 );

    C4 = vMADD( vSplat( A4, 0 ), B1, C4 );

    // Bの第2行とAの第2列を乗算し、

    // 前の結果に加算する（合計はしない）

    C1 = vMADD( vSplat( A1, 1 ), B2, C1 );

    C2 = vMADD( vSplat( A2, 1 ), B2, C2 );

    C3 = vMADD( vSplat( A3, 1 ), B2, C3 );

    C4 = vMADD( vSplat( A4, 1 ), B2, C4 );

    // Bの第3行とAの第3列を乗算し、

    // 前の結果に加算する（合計はしない）

    C1 = vMADD( vSplat( A1, 2 ), B3, C1 );

    C2 = vMADD( vSplat( A2, 2 ), B3, C2 );

    C3 = vMADD( vSplat( A3, 2 ), B3, C3 );

    C4 = vMADD( vSplat( A4, 2 ), B3, C4 );

    // Bの第4行とAの第4列を乗算し、

    // 前の結果に加算する（合計はしない）

    C1 = vMADD( vSplat( A1, 3 ), B4, C1 );

    C2 = vMADD( vSplat( A2, 3 ), B4, C2 );

    C3 = vMADD( vSplat( A3, 3 ), B4, C3 );

    C4 = vMADD( vSplat( A4, 3 ), B4, C4 );

    // 結果を書き込み先に書き込む

    vStore( C1, C );

    vStore( C2, C + 1 );

    vStore( C3, C + 2 );

    vStore( C4, C + 3 );

}