非対称性は real で既に起きている様子 ...こっちには見込みなさそうなので、別のアプローチ

• 
  LU decomp ルーチンを使って inverse を計算 ...negative matrix になっている場合がある
• 
  bug 発見 ; M の self part の matrix が identity tensor の ２倍だった
• 
  修正後は問題無くなる！！-> cholesky にもどす
• 
  しかし、昨日の「直った事件」は、なんだったのか？ (多分 N=1 としてたのでしょう)