2004年04月09日(Fri) いい天気なのに [長年日記]
● 家の中で数値積分
いくつかパラメータのある数値積分をして、観測と合うようにフィッティングをしています。微分できないので最小二乗法じゃなくて極小値探し。
しかしなあ…。プロファイルを取ってみると一番時間を使ってるのはpow( x, 2+n ) / (exp( x ) - 1)
こんなコード。2階の数値積分の内側の被積分関数なのです。powとexpのxもnも極小値探しで実数を動きまわるのです。どうしようもないよねぇ…。
いくつかパラメータのある数値積分をして、観測と合うようにフィッティングをしています。微分できないので最小二乗法じゃなくて極小値探し。
しかしなあ…。プロファイルを取ってみると一番時間を使ってるのはpow( x, 2+n ) / (exp( x ) - 1)
こんなコード。2階の数値積分の内側の被積分関数なのです。powとexpのxもnも極小値探しで実数を動きまわるのです。どうしようもないよねぇ…。
最近のツッコまれどころ