数値解析に関する質問です

質問

以下のプログラムでは、関数pivotは，３つの引数を利用しているが，これを2つに減らす方法はないか？
教えてください。よろしくお願いします。

/*Gaussian Elimination*/
#include <stdio.h>
#include <math.h>

/*--------------*/
void inputdata(double [3][4]);
void calsol(double [3][4],double [3]);
void elimination(double[3][4],double[3][4]);
void pivot(int *, double [3][4],int *);
void outputdata(double [3]);
void checksol(double [3][4], double [3]);

/*--------------*/
main(){
double c[3][4], x[3];
/*-----------------------*/

printf("Program of Gauss Elimination\n");
/*-----------------------*/
inputdata(c);
/*-----------------------*/
calsol(c,x);
/*-----------------------*/
outputdata(x);
/*-----------------------*/
checksol(c,x);
/*-----------------------*/
return(0);
}
/*-----------------------*/
void inputdata(double c[3][4]){
int ie,ic;
printf("------\n");
printf("Input Data\n");
printf("(a0*x0 + a1*x1 +a2*x2 = a3)\n");
printf("\n");

for(ie=0; ie<=2; ie++){
printf("Equation %d\n",ie);
printf("Input a0 a1 a2 a3\n");
for(ic=0;ic<=3;ic++){
scanf("%lf", &c[ie][ic]);
}
}
}
/*-----------------------*/

void calsol(double c[3][4], double x[3]){

　 double cc[3][4];

elimination(c,cc);
x[2]=cc[2][3];
x[1]=cc[1][3]-cc[1][2]*x[2];
x[0]=cc[0][3]-cc[0][2]*x[2]-cc[0][1]*x[1];
}
/*-----------------------*/
void elimination(double c3[3][4], double cc[3][4]){
int ie,imax,eli;
double c2[3][4],c1[3][4];

eli = 0;
pivot (&eli,c3,&imax);

cc[0][0]=1;
cc[0][1]=c3[imax][1]/c3[imax][0];
cc[0][2]=c3[imax][2]/c3[imax][0];
cc[0][3]=c3[imax][3]/c3[imax][0];

for(ie=0;ie<=2;ie++){
if(ie != imax){
c2[ie][0]=0;
c2[ie][1]=c3[ie][1]-cc[0][1]*c3[ie][0];
c2[ie][2]=c3[ie][2]-cc[0][2]*c3[ie][0];
c2[ie][3]=c3[ie][3]-cc[0][3]*c3[ie][0];

}
}

eli = 1;
pivot(&eli,c2,&imax);

cc[1][1]=1;
cc[1][2]=c2[imax][2]/c2[imax][1];
cc[1][3]=c2[imax][3]/c2[imax][1];

　 for(ie=0;ie<=1;ie++){
if(ie !=imax){
c1[0][2]=c2[ie][2]-cc[1][2]*c2[ie][1];
c1[0][3]=c2[ie][3]-cc[1][3]*c2[ie][1];

}
}

cc[2][2]=1;
cc[2][3]=c1[0][3]/c1[0][2];
}
/*-----------------------*/

void pivot(int *eli,double c[3][4],int *imax){
int ie;
double cmax;

if(*eli ==0){
cmax=0;
for(ie=0;ie<=2;ie++){
if(fabs(c[ie][0])>cmax){
cmax=c[ie][0];
*imax=ie;
}
}
}
else{
cmax=0;
for(ie=0;ie<=1;ie++){
if(fabs(c[ie][1])>cmax){
cmax=c[ie][1];
*imax=ie;
}
}
}
}
/*-----------------------*/

void outputdata(double x[3]){
int i;
printf("------\n");
printf("Solutions\n");

for(i=0;i<=2;i++){
printf("x%d=%14.6e\n",i,x[i]);
}
}
/*-----------------------*/

void checksol(double check[3][4],double x[3]){
double crt,clt;
int ie,ic;
printf("------\n");
printf("Check Solutions\n");
printf("Left Term Right Term\n");

for(ie =0;ie<=2;ie++){
clt =0 ;
for(ic=0;ic<=2;ic++){
clt = clt + check[ie][ic]*x[ic];
}
crt = check[ie][3];
printf("%14.6e%14.6e\n",clt,crt);
}
}

回答No.3

yama5140

＞関数pivotは，３つの引数を利用しているが，これを2つに減らす方法はないか？

　その前に、「関数は、void 型」という拘りがあるのでしょうか。
　無い場合、　関数 pivot を、int 型とし、
　呼び出し元で、

　　　imax = pivot( 0, c3 );
　　　　　・
　　　imax = pivot( 1, c2 ); とするのが普通かな。

　関数 pivot は（共通部分を簡略化して）、

　　　int　pivot( int eli, double c[][4] )　　　　　eli は、int * でない
　　　{
　　　　　int ie, imax;
　　　　　double cmax = 0.0;

　　　　　for( ie = 0; ie <= ( 2 - eli ); ie++ ){
　　　　　　　if( fabs( c[ ie ][ eli ] ) > cmax ){
　　　　　　　　　cmax = c[ ie ][ eli ];
　　　　　　　　　imax = ie;
　　　　　　　}
　　　　　}
　　　　　return( imax );
　　　}

回答No.2

Tacosan

なぜ減らしたいんでしょうか?

回答No.1

prophetok

構造体で渡せば２つでも１つでも可能。