- 締切済み
初心者な質問その５
************************************************************************* void mmul(double a[],double b[],double c[],int imax,int jmax,int kmax){ int i,j,k; for(i=0;i<imax;i++){ for(k=0;k<kmax;k++){ c[(1)]=0.0; for(j=0;j<jmax;j++){ c[i*kmax+k] += a[(2)]*b[(3)] } } } } ************************************************************************* 行列の乗算を計算するプログラムです。 ⅰ，(1)～(3)に入る式 ⅱ，この関数mmulを使って下の行列を計算するプログラム についてお願いします。
- rokkaku
- お礼率27% (3/11)
- C・C++・C#
- 回答数1
- ありがとう数2
- みんなの回答 （1）
- 専門家の回答
みんなの回答
- hitomura
- ベストアンサー率48% (325/664)
i (1):i*kmax+k (2):i*jmax+j (3):j*kmax+k ii …って、お～い、下の行列って何だ～？
関連するQ&A
- ファイルの読み込みが上手くいかず困っています。
#include<stdio.h> #include<math.h> #include<time.h> #define i0 2 /*流体部分の最小格子番号*/ #define j0 2 #define k0 2 #define in 41 /*流体部分の最大格子番号*/ #define jn 41 #define kn 41 #define imax 44 /*境界を含めた格子の数*/ #define jmax 44 #define kmax 44 #define N 100000 void sabun(void); int i,j,k,n,pgosa,vgosa,x,y,z; double u[imax][jmax][kmax], v[imax][jmax][kmax], w[imax][jmax][kmax], p[imax][jmax][kmax]; double u2[imax][jmax][kmax], v2[imax][jmax][kmax], w2[imax][jmax][kmax], p2[imax][jmax][kmax]; double Vn=-1.0, a=10.0; /*初期値*/ double dt=0.001, dx=0.025, dy=0.025, dz=0.025, rei=250000.0; double dx2, dy2, dz2, po; main() { FILE *fpin,*fpout; clock_t start,end; start = clock(); /*時間計測開始*/ fpin=fopen("H-clock.dat", "r"); /*入力ファイル*/ if(fpin==NULL) { printf("入力ファイルのオープンに失敗しました"); return 0; } x=i-1; y=j-1; z=k-1; for(i=i0;i<=in+1;i++){ for(j=j0;j<jn+1;j++){ for(k=k0;k<kn+1;k++){ fscanf(fpin,"%d %d %d %.6lf %.6lf %.6lf %.6lf\n",&x,&y,&z,&u[i][j][k],&v[i][j][k],&w[i][j][k],&p[i][j][k]); printf("%d %d %d %.6lf %.6lf %.6lf %.6lf\n",x,y,z,u[i][j][k],v[i][j][k],w[i][j][k],p[i][j][k]); } } } 読み込みの部分のある部分まで載せているのですが読み込みの部分に問題が あると思うのですがどうすれば読み込んでくれるのか分かりません。ですので解決法を教えていただけると助かります、よろしくお願いします。
- 締切済み
- C・C++・C#
- 行列に関して。
以下の行列を対角成分が最も大きくなるような プログラムをつくりたいのですが 0 1 4 3 7 0 1 -5 1 3 0 7 -2 4 4 -8 #include<stdio.h> #include<math.h> #include<stdlib.h> #define NUM 4 int main(){ /* 初期化 */ int i; /* ループ変数 （行）*/ int j; /* ループ変数 （列）*/ int k; /* ループ変数 */ int imax; /* 最大値のｉ成分 */ int jmax; /* 最大値のｊ成分 */ double nmax; /* 最大値を入れるための変数 */ double tmp; /* ソートのための一時的な変数 */ double eps = 1.0e-6; /* 判定値 */ double a[NUM][NUM] = { { 0.0, 1.0, 4.0, 3.0}, { 7.0, 0.0, 1.0, -5.0}, { 1.0, 3.0, 0.0, 7.0}, {-2.0, 4.0, 4.0, -8.0}}; /* 対象とする要素位置 = k */ for( k = 0 ; k < NUM ; k++ ){ /*対象とする要素位置 = k*/ nmax = 0.0;/* 初期化 */ imax = k;/* 初期化 0とは限らない */ jmax = k;/* 初期化 0とは限らない */ /* 1．絶対値の最大値を求める */ for( i = k ; i < NUM ; i++ ){ for( j = k ; j < NUM ; j++ ){ if(fabs(nmax) < fabs(a[i][j])){ nmax = a[i][j]; imax = i; jmax = j; /* 対象としている行列要素の中で，最大となる行列要素のｉ，ｊ成分を探す */ /*ここに絶対値の最大値を求め、nmax、imax、jmaxに代入するプログラムを書く*/ /*絶対値の求め方、数学関数fabs(x)を用いる*/ } } } /* 2．行入れ換え */ if( k != imax ){/* 対象としている行と最大値を持つ行が同じなら，入れ替える必要が無い */ for( j = 0 ; j < NUM ; j++ ){ tmp=a[k][j]; a[k][j]=a[imax][j]; a[imax][j]=tmp;/* 対象としている行と最大値を持つ行との入れ替え */ /*ここに行の入れ替え行うプログラムを書く*/ } } /* 3．列入れ換え */ if( k != jmax ){/* 対象としている列と最大値を持つ列が同じなら，入れ替える必要が無い */ for( i = 0 ; i < NUM ; i++ ){ tmp=a[i][k]; a[i][k]=a[i][imax]; a[i][imax]=tmp;/*対象としている列と最大値を持つ列との入れ替え*/ /*ここに列の入れ替え行うプログラムを書く*/ } } } /* 結果の出力 */ for(i=0; i<NUM; i++){ for(j=0; j<NUM; j++){ printf("%5.1f ",a[i][j]); } printf("\n"); } /*ここに行列の表示プログラムを書く*/ return(0); } 上のぷろぐらむを実行しても、1個目のfor文の１ループ目で終了してしまい 1行目と、4行目を交換、4列目と1列目を交換した行列だけが表示されます。 結果が以下のような行列にするにはどう改善すればよいでしょうか？ -8.0 -2.0 4.0 4.0 -5.0 7.0 1.0 0.0 3.0 0.0 4.0 1.0 7.0 1.0 0.0 3.0 です。
- 締切済み
- C・C++・C#
- 行列の計算について。
こんにちは。 JAVAで行列の計算をするプログラムを作っていますが、なかなかうまくいきません。 自分で考えたのはfor文を使った物です。行列数を入力するとその数だけの行列を作る、という物なのですが、どうもうまくいきません。 「2」を入力すると、 │1 2││5 6│ │ ││ │ │3 4││7 8│ こんな2行2列の行列の計算をする、というものです。 for (int i =0; i < 2; i++){ for (int j = 0; j<2; j++){ for (int k = 0; k<2; k++){ C[i][j]=C[i][j]+a[i][k]*b[k][j]; } } } 上のようなプログラムを考えました。 JAVAには行列の計算をする専用の関数（クラス？）などあるのでしょうか？ また、ソースなど教えていただけると、とても助かります。 CからJAVAに急に移ったので、つまづくことが非常に多いです。 よろしくお願いします。
- ベストアンサー
- Java
- スレッドについて質問があります。
こんばんは。 今、「スレッドを使った並列計算」というテーマを考えています。 並列処理が可能である行列の計算を、スレッドを用いて処理してみる、といったことをやっているんですが、スレッドの使い方がいまいちわかりません。 プログラムなんですが、はじめに、n行n列の配列を定義します（nは自分で入力）。 double A[][] = new double [n][n]; double B[][] = new double [n][n]; double C[][] = new double [n][n]; その後、A行列、B行列として使う配列に数値を代入します（A行列には1～n^2、B行列には(n^2)+1～(n^2)*2までの数値が代入されるようになっています）。 そして、次の計算を行なうわけですが・・・ for (i = 0; i < n; i++){ for (j = 0; j < n; j++){ C[i][j]=0; for (k = 0; k < n; k++){ C[i][j] += A[i][k] * B[k][j]; } } } この計算をどうスレッドに対応させればいいのか、わからなくてずっと悩んでいます。 ぜひみなさんの意見、アドバイスを聞かせてください。 よろしくお願いします。
- ベストアンサー
- Java
- C＋＋についての質問です
プログラミング初心者です 以下の通りに正方行列の積を求めるプログラムを作成したのですが、うまくいきません。 #include<stdio.h> #define DTM 20 void InputMatrix(double[][DTM], int, char); void PrintMatrix(double[][DTM], int, char); void MatrixMulti(double[][DTM], double[][DTM], double[][DTM], int); int main(void) { double matrixA[DTM][DTM]; double matrixB[DTM][DTM]; double matrixC[DTM][DTM]; int n; printf("正方行列の積を求めるプログラムです\n"); printf("正方行列の次元を入れてください(<=20):"); scanf_s("%d", &n); InputMatrix(matrixA, n, 'A'); InputMatrix(matrixB, n, 'B'); MatrixMulti(matrixA, matrixB, matrixC, n); printf("\n行列 C =A×B\n"); PrintMatrix(matrixC, n, 'C'); return 0; } void InputMatrix(double a[][DTM], int n, char ch) { int i, j; printf("行列 %cの入力\n", ch); for (i = 0; i < n;i++) { for (j = 0;j < n;j++) { printf("%c[%d][%d] =", ch, i + 1, j + 1); scanf_s("%lf", &a[i][j]); } } } void PrintMatrix(double a[][DTM], int n, char ch) { int i, j; printf("行列 %c の出力\n", ch); for (i = 0;i < n;i++) { for (j = 0;j < n;j++) { printf("%5.2f\t", a[i][j]); } printf("\n"); } } void MatrixMulti(double a[][DTM], double b[][DTM], double c[][DTM], int n) { int i, j, k; for (i = 0;i < n;i++) { for (j = 0;j < n;j++) { c[i][j] = 0; for (k = 0;k < n;k++) { c[i][j] =a[i][k] * b[k][j]; } printf("%5.2f\t",c[i][j]); } printf("\n"); } }
- ベストアンサー
- C・C++・C#
- ガウスの消去法
逆行列を求めたいのですが どうも上手くできません。 下にソースコードを載せますので どなたかどこが悪いのか ご教授してくれませんか？ よろしくお願いします。 void MINV(double A[][100],double InvAf2[][100],int G,int R){ int i,j,k; double temp;//行列値の一時記憶 double B[100][100],C[100][100];//行列の一時記憶 for(k=0;k<R;k++){//k列目を調整する for(i=0;i<G;i++){//行数を指す for(j=0;j<R;j++){//列数を指す temp=A[k][k]; B[i][j]=A[i][j];//InvAf2を計算するのに必要 C[k][j]=InvAf2[k][j];//InvAf2を計算するのに必要(更新した値を使われると困る) if(i==k){ A[i][j]/=temp; InvAf2[i][j]/=temp; } else{ A[i][j] -= A[k][j] * A[i][k] / temp; InvAf2[i][j] -= C[k][j] * B[i][k] / temp; } } } } }
- ベストアンサー
- C・C++・C#
- c言語 行列の積に関して
＜問＞ 4行３列の行列aと３行４列の行列bの積を、４行４列の行列cに格納する関数を作成せよ。 void mat_mul(const int a[4][3], const int b[3][4], int c[4][4]) 入門レベルのスキルしかありません。 上手く行列の積のプログラムが組めません。 行列の積の計算結果が何も出てきません。 どの様にしたら良いかご指導の程、宜しくお願いします。 ＜プログラム＞ void mat_mul(const int a[4][3], const int b[3][4], int c[4][4]) { int i, j, k; for (i = 0; i < 4; k++) { for (j = 0; j < 4; i++) for (k = 0; k < 3; j++) c[i][j] = c[i][j] + (a[i][k] * b[k][j]); } } void mat_print(const int m[4][4]) { for (int i = 0; i < 4; i++) { for (int j = 0; j < 4; j++) printf("%4d", m[i][j]); putchar('\n'); } } int main(void) { int i, j ,k; int tensu1[4][3]; int tensu2[3][4]; int seki[4][4]; for(i = 0; i < 4; i++) { for (j = 0; j < 3; j++) { scanf("%d", &tensu1[i][j]); } putchar('\n'); } for(i = 0; i < 3; i++) { for(j = 0; j < 4; j++){ scanf("%d", &tensu2[i][j]); } putchar('\n'); } putchar('\n'); mat_mul(tensu1, tensu2, seki); puts("行列の積"); mat_print(seki); return 0; }
- ベストアンサー
- C・C++・C#
- LU分解について
こんにちは。 行列のＬＵ分解のｃによるプログラムを勉強している者です。 ＬＵ分解のｃによるプログラムについてお尋ねしたいことがございます。 ニューメリカルレシピに、ＬＵ分解の方法として 以下のようなコードが記載せれておりました。 void ludcmp(double **a, int n, int *indx, double *d) { int i,imax,j,k; double big,dum,sum,temp; double *vv; // 各行の暗黙のスケーリングを記録する． vv=vector(n); *d=1.0; // まだ行交換していない． for (i=0;i<n;i++) { // 行についてループし，暗黙のスケーリングの情報を得る． big=0.0; for (j=0;j<n;j++) if ((temp=fabs(a[i][j])) > big) big=temp; if (big == 0.0) printf("Singular matrix in routine ludcmp\n"); // 最大要素が０なら特異行列である． vv[i]=1.0/big; // スケーリングを記録する． } for (j=0;j<n;j++) { // Crout法，列についてのループ for (i=0;i<j;i++) { // 方程式(2.3.12)のi=j以外 sum=a[i][j]; for (k=0;k<i;k++) sum -= a[i][k]*a[k][j]; a[i][j]=sum; } big=0.0; for (i=j;i<n;i++) { sum=a[i][j]; for (k=0;k<j;k++) sum -= a[i][k]*a[k][j]; a[i][j]=sum; if ( (dum=vv[i]*fabs(sum)) >= big) { big=dum; imax=i; } } if (j != imax) { for (k=0;k<n;k++) { dum=a[imax][k]; a[imax][k]=a[j][k]; a[j][k]=dum; } *d = -(*d); vv[imax]=vv[j]; } indx[j]=imax; if (a[j][j] == 0.0) a[j][j]=TINY; if (j != n) { dum=1.0/(a[j][j]); for (i=j+1;i<n;i++) a[i][j] *= dum; } } free_vector(vv); } 例えば、ＬＵ分解させたい行列３×３の（｛4,7,6｝｛2,5,9｝｛3,1,8｝）だとしたら、 メイン関数にどのように書いたらＬＵ分解を実行できるでしょうか？ ニューメリカルレシピを買って読んでみたものの、よく理解できず、質問させていただきます。 稚拙な質問かもしれませんが、どうぞよろしくお願いいたします。
- ベストアンサー
- C・C++・C#
- 行列の積を関数を使って求める・・？
２つの行列の行と列を入力し、積を計算するプログラムを関数を使って書きたいのですが、上手く行きません。どこをどのように直したらよいか教えてください！お願いします！！ 以下が私が書いたプログラムです。 #include<stdio.h> #define NUMBER 10 int first(int x1,int x2,int y1,int y2,int i,int j,int k) { int a[NUMBER][NUMBER] = {0}; int b[NUMBER][NUMBER] = {0}; int c[NUMBER][NUMBER] = {0}; do{ printf("２つの行列の行と列を入力してください\n"); scanf("%d", &x1); scanf("%d", &x2); scanf("%d", &y1); scanf("%d", &y2); if(x1 != y2){ printf("行列の積は計算できません\n"); } }while(x1 != y2); printf("行列Ａの要素を入力してください\n"); for(i=0; i<x1; i++){ for(j=0; j<x2; j++) scanf("%d", &a[i][j]); } printf("行列Ｂの要素を入力してください\n"); for(j=0; j<y1; j++){ for(k=0; k<y2; k++) scanf("%d", &b[j][k]); } } int second(int x1,int x2,int y1,int y2,int i,int j,int k) { int a[NUMBER][NUMBER] = {0}; int b[NUMBER][NUMBER] = {0}; int c[NUMBER][NUMBER] = {0}; for(i=0; i<x1; i++){ for(k=0; k<y2; k++){ for(j=0; j<x2; j++) c[i][k] = c[i][k] + a[i][j]*b[j][k]; } } for(i=0; i<x2; i++){ for(k=0; k<y2; k++) printf("%3d", c[i][k]); printf("\n"); } } int main(void) { int a[NUMBER][NUMBER] = {0}; int b[NUMBER][NUMBER] = {0}; int c[NUMBER][NUMBER] = {0}; printf("行列の積を計算します\n %d\n", first(x1,x2,y1,y2,i,j,k)); printf("行列Ａと行Ｂの積は\n %3d",second(x1,x2,y1,y2,i,j,k)); }
- 締切済み
- C・C++・C#
- C言語を用いた45×４５の逆行列の表示について・・・・。。
次のプログラムを動かすと逆行列を求めることができます。 １０×１０の計算までは、スムーズに計算する事ができたのですが、 私が知りたい４５×４５の行列の計算になるとなぜかうまくいきません。 そこでみなさんにお聞きしたいのは、計算がうまくいかない理由がプログラムにあるのか、単にＰＣの限界なのか、それ以外なのか。みなさんからアドバイスをいただけないでしょうか？ ※下のプログラム、見づらくてすいません。 #include <stdio.h> #define N 10 /* 元の数 */ #define N2 2*N /* プロトタイプ宣言 */ void Print(int k); /* 大域変数 */ double a[N][N2]= {/* 行列 Aおよび単位行列 */ { ここに45×45の正則行列＆45×45の単位行列を入れます・・・。。 }; int main(void) { double d; int i,j,k; for (k=0;k<=N-1;k++) { Print(k); for(i=0;i<=N-1;i++) { if(i!=k) { d=a[i][k]/a[k][k]; for (j=k+1;j<=N2-1;j++) a[i][j]=a[i][j]-a[k][j]*d; a[i][k] = 0.0; } } d=a[k][k]; for (j=k;j<=N2-1;j++) a[k][j]=a[k][j]/d; } Print(N); } void Print(int k) { int i,j; printf("\n ステップ %d \n",k); /* 小見出し */ for(i=0;i<=N-1;i++) { /* 第ｉ行を印刷 */ for(j=0;j<=N-1;j++) printf("%6.4f,",a[i][j]); printf(" : "); for(j=N;j<=N2-1;j++) printf("%6.4f,",a[i][j]); printf("\n"); } }
- 締切済み
- C・C++・C#
補足
ああ、すみません！ 「下の行列」は （1 2 3) （3 4） Ａ＝（5 6 7) Ｂ＝ （2 6） Ｃ＝Ａ・Ｂ （4 5） うう、行列の表現がわからない（；；） 括弧はつながってると思って下さい。 Ａは２行３列、Ｂは３行２列の行列です。