スペクトル法で時間発展拡散方程式(1D)

勉強のためやってみたので、Octaveのスクリプトを。
もちろん、間違ってても責任は持ちません。

function laplacian = calc_lap(nDim, dif_coef, cur_func)
   for i=1:1:nDim
      laplacian(i) = -1.0*dif_coef*i^2*cur_func(i);
   end
endfunction
nGrid=128;
nSteps=2000;
deltaT=0.0001;
dif_coef = 1.0;
#Define initial condition
initial_value=zeros(1,nGrid);
initial_value(1,nGrid/2)=1.0;
init_fft=fft(initial_value);
#Begin time integral loop
cur_coef=zeros(1,nGrid);
hlf_coef=zeros(1,nGrid);
nxt_coef=zeros(1,nGrid);
lap_coef=zeros(1,nGrid);
cur_coef = init_fft;
#hold on;
#plot(abs(initial_value));
for iStep=1:1:nSteps
# 2-stage Runge-Kutta time integral
   lap_coef = calc_lap(nGrid, dif_coef, cur_coef);
   for i = 1:1:nGrid
      hlf_coef(i) = cur_coef(i) + deltaT/2.0 * lap_coef(i);
   end
   lap_coef = calc_lap(nGrid, dif_coef, hlf_coef);
   for i = 1:1:nGrid
      nxt_coef(i) = cur_coef(i) + deltaT * lap_coef(i);
   end
   cur_coef= nxt_coef;
   if(mod(iStep,50) == 0)
     cur_val = ifft(cur_coef);
     hold on;
     plot(real(cur_val));
   endif
end
plot(real(cur_val));
hold off;
print -deps time_diffusion

アンドラ＆バルセロナ旅行(5/5日目)

帰りの飛行機が遅かったので、カタルーニャ美術館に行こうとしたら、なんとメーデーで休み！
ピカソ美術館も同じく。
事前情報が足りてなかった。

やることもないので、カタルーニャ美術館からとぼとぼと砦？まであるく。
途中、誰かの居住地跡みたいなところによるが、詳細はわからじ。

お昼は、地球の歩き方に乗ってる、バレンシアパエリャが食べられるお店で、イカスミパエリャを食べた。スープとワイン付き（色々選べるが私はそれを選んだ）で、20ユーロ以下だった。休日はちょっと高い値段付けであると、どうどうとメニューに書いてあった。

アンドラ＆バルセロナ旅行(4/5日目)

モンセラットという、バルセロナから電車を乗り継いで1時間半ぐらいのところに行って来ました。
黒いマリア像がある教会です。

私はスペイン広場の近くに泊まっていたので、スペイン広場のRENFE駅で、
往復＋モンセラットまでの登山電車＋そこから山の頂上までのゴンドラ＝約27ユーロ
のチケットを買いました。
行き方は二通りあるようですが、このチケットを買った場合はMonistrol de Montserratで乗り換えです。
一個手前の駅で降りて、ロープウェイでもいけるようです。

いわゆる教会なのですが、見どころは

(1)平日13時から始まる少年合唱団のコーラスは。ウィーンのものより荘厳で私は好きでした。
(2)黒いマリア像
です。

黒いマリア像を見るために30分ぐらい？並ばないといけませんので、頑張ってください。

この日もTonka barで夕食。Bioワインがおいしゅうございました。

アンドラ＆バルセロナ旅行(3/5日目)

特に書くことも無い気がするけど、市内観光。サグラダ・ファミリアみてきました。
室内は、森をイメージしてるらしいんですが、ナウシカの腐海の底のシーンを思い出します。

お昼ごはんは、La Rita (地球の歩き方に乗ってる)で、安くておいしゅうございました。
夕飯は、地球の歩き方に乗ってるアンダルシアタパス。ガスパチョ、おいしゅうございました。

アンドラ＆バルセロナ旅行(2/5日目)

アンドラからバルセロナへ。
あ、正確にはアンドラ　ラ　ベッリャですね。

行きたかったレストランは、小さいのに予約してなくて入れず。
代わりに、近くにあったバルで、ちょこちょこ食事を。多分そば粉のクレープに、チーズなどなどを包んだのがおいしゅうございました。

アンドラは綺麗で、ゆっくりするのにいいと思うのですが、観光するところが本当になんもないので、バルセロナからレンタカーで来るのが良いと思います。一週間ぐらい、何もしないでぐったりしながら、料理を楽しむところという感じ。

ガッツリ降りすぎだろう、というぐらい雪が降ってたので、そうそうにアンドラ観光を切り上げて帰りました。

夕食はホテルに近かった、tonka barに行きました。通常メニューじゃなくて、日替わりから選ぶべし。
http://tonkabar.com/

spainでprepaid sim (llamaya)

バルセロナの北バス停で買ったSIMカードはllyamayaという会社のものです。
たぶん、OrangeのMVNO.

10ユーロでSIMを買って、5ユーロクレジットが入ってました。
5ユーロあると、500MB/30daysのパッケージが買えるので、小旅行には十分でしょう。
他には、3ユーロ 250MB/7days, 9ユーロ 1GB/30daysのパッケージあり。
サイズは普通のを買ったので、micro SIMとかがあるのかどうかは知りません。

私はアクティベーションがうまく行ってなかったのか、試行錯誤中に少額クレジットを使ってしまい、250MB/7days を買いました。

同じトラブルに見舞われた人用に、多分私が陥った症状。

(1) 電波は拾う。事業者がllamaya movilであることも認識している。
(2) SMSを送れない。ので、パッケージが買えない。
(3) パケット通信もできない。
(4) *113#　にダイヤルするとクレジット額が帰ってくる。だから、完全に使えないわけでもない。

この状態で、あれやこれややっているうちに使えるようになったのだけれども、最終的に使えるようになったトリガーは、
902902909
に電話したこと。スペイン語が読めないのだけど、オペレーターにつなぐ番号みたい。
繋がってしまって、慌てて切ったけど、ここでクレジットを少し使ってしまった。

ここに電話をかけたあと、色々SMSが飛んできて、また、自分でも送れるようになった。
通話するのが儀式みたいなので、お店のお兄さんに、自分宛に電話をかけてもらえばうまくいきそうな気がしている。

APNの設定は、APNにmoreinternet　を入れるだけ。あとはAndoroidの初期値で良いみたい。

しばらくは、10ユーロで20ユーロ、20ユーロで40ユーロのクレジットが買えるキャンペーンらしいです。

ちなみに、アンドラでは電波入りませんでした。

URL:
http://www.llamayamovil.com/

アンドラ＆バルセロナ旅行 （1/5日目）

アンドラまでの行き方は、フランス　トゥールーズかスペイン　バルセロナからバスを使うしかなく、便数が多いバルセロナ経由にしました。

空港から空港バスでカタルーニャ広場まで移動。

地下鉄一号線で、カタルーニャからArc de triomph(スペル違うかも)まで移動。

で、北バス停のほぼ目の前にでます。

北バス停の中に、雑然とした携帯屋があるので、そこでSIMカードゲット。詳細は別ページ。

北バス停からはALSAのバスで、3時間ぐらい。

四時半に家をでて、13時ごろにはアンドラに着けます。

近くはないなぁ。

行ってから気がついたのでありますが、ピレネー山脈にあるわけで、アンドラではこの時期でも普通に雪がふります。