ひよこっ!

BibDeskのbibファイルをMendeleyに取り込む

2012-01-30T14:08:00.000+09:00

Mendeleyはクラウド対応の文献管理アプリケーション。Wordでの引用、文献リスト作成にはMendeleyが便利。
ToolsメニューからInstall Word pluginを実行するとAppleScriptがインストールされる
引用箇所でCommand+Control+iを叩いて引用を作成。文献リスト作成箇所ではCommand+Control+b。
PDFから適宜データを取り込んでくれるというふれこみだが、うまくいかなかった。不完全なデータの手直しは不便に感じられた。
試しに、BibDeskのbibファイルをDrag&DropでMedeleyに取り込んでみる。省略記号（例:MWR）で記入していたjournal名が落ちていた。定義をbibファイルの冒頭に書いてあるものは、置き換えられる。Mendeleyでは、journal名の略記（例: Mon. Wea. Rev.）中の\ 、title中の{}が不要である。そのまま取り込むと\や{}が残ってしまう。そこで、Mendeleyの取り込みに適した形にbibファイルを加工した。加工内容は

Journalの定義とデータベースを結合
Journal名の略記中の\を削除
Title中の{}を削除

である。3はTitle = {を[に},を]に一旦置き換えて{}を削除し、[,]をそれぞれTitle =と},に戻している。
bibdesk2mendeley.sh

#!/bin/sh
bibdir="${HOME}/Documents/Bibdesk/bib"
abb=$bibdir/ametsoc.bib
ref=$bibdir/references.bib
tmp=./tmp.bib
cat $abb $ref > $tmp
sed -e 's/\\\ / /g; s/Title = {/[/; s/},/]/' $tmp | \
awk '/^[ \t]*\[/{gsub(/{/,""); gsub(/}/,""); print}; !/^[ \t]*\[/{print}' | \
sed -e 's/\[/Title = {/;s/\]/},/'
rm -f $tmp

iPadでフリック入力

2012-01-09T18:17:00.001+09:00

iOS5では、iPadでもフリック入力ができる。「日本語かな」で右下のキーボードキーを長押して、キーボードを分割すると現れる。今迄気がつかなかった。iOS5になってからも、よく知らずにアップルにフリック入力を要望してしまった。五十音の並びは「あ」が右上から始まるように改善して欲しい。

Bluetoothキーボード

2012-01-08T18:09:00.001+09:00

ELECOM TK-FBP014BKを購入。9台に接続し、切り替えて使える。

配列で不満な点。右shiftの左がpage up。Macに繋いで「ろ」や長音記号が入力できないときは、キーボードの設定が必要。その際、右shiftの左を聞かれたら、page upでなく、もうひとつ左の「ろ」をタイプする必要がある。
Fn+F12がMacとWindowsのキー配列の切替に割り当てられている．Caps LockとCtrl，AltとWindowsとの入れ替えはMacでできる．
iPadやiPhoneに繋いだときは、英語はUS配列となる。これはiOS側が原因のようだ。日本語を入力するときに、JISが指定されているかな入力になる。かな入力派としては一瞬喜んだが、それも束の間「ろ」や長音記号「ー」の入力ができない。Apple純正だとできるのだろうか。

GNU enscriptをMacPortsでインストール

2012-01-07T14:33:00.000+09:00

新・OS XハッキングでGNU enscriptが紹介されていた．MacPortsで探してみる．


$ port search enscript
enscript @1.6.4 (print)
    Replacement for Adobe's 'enscript' program

variantsを確認．


$ port variants enscript
enscript has the variants:
   mediaA4: use A4
   universal: Build for multiple architectures

+mediaA4を指定した方が良さそうだ．


$ sudo port -d install enscript +mediaA4

デフォルトではF77のみ対応．F90用のファイルを紹介する記事．~/.enscriptの下に配置．

CSSの横幅

2011-12-21T16:08:00.000+09:00

CSSレイアウトの基本が分っていなかった。

多摩美 - メディアリテラシー 2010 HTML5+CSS3 入門にある例を考える。

bodyのwidthが960pxとなっている。ウィンドウとのmarginは20px。

navとsectionのwidthは、それぞれ236pxと620px. コンテンツは、856px。

横にnavとsectionを並べるとmargin, border, paddingは、左右でふたつが2列あるので計4つ。

ひとあたりmarginが5px, borderが1px, paddingが20pxなので、コンテンツ以外の幅は合わせて104px。コンテンツの856pxと合わせてbodyのwidthと一致する。

MathJax

2011-12-21T13:32:00.002+09:00

MathJaxは、多くのブラウザで数式を表示するJavaScript. Bloggerでも使える。 Wordpress.comでは許可されていない。その代わり、$\LaTeX$で書いた数式を絵にして貼るプラグインが導入されている。MathJaxで描いた数式は、文字としてコピーできるのがメリット。

使えるようにするには、ドキュメントを参考にして「デザイン」「HTMLを編集」して<head>タグの後でMathJaxを読み込むにする。

Exner函数$\pi$で表した静力学平衡の式

\[\frac{\partial\pi}{\partial z} \equiv -\frac{\;g\;}{c_p\theta}\]

WM3500R+Time Capsule

2011-05-08T23:37:00.002+09:00

WM3500RにクレードルをつけてTime CapsuleとLANケーブルで接続．

AirMacユーティリティでインターネット>インターネット接続はブリッジモードに，AirMac>ワイヤレスは切にした．無線LANはWM3500Rの方を使う．

AirMacユーティリティで無線LANを切った後に設定を確認・変更するときは，LANケーブルで接続するとよい．

うまくいかないときは，設定だけでなく，ケーブルの接続や無線LANのアクセス制御もよく確認する．

iPadからiPhoneにFaceTime発信

2011-05-06T22:27:00.003+09:00

iPadからiPhoneにFaceTimeで発信する場合，iPhoneの電話番号を利用する．iPhone電話番号は国番号から入力する必要がある．iPad 2ユーザーズガイドにある通り，iPhone電話番号は国番号から入力する必要がある．

DR-150で自炊

2011-05-03T07:56:00.002+09:00

先人の知恵を参考に本を裁断してスキャンした．

Canon ドキュメントスキャナ― imageFORMULA DR-150

オルファロータリーカッター S

シンワアル助ゴム付 300mm

オルファリミテッドNL LTD-07

を購入．

道具は揃ったが，いざ本を裁断しようとすると，ためらいを感じた．モノとして気に入っているもの，本棚に飾っておきたい本は裁断すべきではない．スキャンした画像は解像度を上げたとしても，拡大するとジャギーが目立つ．やはり紙はきれいだ．自炊はあくまでも便宜的なもの．日本の出版社がiBook Storeで安価に販売するようになったら，買えば済むことだ．
裁断は以下のような明確な目的を持ってしないと後悔するだろう．

読めればよい．
検索したい．
スペースを節約したい．
たくさん持ち歩きたい．

ロータリーカッターはよく切れる．

DR-150はコンパクトだが，白黒にするとスキャンのスピード十分実用的だ．重送検知はないので，目視で確認しながらスキャンする．OCRは縦書き非対応のようだ．傾きの修正が十分でない．Acrobatはついていないが，もっておいたほうがよさそうだ．

鏡面加工は美しいが，指紋が目立つ．

Windows用にドライバやソフトウェアをインストールしなくても，スキャナをつなぐとUSBメモリのようにドライバやソフトウェアがマウントされる．買った状態ではこの機能がONになっている．背面のスイッチでこの機能をOFFにしないと，Macではスキャンができない．

地球上で最もお客様を大切にする会社

2011-01-29T21:40:00.002+09:00

Amazonでヤマハ piaggero NP-V80を購入．

ところが右手前が浮いていて，ガタつく．ヤマハに問い合わせたら，不良品なので，販売店を通じて交換するとのこと．携帯からナビダイヤルがつなかった．電話代は自腹... 三洋さんはAQUAWARMの水受けが歪んでるといったら，二人で訪ねてきてくれた．親切な会社は生き残れないのかな...

さて，発泡スチロールを正しい向きにして梱包するのに苦労したが，なんとか箱に収まり，配送業者に持っていってもらった．ところが後でACアダプタをマッサージクッションのものと間違えてしまったことに気づいた．ACアダプタは送り返すことになったが，Amazonは着払いでよろしいという．間違えて送ったACアダプタも返してくれるそうだ．

メールの最後に


Amazon.co.jp は、お客様からのご意見により、地球上で最もお客様を大切にする会社を目指しています。

とあった．一日落ち込んでいたが，親切な対応がうれしかった．担当の方，ありがとうございます．

Time Machineからフォルダが復元できない

2011-01-27T09:27:00.002+09:00

毎日送られてくる迷惑メール一覧をまとめて消そうと思い，Command+Aで選択したら全ての選択されていて，全てのメッセージがゴミ箱に行きなった．整理せず「受信」に全てためているので，25,000件くらい．ゴミ箱への移動も時間がかかる．移動を待って元に戻せば良かったのだが，短気をおこして止めようとしたのがいけなかった．Mailを強制終了することになり，ゴミ箱にある1か月以上前のメッセージは再び立ち上げたときに消されてしまった．

Time Machineがあるじゃないか．

ところが，復元したはずの~/Library/Mailは空．Finderでコピーしても同様．ファイルは復元されるがフォルダ（ディレクトリ）は復元されない．

Time MachineのGUIを使わずに


tar cf - | (cd 保存先; tar xvf -)

しても特殊なハードリンク (Archive Directory Link)がコピーされるだけで実体ではない．

試行錯誤の結果，別のマシンにTime Machineに使っているディスクをマウントしたところ，フォルダを復元できた．Finderでコピーしても，DockのTime Machineアイコンを長押しするか右クリックして「ほかのTime Machineディスクをブラウズ...」を選んでGUIを使ってもできた．

lsbomでreceiptの中身を調べる

2011-01-18T12:17:00.002+09:00

MacPorts-develでMac OS X 10.6.6で多数の/usrのファイルが追加または更新されたことが話題になった．更新されたファイルの確認は，lsbomでできる．忘れそうなのでメモ．


$ lsbom /var/db/receipts/com.apple.pkg.update.os.10.6.6.patch.bom | grep /usr

VimでFortranのsyntax colouring

2011-01-06T16:39:00.003+09:00

このところsyntax colouringを使ってきたが，拡張子がf90なのに固定形式として認識されているようだった．
調べてみると，先頭250行を見て判断していることが分かった．
このソースは，左端から6桁開いているが，72桁より長いところが
無効としてハイライトされていた．

拡張子がfとFの時に固定形式とするには，.vimrcに次のように書く．
filetype plugin indent onはsyntax onよりも前でなくてはならない．

なお，既定では，新しいファイルは固定形式と仮定するようだ．
これまでは，少し書いてから，いったん保存して開き直していた．


filetype plugin indent on
let s:extfname = expand("%:e")
if s:extfname ==? "f" || s:extfname ==? "F"
    let fortran_fixed_source=1
    unlet! fortran_free_source
else
    let fortran_free_source=1
    unlet! fortran_fixed_source
endif
syntax on
autocmd FileType * set formatoptions-=tcq

最後の行は，自動改行を防ぐおまじない．

風呂掃除

2010-12-30T22:52:00.002+09:00

風呂用洗剤をきらした．特に自然派お掃除にこだわるわけでも，節約志向でもないが，買いにいくのも面倒だったので，重曹水で代用．

水 200 cc をメジャーカップに入れ重曹大さじ 1 を加える．そのままではすべて溶けないので，電子レンジで1分．よくかきまぜて，残りも溶かす．霧吹きに移す．

浴槽に吹きかけて軽くこする．泡が立たないので地味だ．でも，実際入浴時にとてもきれいになったと感じた．

Mac OS X ServerでMacPortsを共有

2010-12-28T12:11:00.004+09:00

MacPortsを使うと，フリーウェアのインストール手間が大幅に軽減される．しかし台数が増えてくると，アップデートも手間に感じられるようになってくる．そこでMacPortsを共有することを考えた．

MacPortsのインストールされたXserveで/optを共有しようとしたら，サーバ管理の共有には表示されなかった．表示されないディレクトリに移動するときは，Finderで移動>サーバに接続...を選ぶところだ．サーバ管理には隠しディレクトリを表示するオプションは見当たらない．アップルの記事によると，コマンドラインを使う必要があることが分かった．

以下を実行すると共有ポイントとして表示される．/optが格納されているボリューム（たとえばServer HD）には表示されない．


sudo sharing -a /opt

AFPで思ったところ（/opt/local）にマウントできなかったので，上記の方法は中止．

昔ながらのNFSでエクスポートとマウントすることにした．

サーバ側で/etc/exportsに記述.
NFSを起動 (sudo nfsd enable).
クライアント側でマウント (sudo mount_nfs サーバ:/opt/local マウントポイント)
クライアント側で/etc/fstabに記述.

XserveにSSを追加できない

2010-12-20T22:06:00.004+09:00

SSDを起動ディスクにしたXserveは起動が速いと同様に進めていたからか,
同僚がSSDを搭載したXserveを発注しようとしたところ,
SSDを追加できなくなっていた.

アップルに問い合わせをしたが, カスタマイズでのSSDの追加は終了とのこと.

どうも1週間くらい前に密かに終了していたようだ.

写真はXserveのSSD周辺. iPhone 3Gで撮影した最後の写真.

Octaveで2Dグラフ

2010-12-17T09:59:00.005+09:00

Octaveのマニュアルなどの入門では, sinのような組込函数を描く例が多い.

描きたいのはsinではない. 論文に出てきた式を図にして確認したい.

例として描いてみるのは, CALIPSO搭載CALIOPセンサを使って見積もられた過冷却水の気温に対する依存性 (Hu et al. 2010のFig. 6dの赤い線) である.

pやfを函数定義 (function endfunction) してもよいが, ここでは単にx軸に対応する値を計算するにとどめる.

tの階乗は要素毎なので, スカラの階乗**に.をつけて.**とする. %のあとの文字はコメントである.

MacPortsでインストールしたばあい, AquaTermがデフォルトなので, 描かれた図はPDFかEPSで保存できる.


t = -40:0.1:0; % -40°Cから0°Cまで0.1°C刻み
p = 7.6725 + 1.0118*t + 0.1422*t.**2 + 0.0106*t.**3 + 3.39e-4*t.**4 +3.95e-6*t.**5; % 式(4)
f = 1./(1+exp(-p)); % 式(1)
plot(t,f)

zencoding.vim

2010-12-05T17:33:00.002+09:00

HTMLやCSSを手早く書くためのvimスクリプトzencoding.vim

~/.vim/bundleにpathogenをclone.

~/.vim/autoload/にautoload/pathogen.vimのリンクを張る.

.vimrcの設定.


call pathogen#runtime_append_all_bundles()
call pathogen#helptags()

Xserve 2010年1月31日で販売終了

2010-11-07T14:24:00.002+09:00

Xserveの新製品がでないないので, 予感はしていたが, 2010年1月31日で販売終了が発表された.

Mac ProかMac miniへの移行が示されているので, Mac OS X Serverは継続して販売するようである.

Mac Proの方が拡張性があり, 製品の更新が早いので, パフォーマンスも上である.
Mac OS X Serverは社内の情報共有にカスタマイズされていると思われるが, その用途であればMac miniで十分かもしれない.
したがって, 大して売れていなかったと思われるし, 困る人も少ないのだろう.

私の職場では, Mac OS X Serverを導入したMac Proも使っているが, ここ数年複数のXserveを導入した. Xserveは, RAIDの頭として使っている. ラックに据え付けるには, 1UのXserveが最適なので, 販売終了は残念だ. Mac Proは12Uもあり, ラックには向かない. Mac miniは1Uだそうだが, インターフェースが不足している.

Xserveはラックマウントキットのできがすばらしい. レールは手を傷つけない塗装がしてある. 金具にはマグネットがついているので, 一人でも簡単に設置できる. ついでに箱もすばらしい.

Cybershot DSC-HX5V

2010-11-07T11:29:00.004+09:00

SONYのGPS付デジカメDSC-HX5Vを購入.

位置取得時間を短縮するために，GPSアシストデータを更新する必要がある. マニュアルによると, カメラを専用ケーブルでPCにつなぎ, Windows版しかないPMBから更新することになっている.

Kakaku.comのクチコミ掲示板にメモリカードの所定の場所にファイルを保存しているだけだとあったので, シェルスクリプトを作成した. このスクリプトを使えば, 専用ケーブルもWindows PCも不要だ. 引数にはカードの名前を指定する. DSC-HX5Vはメモリースティックだけでなく, SDカードも使える. 最近のMacBook Proには, SDカードスロットがあるので, カードリーダもいらない.


#!/bin/sh
URL=http://control.d-imaging.sony.co.jp/GPS/assistme.dat
CURL=/usr/bin/curl
CURLOPTS="--create-dirs -o"
if [ $# -lt 1 ]; then
  echo "Usage:: $0 Volume"
  exit
fi
VOL=$1
OUT=/Volumes/$VOL/PRIVATE/SONY/GPS/ASSISTME.DAT
$CURL $CURLOPTS "$OUT" $URL

リソースフォーク

2010-10-21T17:00:00.002+09:00

Safariでデータをダウンロードしたら，リソースフォークがついてしまった．tarしたときに気がついた．過去の遺産, そろそろ廃止してほしい.

tarするときは分離されるが, ゴミが含まれてしまう.
環境変数を設定すると, データだけがアーカイブされる.


export COPYFILE_DISABLE=true

MacWiki: リソースフォーク

atlasなしのoctave

2010-09-12T21:14:00.003+09:00

Accelerate.frameworkには，問題があるために，octaveは行列ライブラリとしてatlasに依存させている．

atlasはコンパイルに時間がかかる上に, g95のvariantがない．gcc44のコンパイルも加えると，かなり長い時間かかる．

そこでarpack, octave用にパッチを作成した．

arpack用

octave用

2010-11-21追記: mkoctfileがうまくいかなかった問題を修正.
octaveのビルド中に一時的なディレクトリに生成されるdotwrp.oへの
パスが残ってしまった.
以下のラッパは独立したportとしたので, 下記にあるように,
filesに入れておく必要はない.

g95を使う場合は, octaveのPortfileパッチをあてるだけでなく，
filesの中にdotwrp.f90を入れておく必要がある．

gfortranの場合は


sudo port -d install octave +no_atlas

g95の場合は


sudo port -d install octave +g95

とする．

Accelerateのcdotu, cdotc, zdotc, zdotuで出るsegmentation faultを回避する

2010-08-09T22:17:00.006+09:00

MacPortsのOctaveがATLASに移行してしまい, Accelerateを使わなくなった.
ATLASはコンパイルに時間がかかる上, g95がサポートされていない.
ATLASがg95サポートするようにしてもよいのだが, OSに添付され最適化されているAccelerate frameworkを使いたい.
そこでAccelerateを復活させようとしている.

問題となるのは, OctaveのconfigureでCDOTU, ZDOTUのテストにSegmentation faultなどが出て失敗することである.

原因は, AccelerateはCBLASなのでCやg77のABIに従い, 返り値が複素数を指すポインタであるのに, g95, gfortranは浮動小数点レジスタファイルであること. gfortranでは-ff2cとすることで解決する. g95ではラッパが必要. Appleが示したラッパをちょっとだけ改変して関数名に_をつけたので, ラッパは-fno-underscoringをつけてコンパイルするが, 呼び出し側ではオプションは不要.

2010-11-20追記: githubにて公開. MacPortsを使ってインストール可．


!g95 -fno-underscoring dotwrp.f90 -c
double complex function zdotc_(n, zx, incx, zy, incy)
double complex zx(*), zy(*), z
integer n, incx, incy

call cblas_zdotc_sub(%val(n), zx, %val(incx), zy, %val(incy), z)

zdotc_ = z
return
end

double complex function zdotu_(n, zx, incx, zy, incy)
double complex zx(*), zy(*), z
integer n, incx, incy

call cblas_zdotu_sub(%val(n), zx, %val(incx), zy, %val(incy), z)

zdotu_ = z
return
end

complex function cdotc_(n, cx, incx, cy, incy)
complex cx(*), cy(*), c
integer n, incx, incy

call cblas_cdotc_sub(%val(n), cx, %val(incx), cy, %val(incy), c)

cdotc_ = c
return
end

complex function cdotu_(n, cx, incx, cy, incy)
complex cx(*), cy(*), c
integer n, incx, incy

call cblas_cdotu_sub(%val(n), cx, %val(incx), cy, %val(incy), c)

cdotu_ = c
return
end

pandocで日本語epub

2010-08-03T14:18:00.003+09:00

pandocでepubを作る時, metadata.xmlに言語を指定できるが, content.opfに書き込まれるだけで, 個々のxhtmlには言語は指定されない. このままだとAdobe Digital Editionsでは文字化けしてしまう. 横浜工文社のbatスクリプトを参考にpandocが生成したepubを解凍, xhtmlを編集, 再びepubにまとめるスクリプトを作成した.


#!/bin/sh
epub=$1
bname=`basename ${epub}`
base=${bname%.*}
epubo=${base}_jp.epub
TMPDIR=/tmp/${base}

if [ -d ${TMPDIR} ]; then
  rm -rf ${TMPDIR}
fi
if [ -f ${TMPDIR}/../${epubo} ]; then
  rm -f ${TMPDIR}/../${epubo}
fi
if [ -f ${epubo} ]; then
  rm -f ${epubo}
fi
mkdir ${TMPDIR}

unzip ${epub} -d ${TMPDIR}

cd ${TMPDIR}
for f in *.xhtml; do
  g=${f}j
  sed -e 's|html xmlns="http://www.w3.org/1999/xhtml"|html xmlns="http://www.w3.org/1999/xhtml" xml:lang="ja" lang="ja"|' ${f} > ${g}
  mv -f ${g} ${f}
done

zip -0 ../${epubo} mimetype
zip -r ../${epubo} * -x mimetype

mv ../${epubo} ${OLDPWD}
rm -rf ${TMPDIR}
rm -f ${TMPDIR}../${epubo}

連番画像をQuickTimeに

2010-08-02T14:23:00.004+09:00

MacWikiにあるものを少し修正. Snow LeopardではQuicktime PlayerではなくQuickTime Player 7とする.


tell application "QuickTime Player"
 activate
 set firstFile to choose file with prompt "連番の最初のファイルを選択"
 display dialog "1秒あたりのフレーム数" buttons {"OK"} default answer "2" default button 1 with icon 1
 -- open image sequence firstFile seconds per frame 2
 open image sequence firstFile frames per second text returned of result
 save document 1 as self contained
end tell

2011/1/17追記. Snow Leopardでは上記スクリプトはコンパイルできず．Leopardで書いたスクリプトをSnow Leopardで開いたら，以下のようになっていた．


tell application "QuickTime Player"
 activate
 set firstFile to choose file with prompt "連番の最初のファイルを選択"
 display dialog "1秒あたりのフレーム数" buttons {"OK"} default answer "2" default button 1 with icon 1
 -- open image sequence firstFile seconds per frame 2
 «event MVWRopis» firstFile given «class fpsc»:text returned of result
 save document 1 as «constant savkflat»
end tell

pandocでepub

2010-08-01T16:35:00.006+09:00

さまざまな形式の文書が作成できるpandocが1.6からepubに対応.

残念ながらMacPortsのpandocは古く, ビルドにも問題があるようだ.

pandocはHaskellで書かれている. Haskellは, コンパイラghcだけでなく, 独自のパッケージ管理システムcabal等を含む包括的なプラットフォーム.

まず, Mac OS X用のHaskellをインストール.

ghcは/usr/bin, cabalは/usr/localにインストールされる.

次にpandocのインストール.

cabal update
cabal install pandoc

バイナリは~/.cabal/binに, マニュアルは~/.cabal/share/manにインストールされる.

MacPortsのptexにはunicodeパッケージは含まれないので，別途入手し~/texmf/texなどに置く（置いたらmktexlsr）．

Mac上ではAdobe Digital Editionsでepub形式のファイルを読むことができる.

多言語ウェブサイト

2010-06-14T16:15:00.005+09:00

同僚に教えてもらったことを糸口に, ウェブサイトの日本語と英語となど言語の切り替えについて調べた．

言語やエンコーディングなどの処理は, mod_mimeが担当. 言語, 文字セットと拡張子との対応は, MacのApache2の場合/etc/httpd/extra/httpd-languages.confに定義されていて，特に設定の必要はない. もしhttpd.confで読み込みがコメントアウトされている場合は, コメントアウトを外す. またOptionsにMultiviewsが必要なのでうまくいかないときは確認する.

順序は何でもよさそうなのだが, MacのApache2では, ja.htmlではうまくいかずhtml.jaだとうまくいった

Safariの場合, ブラウザ自体の言語設定はなく, システム環境設定の「言語とテキスト」で設定. Safariを再起動すると有効になる.

iPad最大の問題点

2010-06-10T21:26:00.002+09:00

かな入力がなくローマ字かな変換だけ．中国語にはあるのに，日本語には手書きをどうして用意してくれないんだろう．

ncarg 5.2.0

2010-05-26T15:53:00.002+09:00

MacPortsのncargを5.2.0に更新. 変更点は，

HDFEOS5ファイルが読めるようになる.
PNGでの書き出しがサポートされた．
サイズの小さいPDFが作られる (newpdf)．
カラーテーブルが充実．
UDUNITSをやめてUDUNITS-2を採用.

などなど.

JRA25

2010-05-26T15:50:00.003+09:00

JRA25のGRIBファイルをMacPortsで入れたGrADS, GrADS2で読もうとしたが, Bus errorやSegmentation faultが出る.
先のTigerでのコンパイルとは別のマシン.

nclではうまく読めるようだ.

alloc: invalid block

2010-05-26T11:25:00.003+09:00

Tigerで更新したportsのチェックをしようとしたところ，

alloc: invalid block

というエラーが出た．Tclから出ているらしい．

export MACOS_DEPLOYMENT_TARGET=10.4

としたらエラーは出なくなった.

しかしコンパイル時にバスエラーが出る．
TigerはもうMacPortsではサポートしていないが，明らかに何かおかしい．

upc

2010-05-19T18:03:00.004+09:00

Unified Parallel Cは，Co-array Fortran同様にメッセージ交換ではなく，PGAS (partitioned global address space)ブログラムモデルに基づく言語. 1ノード用のgccupcをMacPortsに登録. 複数ノードで実行するには, Berkeleyのラインタイムが必要. 連休中にこちらのコンパイルもできたが, パッケージとして調整中.

szip

2010-05-16T17:49:00.004+09:00

szipのdecoderはフリーだが, encoderは営利目的にはライセンスを取得することが必要．MacPortsのszipのデフォルトではencoderは無効. +encoderを付けて有効化する必要がある. HDF4などからszipによる圧縮を利用したい場合は注意.

第110回地球情報館公開セミナー

2010-04-17T21:01:00.005+09:00

公開セミナーには，寒い中たくさんの方が集まってくださいました．話が長くなってしまい，測候所の廃止，アイスランドの火山の噴火，寒い春と北極振動の関係に関する質問をいただきました．お答えになったか分かりませんが，聴衆の皆様の気象に対する強い関心を感じることができました．セミナースライドをアップロードしました．

掃除機の回転ブラシ

2010-03-28T23:49:00.006+09:00

3年ほど前に, SHARPの縦型サイクロン掃除機を購入した. ずいぶん安いが, 母がくれたTupperwareと共同開発したという掃除機と似た構造. 洋服のほこりとりのような部品やスパチュラないしはスクイージーのようなの羽根がついている回転ブラシもよく似ている. いかにもゴミがよくとれそうな構造だ.

ところが, このところヘッドの回転ブラシが回らなくなっていた. 先日ブラシを外して掃除した. ゴミ, とくに糸や髪の毛がからまっていた. きれいにしたのだが, 回転しなかった.

ヘッドの交換が必要かもしれないと思ったが, ベルトは切れていないし, モーターで回しているようにも思えなかった. 気流で回しているのなら, どこか詰まっていない限りブラシは回転するはずだが, 詰まっているところはない. ねじが見えたので思い切って分解してみることにした. ねじは6本. 前の方にツメがあり, ふたを外すときにいやな音がしたが, 折らずにすんだ.

ふたを開けてみると, 内部にごみが溜っていた. 想像通り, 羽根のついた回転軸でブラシは回転する構造. ブラシ回転のスイッチは, 気流の流れを変えるだけの簡単なもの. 外からは見えない回転軸の歯車の部分にも髪の毛が絡まっていた. これでは回るはずがない. きれいにそうじをして元に戻す.

ふたを閉めてねじを元通りに締める. 掃除機のスイッチを入れると, 回転ブラシが勢いよく, 「ウィーン」と音を立てて回った.

Time Capsule導入

2010-03-08T22:17:00.005+09:00

AmazonでAPPLE Time Capsule 2TB MC344J/Aを購入．

自作PCの外付HDDでのTime Machineには成功したが，PCとHDDの電源を入れっぱなしにしておくことには抵抗がある．

サードパーティのNASでもTime Machineは可能だし，NAS対応のケースに安価なバルクHDDを入れても良い．実際PCにつないだのは，ロジテッククレードルタイプ HDDリーダーライタ eSATA&USB2.0接続タイプ LHR-DS01SAU2.

Apple AirMac Extremeベースステーション MC340J/Aに，現在使っているハードディスクをNASケースに入れてつなごうと考えた.

電源内蔵のハードディスクケースは現在製造物責任法のため，手に入りにくいようだ．しかも安いとはいえない.

Time Capsuleは割高だと言われているが，無線LAN, ギガビットの有線ハブ (3つ), 2TBハードディスク，電源内蔵ハードディスクと考えると決して高くないように思われる.

Ethernetケーブルと電源をつなぐとAirMacユーティリティがTime Capsuleを自動的に検出. AirMacユーティリティの更新を求められる. 添付のCDからイントールすると，1) 新たなネットワークを作るか, 2) 既存の無線LANを置き換えるか, 3) 既存の無線LANに接続するか聞かれる. プロパイダから借りているルータに差している無線LANカードを置き換えるので, 2) に該当. 無線LANカードを抜くと, 設定がすべて引き継がれている. すばらしい.

いったんTime Machineの設定は後回しにしたが, 後で有線で接続し, Time CapsuleのDataを選択してバックアップ実行.

なお，比較的最近 (2010/2/28) 付のユーザーレビューで最悪との評がある.
買ったばかりなので分からないが, 旧製品の問題はアップルが対応している. また有線がギガビットでないことで躊躇したとのコメントも目にしたが，現在の製品はWANx1, LANx3ポートすべてギガビットになっている．新旧製品の情報が入り乱れているので注意．

f2py

2010-02-08T22:30:00.001+09:00

f2pyは, numpyに含まれF77/f90のサブルーチンをpythonから利用するために用いるスクリプト. F77の場合どれが出力となるかをコメントに書く必要が生じることもあるが, 簡単に共有ライブラリの生成ができる. Eigen2やmkoctfileくらい簡単だ. NCLには独自のWRAPITがあるが, PyNGLではnumpyのf2pyを使う.

PyNGL

2010-02-08T22:29:00.000+09:00

NetCDFなどのファイルを読み書きするPyNIOとNCLをPythonから使えるようにするPyNGLをMacPortsに追加した (py26-nio．py26-ngl). NCL (ncarg) 同様, 入手には登録が必要なのでEarth System Gridから手動でソースを取得する必要がある. NCLで描ける絵のほとんどに対応しているようだ. Pythonのテキスト処理や辞書などのデータ構造や, numpy, scipyを使った数値処理が利用できるようになる. 依存するNCL等のコンパイル時に用いたコンパイラがg95のときは, インストールするときに+g95をつける.

英語の月名

2009-12-07T14:51:00.002+09:00

数字の月名を英語の月名に置き換えるとき配列があると便利だが，/bin/shには無い. caseで書くのも面倒．

通常使われるテクニック．


set dummy JAN FEB MAR APR MAY JUN JUL AUG SEP OCT NOV DEC
shift $i
moneng=$1

dateコマンドを使ったもの．


moneng=`LANG=C date -j -f "%m" $mon "+%b" | tr "[:lower:]" "[:upper:]"`

install_name

2009-11-17T16:17:00.005+09:00

Macの共有ライブラリには，インストール先の情報が記録される．インストール先は，-install_nameの指定があればそのディレクトリ，なければ-oで指定したものが使われる．ライブラリのinstall nameは，これをリンクしたライブラリやバイナリに転写されていく．

-install_nameの後はスペース．=を入れると正しく動作しないので注意．

オブジェクトファイルに関する情報は，otoolで表示することができる．
ライブラリの名前とバージョン，リンクしているオブジェクトやライブラリを表示するには，

otool -L ライブラリまたはバイナリ

ライブラリのinstall nameを表示するには，

otool -D ライブラリ

を使う．

共有ライブラリのinstall nameの変更は，install_name_toolを使う．

install_name_tool -id 新しいパス

でできる．ライブラリやバイナリに書き込まれた他のライブラリのinstall nameを変更するには，

install_name_tool -change 古いパス 新しいパス

とする．

多くのフリーウェアでは，libtoolがこのあたりの面倒を見てくれるが，シェルスクリプトやGNU Makeを駆使して機種依存の問題を解消しようとしているものもある．

vaporは，configureの変わりにoptions.mk, site.mkを手動で編集したり，機種依存のインクルードファイル（例: Darwin.mk）を用意している．ライブラリには，生成したディレクトリのパスがinstall_nameとして記録されているため，うまく動作しない．DYLD_LIBRARY_PATHを設定してライブラリを見つけようとしている．既定の-two_levelnamespaceだとコマンドラインと共有ライブラリに記録されたパスしか検索しないので，ライブラリは見つからない．

make/config/Darwin.mkのSHARED_LDFLAGSに-install_nameを指定して対処．

Octave-3.2.x

2009-11-09T17:59:00.003+09:00

Octave-3.2以降の更新が滞っていた. doxygenのコンパイルが通らないため. doxygenをリンクするときに, srcからみて../libにあるlibmd5.aがなぜか見つからないためだった. 直接../lib/libmd5.aを指定したところ, リンクがうまくいくようになった (#22426).

doxygenは, octaveが依存しているftglがHTMLのドキュメント作成のために要求していたが, ドキュメントは必ずしも必要ない. またPDFのドキュメントを作るためにconfigureはpdflatexなども探しているが依存関係に指定されていない. そこで, docをvariantととして分離するように提案した (#22427).

Leopardでx86_64

2009-11-01T17:30:00.002+09:00

MacPorts 1.8では, ${prefix}/etc/macports/macports.conf (${prefix}はたいてい/opt/local) にコンパイル時のアーキテクチャbuild_archを設定できる. Snow Leopardでは, x86_64がデフォルトなのですべて64-bitでコンパイルされる. Tiger, Leopardでは, Intel Macの場合はi386, PowerPCの場合はppcなので32-bitになる.

Leopardでbuild_archにx86_64を設定して, さらからコンパイルしてみた. p5-locale-gettextがデフォルトのアーキテクチャ以外でコンパイルできない (19381). autoconfなどが依存しているhelp2manがp5-locale-gettextに依存しているので, 先に進めない. nomaintainerなので困ったものだ.

Tclパッケージのパス

2009-10-27T12:28:00.003+09:00

Mac OS X添付の/usr/bin/tclshは, ~/Library/Tclにパッケージをおけば見つけてくれる.
MacPortsでインストールした/opt/local/bin/tclshは,


$ /opt/local/bin/tclsh
% set auto_path
/opt/local/lib/tcl8.5 /opt/local/lib

しか見ていない.

環境変数TCLLIBPATHに追加するか，スクリプトの冒頭で


lappend auto_path /Users/foo/Library/Tcl

と書けば追加できる. Tclで完結する後者の方が良さそう.

Tcl httpでssl

2009-10-27T12:14:00.002+09:00

SSLでログインしたあと, ファイルをダウンロードするためのTclスクリプトを書いてみる. wgetやcurlを使ってもよいが, tlsパッケージを併用するとhttpでSSL (https) が扱える.


package require http
package require tls
::http::register https 443 ::tls::socket

SSLサイトでのログイン. まず::http::formatQueryでPOSTするデータを作る. 送信すべきデータはサイトにより異なる.


set login [::http::formatQuery email $email passwd $passwd action login]

メタデータからクッキーをリストに保存.


set tok [::http::geturl $loginurl -query $login]
upvar #0 $tok state
set cookies [list]
foreach {name value} $state(meta) {
if {$name eq "Set-Cookie"} {
  lappend cookies [lindex [split $value {;}] 0]
}
}
::http::cleanup $tok

クッキーは, ヘッダにつけて送信.


set tok2 [::http::geturl $dataurl$fname -headers [list Cookie [join $cookies {;}]]]

リダイレクトされる場合は, メタデータのLocationを使う.


upvar #0 $tok2 state
array set meta $state(meta)

$meta(Location)にURLが入っているはず.
リダイレクト先からのダウンロードにもクッキーが必要.


set tok3 [::http::geturl $meta(Location) -channel $f -headers [list Cookie [join $cookies {;}]]]
::http::cleanup $tok2
::http::cleanup $tok3

Tcl httpバッケージ

2009-10-26T22:39:00.002+09:00

Tclのhttpパッケージは，httpクライアントの動作を実装したもの. ウェププラウザやwgetやcurlがサーバとのやりとりでするようなものがTclのコマンドで簡単に書ける．

::http:geturlの-validateオプションは, 0でない引数を与えると, ヘッダだけを取得する．結果は返り値に保存されるが, upvarを使ってアクセスする必要がある．ここでは, ヘッダにContent-Lengthが含まれているか調べて, ファイルが存在するかどうかを確認している.


  set token [::http::geturl $url -validate 1]
# obtain results of geturl
  upvar  #0 $token state
  set meta $state(meta)
# look for Content-Length
  return [lsearch $meta Content-Length]

g95 on Snow Leopard

2009-10-26T21:41:00.005+09:00

Snow Leopardでg95がコンパイルできずに苦労している．

gcc-4.0.4はi386では一応コンパイルできるが, g95のコンパイル途中で失敗する．

gcc-4.0.4とgcc-4.1.2では, configureで--buildにx86_64-apple-darwin10を与えるとうまくいかない．
x86_64はサポートされていないようだ．

Appleが修正したgcc-5465は，i386とx86_64両方のライブラリを作る．g95のコンパイルのときに，未定義のシンボルが多数出てしまいどうしようもなくなった．

gcc-4.3.4はいろいろと細工をして最後までコンパイルすることはできたが，g95の実行時にエラーが出る．gcc-4.3では，仕組みが大きく変わったのだろう．そこまで直せそうにない．

それならばとgcc-4.2.4を試したみた．gcc-4.3.4でしたような細工がいくつか必要だったが，gcc-4.3.4よりは難しくなかった．一応コンパイルに成功し，簡単なプログラムが動作することを確認した．細工にはいい加減なものもあるので，あとで問題が発生するかもしれない．それでもともかくめでたい．fileで確認したところ，確かにx86_64のコードが生成されている．

次はこれをPortfileにまとめなくてはならない．

11/9追記: gcc-4.2.4にリンクしたg95は, Snow LeopardだけでなくIntel MacであればTigerでもLeopardでも動いたが, PowerPCのMacでは動かなかった. そこで, 従来の4.0.4をベースにしてgcc42はvariantととした. r60290で変更をコミットした.
さらにr60340でSnow Leopardの検出を簡素化した.

g95でLAPACK

2009-10-14T18:20:00.006+09:00

g95で実対称行列の固有値を求めるサブルーチンには，

dsyev (driver),

dsyevd (divide and conquer)

dsyevx (expert)

dsyevr (relatively robust representations)

といくつかある（LAPACKのユーザガイド及びmanページ参照）．以下の簡単な例では，最初の3つの固有値は一致した．dsyevrはe-14〜e-15程度の差があった．

Macには，CBLAS, CLAPACKがAccelerate FrameworkとしてOSに添付されている．g95がインストールしてあれば，


g95 foo.f90 -framework accelerate

でライブラリをリンクできる．とっても簡単．


program eigen
  implicit none

  character, parameter :: jobz = "V", uplo = "L"
  integer, parameter :: n = 3, lda = n, lwork = 3*n-1
  integer :: info, i
  real*8, dimension(lda,n) :: a
  real*8, dimension(n) :: w
  real*8, dimension(lwork) :: work

  a(:,1) = (/1.0d0,  2.0d0,  3.0d0/)
  a(:,2) = (/2.0d0,  3.0d0, -4.0d0/)
  a(:,3) = (/3.0d0, -4.0d0,  5.0d0/)

  call dsyev(jobz, uplo, n, a, lda, w, work, lwork, info)
  print *, "info=", info
  do i=1, n
    print *, "mode:", i, " eigenvalue=", w(i)
    print *, "         eigenvector", a(:,i)
  end do
end program eigen

MacPortsでphp5をインストール

2009-10-03T18:16:00.006+09:00

Snow LeopardのPHP5はGDが有効化されているが，freetypeがサボートされていない．そこでMacPortsからインストール
MacPortsのHowToも参考に．


$ sudo port -d install php5 +pear
$ sudo port -d install php5-mysql
$ sudo port -d install php5-gd

標準のapache2にphp5のモジュールをロードするために，/etc/apache2/httpd.confを編集.


LoadModule php5_module /opt/local/apache2/modules/libphp5.so

MySQL5との接続.


$ cd /opt/local/etc/php5/
$ sudo cp php.ini-development php.ini

MySQL5に接続するために，php.iniを編集


pdo_mysql.default_socket= /opt/local/var/run/mysql5/mysqld.sock
mysql.default_socket = /opt/local/var/run/mysql5/mysqld.sock
mysqli.default_socket = /opt/local/var/run/mysql5/mysqld.sock

Apache2を再起動．


$ sudo apachectl graceful

MacPortsでmysql5をインストール

2009-09-25T23:48:00.001+09:00

mysql5のインストール

sudo port install mysql5

サーバとして使うにはmysql5-serverもインストールする.

データベース格納用ディレクトリの作成

sudo mkdir -p /opt/local/var/db/mysql5
sudo chown mysql /opt/local/var/db/mysql5

プロセスIDファイル用ディレクトリの作成

sudo mkdir -p /opt/local/var/run/mysql5
sudo chown mysql /opt/local/var/run/mysql5

初期化

sudo -u mysql mysql_install_db5

起動

sudo /opt/local/share/mysql5/mysql/mysql.server start

起動の確認

mysql5 -u root

パスワード変更

/opt/local/lib/mysql5/bin/mysqladmin -u root password 'password'

Snow LeopardとMacPorts 1.8

2009-08-29T19:48:00.004+09:00

Snow Leopardを上書きインストールすると, 一見/optが消去されたように見える. ディスクの空き領域も増えている. 実は, /optはFinderで見えなくなっているだけで, 消去されていない. ディスクに余裕ができたのは, システムのPowerPC用の部分がSnow Leopardではなくなったから.

ただし, Snow Leopardと相前後してリリースされたMacPorts 1.8は以前のバージョンと仕様が変わり, デフォルトではCPUに応じたアーキテクチャが選択される. 例えば, Core 2 Duoを使ったMacBook ProやXeonを搭載したMac Proなら, 既定ではx86_64用にコンパイルされる. 1.7以前はi386用にコンパイルされていた. そのため, この例の場合, これまでのものを引き継ぐとアーキテクチャの違いから, リンクがうまくいかなくなる.

64-bit CPUを搭載したマシンで64-bitのコード生成は, 科学技術計算にはメリットがあるので, 歓迎すべき仕様変更である. 手間となるが, コンパイルしなおすのがよさそうだ. Leopardでコンパイルしたものは32-bitだが, ディレクトリの名前を変更すれば使い続けることができる (例えば/opt/localを/opt/local_i386).

MacPortsはボランティアにより支えられているので, 一朝一夕にすべてが移行できるわけではなく, 中にはSnow Leopardでうまくコンパイルできなくなるものもある. しばらく時間がかかるかもしれないが, 徐々に移行されていくはずだ. MacPortsは現行OSとその前までをサボートするので, Tigerはlegacyとなり積極的なサポートはなくなる. TigerユーザはLeopard, Snow Leopardへの移行を迫られるか, コンパイルできなくなるパッケージが出てくるという心づもりが必要だ.

Snow Leopardでは, XCode 3.2となった. XCode 3.2は, Snow LeopardのDVDのOptional Installsにある. デフォルトのコンパイラは, gcc-4.2となったようである.

gnudatalanguage +universal

2009-08-14T15:21:00.003+09:00

gnudatalanguageのPortfileには問題はないが, +universalが失敗することがある.

これは, python25+universalがuniversalとは言っても, i386とppcの2種類のアーキテクチャ用のバイナリをつくるため (r17501). Intel Macでsources.confのuniversal_archsにx86_64が入っていてもi386とppc用に作られてしまう.

octave-3.2

2009-08-12T10:09:00.004+09:00

グラフィックが充実. FLTK/OpenGLを使ったものも実験的に追加されている.

MacPortsのoctaveは, ftgl, fltk, GraphicMagickのほか, fftw-3-single, arpackへの依存関係がなく, 単にバージョンアップしただけの様子. qrupdateはOctaveから分離したようなので, 別途パッケージを作った.

チケット#20648.

netcdfのmaintainer

2009-08-12T10:01:00.002+09:00

MacPortsのnetcdfのmaintainerをすることにした. netcdfに依存するパッケージの大半は自分が管理しているので, netcdfを自由にできた方が都合が良い.

--enable-netcdf-4はvariantにまわして, netcdf3にリンクしていたgrads, grads2, gmtもnetcdfにリンクしてコンパイルできるようにした. netcdf 3は削除.

nco, cdo, wgrib2は, netcdf4にvariantで対応.

TigerでMacPorts

2009-08-10T19:52:00.020+09:00

Mac OS X 10.4 でoctave, gnudatalanguage, ncargのコンパイルに問題があるようだ. Tigerにさらからインストールしてみる.

Ports treeはrsyncでもよいが, 更新が遅いのでここではsubversionを使うことにする. 通常はそのようにする必要はない. rsyncのports treeもしばらくすればアップデートされる.

Tigerにはsubversionは入っていないので, まずsubversionをインストールする.


sudo port -d sync
sudo port -d install subversion

/opt/local/etc/macports/sources.confを編集.


file:///Library/MacPorts/ports [default]

を追加し, rsyncのエントリをコメントアウト.
Ports treeを取得する.


mkdir /Library/MacPorts/
cd /Library/MacPorts
svn checkout http://svn.macports.org/repository/macports/trunk/dports ports

パッケージのインストール. まずはgrads.


sudo port -d install grads

途中でgcc43がビルドされるので時間がかかる. でも, port dependentsには出てこない. どこかでdepends_buildになっているのだろうか.
gradsの動作確認. TigerではX11は手動で起動. 環境変数の設定も必要.


export DISPLAY=local:0.0
export GADDIR=/opt/local/share/grads

LeopardではDISPLAYは設定しない. Xが必要になると, 自動的に起動する. ドックにX11.appをおいておく必要はないし, おくべきではない.

インストーラーでMacPortsをインストールすると, .bash_profileが書き換えられ, PATHに/opt/local/binが追加される. にもかかわらず, Tigerのxtermでは, portコマンドが見つからない. これはX11起動時も, xterm起動時も.bash_profileが読み込まれないためである. X11の「アプリケーション>ターミナル>メニューをカスタマイズ」の「コマンド」にオプション-lsを与える. このようにすると, xtermを起動すると, .bash_profileが読み込まれるようになる. ターミナル (Terminal.app) にはこの問題はない.

ncargは, OPeNDAP対応にするとlibmfhdf.aとlibnc-dap.dylibの_cdf_routine_nameが衝突してしまう. Nclのソースでは使っていないようだが. OPeNDAP対応のTigerバイナリはなさそうなので, うまくいかないのかもしれない. Tigerでは, OPeNDAPはあきらめ, Portfileを変更してコミットした.

今度はnclのビルドが成功し, うまくインストールできるはずである. 遠隔サーバへの接続はできないが, OPeNDAPサーバのデータが全く使えなくなる訳ではない. クイックルックにはgradsのgradsdods, あるいはgrads2のgradsdapが使える. 解析のためには, ncoのncksを使って必要な部分を切り出してダウンロードすればよい.


sudo port -d install ncarg +g95

ncoをインストール.


sudo port -d install nco

GSLをビルドするので時間がかかる.
cdoのコンパイルはすぐに終わるはずである.


sudo port -d install cdo

Octaveは少しPortfileに手を入れた (関連記事). コミットされてからインストールした方がいいだろう. コンパイルにはかなり時間がかかる.


sudo port -d install octave +g95

gnudatalanguageの描画には, plplotが使われている. plplotはswigに依存しているのだが, デフォルトでたくさんのスクリプト言語に対応する. 不要な場合は, 取り除くとコンパイルが早くすむ. 特にphp5はapacheもインストールするので時間がかかる. 既にphp5などがインストールされている場合は, 取り除く必要はない.

gnudatalanguageが依存しているpython25, py25-numarrayをビルド, インストールするがあまり時間はかからない.


sudo port -d install swig -python -perl -ruby -php5
sudo port -d install gnudatalanguage +proj +g95

gradsの代わりにgrads2でもよいかもしれない. アンサンブルに対応しているほか, Google Earthに貼るためのKMLファイルを出力できる. ncarg同様にlibnc-dapとlibmfhdfとで_cdf_routine_nameが衝突していたので, TigerではOPeNDAP非対応とした.


sudo port -d install grads2

grads2のコマンドには, grads-2のように-2がついている. GADDIRはgrads2用にする.


export GADDIR=/opt/local/share/grads2

Exif情報の消去

2009-08-02T15:59:00.002+09:00

Exif情報も一種の個人情報. ImageMagickでExifなどを消すには,

-strip strip image of all profiles and comments

オプションを使う.

convert -strip foo.jpg foo_strip.jpg

みんなの体操祭

2009-08-02T15:06:00.003+09:00

開港150周年記念イベントの一環で横浜赤レンガ倉庫広場で開催．中田市長もあいさつしたが，本番では名前がでなかった．長野先生が忘れたのか，それとも体操しないでかえったのか．ラジオ体操連盟の皆さんは，15分のために全国から集まる.

集中講義予告

2009-07-29T22:50:00.003+09:00

8/5から奈良女子大で「グローバル気象予測入門」と題し, 気象学の基礎やモデルの仕組み，現象論について集中講義をする予定．3日間連続で一日中立っているだけでも, 軟弱な研究者には奈良での修行.

万が一受講予定の学生さんが読んでいたら, 関数電卓 (iPhone等添付のものでも可) を持参してきてください.

ncarg on MacPorts

2009-07-29T22:45:00.003+09:00

nclのソースが公開されて随分経つが, ソースのコンパイルは先延ばしにしてきた. ついに, Portfileを作成してコミットした. ビルドシステムは, configureやCMakeではなく, 独自cshスクリプトとymakeだが案外順調に進んだ.

HDFのヘッダをhdf/mfhdf.hで参照しているので, hdf/を削除するというパッチは必要だった. g2clibはwgrib2のものを使った. hdfeosにはPortfileを書いた. gcc43が既定のfortranコンパイラとなったが, +g95でg95でもコンパイルできる.

BLASのコンパイルに時間をかけているが, Macにはaccelerateフレームワークがあるから無駄なような気がする. でも, 最初は安全重視でncarg添付のものを使った.

花はきりしま

2009-07-29T22:25:00.006+09:00

イージス艦「きりしま」を見に横浜大棧橋へ．見学は30分待ち，と思ったらもれなく砕氷船「しらせ」もついてきた．「しらせ」では飲み物を売っていたし，船内をくまなく見られたのでよかっただが，なかなか「きりしま」に到達できず．「きりしま」はデッキのみの公開．「えんしゅう」の海難救助訓練が見られたのでよしとする．クラゲがうじゃうじゃいたのも印象に残る．

官僚たちの夏

2009-07-29T22:15:00.003+09:00

2009/7/26横浜関内．神奈川県庁を「通商産業省」に見立てて，撮影をしていました．夏なのにコートを着た人が出たり入ったり．入っていく人は，杉本哲太さんらしき俳優さん．合間にはシャツになってスタッフに団扇であおいでもらっています．エキストラさんは年長なのにコートを着たまま．中華街もそばですが，儒教道徳的にいかがなものでしょう．怒ったついでに，公共物である神奈川県の建物を営利会社に撮影させて，通りがかりの県民や大事な観光客のみなさんには，何の権限か警備員が撮影を妨げるのは解せぬ．

KML対応grads2

2009-07-19T16:38:00.003+09:00

MacPortsのgrads2を2.0.a5.oga.5にアップグレード.

このバージョンの特長は, Google Earthに結果を表示するKMLに対応したこと.

MacPortsのgrads2では+geotiffをつける.
ラスター形式のファイル出力に便利な+printimもついでにつけてインストールする例.


sudo port -d install grads2 +geotiff +printim

netcdf4に対応するとのことで, しばらく格闘していた. configure.acなどを解読して, --enable-dyn-supplibsをつけないで静的ライブラリをリンクするようにすると, netcdf4を探してくれることをようやく理解. 現状では, hdf5-18やlibgeotiffは静的ライブラリをインストールしない. しない理由はなさそうだが. 一応コンパイルはできたが, grads-2を起動したディレクトリに変なファイルができてしまうこと, x軸を探せずにsdfopenができなかったことから, 当面netcdf4対応は見送り, 復活させたnetcdf3にリンクすることにした.

これでnetcdf4騒動は一段落.

grads2に含まれるコマンドは-2がつくので注意 (grads-2, gradsdap-2など). grads2が作るwgribはwgrib2ではないwgribだが, バイナリ名はwgrib-2. gradsのwgribと同等だと思われる. wgrib2がインストールするバイナリはwgrib2.

autoconf

2009-07-19T14:21:00.002+09:00

パッケージを管理していると, たまにaclocal.m4が古すぎるから作り直せと言われることがある. また, configureがなぜ期待通りに動作しないことがある. このようなときにautoconfの知識が少し必要になる.

configureはconfigure.acからautoconfを使って作成される. configureは長大なシェルスクリプトで読み解くのが大変だが, configure.acはm4マクロ置換をする前のものなので, 短めである.

aclocalは, configure.acに取り込まれるaclocal.m4を作成する. 同一ディレクトリにあるacinclude.m4は読み込まれるようだ. その他のマクロはaclocalに-Iオプションでディレクトリを渡す. そのアプリケーション独自のマクロは, acinclude.m4やconfigure.acと同一または別のディレクトリ (例: m4）にAC_DEFUNで定義されている.

netcdf 4騒動

2009-07-18T22:24:00.003+09:00

有志でパッケージ管理をしていると, maintainerの中には連絡が取れなくなる人もいる. netcdfはいつの間にかnomaintainerとなっていた. そのまま放置されていただけならまだましだったのだが, どなたたかがnetcdf 4にアップグレードしてくれていた. netcdfに依存するライブラリが軒並みコンパイルできなくなって, 大迷惑だった. とりあえずnetcdf4の機能を分離してdefault variantsに指定, -netcdf4でオフにできるようにした.

netcdf 4は上位互換であるが, HDF5と統合されたので，リンクのときにHDF5のライブラリも指定する必要がある. 必要なライブラリは, pkg-configを使って,


$ nc-config --libs

などととして取得する.

cdo, ncoはnetcdf4に対応していたので, 依存するライブラリの指定やconfigure.argsの変更で済んだ. ただし, ncoはHAVE_NETCDF4_Hがconfig.hに書き込まれないので，陽に指定する必要があった. wgrib2, gnudatalanguage, vis5d+は-lnetcdf以外のリンクするライブラリを指定する必要があった. gmtとgradsは, netcdf3を復活させて, そちらにリンクするようにした. netcdf3は${prefix}/local/lib/netcdf3に「隠して」ある.

grads2はnetcdf4に対応しているはずだが, configureがnetcdf4を検知する部分をスキップしてしまう. Google EarthのKMLファイルが作れるので, 早く提供したいのだがまだ作業が必要なようだ.

netcdfでないけど, ついでに. gnudatalanguageは, aclocal.m4をconfigureの前に再作成して, use_autoconf yesとして, configureを作り直す必要があった. HDF5の古いAPIを使うために, -DH5_USE_16_APIをconfigure.cxxflagsに追加.

GrADSで書いた絵の縞

2009-06-05T12:08:00.002+09:00

GrADSのメタファイルから作成したEPSには縞が出る．

GrADS:


ga-> set gxout shaded

ga-> enable print foo.gx
ga-> print
ga-> dispable print

コマンドライン:


$ gxeps -c -i foo.gx -o foo.eps

これは，アンチエイリアスが原因．

ImageMagickで縞のないラスタ画像を作るには，


convert +antialias foo.eps foo.png

とする．

Mac OS X Server 10.5のMySQL

2009-04-13T13:41:00.002+09:00

Mac OS X Server Leopardには，MySQLのライブラリ等がイントールされていない．

下記から入手できる．Mac OS X Server 10.5.6に対応するのはMySQL-45.binaries.tar.gzで，45はMySQLのバージョンではなくAppleが付けた番号のようである．

http://www.opensource.apple.com/darwinsource/other/

aolserverのインストール

2009-04-11T11:38:00.003+09:00

MacPortsからaolserver-4.5.1をインストール．
aolserverはtclで書かれているので，tclをインストールしておく．
tclは+universalなし+threadsありでインストールする．

ローカル用のソースツリー（例えば/Library/MacPorts/ports_local/www/aolserver）に
aolserverのPortfileをコピー．
ドキュメントをインストールするように修正．


variant doc description {install documentation} {
    destroot.target-append install-docs
}

インデックスを更新し，インストール．


$ cd /Library/MacPorts/ports_local
$ portindex
$ sudo port -d install aolserver +doc

/opt/local/aolserverの下にhtmlとmanができる．
shareの下に置くべきかは分からないので，maintainerに任せる．

aolserverのデーモンはnsdスーパーユーザでは起動できないので，wwwで起動することにする．wwwがログを書き込めるように，所有者を変更する．ログファイルとプロセスIDがlogに作られるので，ディレクトリを作成する．servers以下にサイトが作られる．サイトは複数置ける．あとで，サンプルservers/server1を使ってnsdを起動するが，ここにもログが出力されるので，所有者を変更しておく．


$ sudo mkdir /opt/local/aolserver/log
$ sudo chown /opt/local/aolserver/log www
$ sudo chown -R /opt/local/aolserver/servers/server1

Foreground (-f) でnsdを起動する．


$ sudo bin/nsd -u www -ft base.tcl

IPアドレス:8000にアクセスすると，設定情報が表示される．

CocoaなTk

2009-03-29T21:34:00.004+09:00

Mac版のTkはAqua化されているものと，X11と二通りある．Aquaは，Carbonを使っている．以前のMac OSでもTcl/Tkは動いていたので，これを移植したものだろう．Cocoaを使ったものも開発されているようである．

Tcl and the Tk toolkit second edition

2009-03-29T12:51:00.003+09:00

Amazon.co.jpでTcl and the Tk Toolkit (Addison-Wesley Professional Computing Series)の第2版の予約を受け付けている．2009/7/31発売予定．この本は，Tclの著者によるTcl/Tkの解説本で，Tcl/Tk 8.5に対応する．C言語は一通り知っていることを前提としており，TclとCについて解説していることが特徴．

今の版より安くなるようだ．

C言語でTclを拡張

2009-03-29T12:26:00.004+09:00

Welch, Jones, and HobbsのPractical Programming in Tcl and Tkの第47章には，C言語でTclを拡張する方法が説明されている．

Example 47-2にはバグがあり，古い書き方が残っている．


interp->result = "Usage: random ?range?";

のように直接Tcl_interpのフィールドを操作するのではなく，本文中にあるように


Tcl_SetResult(interp, "Usage: random ?range?", TCL_STATIC);

とした方が安全．

文字列はbufferという名前で宣言しているので，bufは誤り．

リンクオプションは，-bundleでも-dylibでもよい．tclのソースのmacosx/READMEには-dylibのみunload可能とあるが，どうしたらよいかは分からない．

MacPortsのTclでは，libtcl.dylibにリンクするが，Mac OS X添付のTclでは，-framework Tcl -framework Systemが使える．自分でソースからコンパイルするときは，ソースのmacosxで./configureのオプションに--enable-frameworkを付け加える．

ところで，Leopard添付のtclshでは共有ライブラリの読み込みができたが，MacPortsのtcl-8.5.6ではRandom_SafeInit, Random_SafeUnload, Random_Unloadをソースに定義していないとエラーとなる．


Random_SafeInit(Tcl_Interp *interp) {
    return Random_Init(interp);
}
int
Random_Unload(Tcl_Interp *interp, int flags) {
    return TCL_OK;
}
int
Random_SafeUnload(Tcl_Interp *interp, int flags) {
    return TCL_OK;
}

Ubuntuでやってみたが，これらは必要ないので，MacでTcl-8.5を使うときに発生する問題のようだ．

MacPortsのtclsh（tclを+universalでインストール）では，ロード後コマンドを実行したら，mallocでエラーが出た．ActiveTclや/usr/localにインストールしたものは大丈夫だった．MacPortsでも+universal +threadsとしたら，正常に動作するようになった．

FLTK

2009-03-15T20:05:00.007+09:00

FLTKは，高速，軽量なツールキット．クロスプラットフォームで，Macでは，AquaなGUIアプリケーションが作れる．Carbonを使っていて，ウィジェットは独自のものなので，X11でないという程度で見た目はよくない．MacPortsからインストールできる．

ヘッダ二つ，本体3行でウィンドウが出る．


#include <FL/Fl.H>
#include <FL/Fl_Window.H>
int main()
{
 Fl_Window window(300,180);
 window.show();
 return(Fl::run());
}

fltk-config --cxxflagsやfltk-config --ldflagsでフラグを調べて，Makefileに書いておくと楽．


$ fltk-config --cxxflags
-I/opt/local/include -O2 -D_THREAD_SAFE -D_REENTRANT
$ fltk-config --ldflags
-L/opt/local/lib -lfltk -lpthread -framework Carbon -framework ApplicationServices
$ cat Makefile
CXX=c++
CXXFLAGS = -I/opt/local/include -O2 -D_THREAD_SAFE -D_REENTRANT
LDFLAGS = -L/opt/local/lib -lfltk -lpthread -framework Carbon -framework ApplicationServices
TARGET = hello

$(TARGET) : $(TARGET).o
$(CXX) $(LDFLAGS) -o $@ $^

clean :
rm -f *.o $(TARGET)

%.o : %.cxx
$(CXX) $(CXXFLAGS) $< -c

Ercoさんのビデオのチュートリアルでは，この簡単な例から始めて，ウィンドウのタイトルやボタンをつけ，最後にはクラスを使ったオブジェクト指向のプログラムにする．自然な流れが見事．

Interface Builderのように，FLUIDというGUIを使ったプログラミングもできる．

VTK5

2009-03-14T18:40:00.003+09:00

+pythonでビルドし直そうとしたら，PythonのWrappingは構築されていない．

まず，DBUILD_SHARED_LIBS:BOOL=ONが必要だった．そのままだとdyldがライブラリを見つけられない．DYLD_LIBRARY_PATHを指定すると，vtkやvtkpythonコマンドが動作した．Tkのウィンドウがうまく動かないことがあったのは，-DVTK_USE_COCOA:BOOL=ONなのに-DBUILD_SHARED_LIBS:BOOL=OFFだったからかもと思ったが，これが原因ではなかった．

RPATHの指定．otool -Lでリンクされているライブラリのパスを確認できる．Mac OS Xでは，リンク時に-install_nameにライブラリのパスを指定することで，ライブラリに書き込める．vtkのライブラリは，/opt/local/lib/vtk-5.2にあり，install nameがフルパスでなかった．フルパスでなければ，DYLD_LIBRARY_PATHなどから探すことになる．

CMake 2.4の途中からか，CMAKE_INSTALL_NAME_DIRでパスを指定できるようになった．その他RPATHに関する記述やFinkのパッケージを参考にした．

-DCMAKE_BUILD_WITH_INSTALL_RPATH:BOOL=ONはlibvtkCommon.dylibのリンクに失敗する．-DVTK_INSTALL_NO_PYTHON:BOOL=ONとして，post-destrootで手動でsetup.py installするようにした．CMakeを介すると，クウォートのためかうまくsetup.pyの引数が渡らなかったから．ここはFinkとは異なるところ．

+x11もうまく動作していなかった．CMakeはOS添付のTclやTkを見つけてしまう．こちらは，Aqua (Carbon)．x11に限らず，いつでもMacPortsのTcl/Tkを指定するようにした．


-DTCL_INCLUDE_PATH=${prefix}/include \
-DTCL_LIBRARY=${prefix}/lib/libtcl.dylib \
-DTK_INCLUDE_PATH=${prefix}/include \
     -DTK_LIBRARY=${prefix}/lib/libtk.dylib

C++でVTK

2009-03-09T20:03:00.003+09:00

C++にするとTclに比べて，ヘッダのインクルード，変数の宣言，余計な記号たちに加えて，インクルードパスやリンクするライブラリの指定が面倒．

CMakeを使ってみる．CMakeLists.txt


CMAKE_MINIMUM_REQUIRED(VERSION 2.6)
PROJECT(foo)

FIND_PACKAGE(VTK REQUIRED)

INCLUDE(${VTK_USE_FILE})

ADD_EXECUTABLE(foo foo.cxx)
TARGET_LINK_LIBRARIES(foo vtkCommon)

コマンドラインで

cmake .

VTKを使わない場合でも，PROJECTとADD_EXECUTABLEを指定してMakefileが自動生成されるのは便利．他にライブラリvtkRenderingなどが必要なら追加．

C++で動的配列

2009-03-09T19:59:00.004+09:00

Cのmallocを使ってもいいが，C++ではnew演算子を使う．

float* x = new float[n];

配列は，やっぱりFortranが便利．

C++でバイナリファイルを読む

2009-03-09T19:54:00.003+09:00


#include <fstream>

std::ifstream f;

f.open(filename, std::ios::in | std::ios::binary);
f.seekg(pos);
f.read((char*)data, size);
f.close();

openするときの二つ目の引数がポイント．

C++でHello, World!

2009-03-09T19:51:00.006+09:00


#include <iostream>
using namespace std;

int main()
{
    std::cout << "Hello, World!" << std::endl;
}

名前空間std::が必要な今日この頃．

C++の標準数学関数

2009-03-09T19:48:00.006+09:00

C同様


#include <math.h>

でもよいが，C++には

#include <cmath>

がある．

地球を描く

2009-03-08T22:33:00.003+09:00

vtkEarthSourceで海岸線付きの球が描ける．SetOnRatioの既定は10で1が一番細かい．日本の形状は少しまともになるが，北海道，四国，九州が出てこない...


package require vtk

vtkEarthSource earth
    earth SetOnRatio 1

vtkPolyDataWriter writer
    writer SetInput [earth GetOutput]
    writer SetFileName earth.vtk
    writer Write

writer Delete
earth Delete

VTKで球を描く (手作り)

2009-03-08T22:16:00.005+09:00

vtkSphereSourceは便利だが，勉強にならないので，球をセルから構成してみる．

定数


package require vtk

set nlon 100
set nlat [expr {$nlon/2 + 1}]
set pi [expr {acos(-1.)}]
set pi2 [expr {$pi * 2.}]
set npts [expr {$nlon * $nlat}]
#puts stdout "pi=$pi pi2=$pi2 npts=$npts"

経度φと余緯度θからデカルト座標に変換するコマンドを作る．
VTKでは，経度がθ余緯度がφのようだが気にしない．


proc phitheta2xyz args {
         set phi [lindex $args 0]
         set theta [lindex $args 1]
         set xyz {}
         set sin_theta [expr {sin($theta)}]
   lappend xyz [expr {cos($phi) * $sin_theta}]
   lappend xyz [expr {sin($phi) * $sin_theta}]
         lappend xyz [expr {cos($theta)}]
}

経度，余緯度を生成．


set lon {}
set dlon [expr {$pi/$nlon}]
for {set i 0} {$i<$nlon} {incr i} {
    lappend lon [expr {$pi2 * $i / $nlon}]
}
#puts stdout $lon
set lat {}
set dlat [expr {$pi/($nlat-1)}]
for {set j 0} {$j<$nlat} {incr j} {
    lappend lat [expr {$dlat * $j}]
}
#puts stdout $lat

座標点のデカルト座標を生成してvtkPointsに登録．


vtkPoints pts
pts SetNumberOfPoints $npts
set i 0
foreach theta $lat {
    foreach phi $lon {
# puts stdout "$phi $theta"
# puts [phitheta2xyz $phi $theta]
      eval pts SetPoint [concat $i [phitheta2xyz $phi $theta]]
      incr i
    }
}

4点からセルを作る．経度方向に周期境界．$nlon-1番目の隣は，0番目に戻る．


vtkCellArray clls
set nlonm1 [expr {$nlon - 1}]
for {set j 0} {$j<$nlat-1} {incr j} {
    for {set i 0} {$i<$nlon} {incr i} {
        set k0 [expr {$nlon*$j + $i}]
        if {$i!=$nlonm1} {
            set k1 [expr {$k0 + 1}]
        } else {
            set k1 [expr {$k0 - $nlonm1}]
        }
        set k2 [expr {$k1 + $nlon}]
        set k3 [expr {$k0 + $nlon}]
        clls InsertNextCell 4
        clls InsertCellPoint $k0
        clls InsertCellPoint $k1
        clls InsertCellPoint $k2
        clls InsertCellPoint $k3
    }
}

点とセルをvtkPolyDataに登録．


vtkPolyData sfc
    sfc SetPoints pts
    sfc SetPolys clls
    pts Delete
    clls Delete

ファイルに保存して完了．


vtkPolyDataWriter writer
   writer SetInput sfc
   writer SetFileName sphere2.vtk
   writer Write
   writer Delete
sfc Delete

VTK形式のファイルを描画

2009-03-08T22:12:00.004+09:00

vtkinteractionを使う一部のサンプルでBus errorやSegmentation Faultが出る．VTK-5.2.1に更新したがだめだった．仕方がないので，単に表示してしばらくしたら終了するようにした．仮に動作しても，表示はParaViewの方がおそらく便利．


#!/bin/sh
# ¥
exec vtk "$0" ${1+"$@"}
package require vtk

if {$argc<1} {
    puts stderr "vtk draw.tcl filename"
    exit
}

vtkPolyDataReader reader
    reader SetFileName [lindex $argv 0]

vtkPolyDataMapper sphereMapper
    sphereMapper SetInputConnection [reader GetOutputPort]
    reader Delete

vtkActor sphereActor
    sphereActor SetMapper sphereMapper

vtkRenderer ren1
    ren1 AddActor sphereActor

vtkRenderWindow renWin
    renWin AddRenderer ren1
    renWin SetSize 300 300
    renWin Render

after 1000

vtkCommand DeleteAllObjects

VTKで球を描く (vtkSphereSource)

2009-03-08T22:07:00.003+09:00

MacPortsでvtk5をインストール．デフォルトは，Cocoa．X11ではコンパイルに失敗した．まずは，お手軽にTclで球を描いてファイルに保存．次のスクリプトを適当な名前（例えばsphere1.tcl）で保存して，実行（vtk sphere1.tcl）．


package require vtk

vtkSphereSource sphere
        sphere SetRadius 1.0
        sphere SetPhiResolution 100
        sphere SetThetaResolution 50

vtkPolyDataWriter writer
        writer SetInput [sphere GetOutput]
        writer SetFileName sphere1.vtk
        writer Write

sphere Delete
writer Delete

SPHEREPACK

2009-02-23T23:40:00.001+09:00

07/10/06 19:51

SPHEREPACKはLegendre変換をするFortranライブラリ.
Mac OS X上でg95を使ってコンパイルして使っている。

解像度変換のユーティリティルーチンtrssph()が配列サイズが大きいときにNaNを返すことがある。試した範囲では、288×145を1200×600にするのは大丈夫で、1440×720にするのはだめだった。

メモリ節約のため、trssph()は、Legendre陪関数を逐次計算するサブルーチンsphcom.fにあるlegin()でルジャンドル陪関数が正しく計算されないようである。アルゴリズムとしては、robustな方法を使っているはずだし、ベクトルのLegendre変換をするルーチンはうまく動作する。また、SPHEREPACKを呼んでいるNCLでは、正しく動作する。

ソースを見たところ、ルジャンドル陪関数を計算するときに必要な係数abel, bbel, cbelが大きな整数のかけ算になっていたためであることが分かった。shagc.f, shsgc.f, shags.f, shsgs.f, shigc.f, shigs.fでは


-      abel(imn)=sqrt(float((2*n+1)*(m+n-2)*(m+n-3))/
-     1               float(((2*n-3)*(m+n-1)*(m+n))))
-      bbel(imn)=sqrt(float((2*n+1)*(n-m-1)*(n-m))/
-     1               float(((2*n-3)*(m+n-1)*(m+n))))
-      cbel(imn)=sqrt(float((n-m+1)*(n-m+2))/
-     1               float(((n+m-1)*(n+m))))
+      abel(imn)=sqrt((2.*n+1.)*(m+n-2.)*(m+n-3.)/
+     1               ((2.*n-3.)*(m+n-1.)*(m+n)))
+      bbel(imn)=sqrt((2.*n+1.)*(n-m-1.)*(n-m)/
+     1               ((2.*n-3.)*(m+n-1.)*(m+n)))
+      cbel(imn)=sqrt((n-m+1.)*(n-m+2.)/
+     1               ((n+m-1.)*(n+m)))

と修正 (diff -uの出力)したところ、問題なく動作した。またvhags.fとvhsgs.fは以下の通り修正した。


-      abel = dsqrt(dble(float((2*n+1)*(m+n-2)*(m+n-3)))/
-     1                dble(float((2*n-3)*(m+n-1)*(m+n))))
-      bbel = dsqrt(dble(float((2*n+1)*(n-m-1)*(n-m)))/
-     1                dble(float((2*n-3)*(m+n-1)*(m+n))))
-      cbel = dsqrt(dble(float((n-m+1)*(n-m+2)))/
-     1                dble(float((m+n-1)*(m+n))))
+      abel = dsqrt((2.d0*n+1.d0)*(m+n-2.d0)*(m+n-3.d0)/
+     1                ((2.d0*n-3.d0)*(m+n-1.d0)*(m+n)))
+      bbel = dsqrt((2.d0*n+1.d0)*(n-m-1.d0)*dble(n-m)/
+     1                ((2.d0*n-3.d0)*(m+n-1.d0)*(m+n)))
+      cbel = dsqrt((n-m+1.d0)*(n-m+2.d0)/
+     1 ((m+n-1.d0)*(m+n)))


-      abel = dsqrt(dble(float((n+m)*(n-m+1))))/dcf
-      bbel = dsqrt(dble(float((n-m)*(n+m+1))))/dcf
+      abel = dsqrt((n+m)*(n-m+1.d0))/dcf
+      bbel = dsqrt((n-m)*(n+m+1.d0))/dcf

2009/2/23 追記

パッチにバグがありました．vhags.fとvhsgs.fの修正の前半，abel, bbel, cbelの分母全体の括弧がありませんでしたので訂正しました．

D 2.025 for Mac OS X

2009-02-22T20:58:00.004+09:00

D 2.025には，Mac OS X版が含まれている．D言語を開発したDigital Mars社純正で，Intel Macで動作する．writeln()等D 2.0で追加された機能が使える．ChageLogを見ると2009/2/14に公開された2.025からMac OS X版が追加されたようだ．マニュアルrdmd.1はあるが，rdmdはない．

2009/2/24 追記

#!/usr/bin/dmd -run

と書けば良いようだ．スクリプト風に動くというのは面白い．

私の場合は，/usr/local/dmdというフォルダを作り，bin, libをコピー．shareにhtmlとmanをコピー．ランタイムはsrc以下にあるので，これもコピー．dmdは/etc/dmd.confの設定を読む．/usr/local/dmd/etc/dmd.confからシンボリックリンクを張った． dmd.confには，次のように書いた．


DFLAGS=-I/usr/local/dmd/src/phobos -I/usr/local/dmd/src/druntime/import -L-L/usr/local/dmd/lib

gdcで変なメッセージが出ることに対応するパッチは，MacPortsのリポジトリに登録されたようだ．

plplot-5.9.2 revision 1

2009-02-22T19:06:00.003+09:00

D言語のbindingとcairo/pangoドライバをvariantsとして追加．cairoドライバによりPDFが作成できるようになる．g95, gdc, cairoを有効にしてインストールするには，

sudo port -d install plplot +g95 +gdc +cairo

とする．

D言語サンプルは，bindingがexperimentalなので少ないが，/opt/local/share/plplot-5.9.2/examples/dにいくつかある．コンパイルは，ヘッダ/opt/local/include/plplot/plplot.d，モジュール本体は/opt/local/lib/libplplotd.dylibなので，

gdmd -I/opt/local/include/plplot x03d.d /opt/local/lib/libplplotd.dylib

とする．gdcでなくgdmdを使うとソースコードから.dをとった名前のバイナリができる．簡単なプログラムでは便利．-lオプションはないようだ．

D言語

2009-02-22T16:44:00.004+09:00

Windowsの.NETやLinuxやMacでも.NET互換の環境Monoの開発言語として，C#が注目されている．言語仕様は，よりよいC, Javaになっている．多次元配列もある．しかし，中間言語に翻訳するためか，Cの函数を直接呼べず，動的ライブラリをロードする必要がある．このあたりは不便であるように思われた．

plplotに対応している言語を見ていて，D言語があるのを見つけた．聞いたことはあったが，仕様はよく知らなかったので，調べてみた．「わかったつもりになるD言語」にはドキュメントの翻訳を含めて，詳しく紹介されていた．コンパイル型言語で，Cの函数をラッパー無しで直接呼べる．複素数型，動的に伸張配列可能な配列，文字列がUnicode文字の配列あること，標準ライブラリに正規表現が含まれることなどが大変便利に思えた．

MacPortsにもGNUのDコンパイラgdcがあったので，インストールしてみた．Leopardでは変なメッセージが出るが問題なく使える．かつてg95にもあり既に解決済なので，Portfileのパッチを作り，Ticketを発行して添付しておいた．

GNU Dも仕様通りUnicodeに対応しているので，Hello, worldは日本語で試せた．gdcはDMD 1.0仕様なので，writelnはないことに注意．

ちなみに，NHKの地上波の朝のニュースは「おはよう日本」，BSは「おはよう世界」である．

ソースはHTMLでもよいというのも面白い．の間のソースをコンパイルする．の間にタグが含まれていても無視するので問題ない．

Vis5D+

2009-02-15T14:35:00.001+09:00

MacPortsにVis5D+を追加．以前にインストールメモを見てコンパイルしたことがあった．Fortranのコンパイルがうまくいかなくて困っていたが，configureスクリプトを改変すればg95でコンパイルが通るようになった．gcc43では，うまくいかない．Fortranコンパイラは使われておらず，Fortran用のインターフェースをコンパイルするかどうかだけ．+g95をつければ，-lvis5dをリンクしてFortranでvis5d+形式のデータを作成するツールが作れるはず．GTK1は，
うまくいかなかった．

モードライン

2009-02-07T09:15:00.003+09:00

Portfileの最初の行におまじないのような一行がある．

# -*- coding: utf-8; mode: tcl; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*- vim:fenc=utf-8:ft=tcl:et:sw=4:ts=4:sts=4

EmacsとVim用にソースの種類やタブの取り扱いを設定するためのもので，モードラインと
呼ばれている．

Portfileの書き方として推奨されているので，よく考えずにこの行をつけていた．

モードラインの設定に従うはずなのに，この行があるにもかかわらず，~/.vimrcの設定が使われるので，モードラインが有効になっていないことに気づいた．

モードラインを有効にするには，

set modeline
set modelines=1

とする．いつも使うのなら，~/.vimrcに書いておけば良い．モードラインはset modelineで有効になる．モードラインの設定を検索する行数は，set modelinesで指定する．既定では，モードラインは無効でmodelines=0のようであるので，モードラインを使うにはこの設定が必要である．

後半にあるvimの部分では，テキストのエンコーディングをutf-8に（fenc=utf-8），ファイルタイプをtclに（ft=tcl），タブを展開（et），インデント，タブ，ソフトタブをそれぞれ4文字に設定している．ソフトタブに0でない値を設定すると，タブを入力したとき設定した文字数の空白に置換する．

ソースの種類により，スタイルは異なっている．Makefileはタブが必要なので，いつもタブの展開（set expandtab）しておくわけにはいかない．Cなど通常のソースでは，タブが4つの空白にすることが推奨されることが多い．長い数式を書き，doループが多重になることが多いFortranのソースでは，タブは2つくらいにしておきたい．このような場合でも，ソースの最初にモードラインを書いておけば，様々なスタイルに対応することができる．

Portfileは拡張子がないが，モードラインの設定のおかげで，文法に基づいた色づけ（syntax enable）を~/.vimrc等に設定しておけば，キーワード等に色がつく．インデント等の設定があるために，文字を揃えるのが容易になった．大変便利だ．Fortranやnclのソースにも適宜モードラインを設定して，見やすいソースを書くようにしたい．

gdl-0.9r2, GMT-4.3.1

2009-02-06T19:54:00.003+09:00

GDLとGMTのMacPortsパッケージを更新．universalにも対応．GDLが依存しているhdf5のuniversal対応は，管理者に連絡中．GDLのlibtoolが古かったので，/opt/local/share/libtool/config/ltmain.shをコピー．

Eigen2

2009-02-01T18:27:00.004+09:00

Eigen2は，配列を扱うためのC++テンプレートライブラリ．MacResearchに掲載されたインタビューによると，SSE2やAltivecにも対応していて，ATLASに劣らない性能が出るそうだ．データ解析とお絵描きには使えるかもしれない．MacPortsで簡単にインストールすることができる．MacPortsで入れた場合，コンパイルは

c++ -I/opt/local/include/eigen2 foo.cc -o foo

とすればよい．

GrADS 2

2009-01-18T14:05:00.002+09:00

OpenGrADSプロジェクトのArlindo da Silvaさんのすすめがあり，OPeNDAP，GRIB2に対応したGrADS 2のパッケージを作成．GRIB2のライブラリ・ユーティリティとOPeNDAPをステーションデータに対応させるためのgadapが必要なので，これも別のパッケージとして作成．

インストールするときに+universalとすると，複数のアーキテクチャに対応するバイナリが生成される．既定ではi386とppc，つまりIntel MacとPower Macで動作する．/opt/local/etc/macports/macports.confのuniversal_archsを編集すると，ほかのアーキテクチャ用にもコンパイルできる．例えば，ppc64とすればPowerPC G5用，x86_64とすればAMD64用（AMDのプロセッサを搭載したMacはないので，これに互換のIntelのプロセッサ）の64ビットのバイナリになる．

GrADS2自体のコンパイルは簡単だったが，上記のふたつのライブラリ用のパッケージやlibwww, hdf4, libdap, libnc-dapのPortfileをいじってuniversalに対応させたるのに手間がかかった．

Windows 7 beta

2009-01-12T19:09:00.002+09:00

自作PCは64ビットなので，できれば64ビットOSがよいのですが，32ビットでないといろいろと問題がおきるようなので，XPは32ビットでした．64ビットOSで遊んでみるためにWindows 7はもってこいです．

2つ目の内蔵ディスクにインストール．XPでダウンロードしてImgBurnでDVDに焼きます．BIOSではDVDドライブがHDDに優先しているのですが，DVDから起動したらすぐになにかキーを押さないとインストーラは起動しません．Windowsユーザには常識なのかもしれませんが．

ファイルのロードとSetupの起動はハングしたのではないかと思うほど待たされますが，気長に待ちましょう．Setupの途中でうっかり「ドライバのロード」を選ぶと戻れなくなり，振り出しに．ファイルの展開とコピーが終わると再起動になります．LinuxとXPのデュアルブートになっていたので，心配してF8を押してしまいましたが，これは誤りでインストールやり直しになりました．インストールしたHDDから起動してもだめでした．

Windows 7の画面はXPよりもきれい．問題なく解像度も最大に．Windows UpdateにATI X1200シリーズのアップデータがあり早速インストール．このあたりはXPより改善されています．

Bootcampまたは仮想化環境でWindowsを使いたいMacユーザでWindows 7を試してみようというという方もいるかと思います．2009/8/1まで無料で使えるのでよいのですが，Yahoo!動画はだめでした．そのうち対応してくれるとよいのですが．

MacユーザにやさしいWindows XP

2009-01-12T18:43:00.004+09:00

Windows 7が出たところですが，Macユーザにも使いやすいようにWindows XPを整備してみました．Yahoo!動画だけではもったいないので．インストールも何度も再起動して大変ですが，使える環境にするのも結構大変です．UIを似せることにはこだわらず，シンプルにいきます．

Microsoft

Windows update
Windows XP Service Pack 3
DirectX 9.0c
.NET Framework 2.0, Service Pack 1, 言語パック
Internet Explorer 7
Windows Media Player 11
メイリオ
MSゴシック，MS明朝JIS2004
Windowsサーチ4.0

メイリオを入れてシステムやブラウザで指定すると少しは見た目がきれいになります．

設定はお好みで．

サウンドなし
ウィンドウのデザイン，スタートメニュー　クラッシック表示
背景単色

ビデオとチップセットドライバ

Catalyst 8.12

環境に合わせて入れます．プリンタドライバなども．

キーボード

xkeymacs ctrlとcaps lockとを入れ替え．Emacs風のキーにする機能はオフ
秀Caps 半角/全角をEscに，無変換をIMEオフに変換をIMEオンに，Shiftを押さずにアンダースコア
MS IMEはVJE風のキーバインド

キー配列が違うとかなりストレスがたまりますので，このあたりは重要．

解凍ユーティリティ

Lhaplus

何が定番かはわかりません．ところで，Macだからといってsitはそろそろやめてほしい．

Apple

QuickTime
iTunes
Safari
Bonjour for Windows

Safariが入ると，かなりMacらしくなります．

Adobe

Acrobat Reader 9
Flash Player

Macにはプレビューが入っていますが，Windows XPにはPDF閲覧ソフトウェアはありません．

アプリケーションをいくつか

OpenOffice.org 3.0
Picasa2
ImgBurn

OpenOffice 3はよくできています．でも数式はやっぱりTeX．

TeX環境

TeXインストーラ3
ImageMagick
gvim-7.1
TeX2img

開発環境，科学技術アプリケーションは今なら仮想環境かデュアルブートでLinuxが多いのでしょうか．まだCygWinやMinGWも便利かも．

やっぱりMacのほうが楽です．

やっぱりFortran

2008-12-23T17:54:00.003+09:00

HIMに触発されて，C言語を数値計算に使えないかしばらく考えていた．

C言語の文法をおさらい．さまざまな入門書があり，それとは別にポインタの解説が売られている．代入やら繰り返し，変数のような基礎概念を解説しなければならないので，ポインタについては簡単に済ますしかないのかもしれないが，肝心なところをごまかして，2冊目を売ろうというしているのではないかとも思える．プログラミングの基礎は，スクリプト言語でしておくとよいかもしれない．やはり，K&Rのプログラミング言語C ANSI規格準拠が要領よくまとまっているように思える．Cの基本的な文法自体は難しくなくい．むしろBetter CとされるC++やJavaの方が覚えることが多い．

現在のFortranプログラミングで重要なのは，モジュールと行列．Cはこのあたりが弱い．Cは大域変数の名前空間は基本的にひとつ．ヘッダファイルにマクロ定義をしたり，函数を定義したりして，必要なときにインクルードするというスタイルになる．Fortranだと，useのように陽に指定できる．配列については，多重配列は配列の配列しかないし，添字は0からに限られ，行列式での代入のような便利な機能はない．ポインタとの関係や値渡しの引数には注意が必要．スパコンのコンパイラでは，ポインタは最適化を阻害する要因になりうる．

UNIXワークステーションによる科学技術計算ハンドブック―基礎篇C言語版のような数値計算の本は，数値計算のアルゴリズムが中心に解説されている．それはそれでよいのだが，よくできているものはライブラリを使うのが普通で，データをいかに処理するかが数値計算の中心．配列のことやモジュール化をどうするか，I/Oについてもっとページを割くべきであるように思う．

Cは簡単な函数を書くにはよさそうだ．Octaveなどで時間がかかる部分を高速化するのに役立つ．グラフィックやシステムの機能を使うプログラムにもよいだろう．しかし，HIMのような大きなものを書くにはCに慣れておかしなことをしないように注意しないといけない．通常は，文法がやさしく，数値計算，データ処理に便利な機能がそろっているFortranを使うのがよいように思う．

Xcodeでmkoctfile

2008-12-21T15:57:00.010+09:00

XcodeのビルドシステムはGNU Makeではないが，少しカスタマイズすれば，任意のコンパイラやスクリプトを実行することが可能である．ここでは，Octaveのライブラリにリンクするスタンドアロンのバイナリを作成してみる．

まず，Xcodeを起動しプロジェクトを作成する．プロジェクトの種類が並んでいる左側からCommand Line Utilityを選ぶ．右側からStandard Toolを選ぶ．プロジェクト名はhellooとして適当な場所に保存する．プロジェクト名と同じ名前のディレクトリが自動的に作成されて，プロジェクトファイルhelloo.xcodeprojやサンプルソースがこのディレクトリに作成される．

main.cやhelloo.1は不要なので，プロジェクトのウィンドウで選択して削除する．ファイルメニューから新規ファイルを選択してhelloo.ccというファイルを作成する．テンプレートはその他から空のファイルを選ぶ．このファイルにOctaveのマニュアルにある例をコピーして保存．この例にはバグがあるのだが，ここではそのままにしておく．

プロジェクトウィンドウ左側のターゲット，さらにhellooも開く．3段階のフェーズのうち，ソースのコンパイル以外は不要なので削除する．次に，ソースのコンパイルをカスタマイズしてmkoctfileを呼ぶようにする．ターゲット>hellooを選択した状態でボタンバーの情報をクリック．Command+iでもよい．ルールタブをクリックして表示したら，左下の+ボタンをクリックしてルールを追加する．プロセスは同名のファイルとして，*.ccと入力する．使用はカスタムスクリプトとし

/opt/local/bin/mkoctfile --link-stand-alone ${INPUT_FILE_DIR}/${INPUT_FILE_NAME} -o ${TARGET_BUILD_DIR}/${INPUT_FILE_BASE}

と入力する．出力ファイルを追加し

${TARGET_BUILD_DIR}/${INPUT_FILE_BASE}

と入力する．ここで用いた環境変数は，右下のボタン?をクリックして表示されるマニュアルに説明されている．

ビルドをクリックしてビルドしてみる．二重forループの中に

a_matrix(row,column)

となっていればエラー箇所として指摘されているはずである．row, columnは定義されておらず，おそらくi, jの誤りだろう．ここを修正すればコンパイルが通るはずである．

ビルドして実行ボタンをクリックしても，コマンドラインツールなので，何も表示されない．結果はツールからコンソールを選べば表示できる．プロジェクトウィンドウのProductの下に表示されるhellooをダブルクリックするとターミナルが開いてバイナリが実行される．

mkoctfileは，シェルスクリプトで，コンパイルにはg++が呼ばれている．Xcodeでは，g++のエラーメッセージを解釈して，エディタに該当箇所を表示できるので，効率的な開発役立ちそうである．MPIやOpenGLのプログラムは，同様な方法でmpiccやglccを呼べばよさそうである．シェルスクリプトが面倒を見てくれるので，ヘッダやライブラリをプロジェクトに追加する手間が不要だ．C, C++のラッパ以外でも同様の方法で良い．ただし，g95, gfortranだと，エラーメッセージは正しく解釈されない．エラーメッセージは，ソースの上に表示される．

MacPortsの大掃除

2008-12-21T09:08:00.007+09:00

MacPortsのパッケージは，コミッタ，パッケージ管理者そしてユーザ有志の努力により，絶えず更新が続いている．インストールされているものがパッケージよりも古くなったものの一覧は

port outdated

とすると表示される．

MacPortsのパッケージを更新するときのコマンドはport upgradeである．特定のパッケージfooを更新するには，

sudo port upgrade foo

とする．すべてを更新するには，

sudo port upgrade outdated

とする．upgradeをする前に

sudo port sync

をして，ソースツリーを更新しておく．

更新したパッケージやそれに依存するパッケージなど旧バージョンはdeactivateされ消されずに残る．更新をするたびに「ごみ」がたまっていく．既定で古いものを残すようにしているのは，新しいパッケージが動くことを確認してから古いものを消せばよいからである．

port -u uninstall

とすると「ごみ」をすべて清掃してくれる．古いものを残す必要がないときは，

port -u upgrade foo

とすると旧バージョンは消去される．更新対象のfooは，最上位に位置する「大物」だと依存するライブラリも更新され，-u付きの場合は古いものが消去されるので楽である．

ACML

2008-12-20T18:18:00.003+09:00

Athlon X2 5400+を載せた自作PCのベンチマークを改善しようと，AMD Core Math Libraryをインストールし，Octaveをコンパイルしてみた．Linux用でgfortranでコンパイルしたものは，intが32-bitのものと，64-bitのものの2種類あった．どちらも取得してインストールした．取得には登録が必要だった．インストールは添付のシェルスクリプトで行う．OpenMPに対応したライブラリが含まれている．

PowerPC G4のMac mini以下の性能だったが，atlas-3gf-baseを使うと3.5Gflopsくらい出た．ACMLはもう少し性能がよく5Gflopsでたこともあった．CFLAGSやCXXFLAGSに-fopenmpを指定し，マルチプロセッサ用のライブラリとリンクしてみたが，Mac mini以下に逆戻りしてしまった．gccのOpenMPのできが問題なのだろうか．Athlon X2 5400+は，2コア使うとPowerPC G5やCore Duoと勝負できそうだが，クロックの近いCore 2 Duoには劣る．

OctaveをMacPortsから普通に入れてもLeopard添付のAccelerate.frameworkはスレッドを立ち上げてコアをすべて使ってくれているようである．手軽に性能を出すならMac Proが楽なのではないだろうか．

HIM

2008-12-14T22:06:00.006+09:00

C言語で書かれた海洋モデルHIMを試してみたところ，3つの例をうまく動かすことができた．アメリカ気象学会のページで検索してみたところ，100を超える引用がある．鉛直座標に密度を用いているisopycnalモデルだけあって，とくに浮力により生ずる流れに関する研究に使われている．Park and Bryan (2001)は密度座標のモデルと高度座標とを比較して，亜寒帯ジャイアに差異が見られ，大気モデルと結合した場合に大気海洋相互作用に影響が表れることが示唆されると述べている．

Mac OS Xではnetcdfさえインストールすれば，LeopardのOpen MPIを使って並列計算をすることができる．LeopardのOpen MPIはFortranに対応していないので，Fortranでかかれたプログラムの場合は，MacPortsなどからOpen MPIをインストールする必要がある．OSのOpen MPIは32/64ビットの両方，PPC/Intelの両方，計4つのアーキテクチャに対応し，universalである．

netcdfを64ビットでインストールし，HIMをコンパイルするときのオプションに-m64をつければ64ビットの実行ファイルができる．64ビットなら，4GBを超えるメモリを必要とする場合でも動作するはずである．MacPortsを64-bitでコンパイルする方法はここに書いてある．MacPorts 1.6まではuniversalとはi386とppcを差していたが，MacPorts 1.7からはi386, x86_64, ppc, ppc64のいずれをユーザが指定するのだそうである．MPIを使うときは，CC=/usr/bin/mpiccをMakefileに書く．Mac Proで8コアをフルに使うと，シングルでは少し待った計算がすぐに終わり，並列計算の威力を体感できた．

C言語で実用になるモデルがあることに感心した．数値計算の分野，とくにスパコンではFortranが有利だと言われているが，パソコンではCの方が便利かもしれない．最近はg95やgfortranがでてきているが，可視化など数値計算以外のものとのリンクはCの方が便利である．数値計算のライブラリについても，GSLやFFTWはCで書かれている．Octaveを使えば自分でベクトルや行列及びその演算を定義しなくても演算や様々な函数が使える．

Octaveのライブラリを使ってスタンドアロンのプログラムを作る方法は，ドキュメントに書いてあるが，サンプルプログラム中a_matrix(row,column)はa_matrix(i,j)の誤りであるようだ．C++のためか，Octaveの初期化のためか，起動のオーバヘッドがあるように感じた．Octaveから使える函数を書く方法も簡単であるようだ．

しばらく，CやOctaveのライブラリを使ったプログラムで遊んでみてもいいかもしれない．

uptex

2008-12-14T21:59:00.005+09:00

JISで行われている内部処理をUTF-8にしたuptexを試してみた．JISのmapファイルはインストールしていないが，UTF-8のテキストはuptex/uplatexとdvpidmxでうまくタイプセットできているようだ．まだ自信はない．以下インストール時のメモ．

cbfonts.zip のディレクトリ構造が変わった．
ptetex3-cmapをインストールする必要がある．
my_option中export CFLAGS="-O2 -DSTDC_HEADERS"とする必要があった．

--- uptex-0.24/patch/0uptex.sh    2008-10-17 21:29:24.000000000 +0900
+++ ptetex3-20080616/0uptex.sh    2008-12-14 11:45:49.000000000 +0900
@@ -114,20 +114,20 @@
  FONT_DIR=$UPTEX_TEXMF_DIR/fonts
  unzip $SRC_DIR/cbfonts.zip -d $PTETEX_TMP
  $MKDIR $FONT_DIR/type1/cbfonts
-  mv $PTETEX_TMP/cbfonts/type1/*.pfb $FONT_DIR/type1/cbfonts
+  mv $PTETEX_TMP/cbfonts/fonts/type1/cbgreek/*.pfb $FONT_DIR/type1/cbfonts
  $MKDIR $FONT_DIR/tfm/cbfonts
-  mv $PTETEX_TMP/cbfonts/tfm/*.tfm $FONT_DIR/tfm/cbfonts
+  mv $PTETEX_TMP/cbfonts/fonts/tfm/cbgreek/*.tfm $FONT_DIR/tfm/cbfonts
  $MKDIR $FONT_DIR/source/cbfonts
-  mv $PTETEX_TMP/cbfonts/mf/*.mf $FONT_DIR/source/cbfonts
+  mv $PTETEX_TMP/cbfonts/fonts/source/cbgreek/*.mf $FONT_DIR/source/cbfonts
  $MKDIR $FONT_DIR/enc/cbfonts
-  mv $PTETEX_TMP/cbfonts/encodings/*.enc $FONT_DIR/enc/cbfonts
+  mv $PTETEX_TMP/cbfonts/fonts/enc/*.enc $FONT_DIR/enc/cbfonts
  $MKDIR $FONT_DIR/map/dvips/cbfonts
-  mv $PTETEX_TMP/cbfonts/dvips/*.map $FONT_DIR/map/dvips/cbfonts
+  mv $PTETEX_TMP/cbfonts/fonts/map/*.map $FONT_DIR/map/dvips/cbfonts/
  $RM -r $PTETEX_TMP/cbfonts
cat $TMP_PREFIX/share/texmf/web2c/updmap.cfg || exit

# cbfonts
-MixedMap cbgreek.map
+MixedMap cbgreek-full.map
EOF
fi

加湿器+温風

2008-11-30T16:02:00.003+09:00

SANYOのAQUAWARMは加湿器とセラミックヒータを兼ねた商品. セラミックヒータはスイッチを入れてすぐに温風が出ることが利点.同時に加湿できるので，体感温度が上がる.加湿は気化式なので，電気代は少なくてすむ.空気清浄器ではないが，空気を除菌する機能もある.ポイント暖房として購入したが，加湿機能のためか意外と部屋もよく暖まる.

2010-11-20追記: 2010年バージョンは, RSF-VW13Dのようだ. 前から水の残量が見えやすくなっている. シルバーに加えて, ブルーが追加されている.

小売業界下克上

2008-11-16T23:32:00.003+09:00

かつては百貨店はカメラ量販店より格上だった．一昔前から，百貨店が量販店を上の方の階に入れて集客力を高めることが始まった．例えば，新宿のタカシマヤはベスト電器を上大岡の京急百貨店はヨドバシカメラを店子としてレストラン街の下あたりに入れている．最近は，量販店が百貨店と入れ替わるということが起きている．横浜のヨドバシカメラは三越の跡に入った．
その結果何が起きたのか．今度は，ユニクロが客を上階にあげる役割を担ってヨドバシカメラの店子となっている．ユニクロはかつての百貨店の商品である洋服を売る店．大家と店子が入れ替わっている．そのほか百貨店の商品として思い浮かぶ，スーツのビッグビジョンや靴の店歩人館もある．それぞれ既存の百貨店の売り方とは異なっている．小売りの仕方が10年規模で絶えず転換していることが感じられる．
そういえば三越に来る前のヨドバシ跡地はドンキホーテだ．整然と商品が売られている家電量販店とはドンキホーテは対照的で，数十年前日本でも見られ，ソウルの市場や中東のバサールのような宝探しのような雰囲気がある．今後はここに回帰していくのかもしれない．

浄水器を新調

2008-11-16T23:03:00.005+09:00

某量販店で駐車券をもらいたいばっかりに約4,500円のトレビーノのカートリッジを買おうとしたところにメーカーから来たと思われる販売員に言葉巧みにクリンスイを勧められた．「浄水器は何年くらいお使いでしょうか．」本当は6年以上使っているのだか，控えめに「5年くらいです」と答えたのだが，そろそろかえた方が良いような心理状態に追い込まれた．すかさず，クリンスイのカートリッジが二つで約7000円と割安で性能も同等だと言ってきた．鉛，トリハロメタン，農薬も濾過するタイプ (スーパーハイグレード) が付属し，浄水能力もカートリッジのもちも同じだった．「ようやく追いついたということですか」と切り返すのが精一杯で，立ち去るときには，「交換が面倒だなぁ」と思いつつクリンスイを手にしていた．プリンタのような消耗品ビジネスなので，本体は安いので乗り換えの障壁は交換の手間だけである．
トレビーノを外してみると，ネジ山などが汚れていた．掃除すればよいのだが，やはり交換するといった機会でないと気がつかない．トレビーノは専用の工具がついていたが，クリンスイは硬貨を使うように説明書にある．不親切なようだが滅多に使わないのだから，合理的とも言える．作業は5分程度で完了．
使い勝手は微妙に異なる．トレビーノはレバー式でクリンスイはダイヤル式．後者は切り替えにはつまむ必要があるが，まわしやすい．トレビーノはフタをあけないと汚れの状態が分からないものだったが，クリンスイは窓がついている．交換を促す意味ではメーカーにとって消費者が気がつきやすい方がよい．フタはよっぽどのことがない限りあけないもので，トレビーノは水がなかなかでなくなるまで使い続けたことがあった．トレビーノの浄水はシャワーだったが，クリンスイはシャワーにならない．それでもいいと思った．水道水にあるにおいや味が味気ないくらいに除去されている．BRITAの方がまろやかな感じがする．まぁ良い買い物だったと思うが，amazonではさらに安かったのであまり得をした気分にはなれない．