Vector Norm, Matrix Norm 을 구하는 수치해석 알고리즘 소스
- 최초 등록일
- 2005.01.01
- 최종 저작일
- 2005.01
- 4페이지/
한컴오피스
- 가격 1,000원
소개글
Vector Norm, Matrix Norm 을 구하는 수치해석 알고리즘 소스입니다. 과제수행시 참고하시기 바랍니다.
컴파일 실행환경
Visual C++ 6.0
1.5G 노트북
본문내용
/*
Norm
*/
#include<stdio.h>
#include<math.h>
float absol(float x);
float MatrixMulti(float *a, float *b, int c, int d, int e);
float result[3][1];
main(){
int i;
float x[3]={0,-7,5}, xh[3][1]={{-0.33},{-7.9},{5.8}},b[3]={1,-1,2};
float A[3][3]={{1,2,3},{2,3,4},{3,4,6}}, norminf=0,norminf2=0,t=0;
for(i=0;i<3;i++){
t=absol(x[i]-xh[i][0]);
if(t>norminf)norminf=t;
}
MatrixMulti(A,xh,3,3,1);
for(i=0;i<3;i++){
t=absol(result[i][0]-b[i]);
if(t>norminf2)norminf2=t;
}
printf("infinite norm of |x-xh| = %f infinite norm of |Axh-b| = %f\n",norminf,norminf2);
}
float absol(float x)
{
if(x>=0)return x;
else return (-1) * x;
}
float MatrixMulti(float *a, float *b, int c, int d, int e){
int i,j,k;
for(i=0;i<c;i++){
for(j=0;j<e;j++){
for(k=0;k<d;k++){
result[i][j] += *(a+i*d+k) * *(b+k*e+j);
참고 자료
없음
