Enjoy our world with numerical computations!
薄寒く，驟雨の一日。コートを着るほどではないが，日中での厚めの上着が必要になる気温。今週前半はこの調子で，後半になると温かくなるようである。この調子で梅雨に向かってジメジメする日と五月晴れが交互にやってくるのであろう。
本日は何をやっても「平成最後の」という枕詞が付く日。警備が厳しくなるということで，昨日中に東京日帰り旅行を敢行してきた。先ずは上野で困っしゃくれた後ろ姿とバリボリとポッキーのように竹を食むパンダの親子を堪能。
その後は皇居の様子を伺い，そそくさと帰宅した。
明日の令和元年初日はN先生のお通夜に出席予定。それまでにあれこれGW中にやらなきゃならん仕事に手を付けなきゃいかんのだが，別テーマの本に手を出して進展せず。まぁ新しいこと（そーでもないかな）を勉強するのは楽しいのでつい横道に反れてしまう。まぁあと10年は少しでも面白いと思えることをしていきたいからなぁ。
寝ます。
夏も近づく八十八夜がすぐそこに迫っている。日中の気温は２０度を超えることも珍しくなくなり，朝晩もコートが不要になった。ってことは，戦後日本初の十連休を過ぎた辺りでジメジメベトベトな毎日がやってくることになるわけで，それが終わると灼熱の日々が待っている。ハンコ押しをこなしているうちに夏になりそうな予感。
今年１月に頼んだSIAMの本が，クレームを送ってようやく届いた。なんとなく，○○まで届いて消えちゃった感じがあるんだけど，真相は藪の中。届いていないモノは届いていないわけで，「とっとと届けて，届かないなら返金して」と催促したら返事があった後，最速便で三日で届いた。
Strangの本って，理論も概説もコンパクトにまとまっているんだけど，ワシみたいなバカにはスンナリ染みこんでこないんだよなぁ。機械学習とDeep Learningの解説は参考にしつつ，ワシなりのモノにしていかないと。
先週は予定通り学内教育プロジェクトの報告書がてらの紀要原稿を上げ，明日の新入生向けガイダンス（何度目だ）で使う資料を，前任者制作物に手を加えて仕上げて印刷し，6/9～11に開催されるIEEE国際研究集会に参加申し込みし($480！高ぇ)，Python数値計算入門資料をちょろっと整えた。丸めモード変更手続きについてはイロイロ調べた結果，こういう風にすればよさげなことが分かった。
|
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 |
# 丸めモード変更 # [Windows&Linux] https://rafaelbarreto.wordpress.com/2009/03/30/controlling-fpu-rounding-modes-with-python/ # [macOS] https://stackoverflow.com/questions/16000574/tie-breaking-of-round-with-numpy # OS種別 # https://stackoverflow.com/questions/1854/python-what-os-am-i-running-on import platform # 丸めモード定義 global FE_TOZERO, FE_DOWNWARD, FE_UPWARD, FE_TONEAREST # set_rmode(丸めモード), 丸めモード = get_rmode() global set_rmode, get_rmode # DLL呼び出し # from ctypes import cdll from ctypes import * # OS名 os_name = platform.system() # Windows if os_name == 'Windows': msvcrt = cdll.msvcrt MW_RC = 0x00000300 set_rmode = lambda mode: msvcrt._controlfp_s(0, mode, MW_RC) # get_rmode = lambda : msvcrt._controlfp_s(pointer(current_mode), 0, 0) FE_TOZERO = 0x00000300 FE_DOWNWARD = 0x00000100 FE_UPWARD = 0x00000200 FE_TONEAREST = 0x00000000 def get_rmode(): current_rmode = c_uint() msvcrt._controlfp_s(byref(current_rmode), 0, 0) return current_rmode.value & MW_RC # Linux elif os_name == 'Linux': from ctypes.util import find_library libm = cdll.LoadLibrary(find_library('m')) set_rmode, get_rmode = libm.fesetround, libm.fegetround # x86 FE_TOZERO = 0xc00 FE_DOWNWARD = 0x400 FE_UPWARD = 0x800 FE_TONEAREST = 0 # macOS elif os_name == "Darwin": from ctypes.util import find_library libc = cdll.LoadLibrary(find_library('c')) set_rmode, get_rmode = libc.fesetround, libc.fegetround # x86 FE_TOZERO = 0xc00 FE_DOWNWARD = 0x400 FE_UPWARD = 0x800 FE_TONEAREST = 0 # 丸めモード確認 def print_rmode(): rmode = get_rmode() if rmode == FE_TOZERO: print('切り捨て') elif rmode == FE_DOWNWARD: print('-Infへの丸め') elif rmode == FE_UPWARD: print('+Infへの丸め') elif rmode == FE_TONEAREST: print('最近偶数値丸め') else: print('不明 ', rmode) return rmode # テスト default_rmode = print_rmode() # 切り捨て set_rmode(FE_TOZERO) print_rmode() # 元に戻す set_rmode(default_rmode) print_rmode() |
さて，今週中に「Python数値計算入門（仮）」第２章を上げなければ。上げてGWはのんびりテキストのネタを構築するのだ。
桜満開の頃を過ぎると急激に気温が下がり，朝晩はコートが必要なほど。それでも日中は徐々に春めいてきており，次週は夏も近づく八十八夜を体感できるようである。花粉の飛散は最高潮に達しているようだが，それを除けばうららかな陽気を満喫できる季節になった。良きかな。
つーことで，風疹の予防注射をしてきた。
ワシは昭和44年生まれだが，この近辺の中高年男性は風疹の予防接種を受けておらず，感染するとたちどころに発症してはた迷惑なウィルス拡散生物と化してしまうとのことで，全国的に無料クーポン券を自治体が配布し，少しでも風疹の拡散を防ごうということになっているようだ。ワシについては，静岡市の第一弾クーポン対象者にはなっておらず，今年の後半あたりに送られてくる予定だったようだが，出張の機会が多いワシとしてはいてもたってもいられず，抗体検査の結果，抗体がないということが分かった時点で早々に予防接種を希望してきたのである。もう少し遅れるのかなと思ってたら，需要に生産が追い付いていないワクチンが届き，クーポンも今週半ばに届いたので，本日，近場の内科クリニックで注射を打ってもらったのである。風疹と麻疹の混合ワクチンだそうだが，まぁ目的が果たされれば何でもよろしい。副作用は１０万人に一人ぐらい出るとのことで，担当医師から詳細を伺ったが，なったらなったであきらめるしかないし，それより風疹患者になるリスクの方がはるかに高いから，何の問題もないのである。ま，これで一安心。
新学期が始まって，ハンコ押しマシン，決済処理マシンと化している。いろいろな情報が見られるのはメリットではあるが，知ったからといって何ができるっていうモンでもないことも理解できた。空き時間があれば会議が押し込まれるし，今週は特に会議やらメール処理やらで忙殺されつつ，refereesにボロカスに言われながらも通っちゃった三度目の正直である国際研究集会の最終原稿を金曜日までにsubmitしなければならず，あまつさえ，Python講義資料も作らなきゃならんという４重苦。昨日，PDFチェック済みの原稿を投稿し終えた時には心底ホッとした。講義資料も，とりあえず１週分はクリアしたから，明日までに残りの分を仕上げておしまいになる予定。これをベースに，今月中に「Python数値計算入門（仮）」第２章を完成させたいものである。Yさん，少々お待ち下さいまし。
・・・と，その前に昨年度，学内補助金貰った教育プロジェクトの報告を紀要原稿としてまとめねばならず，締切は次週末。つーことで，各週毎に原稿の締め切りがあるという充実しつつも，疲労困憊50歳の加齢臭をまき散らしながらの過酷な日々がGW前まで続くのである。今年中にあと2本ぐらいは論文投稿もしたいし，「多倍長数値計算入門」も初稿が出るから後期講義までには出てほしいし（他力本願），GW後の方が結構キツイかも。
死なないように頑張りたいと思います。
駿府城公園の桜は満開には至らず，半分程度というところか。今週末の静岡まつりには満開を過ぎたあたりになりそう。ゆっくり散策するのは再来週，その頃には散り際のまばらな桜になっているだろうが，これもなかなか風情があってよろしいモンである。
ということで２０１８年度が本日で終了，明日から２０１９年度で，初日から会議やら辞令交付式やらガイダンスやら卒研ゼミやら，一日のスケジュールがパンパンである。おまけに勤続２０周年というから我ながら長く務めたなぁと嘆息するばかり。節目の年だから役職が振ってきたわけではないけど，まぁ長く務めたモンの義務というものであろう。せいぜい義務を履行した後の隠居生活を楽しみに日々生きていくことにしよう。
「Python3数値計算入門」は，とりあえず第１章をちょろっとPython向けに改変して編集さんに投げておく。若書きの部分は今読むとこっぱづかしいが，何か言いたくて仕方ないという妙な情熱だけは褒めておこう。つーことでこの部分はあまり変える必要がないと判断，饒舌すぎる部分を削ってPythonやJuliaという用語を追加するぐらいで留めた。問題は線型計算の部分，直接法をダラダラ書いても仕方ないし，ScipyとLAPACK/BLASの解説に入れ替えちゃおうかなぁと画策しているが，どの程度になるやら。まぁ，GWまでにはこの山を越えたいところ。
明日に備えて寝ます。
寒の戻りとかで，肌寒い一日。ボチボチ静岡でもソメイヨシノの開花宣言が出るかなと思っていたが，ちょっと遅れそうである。明日は駿府城公園の基準木をもう一度見に行こう。
「多倍長数値計算入門（仮）」，一週間遅れで何とか現担当・Fさんに送付。多重精度計算を何とか最終章に突っ込みたかったのだが，「初期誤差ゼロなら計算中の丸め誤差を消しながら精度維持します」っつー流儀と，「初期誤差は桁を伸ばして末尾に追いやって悪条件問題に対処します」っつー多倍長計算流儀では食い合わせが悪すぎて，あの方法，そのまんまじゃワシら多倍長屋の土台には載らないことが判明。やっぱ正規化過程が重要っつーことで，これはワシの今年の研究テーマにすることにして，書籍からは削除決定。これを決断するまで一週間浪費して，3月半ばには脱稿する予定が狂ってしまった。それでも金曜日には全部終わらせたので，あとは新担当の方とのやりとりだけになる・・・が，さてどーなるやら。まぁ，本来なら昨年には出ている予定がもう一年後ろ倒しになっちゃったからなぁ。ぽしゃらないことを祈るのみ。Webには予定稿を上げたので，プログラム整理してGW前には皆様にURL送りつけますのでよろしゅうに。Twitterにあげたこれとかとこれとか，他にも大量にあるので収拾つけるのが大変でねぇ。GitHubに上げる前に整理しとかんと訳が分からん。4月中は雑務の合間を見てチマチマやります。
つーことで引き継ぎだの見学だのセレモニーだのに忙殺されそうな新年度前のこの時期，そろそろ腹くくって無能な上司役を勤めましょうぞ。桜田五輪相見ていると，人ごととは思えん今日この頃。
今度こそPython3!，金曜までに教材用スクリプト作っておかんと間に合わん！ やるぞ。