月別アーカイブ: 2014年1月

[F#] F#らしくLU分解を書き直してみる

LU分解を F# らしく書き直してみます。 せっかくの F# なのに for ループを使っているのがもったいないし(?)、数式のΣから外れています。ここを sum 関数を使って書き直したいのですが… matr … 続きを読む

カテゴリー: F# | [F#] F#らしくLU分解を書き直してみる はコメントを受け付けていません

[F#] LU分解を作る

旧タイプなので、連立一次方程式を解くときに逆行列を使って解けばOK…と頭から思い込んでいたのですが、Gaussの消去法を使うと逆行列は出てこなくて「?」となるわけです。上三角化を使うのだから、LU分解でいいん … 続きを読む

カテゴリー: F# | [F#] LU分解を作る はコメントを受け付けていません

[F#] 連立一次方程式を解く

逆行列、行列式とできたので、連立一次方程式を解かせる。有限要素法のいわゆる「solve」というやつで、解析の手順として、 pre : 構造を設定、要素に分割 solve : 連立一次方程式を解く(それだけじゃない?) p … 続きを読む

カテゴリー: F# | 1件のコメント

[F#] 行列式を計算する

F# MathProvider だと、あっさりと MathProvider.LinearAlgebra.inv なる関数があるのだが、自作してみる。 let det ( A' : matrix ) = let … 続きを読む

カテゴリー: F# | [F#] 行列式を計算する はコメントを受け付けていません

[F#] 逆行列を計算する

F#で数値・線形代数計算をするためのライブラリ紹介(F# PowerPack, F# MathProvider) http://d.hatena.ne.jp/teramonagi/20111215/1323874810# … 続きを読む

カテゴリー: F# | [F#] 逆行列を計算する はコメントを受け付けていません

年頭なので初心に帰って

去年と一昨年は何気に正月に仕事をしていたのですが、今年は年末に原稿の殴り書きをした後にプログラムを殴り書きして正月休みを作りました。妻の実家にノートパソコンを持って行ってプログラムの下書きなどをすることが多いのですが、今 … 続きを読む

カテゴリー: 雑談 | 年頭なので初心に帰って はコメントを受け付けていません