なかのにっき

素晴らしい。

独立数式で数式番号をつけるのはどうすんじゃ、 という問いへの解答。しかしこれもアッタマ悪い仕様だなあ。 2010 ではどうなんでしょう。

MOS 構造のゲートに 40nm 厚の Si リッチ SiO2 を付け、 17% excess、N2 中 1250℃ で 1h 焼いたところ、 交番電界による電荷注入でルミネッセンスを示したというお話。 Si 基板の上に付けてるのだが、バンド図に SiO2 層があるのはなぜだろう？

Schneider さんのところの。こないだ JAP に載ってた話と基本的には同じ。 Al ターゲットの場合に、 pulse DC のターゲット電位の時間で均した分布と EDF を比較。

Si/SiO2 をリモートプラズマ CVD 交互に積んだ構造の Raman 散乱による研究。 これは基板は quartz、1100℃ 30s の RTA。 Si 層を薄くすると nc-Si 起源の Raman ピークが計算とは逆にブルーシフトする (というか bulk 側に戻る)。 んでこれを圧縮応力由来と仮定すると 5GPa くらいだけど、ということ。

SiO2 中に分散させた CdSe クラスタの PL。 温度・電場依存性。電場が入ると光らなくなる。 10K, 77K では電場の向きに対して対称だけど、300K では変なかたち。 原因として、電界による電荷注入と、電場によって光励起した exciton が壊れる効果を議論。

MOCVD みたいな rod がスパッタでもできるで、という話。 XPS、SEM/TEM、field emission 特性。 3N 4インチの Ru ターゲット、working pressure (全圧か？) は 0.01Torr、 rf 電力 20W、んで O2/Ar を変えて成長。5sccm/10sccm まで行くと rod が見えはじめ、 さらに基板を200→450℃まで上げると faceted な rod 成長になるとのこと。

250〜1000 hPa の Ar/H2O mixture に 15keV 10uA の電子ビームを突込んで 310 nm あたりの励起光を出させる話。効率 3% ってすごく高いな。 しかも実験の範囲だと低圧の下限に max があったらしい。 Ar の purification として SAES getters の Monotorr と言うのを使っている。 日本代理店 もある。今度聞いてみよう。

glass などに stress をかけていったときの塑性変形をどう捉えるか、という話。 確かに結構難しいのかも。二種類の原子を含む (通常の L-J potential) のモデルが amorphous になりやすい傾向を持つ、というのはおもしろいかも。 local な準結晶が重要で、この部分は安定、amorph なところから変形が始まる、 ということらしい。

柱からなる 2D の音響レンズを GA で設計する話。 結果の絵がかなり面白い。

ときどき出てくるけど、hot electron の scattering による lifetime の計算。 考えてるのは背景 electron gas との scattering。このあたり phonon なんかの cross section は全然小さかったりするんだろうか。しらんけど。 一様な electron gas の結果に比べると、 d バンド伝導体では重要になる、という結果。

スパッタで ripple 構造作っといて、431〜466K でその緩和をさせたかたち。 測定は ex situ AFM で。あと再現性確認用に光散乱。 理論的には、Free energy を傾斜の関数として表わせれば、 その (傾斜に対する？) 変化率から chemical potential がわかり、 あとは mobility がわかれば flux が出て、形状の時間変化が導ける、 ということのようだ。それらをどう定義するかは結構難しそうだけど。 離散モデルも出ているから、いじってみると面白いかも。

*CLS
*RST
:SYST:REM

:VOLT:DC:RANG 10
:VOLT:DC:RES 0.0001
:VOLT:DC:NPLC 1
:TRIG:SOUR IMM
:TRIG:DEL MIN
:SAMP:COUN 1

無限ループ{
  (時刻取得)
  READ?
  read();
  printf(時刻・データ)
  fflush(stdout)
  usleep(9e5)
}

という感じで。

測定の無限ループを Ctrl-C で止めたときに 34401A とシリアル線の終了処理をするため、

void intHandler(int signum){
  xtc_send(":SYST:LOC");
  xtc_close();
  exit(0);
}

int main(int argc, char * argv[])
{
    :
  struct sigaction sa;
    :
  memset(&sa, 0, sizeof(sa));
  sa.sa_handler = intHandler;
  if (sigaction(SIGINT, &sa, NULL)) perror("sigaction");
    :

という感じで。 『プログラミング Linux』 の 13 章あたり。

    struct timeval tv;
    struct timezone tv;
       :
    gettimeofday(&tv, &tz);
    printf("%ld.%06ld %s\n", tv.tv_sec, tv.tv_usec, rbuf);

という感じで。

ここが一番はまった。 どうもタイミングによって "READ?" を送った後の read() コールで刺さることがある。 いろいろ試行錯誤していたのだが、 どうも制御文字関連ではまっているっぽい。 34401A は シーケンスの最後に \r\n を返してくる のだけど、 今回 m102 から引っぱってきたルーチンを叩き台にしてたので、 termios 構造体メンバの c_cc[VEOL] が '\r' になっており、 このせいで末尾の \n を 34401A 側のバッファに残しちゃってたようだ。あほだ。

DTR/DSR フローコントロールをしているので、EOL もろもろの制御文字は不要のような気もするのだが、 c_lflag ICANON フラグを落としちゃうと、 しょっぱなの read() 呼び出しの際に文字列が切れる。 ということで ICANON は立てておき、 c_cc[VEOL] を '\n' にしたらうまくいった…っぽい。うーむ。 '\r' は c_iflag に IGNCR を立てて捨ててしまうことにした。 んで昨日のやり方で \n を \0 にして返すことにした。

…それでも、どうしてもときどき刺さることがあるなあ。 つづけて何回か 起動-Ctrl-C を繰り返し、 一度回りだせば続けて回るのだが。

ん、*RST 送ると 34401A で「カッチン」というリレーっぽい音がしてたので、 そのあとに usleep(5e5) 噛ますと大丈夫になったっぽい。 うーむ、このへん 34401A の素性かな。どうも良くわからんな。

XTC のテストモード、こいつは 16.2Angstrom/sec で擬似的な積算膜厚に対応した 出力を出すことになっているのだけど、 実際に取ってみた電圧の増加率は 0.1574 V/s。 これはループ中の usleep() を変えても (=測定間隔を変えても) きっちり同じなので、どうやら XTC のタイマが 3% くらいずれている…のかな？ まあ表裏の ratio が重要なので、 これは覚えといて、とりあえずこのまま測ってしまうことにしよう。

miniterm.c は短いのだけど結構おもしろい。 fork() して親が serial input→stdout、子が stdin→serial output を担うようになっているのだが、 fork() の返り値で switch〜case するやり方とか、 特定の文字が入ると子が親にシグナルを送って終了する様子とか、 親がシグナルハンドラを登録して処理に入り、子を wait してから終了する様子とか。 おそらくこの手のプログラムの定跡なのだろうが、勉強になる。

 const char QCMD[] = "Q:\r";
  :
 write(fd, QCMD, sizeof(QCMD));

のようにしたところ、結果はちゃんと返ってくるのだけど、 終了後にコントローラが 例の HOST/MANUAL が効かない状態 になった。sizeof(QCMD) は (あたりまえだけど) 4 なので、 末端の \0 も write しちゃってるためらしい。

 write(fd, QCMD, sizeof(QCMD)-1);

にしたら問題なくいくようになったけど、頭悪そうだなあ。 どうするのが普通なんでしょうかね。 \0 のない 3 文字の定数って簡単に宣言できないのかな。

むつみさん に

  const char QCMD[3] = { 'Q', ':', '\r'};

と教えていただいた。うまいこと動作しました。ありがとうございます。

うーん、数値入りのコマンド文字列を strcat() でつくるときとかを考えると、 結局 sizeof()-1 しておく方が見通しがいいのかなあ。うーんうーん。

あ、

 const char QCMD[] = "Q:\r";
  :
 write(fd, QCMD, strlen(QCMD));

がいいかな。こうしとくか。

-      900,         0,K,K,R

のような結果が返ってくるのを、 掟破りの strtok() で ad-hoc に実装。はっきり言ってヒドイ代物だが以下のような感じ。

  res = read(fd,buf,255);
  buf[res] = 0;

  token = strtok(buf, ",");
  if (token != NULL){
    *xpos = atol(token+1);
    if (*token == '-') {*xpos *= -1;}
  } else { rcode = 1;}

  token = strtok(NULL, ","); /* ypos */
  if (token != NULL){
    *ypos = atol(token+1);
    if (*token == '-') {*ypos *= -1;}
  } else { rcode = 1;}

  token = strtok(NULL, ","); /* ack1 */
  switch (*token){
    case 'X':
      *ack1 = ACK1_NG;
      break;
    case 'K':
      *ack1 = ACK1_OK;
      break;
    default:
      *ack1 = ACK_ERROR;
      rcode = 1;
      fprintf(stderr, "(m102_q): parse error for ACK1\n");
  }
   :

以下 ack2, ack3 も同様。

Q: !: C: M: G: R: の 6 コマンドを実装した。 cvs モジュール m102 として import。 checkout した場所は ~/project/stage。

あらゆる 前回 のコピペだが、ひとまずまとめた。

日本語ファイル名の扱いで悩む。 やっぱ このへん を使わないとだめかな。( 2001/02 の日記 より)。

要するに locate で取得したパス名の EUC 文字の部分を URL エンコードしたいのだが、

$arg =~ s/(\W)/sprintf("%%%02X", ord($1))/ego;

だと / も置換されちゃうのでちとまずげ。結局アドホックに

$arg =~ s/([^-\/\.a-zA-Z_0-9])/sprintf("%%%02X", ord($1))/ego;

としたが、いかにも「ダセエ」。 小松さんの urlencode.el では例外を指定できるようになってるのだが、こんな感じが良いなあ。

ところで、逆に URL に生で指定するとまずい ASCII 文字って何があるのだろうか? とりあえずスペースはまずそうだな。

馬場さんの 『NAMAZU システムの構築と活用』 を参考にしつつインデックス作り。

まず /etc/namazu/mknamazurc を

$ADDRESS = 'webmaster@surf.ap.seikei.ac.jp';

だけ書き変えて

sudo -u www-data mknmz -O /var/lib/namazu/index /var/www/EXPDB

を実行し、インデックスができたことを確認。 コマンドラインの namazu でのテストもパス。

続いて /etc/namazu/namazurc の

Replace       /var/www/  /
Lang          ja

の両エントリを変更し、cgi の動作を確認。 ツリー構築スクリプトと mknmz とを cron に仕込む。

将来は /home 以下を嘗める別のインデックスも作りたいけど、 まあとりあえずこんなもんで。メンバーにメールする。

いい天気だったのでチャリ乗りついでに。

今日から開けたらしく、10:30 ですでに長蛇の列。

二郎以外の飲食店は学校の周りではほとんど休みだったので、 吉祥寺に出てココイチ。 この自転車だと割に気楽に吉祥寺まで行こうという気になるな。 サンロードあたりは今日でもむちゃむちゃ混んでいたので押して歩く。

学生さん達が研究室に次々と登場。

正月の時差ボケで眠くなってしまったので早々に撤退。 非常に快調に飛ばす。 これで筋肉痛さえ出なければ、 毎日チャリで通えるのだがのお。

<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01//EN"
  "http://www.w3.org/TR/html4/strict.dtd">

<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN"
  "http://www.w3.org/TR/REC-html40/loose.dtd">

ついでにいじっていたら、なぜか 前につくったやつ は 31x31 とか 15x15 で保存されていたらしい。 あきらさんが deb を作ってくれたので、 そいつでもう一回作り直したら mozilla でもちゃんと見えるようになったみたい。 やり方は変えてない。なんでじゃろ。

linux-users での optimization and accuracy with gcc のスレッドで出てきた fpu_control.h を使うのが、 potato (glibc 2.1.3) 〜 woody (ibid 2.2) の時点では portable のようだ。

#include <fpu_control.h>

void setfpu_()
{
	fpu_control_t cw;

	cw = _FPU_DEFAULT & ~(_FPU_MASK_IM | _FPU_MASK_ZM | _FPU_MASK_OM);
	_FPU_SETCW(cw);
}

こんな感じかな。 こいつを Fortran ソースから

      call setfpu

と呼べば良い。

本当はこっちが本命っぽい。 多分 fesetenv() 辺りを使うんだろうけど、 現時点ではまだ「例外を trap する」環境がまだ定義されていないようだ。 #define __USE_GNU すれば、woody の <fenv.h> では

#ifdef __USE_GNU

/* Enable individual exceptions.  Will not enable more exceptions than
   EXCEPTS specifies.  Returns the previous enabled exceptions if all
   exceptions are successfully set, otherwise returns -1.  */
extern int feenableexcept (int __excepts) __THROW;

/* Disable individual exceptions.  Will not disable more exceptions than
   EXCEPTS specifies.  Returns the previous enabled exceptions if all
   exceptions are successfully disabled, otherwise returns -1.  */
extern int fedisableexcept (int __excepts) __THROW;

/* Return enabled exceptions.  */
extern int fegetexcept (void) __THROW;
#endif

だったり、<bits/fenv.h> では

#ifdef __USE_GNU
/* Floating-point environment where none of the exception is masked.  */
# define FE_NOMASK_ENV  ((__const fenv_t *) -2)
#endif

となってるようなのだけど *1 。そのうち C の標準に入るのかしら。

NOTE: This will be the FINAL announcement sent to the old
      list (ldp-announce@lists.linuxdoc.org). All announcements
      will now go ONLY to the new list server (announce@linuxdoc.org).

      If you have not yet subscribed, I suggest you take a look at
      http://www.linuxdoc.org/mailinfo.html  for details.

ということらしい。