Conjugate Gradient Methodで計算しているがソルバとソルバ前段階で工夫したところ一気に３０倍高速に処理できるようになった。 ３００ｘ３００画像で１分だったのが２秒。 ５５０ｘ４５０画像で３０分だったのが１分弱。 後はOpenMPで並列化できるところは…

一応、精度を満たす計算がまともに出来るようになった。 スパース行列なのでNR（Numerical Recipes）に書いてあるConjugate Gradient Method（2.7 Sparse Linear Systems）を使って計算。 ソルバ自体は割りと高速に計算出来ているがその前処理 Indexed Stora…

まだまだ正しい結果にならない。 誤差なのか拘束条件で指定したピクセルか分からなくなる。拘束条件と判定されるとその周囲の色が動かなくなってしまう。 沢山の考慮漏れとバグがあるがちゃんとそれっぽく計算できたやつが下の絵。Photoshop職人達(?)が歴史…

理論的に不可能 TV等で自動的に白黒カラー化が可能になったような放送を見たという人も多いよで私も凄い技術だと思った。しかし、良く考えると不可能な話だ。カラー(24ビットとして）、白黒８ビット、どうやっても失われた情報は復元できない。 視聴者にイン…