Tk

From Ecal

Jump to: navigation, search
  • personal memo by tk.


Contents

bug in ubuntu csh

set exinfo = $?

did show

Newline in variable name.

Use tcsh instead. or "$stasus" in bash. --Tkotani2 21:02, 3 September 2009 (JST)


memo


fortranコードの解析(字句、構文の解析をしてデータベース化;スタンダードはないのか?) それのデータベース的な処理による可視化(粗視化あるいは必要な情報をみせる方法)、 コードへのリンクを張る、graphvizあるいはUMLのviewer? f2pyはどうなるのか?あるいはもっといい方法は?



将来の数値計算プログラミング:おそらく、すごく簡単に書けるようになるだろう。

  • 「indexセット」という型、合成則、順引きと逆引き、 for x in x_index:の様な書式。およそパイソンではだいたいできてる?
  • 「データセット」という型(indexセットの上で定義される)。
  • ユーザーインターフェイスとしてのコーディングと「計算の形式化」としてのコーディングの切り分け。いろいろな可視化。
  • 変数の使いまわしはしない(破壊的な代入はしない)。ただし「積算変数宣言」したものは加算(=+)できる。またループ内での局所変数。
  • 関数型言語のような原則的に計算手順を指定しないコーディング(積算変数などはのぞく)、ただし「計算手順」を適切に並べてやらないとエラーを吐く。
  • 変数の型宣言はしない(型推論でやる)。変数は生成されるときに値と型が与えられる。

bug in gfortran

bug for gfortran ------hqpe.sc.m.F L545:

ctakao May2009
c    I don't know, but gfortran did not allow next line.
c    So I made a replacement. to 
c        evec_invt= dconjg(transpose(evec_inv))  !this did not work for gfortran
        evec_invt= transpose(dconjg(evec_inv))  !this works for gfortran

--Tkotani 17:30, 19 May 2009 (UTC)

Try to make better interpolation for sigm.

 Test. Write LDA overlap matrix in bndfp.F
  Look for ctakaox. 
  hambls.F pass ss (overlap.gas) from bndfp.F

  After "call bloch", we will multiple s(1,1,ispec) by ss.
  Thus we need to devide w(ohrs) (in rdsigm before it is transformed to real space).--Tkotani 16:11, 2 June 2009 (UTC)
Personal tools