Cのプログラムに変換

質問

どなたかこのプログラムをCに変化してくれる方いませんか？
なんで変換するのとかの質問はいらないです。
ただ純粋に変換してくれる方いませんか？

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <complex.h>
#define NREPEAT 1000

int sweepout(float complex **a, int size) {
int n,m,k;
// 前進消去
for(n=0; n<size; n++) {
// 軸の選択
{
float amax = cabs(a[n][n]);;
float aabs;
int npivot = n;
for(m=n+1; m<size; m++) {
if((aabs = cabs(a[m][n]))>amax) {
npivot = m;
amax = aabs;
}
}
if(aabs == 0.0) return n;
if(npivot != n) { // 軸の交換
float complex t;
for(k=n; k<=size; k++) {
t = a[n][k];
a[n][k] = a[npivot][k];
a[npivot][k] = t;
}
}
}

// 消去
{
float complex t;
t = 1.0/a[n][n];
// 対角要素 a[n][n] <-- 1
for(k=n+1; k<=size; k++) a[n][k] *= t;
// 非対角要素 a[m][n] <-- 0
for(m=n+1; m<size; m++) {
for(k=n+1; k<=size; k++) a[m][k] -= a[m][n]*a[n][k];
}
}
}

// 後退代入
{
float t;
for(n=size-1; n>0; n--) {
for(m=0; m<n; m++) a[m][size] -= a[n][size]*a[m][n];
}
}
return -1;
}

// 乱数を発生する
double drand() {
double rn;
rn = (2*((double)random()))/((double)(RAND_MAX))-1.0;
return rn;
}

main(int argc, char **argv) {
int n, m, k, ncount;
double complex **a0;
float complex **a;
double complex t;
double d, r;

if(argc>1) n = atoi(argv[1]);
if(n <=0) n = 3;
srandom((argc>2)?atoi(argv[2]):0);

a0 = (double complex **)malloc(n*sizeof(double complex *));
for(m=0; m<n; m++)
a0[m] = (double complex *)malloc((n+1)*sizeof(double complex));
a = (float complex **)malloc(n*sizeof(float complex *));
for(m=0; m<n; m++)
a[m] = (float complex *)malloc((n+1)*sizeof(float complex));

for(ncount=0; ncount<NREPEAT; ncount++) {
// 乱数により、配列の係数を決める(double型)
for(m=0; m<n; m++)
for(k=0; k<=n; k++)
a0[m][k] = drand() + 1.0I*drand();
// 係数をfloat型に変換し、配列 a に代入する
for(m=0; m<n; m++) {
for(k=0; k<n; k++)
a[m][k] = (float complex)a0[m][k];
t = 0;
for(k=0; k<n; k++) t += a0[m][k]*a0[k][n];
a[m][n] = (float complex)t;
}
// 連立方程式を解き、数値解を a[][n] に求める
sweepout(a,n);
// 2乗誤差を求める
d = r = 0;
for(m=0; m<n; m++) {
d += pow(cabs(a[m][n]-a0[m][n]),2);
r += pow(cabs(a0[m][n]),2);
}
printf("%lf\t\n", 0.5*log10(d/r));
}

for(m=0; m<n; m++) {
free(a0[m]);
free(a[m]);
}
free(a0);
free(a);
}

回答No.8

Tacosan

もはやどうでもいいような気はしますが....

C++ だったら複素数はテンプレートになっているので, むしろ float complex とか書いたらエラーになりますです.

ところで, なんで float complex 使ってるんだろう. 全部 double complex にしちゃダメなのかなぁ?

回答No.7

Wr5

>たぶんタブを利用しているとインデントが消えるようです。スペースに変換してから
>アップすれば良いのかな？

半角スペースでは結局インデントが消えます。
コピペで処理できませんが、TAB1つを全角空白で…というパターンがありますかね。

本題と関係ないですが。

回答No.6

kmee

あえて言えば、 // から行末まで、というコメントはもともとはC++のものなので、エラーになる場合がある、ということでしょうか。

もっとも、新しい規格であるC99では正式にCにも導入されましたし、それ以前でも、コンパイラの独自機能として対応してたりもするので、実際コンパイルしてエラーになったら修正する、くらいでいいかと。

あとは、randomとsrandomは無い場合がある、とか mainの型を指定していないとか(intとみなされるので大丈夫だけど)とか、中身までは確認してないのでロジックが間違ってるかもしれない、とかはあるけど、いずれもC言語の範囲の問題で、このプログラムを「変換」という話しではありません。

回答No.5

jacta

すでに指摘されているように、このソースコードはCのものです。
処理系によっては可能かもしれませんが、原則としてC++ではコンパイルできません。
あと、randomやsrandomに関しては非標準関数ですので、特定の処理系に依存することになると思います。

回答No.4

annys2000

>読みにくくてすいません。

インデントついて無いと読む気が失せます。普段はスルーしてます
たぶんタブを利用しているとインデントが消えるようです。スペースに変換してから
アップすれば良いのかな？

>C++だと思ったんですけど違うんですか？

出てくる関数はCに存在します。C++依存のところは無いようです。
フリーで得られるコンパイラは殆どC/C++ 対応となってますのでどっちにしろコンパイル
できますが。

>ただ純粋に変換してくれる方いませんか？

そういうサイトではないというのは既出ですが、奇特な人はいるかもしれません
このくらい自分で調べられなければ、プログラミングできませんよ

回答No.3

Tacosan

#1 で言われている通り C ですね.

微妙に「ん?」と思うところはありますが.

回答No.2

noname#158634

＞変換してくれる方いませんか？
ここそういうサイトじゃないから(笑)

回答No.1

annys2000

読みにくいので細かく見てませんが、これってC言語に既になってると思うのですが・・
デバッグしてということ？

補足

読みにくくてすいません。

C++だと思ったんですけど違うんですか？