行列方程式を解く（１）

• 問題設定
  • Ax=bという線形方程式を解きたい（Aは行列、xとbはベクトル）
  • Aは対称行列である
  • Aはとりあえず48x48の大きさ程度のものを扱おう
  • Aは対角成分、上三角成分、下三角成分に分けて記録しておく
  • Aの対角成分はファイルDDに記録する。ファイル形式は、"行番号　値”という形で縦に並べたものとする
  • Aの上三角成分はファイルAUに記録する。ファイル形式は、"行番号　列番号　値”という形で縦に並べたものとする。ただし0でない値だけが記録させているとする。
  • Aの下三角成分はファイルALに記録する。ファイル形式は、"行番号　列番号　値”という形で縦に並べたものとする。ただし0でない値だけが記録させているとする。
  • bの成分はファイルBBに記録する。ファイル形式は、"行番号　値”という形で縦に並べたものとする
  • Aは疎な(Sparse)行列とする

このような問題を扱いたい。解がない問題は扱いようがないので例題を選ぶのに注意せねば。
（とりあえずここまで）