やっぱりC++の挙動はわからない

powerpcとx86の問題は解決しましたが、また問題が出てきた。

init()を作って値を初期化しているのですが、いざ初期化した配列を使おうとするとなぜかされていません。変だなと思ってアドレスで無理やりしていしてみましたが、同じアドレスなのに初期化した後も値が入っていな状態になっています。現在原因特定中。

大雪

学校のまわりが雪降りすぎ。買い物にいけない。

更新情報

ホームページ更新すた。みてけろ

前処理付きの安定化双共役勾配法、PBiCGStab ( preconditioned Bi-Conjugate Gradient Stabilised )を追加。一応OpenMPで並列化してある。カリカリにはチューニングしてないけどそれなりに早いはず。ソースコードもあるよ。

http://www.geocities.jp/yosuke_3141592/numerics/bicgstab.html

c++の挙動がよくわからない件について

メッシュデータを読み込ませるときにかなり不可解な挙動をしめす。

メッシュのノードの座標位置を読みたい時、個数が事前にわからないので二回読み込ませている。これをなぜか同じ一回目の読み込み関数から、二回目のファイルアクセスをとエラーになる

もちろん一回目のファイルアクセスのあとifstreamを閉じている。この問題は、関数を二つに分け、mainから呼んだら正常に動きました

あとこの情況で、なぜか新しくポインタをつくるだけ、alloactionもアクセスもしていない、でallocation errorになる。

x86では正常に動くけど、power系では動かないといった謎なコードになってしまった。

ネットさまよっていたらすごい文字見つけた。

これ

☭

どうやってだすの？文字コードは？ソビエト限定？

七重対角行列をBiCGStabでといてみた

画像のような七重対角行列をBiCGStab法で解いてみました
iteration = 0 residual = 0.838021
iteration=59

x=
-0.329031
-0.390997
-0.385867
-0.459732
-0.452909
-0.413003
-0.465395
-0.436181
-0.359735
-0.325728
-0.271722
-0.188105
-0.137904
-0.083966
-0.010057
0.031578

A*x =
0.000000
1.000000
2.000000
3.000000
4.000000
5.000000
6.000000
7.000000
8.000000
9.000000
10.000000
11.000000
12.000000
13.000000
14.000000
15.000000

ちゃんと解けているようです。

いぬ

実家に先日帰ったとき、新しい犬がいました。前の子は自分が下宿先に移った矢先に亡くなってしまいました。とても可愛かったのですが。

新しい娘もなんだか不憫です。どこかのブリーダーに育てられたらしいのですが、商品として不適格と判断されたらしく声帯を手術されたらしいのです。おかげで吠えても力のない声しか出ません。ブリーダーの方に室内でずっと買われていたらしいのですが..ちょっとかわいそう。

姿は、たしかに目が小さく本来の犬種のイメージとは少し違います。でも人間の都合で飼われたり手術されたりするのはかわいそうだと思う。

本人はいつも日向でのんびりしているけど、幸せなのだろうか？