昔 Delphi で作成したものを C で書き直しました。Delphi ではまったく問題なかったのですが、C の場合、方程式の係数によって、掃き出し処理をした後 0 にマイナスが付加されることがあります。書式指定しない printf でも同じです。
なぜでしょうか?
※解は正しいです。
#include <stdio.h>
#include <math.h>
#define M 5
/* これは掃き出し処理、後単位行列になった 0 成分にマイナスがついてしまう
double A[M][M+1] = {
{ 5.0, -3.0, -3.0, 3.0, 2.0, -8.0 },
{ 9.0, -7.0, 1.0, -6.0, -5.0, -3.0 },
{ -1.0, -1.0, -8.0, 7.0, 4.0, -4.0 },
{ -7.0, 2.0, -6.0, 0.0, -1.0, -1.0 },
{ 8.0, 5.0, -3.0, 9.0, -7.0, 5.0 }
};
*/
double A[M][M+1] = { // こちらはOK
{ 2.0, 0.0, 4.0, 2.0, 0.0, 0.0 },
{ 2.0, 2.0, 0.0, 4.0, 8.0,12.0 },
{ 2.0, 4.0, 4.0, 0.0, 2.0, 0.0 },
{ 4.0, 2.0, 0.0, 0.0, 4.0, 0.0 },
{ 4.0, 0.0, 2.0, 2.0, 4.0, 8.0 }
};
//Gauss Jordan
int GaussJordan(int N)
{
double mMax, R_pivot, pivot, temp;
int pRow = -1;
for (int pv = 0; pv < N; pv++) //行ループ(一番外側のループ)
{
mMax = 0.000000001;
for (int k = pv; k < N; k++) {//行ループ ピボット行以下の最大値探索
if (abs(A[k][pv]) > mMax){
mMax = abs(A[k][pv]);
pRow = k;
}
}
if (mMax <= 0.000000001){ //0対策
printf("解が存在しないか、または不定です!");
return (0);
}
if (pv != pRow) //行の入れ替え
{
for (int k = 0; k < N + 1; k++) { //列ループ
temp = A[pv][k];
A[pv][k] = A[pRow][k];
A[pRow][k] = temp;
}
}
//ピボット行 A[pv][j] の処理
pivot = A[pv][pv];
R_pivot = 1.0 / pivot;
for (int j = 0; j < N + 1; j++) //列ループ
A[pv][j] = A[pv][j] * R_pivot; //A[pv][pv] = 1 となる
//ピボット行以外の処理 ⇒ ピボット列を 0 にする。
for (int k = 0; k < N; k++){ //行ループ
temp = A[k][pv]; //各行のピボット列成分
for (int j = pv; j < N + 1; j++) //ピボット列以降を処理
if (k != pv)
A[k][j] = A[k][j] - temp * A[pv][j];
}
}
return (1);
}
int main(void) {
printf("連立一次方程式\n");
for (int i = 0; i < M; i++) {
for (int j = 0; j < M+1; j++)
printf("%8.3g",A[i][j]);
printf("\n");
}
GaussJordan(M);
printf("\n解\n");
for (int i = 0; i < M; i++) {
for (int j = 0; j < M+1; j++)
printf("%8.3g",A[i][j]);
printf("\n");
}
}
No.1
- 回答日時:
・負で 0 に近い数を表示するときに出力指定の影響で「-0」になっている
・実際に値が -0 である (double の形式は処理系依存だが IEEE 60559 では +0 と -0 のどちらも存在する)
のどちらか, かなぁ.
%.20f みたいにぎりぎりまで表示するように出力を指定すれば, どっちかは区別がつく.
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
関連するカテゴリからQ&Aを探す
おすすめ情報
- ・漫画をレンタルでお得に読める!
- ・人生のプチ美学を教えてください!!
- ・10秒目をつむったら…
- ・あなたの習慣について教えてください!!
- ・牛、豚、鶏、どれか一つ食べられなくなるとしたら?
- ・【大喜利】【投稿~9/18】 おとぎ話『桃太郎』の知られざるエピソード
- ・街中で見かけて「グッときた人」の思い出
- ・「一気に最後まで読んだ」本、教えて下さい!
- ・幼稚園時代「何組」でしたか?
- ・激凹みから立ち直る方法
- ・1つだけ過去を変えられるとしたら?
- ・【あるあるbot連動企画】あるあるbotに投稿したけど採用されなかったあるある募集
- ・【あるあるbot連動企画】フォロワー20万人のアカウントであなたのあるあるを披露してみませんか?
- ・映画のエンドロール観る派?観ない派?
- ・海外旅行から帰ってきたら、まず何を食べる?
- ・誕生日にもらった意外なもの
- ・天使と悪魔選手権
- ・ちょっと先の未来クイズ第2問
- ・【大喜利】【投稿~9/7】 ロボットの住む世界で流行ってる罰ゲームとは?
- ・推しミネラルウォーターはありますか?
- ・都道府県穴埋めゲーム
- ・この人頭いいなと思ったエピソード
- ・準・究極の選択
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
std::set<int> で、ある値が何...
-
#define _CRT_SECURE_NO_WARNIN...
-
C言語での引数の省略方法
-
return 1L
-
複数桁10進数の*桁目だけを抽出...
-
C#のコンパイルエラーCS0120に...
-
「指定されたキャストは有効で...
-
C言語 エラーの原因がわからな...
-
C言語のサイコロシミュレート
-
ラップ関数とはどんなものですか?
-
困ってます…nCrを求めるC言語...
-
C言語の配列をC++のvectorに高...
-
int16_t の _t は何?
-
構造体の勉強中です 合計点の高...
-
Arduinoのプログラムにエラーが...
-
C言語の関数で戻り値を返す必要...
-
C言語で分からないところがあり...
-
文字列の構造体キャスト
-
素数問題
-
引数 戻り値 return文について
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
「指定されたキャストは有効で...
-
C言語での引数の省略方法
-
#define _CRT_SECURE_NO_WARNIN...
-
AtCoderABC135の問題Cについて
-
C言語 エラーの原因がわからな...
-
複数桁10進数の*桁目だけを抽出...
-
【C++】関数ポインタの使い方
-
実数の整数部,小数部の取得
-
ラップ関数とはどんなものですか?
-
if と配列の組み合わせ
-
return 1L
-
read関数をノンブロッキングで...
-
(int *)の意味
-
std::set<int> で、ある値が何...
-
Win32APIで作るコンボボックス...
-
C++でvectorにテキストファイル...
-
「{ } で囲むだけ」は正しい?
-
足して100になるような乱数のア...
-
Arduinoのプログラムにエラーが...
-
課題でつまってます・・・
おすすめ情報
> %.20f みたいにぎりぎりまで表示するように出力を指定すれば, どっちかは区別がつく.
やはり、マイナスがつきますね。
ま、これを表示させる必要はほとんどないのでとくに問題はないのですけど、何か引っかかる(笑)。
ちなみにC#ではこんな症状は生じませんでした。