수다닷컴

질문 제목 : 파일 입출력으로 짠 코드를 그냥 입력 받는 방식으로 어떻게 바꿔야되요?질문 요약 :파일 입출력으로 짠 코드를 그냥 입력 받는 방식으로 어떻게 바꿔야되요? 질문 내용 :
#include stdio.h
#include stdlib.h
#include math.h
typedef struct matrix_t
{
int m, n;
double **a;
} matrix_t;

void create(matrix_t *a, int m, int n);
void destroy(matrix_t *a);
void print(matrix_t *a); // 연산 결과 출력함수
void gauss(matrix_t *a); // 가우스 요르단
void make_ai(matrix_t *a, matrix_t *ai); // [a|i] 만들기
void extract_i(matrix_t *ai, matrix_t *i); // [a|i] 에서 i 부분 추출
void product(matrix_t *a, matrix_t *b, matrix_t *p); // p = a * b
void create(matrix_t *a, int m, int n)
{
int i;

a-m = m;
a-n = n;
a-a = (double **)malloc( sizeof(double*) * m);
for(i=0; i m; i++)
{
a-a[i] = (double *)malloc( sizeof(double) * n);
}
}
void destroy(matrix_t *a)
{
int i;

for(i=0; i a-m; i++)
{
free(a-a[i]);
}
free(a-a);
a-m = a-n = 0;
a-a = 0;
}
void load(matrix_t *a, file *fptr)
{
int m, n;
int i, j;
fscanf(fptr,%d %d,&m,&n);
create(a,m,n);
for(i = 0; i a-m; i++)
{
for(j = 0; j a-n; j++)
{
fscanf(fptr,%lf, &a-a[i][j]);
}
}
}
void print(matrix_t *a)
{
int i, j;
for(i = 0; i a-m; i++)
{
for(j = 0; j a-n; j++)
{
printf( % 6.2f\t, a-a[i][j]);
}
puts();
}
puts();
}
void gauss(matrix_t *a)
{
int i, j, k;
double pivot;
printf(\n gauss-jordan elimination 풀이 과정 : \n);
printf( ***************************************************** \n\n);
for(i=0; i a-m; i++)
{
pivot = a-a[i][i];
if(pivot == 0)
{
return;
}
for(j=i+1; j a-n; j++)
{
a-a[i][j] /= pivot;
}
a-a[i][i] = 1;
print(a);
for(k=0; k a-m; k++)
{
if(i == k) continue;
&n pivot = a-a[k][i];
for(j=i+1; j a-n; j++)
{
a-a[k][j] -= a-a[i][j] * pivot;
}
a-a[k][i] = 0;
}
print(a);
}
return;
}
// a 에서 [a|i] 를 만듦.
void make_ai(matrix_t *a, matrix_t *ai)
{
int i, j;
create(ai, a-m, a-n+a-m); // 행 수 만큼 더 긴 행렬
for(i=0; i a-m; i++)
{
// a 복사
for(j=0; j a-n; j++)
{
ai-a[i][j] = a-a[i][j];
}
// i 부분 단위행렬 만들기
for(j=0; j a-m; j++)
{
ai-a[i][a-n+j] = 0.;
}
ai-a[i][a-n+i] = 1.;
}
}
// [a|i]에서 i 부분 추출
void extract_i(matrix_t *ai, matrix_t *i)
{
int i, j;
int m = ai-m;
int r = ai-n - ai-m;
create(i, m, m);
for(i=0; i m; i++)
{
for(j=0; j m; j++)
{
i-a[i][j] = ai-a[i][r+j];
}
}
}
// p = a * b
void product(matrix_t *a, matrix_t *b, matrix_t *p)
{
int i, j, t;
int m = a-m;
int n = b-n;
int k = b-m;
create(p,m,n);
for(i=0; i m; i++)
{
for(j=0; j n; j++)
{
p-a[i][j] = 0;
for(t=0; t k; t++)
{
p-a[i][j] += a-a[i][t] * b-a[t][j];
}
}
}
}
int main()
{
int i, j;
matrix_t a, ai, i, x;
file *fptr;
fptr = fopen(matrix.txt,rt);
load(&a,fptr);
print(&a);
make_ai(&a,&ai);
gauss(&ai);
extract_i(&ai,&i); // 1/a
print(&i);
product(&a,&i,&x); // a * 1/a
print(&x);
destroy(&a);
destroy(&ai);
destroy(&i);
destroy(&x);
fclose(fptr);
}

신청하기