コグノスケ


2022年 1月 4日

Windows 電卓

Windows 電卓は浮動小数点数を扱えます。10^9999 から 10^-9999 という非常に大きい数から、小さい数まで扱えます。

単精度浮動小数点数の最小値は 10^-38(正規化数)もしくは 10^-45 程度、倍精度浮動小数点数でも最小値は 10^-308 (正規化数)もしくは 10^-324(非正規化数)程度ですから、IEEE 754 準拠の浮動小数点数型ではなさそうです。


Windows 電卓で 1 - 1.e-204 - 1 を計算した結果

そのためか少し変わった桁落ちの挙動を示します。同じ程度の大きさの数同士の演算(例: 1.e-9999 + 1.e-9999 = 2.e-9999)であれば桁落ちしませんが、大きい数と小さい数を演算(例: 1 - 1.e-204 - 1 = -9.807971461541689e-149)すると激しく桁落ちします。

Windows 電卓はオープンソース

実は Windows 10 の電卓のソースコードは GitHub に公開されていてGitHub のプロジェクトへのリンク)、誰でも見ることができます。MIT ライセンスで、2019年に公開されたそうです。Windows の標準アプリがオープンソースになっているなんて、調べるまで知りませんでした。

利用方法は簡単で Git clone したあと、Visual Studio で src ディレクトリの下にある Calculator.sln を開くだけです。私のノート PC が非力なせいか、ビルドは非常に遅いですね……。


デバッグの設定

デバッグする際は、設定を「マネージド+ネイティブ」にしておかないと、ネイティブコードを実行している部分のデバッグが全くできず、ブレークポイントなどが無視されるので注意が必要です。

桁落ちしている箇所

コードのありかとデバッグ方法がわかったところで、Windows 電卓のコードを調べます。桁落ちはどこで発生しているでしょうか?

先ほど示した 1 - 1.e-204 - 1 という演算の場合、CalcManager プロジェクトの RatPack/rat.cpp の addrat() 関数にて trimit() を呼んでいる箇所で桁落ちが発生します。

数値表現と加算の実装

//calculator/src/CalcManager/Ratpack/ratpak.h

typedef uint32_t MANTTYPE;

...

//-----------------------------------------------------------------------------
//
//  NUMBER type is a representation of a generic sized generic radix number
//
//-----------------------------------------------------------------------------

#pragma warning(push)
#pragma warning(disable : 4200) // nonstandard extension used : zero-sized array in struct/union
typedef struct _number
{
    int32_t sign;   // The sign of the mantissa, +1, or -1
    int32_t cdigit; // The number of digits, or what passes for digits in the
                    // radix being used.
    int32_t exp;    // The offset of digits from the radix point
                    // (decimal point in radix 10)
    MANTTYPE mant[];
    // This is actually allocated as a continuation of the
    // NUMBER structure.
} NUMBER, *PNUMBER, **PPNUMBER;
#pragma warning(pop)

//-----------------------------------------------------------------------------
//
//  RAT type is a representation radix  on 2 NUMBER types.
//  pp/pq, where pp and pq are pointers to integral NUMBER types.
//
//-----------------------------------------------------------------------------

typedef struct _rat
{
    PNUMBER pp;
    PNUMBER pq;
} RAT, *PRAT;


//calculator/src/CalcManager/Ratpack/rat.cpp

void addrat(_Inout_ PRAT* pa, _In_ PRAT b, int32_t precision)
{
    PNUMBER bot = nullptr;

    if (equnum((*pa)->pq, b->pq))
    {
        // Very special case, q's match.,
        // make sure signs are involved in the calculation
        // we have to do this since the optimization here is only
        // working with the top half of the rationals.
        (*pa)->pp->sign *= (*pa)->pq->sign;
        (*pa)->pq->sign = 1;
        b->pp->sign *= b->pq->sign;
        b->pq->sign = 1;
        addnum(&((*pa)->pp), b->pp, BASEX);
    }
    else
    {
        // Usual case q's aren't the same.
        //★b を変更できないので見づらい書き方になっているが、通分して加算している
        //★pa = pa->pp / pa->pq = A / B, b = b->pp / b->pq = C / D とおくと
        DUPNUM(bot, (*pa)->pq);          //★bot = B
        mulnumx(&bot, b->pq);            //★bot = BD
        mulnumx(&((*pa)->pp), b->pq);    //★pa = AD / B
        mulnumx(&((*pa)->pq), b->pp);    //★pa = AD / BC
        addnum(&((*pa)->pp), (*pa)->pq, BASEX);    //★pa = (AD + BC) / BC
        destroynum((*pa)->pq);
        (*pa)->pq = bot;                 //★pa = (AD + BC) / BD
        trimit(pa, precision);    //★★この呼び出しで桁落ち★★

        // Get rid of negative zeros here.
        (*pa)->pp->sign *= (*pa)->pq->sign;
        (*pa)->pq->sign = 1;
    }

#ifdef ADDGCD
    gcdrat(pa);
#endif
}

Windows 電卓は内部では RAT 型(struct rat 型)で数値を保持しています。RAT は分子 pp と分母 pq の 2つの NUMBER 型で構成された有理数を表す型です。

有理数の加算は小学校で習った通りで、通分して足します。大きな数+小さな数で通分すると、非常に仮数部(構造体の mant、mantissa: 仮数の意味)が長くなってしまいます。例えば 1 + 1.e-9999 だと仮数部が 9511 要素もある int32_t の配列になります。trimit() は一定以上の仮数を切り落とす役目を果たします。


加算の桁落ちを削除した Windows 電卓で 1 - 1.e-204 - 1 を計算した結果

桁落ちを防止したければ trimit() の呼び出しを削れば良いです。しかし不思議なことに完全に trimit() を削除すると Windows 電卓が起動しなくなります。

桁落ちは必須

桁落ちを防止すると起動しなくなる理由は、Windows 電卓はなぜか起動時に CCalcEngine のコンストラクタにて 10^100 を計算しているからです。

Windows 電卓は x^y を x^y = exp(y * ln(x)) で求めているので、log の実装を見ましょう。ソースコードは exp.cpp の _lograt() です。見たところ、十分に小さい項までテイラー展開を続けるアルゴリズムのようです。

log の実装

//calculator/src/CalcManager/Ratpack/exp.cpp

void _lograt(PRAT* px, int32_t precision)
{
    CREATETAYLOR();

    createrat(thisterm);

    // sub one from x
    (*px)->pq->sign *= -1;
    addnum(&((*px)->pp), (*px)->pq, BASEX);
    (*px)->pq->sign *= -1;

    DUPRAT(pret, *px);
    DUPRAT(thisterm, *px);

    n2 = i32tonum(1L, BASEX);
    (*px)->pp->sign *= -1;

    do
    {
        NEXTTERM(*px, MULNUM(n2) INC(n2) DIVNUM(n2), precision);
        TRIMTOP(*px, precision);
    } while (!SMALL_ENOUGH_RAT(thisterm, precision));

    DESTROYTAYLOR();
}

桁落ちを一切なくすとループが終わらなくなってしまいます。変数 thisterm が 0 にならないとループを抜けないところを見ると、ループ条件の SMALL_ENOUGH_RAT がバグっているような気もするんですが……、これ以上の深追いはやめておきます。

他にも √ の計算も非常に遅くなります。パフォーマンスと実用性のバランスから trimit() は必須といえるでしょう。

編集者: すずき(更新: 2022年 1月 6日 03:58)

コメント一覧

  • コメントはありません。
open/close この記事にコメントする



2022年 1月 9日

高速道路の合流と、自動運転の微妙な関係、その解決策

結論から先に書いておくと、高速道路合流の速度問題が解決されていたので、そのメモです。

問題の説明の記事は、プロパイロット2.0の実現で浮かび上がった、高速道路の制限速度問題【岩貞るみこの人道車医】 - レスポンス(Response.jp)、などいろいろ出ているので、詳細はお任せするとして。

簡単に言うと、SAE レベル 3 以上、つまり機械が人の代わりに運転するタイプの自動運転では、法規上の最高速度に問題があります。普段の道路でも最高速度 60km/h で走っている人はあまりいなくて、それも問題なんですけど、そこはさておき。当初、問題として取り上げられていたのが、高速道路への合流でした。

最高速度は道路交通法施行令という法律で定められていて、高速道路の加速車線の制限速度は一般道同様に 60km/h でした。ご存じのとおり、高速道路の本線は 100km/h なので、法律を守って運転すると、次のようになります。

  • 加速車線 60km/h
  • 本線合流(★問題点1: 後ろの車と速度差がありすぎて危ない)
  • 本線 100km/h に加速(★問題点2: 急加速すぎる)

追突の危険性があるのと、本線で急加速せざるを得ない、ギクシャクした運転になってしまいます。

これは昔から存在していた法律上のおかしな点なんですが、人間が運転する場合はあえて法律(=加速度車線の最高速度制限)を無視し、取り締まり側も厳密なことは言わず見逃す、というゆるい運用で問題を避けてきました。しかし自動運転車まで法律を無視する?それは本当に正しいのか??と問題が再燃したわけです。

法律の方がおかしい

引き続き運用でごまかすのも辛いでしょうし、この問題はどうやって直すのかな?と思っていたんですが、どうやら法律の方を直したみたいです。

道路交通法施行令(令和 2年 4月 1日から施行)
第四章の二 高速自動車国道等における自動車の交通方法等の特例

(最高速度)
第二十七条 最高速度のうち、自動車が高速自動車国道の本線車道又はこれに接する加速車線若しくは減速車線を通行する場合の最高速度は、次の各号に掲げる自動車の区分に従い、それぞれ当該各号に定めるとおりとする。
 一 次に掲げる自動車 百キロメートル毎時
   ★★↑60km/h から 100km/h に変わった★★
  イ 大型自動車(三輪のもの並びに牽引するための構造及び装置を有し、かつ、牽引されるための構造及び装置を有する車両を牽引するものを除く。)のうち専ら人を運搬する構造のもの
  ロ 中型自動車(三輪のもの並びに牽引するための構造及び装置を有し、かつ、牽引されるための構造及び装置を有する車両を牽引するものを除く。)のうち、専ら人を運搬する構造のもの又は車両総重量が八千キログラム未満、最大積載重量が五千キログラム未満及び乗車定員が十人以下のもの
  ハ 準中型自動車(三輪のもの並びに牽引するための構造及び装置を有し、かつ、牽引されるための構造及び装置を有する車両を牽引するものを除く。)
  ニ 普通自動車(三輪のもの並びに牽引するための構造及び装置を有し、かつ、牽引されるための構造及び装置を有する車両を牽引するものを除く。)
  ホ 大型自動二輪車
  ヘ 普通自動二輪車
  二 前号イからヘまでに掲げる自動車以外の自動車 八十キロメートル毎時

自動運転の実現に向けて、継続して警察庁にて話し合いが持たれている(自動運転・自動走行 各種有識者会議等 - 警察庁)ので、今後も法律改正されていくことでしょう。

編集者: すずき(更新: 2022年 1月 14日 12:12)

コメント一覧

  • コメントはありません。
open/close この記事にコメントする



2022年 1月 11日

デスクトップ PC の ATX 電源を交換

デスクトップ PC の部品のいくつかは昔の PC から引き継いで使っています。なかでも ATX 電源は交換を怠りがちですが、いざ壊れると起動しなくなるだけでなく、巻き添えで CPU やマザーボードまで故障する可能性もあって故障が怖い部品です。トラブルに遭う前に、予防的に交換します。

買ったのは CoolerMaster V650 GOLD V2(MPY-650V-AFBAG-JP, 650W)です。ヨドバシで 14,000円くらいでした。ちらっと他店の値段を見たら 11,000円くらいでした。ヨドバシは PC パーツがやや高いのかも?

今までは CoolerMaster SilentPro M600(RS-600-AMBA-D3, 600W)を使っていました。Core2 Quad マシンの時代(2009年〜)から使っていましたから、ほぼ付けっぱなしで 13年使ったようです。CoolerMaster の電源は優秀ですね。これからも応援してます。

編集者: すずき(更新: 2022年 1月 12日 17:28)

コメント一覧

  • コメントはありません。
open/close この記事にコメントする



2022年 1月 15日

ハチミツ入りの飴とボツリヌス菌

ハチミツ入りの飴を食べていたら、パッケージに
「はちみつを使用していますので1歳未満の乳児には食べさせないでください。」
と警告があることに気づきました。ハチミツを乳児に与えてはいけないのは、比較的有名な話ハチミツによる乳児のボツリヌス症 - 消費者庁)だと思います。飴の形に加工されていてもやはりダメなのでしょうか?

乳児ボツリヌス症の原因はハチミツに含まれるボツリヌス菌の芽胞です。ボツリヌス菌「食品衛生の窓」 - 東京都福祉保健局によると、ボツリヌス菌の芽胞は熱に強く、殺菌には 120℃ 4分間の加熱処理が必要です。

ハチミツ入りの飴の話に戻ると、

  • ボツリヌス菌の芽胞は土壌に広く分布している(ハチミツにも入る可能性はある)
  • 市販のハチミツは 120℃ 4分間の加熱処理をしない
  • 飴も 120℃ 4分間の加熱処理をしない

なるほど、殺菌するタイミングがなさそうです……。

一般のご家庭でボツリヌス菌を倒せるか?

通常の鍋では、水の沸点(100℃)を超える加熱処理は不可能ですが、圧力鍋を使った場合はどうでしょう?我が家の圧力鍋、パール金属 H-3551(メーカーサイトへのリンク)をみると、高圧側 100kPa、低圧側 60kPa とあります。


飽和水蒸気圧曲線(Wagner の式を使用)

120℃ の飽和水蒸気圧は 198.7kPa のため、大気圧+100kPa = 約 200kPa とすると、沸点は 120℃まで達します。したがって圧力鍋のロックピンが上がり、圧力が規定値に達したときから、4分間以上加熱することで「一般のご家庭でもボツリヌス菌の芽胞は倒せる」と思われます。やるじゃないか、圧力鍋さん。

とはいえ圧力が必ず 200kPa まで達する保証はありませんし、そもそも圧力鍋は殺菌装置ではないので、過信は禁物です。加熱後はちゃんと冷蔵して早めに食べましょう。

飴は……?

圧力鍋が凄いことはわかりましたけど、飴を圧力鍋で煮込むわけにはいきませんから、やっぱりハチミツ入りの飴を乳児にあげてはいけません。

編集者: すずき(更新: 2022年 1月 17日 01:06)

コメント一覧

  • コメントはありません。
open/close この記事にコメントする



2022年 1月 16日

圧力鍋の圧力設定の意味

我が家の圧力鍋は圧力切り替え式でゲージ圧(大気圧 = 0kPa とする記法)で「低圧 60kPa」「高圧 100kPa」となっています。なぜこの数字なんでしょう?

正直に言って設定の意味を理解していなかったんですが、昨日の日記で飽和水蒸気圧曲線(添付の写真、Wagner の式から導出)を見ていて、設定の意味に気付きました。


飽和水蒸気圧曲線(Wagner の式を使用)

グラフの圧力軸はゲージ圧ではなく、絶対圧です。沸点は大気圧を 1kPa として、ゲージ圧+1kPa で概算しました。

  • 通常: 100kPa → 沸点約 100℃
  • 低圧: 160kPa → 沸点約 110℃(実際は 113℃くらいか)
  • 高圧: 200kPa → 沸点約 120℃

これはもう見たままですね。10℃刻みです。とてもわかりやすいですね。ゲージ圧 150〜160kPa(絶対圧 250〜260kPa)の鍋もありますが、さらに上の約 130℃設定(実際は 127℃くらいか)を意味します。

調理器具ですから、温度は切り良く、覚えやすく、メーカーが作りやすい設定値を選んでいるはずです。当然と言えましょう。

全く理解していませんでした

こんなにわかりやすく考慮してくれているにも関わらず、当のユーザーたる俺ときたら……。全く設定の意味を理解せずに「常に高圧の一択」ですからね。メーカーの設計者は泣いてしまいますね。

メモ: 技術系の話は Facebook から転記しておくことにした。一部加筆。

編集者: すずき(更新: 2022年 1月 17日 02:03)

コメント一覧

  • hdk 
    圧力鍋で米を炊く時は、100kPaで炊くと米が変色するのと、加圧に時間がかかるので、60kPaにしています。毎日のことなので常に低圧一択みたいなものです :)
    以前使っていた、カワチで2079円で買った圧力鍋は、単純なおもり式で80kPa固定で、気になるほどの米の変色はなかったように思います。 
    (2022年01月17日 23:59:52)
open/close この記事にコメントする



こんてんつ

open/close wiki
open/close Java API

過去の日記

open/close 2002年
open/close 2003年
open/close 2004年
open/close 2005年
open/close 2006年
open/close 2007年
open/close 2008年
open/close 2009年
open/close 2010年
open/close 2011年
open/close 2012年
open/close 2013年
open/close 2014年
open/close 2015年
open/close 2016年
open/close 2017年
open/close 2018年
open/close 2019年
open/close 2020年
open/close 2021年
open/close 2022年
open/close 過去日記について

その他の情報

open/close アクセス統計
open/close サーバ一覧
open/close サイトの情報