구조체와 배열 함수로 인자전달
쇼콜라데
질문 제목 : 구조체와 배열 함수로 인자전달질문 요약: 제가 c++을 공부하지않은 상태라서 그런지잘 모르겟네요ㅠ
구조체와 배열을 함수로 인자전달할 떄 인자를 어떻게 선언해야할까요ㅠ
아예 방향을 잘못잡고 있는것같기도 하구요...
소스가 긴대요;; 빨간색 칠한부분을 어떻게 처리해야할지 모르겠어요ㅠ
녹색부분참고하셔서 가르쳐주시면 감사하겠습니다ㅠ
질문 내용 :#include stdio.h
#include math.h
const float epsilon = 1.0e-10f;
const int maxsweeps = 32;
void calculateeigensystem(const matrix3d& m, float *lambda, matrix3d& r);
struct arr getdata(void);
struct matrix3d
{
float n[3][3];
float& operator()(int i,int j)
{
return(n[j][i]);
}
const float& operator()(int i,int j)const
{
return(n[j][i]);
}
void setidentity(void)
{
n[0][0]=n[1][1]=n[2][2]=1.0f;
n[0][1]=n[0][2]=n[1][0]=n[1][2]=
n[2][0]=n[2][1]=0.0f;
}
};
struct arr{
float a[4][3];//벡터 성분변수 선언
};
int main()
{
struct arr s=getdata();float p[3];//평균위치 m 연산
p[0]=(s.a[0][0]+s.a[1][0]+s.a[2][0]+s.a[3][0])/4;
p[1]=(s.a[0][1]+s.a[1][1]+s.a[2][1]+s.a[3][1])/4;
p[2]=(s.a[0][2]+s.a[1][2]+s.a[2][2]+s.a[3][2])/4;
printf(평균 위치 m=(%.2lf %.2lf %.2lf)\n,p[0],p[1],p[2]);int i,j;
float c[3][3];
c[0][0]=(pow(s.a[0][0]-p[0],2)+pow(s.a[1][0]-p[0],2)+pow(s.a[2][0]-p[0],2)+pow(s.a[3][0]-p[0],2))/4;//공분산 행렬 연산
c[1][1]=(pow(s.a[0][1]-p[1],2)+pow(s.a[1][1]-p[1],2)+pow(s.a[2][1]-p[1],2)+pow(s.a[3][1]-p[1],2))/4;
c[2][2]=(pow(s.a[0][2]-p[2],2)+pow(s.a[1][2]-p[2],2)+pow(s.a[2][2]-p[2],2)+pow(s.a[3][2]-p[2],2))/4;
c[0][1]=c[1][0]=((s.a[0][0]-p[0])*(s.a[0][1]-p[1])+(s.a[1][0]-p[0])*(s.a[1][1]-p[1])
+(s.a[2][0]-p[0])*(s.a[2][1]-p[1])+(s.a[3][0]-p[0])*(s.a[3][1]-p[1]))/4;
c[0][2]=c[2][0]=((s.a[0][0]-p[0])*(s.a[0][2]-p[2])+(s.a[1][0]-p[0])*(s.a[1][2]-p[2])
+(s.a[2][0]-p[0])*(s.a[2][2]-p[2])+(s.a[3][0]-p[0])*(s.a[3][2]-p[2]))/4;
c[1][2]=c[2][1]=((s.a[0][1]-p[1])*(s.a[0][2]-p[2])+(s.a[1][1]-p[1])*(s.a[1][2]-p[2])
+(s.a[2][1]-p[1])*(s.a[2][2]-p[2])+(s.a[3][1]-p[1])*(s.a[3][2]-p[2]))/4;
printf(공분산 행렬 c는 다음과 같습니다\n);
for(i=0;i3;i++)//공분산 행렬 표시
{
for(j=0;j3;j++)
{
printf(%8.2lf,c[i][j]);
}printf(\n);
}
float lamb[3];
float te[3][3];//struct matrix3d power;
calculateeigensystem(c, lamb, te);------------------이 부분을 모르겠습니다ㅠ
return 0;
}
void calculateeigensystem(const matrix3d& m, float *lambda , matrix3d& r)
{
float m11=m(0,0);
float m12=m(0,1);
float m13=m(0,2);
float m22=m(1,1);
float m23=m(1,2);
float m33=m(2,233=m(2,2);
r.setidentity();
for(int a=0;amaxsweeps;a++)
{
if((fabs(m12)epsilon)&&(fabs(m13)epsilon)&&
(fabs(m23)epsilon)) break;
if(m12 !=0.0f)
{
float u = (m22-m11)*0.5f/m12;
float u2=u*u;
float u2p1=u2+1.0f;
float t= (u2p1 != u2) ?
((u0.0f)?-1.0f:1.0f)*(sqrt(u2p1)-fabs(u))
:0.5f/u;
float c=1.0f/sqrt(t*t+1.0f);
float s=c*t;
m11 -= t*m12;
m22 += t*m12;
m12=0.0f;
float temp = c*m13-s*m23;
m23=s*m13+c*m23;
m13=temp;
for(int i=0; i3; i++)
{
float temp = c*r(i,0)-s*r(i,1);
r(i,1) = s*r(i,0)+c*r(i,1);
r(i,0)=temp;
}
}
if(m13 !=0.0f)
{
float u = (m33-m11)*0.5f/m13;
float u2=u*u;
float u2p1=u2+1.0f;
float t= (u2p1 != u2) ?
((u0.0f)?-1.0f:1.0f)*(sqrt(u2p1)-fabs(u))
:0.5f/u;
float c=1.0f/sqrt(t*t+1.0f);
float s=c*t;
m11 -= t*m13;
m33 += t*m13;
m13=0.0f;
float temp = c*m12-s*m23;
m23=s*m12+c*m23;
m12=temp;
for(int i=0; i3; i++)
{
float temp = c*r(i,0)-s*r(i,2);
r(i,2) = s*r(i,0)+c*r(i,2);
r(i,0)=temp;
}
}
if(m23 !=0.0f)
{
float u = (m33-m22)*0.5f/m23;
float u2=u*u;
float u2p1=u2+1.0f;
float t= (u2p1 != u2) ?
((u0.0f)?-1.0f:1.0f)*(sqrt(u2p1)-fabs(u))
:0.5f/u;
float c=1.0f/sqrt(t*t+1.0f);
float s=c*t;
m22 -= t*m23;
m33 += t*m23;
m23=0.0f;
float temp = c*m12-s*m13;
m13=s*m12+c*m13;
m12=temp;
for(int i=0; i3; i++)
{
float temp = c*r(i,1)-s*r(i,2);
r(i,2) = s*r(i,1)+c*r(i,2);
r(i,1)=temp;
}
}
}
lambda[0]=m11;
lambda[1]=m22;
lambda[2]=m33;
int i;
for(i=0;i3;i++)
{
printf(%f,lambda[i]);
}
}
struct arr getdata(void)
{
struct arr t;
int i,j;
for(i=0;i4;i++)
{
printf(%d번째 벡터 \n,i+1);
for(j=0;j3;j++)
{
printf( %d번째 성분입력:,j+1);
scanf_s(%f,&t.a[i][j]);
}
}
return t;
}
