数値流体

固有値方程式まで戻ったらRiemann不変量の選び方にミスが見つかった。で、計算してみたらニュートン法でちゃんと収束した。やれやれ。テスト計算の解析結果でこんなに苦労してたら先行きoptically thick。