행렬곱셈을 수행하는 함수를 작성하려고 하는데 타입에 대한 이해가 되지 않습니다.
글쓴이: cjscokdd / 작성시간: 금, 2015/03/27 - 12:00오후
#define a11 a->p[0]
#define a12 a->p[1]
#define a21 a->p[2]
#define a22 a->p[3]
#define b11 b->p[0]
#define b12 b->p[1]
#define b21 b->p[2]
#define b22 b->p[3]
#define c11 c->p[0]
#define c12 c->p[1]
#define c21 c->p[2]
#define c22 c->p[3]
#define d11 d->p[0]
#define d12 d->p[1]
#define d21 d->p[2]
#define d22 d->p[3]
typedef union _matrix {
double **d;
union _matrix **p;
} *matrix;Strassen행렬곱을 표현하는 알고리즘을 작성하려고 하는데요, 이미 작성된 다른 소스를 검토하는 중에
각 행렬의 부분행렬 a11,a12, a21, a22등을 a->p[0], a->p[1], a->p[2], a->p[3]와 같이 나타내며,
행렬 타입을 선언하는데 왜 union을 사용하는지, 타입이름이 *matrix와 같이 나타내었는지 이해가 잘 되지 않아 질문드립니다.
Forums:
'왜'냐는 질문에 대답할 수 있는건 만든 사람뿐일
'왜'냐는 질문에 대답할 수 있는건 만든 사람뿐일 듯합니다.
그걸 구현하는 방법은 하나가 아니니까요.
그리고 여기서 typedef한 타입이름은 *matrix가 아니라 matrix입니다.
union _matrix* 를 matrix로 정의한 것입니다.
http://math.ewha.ac.kr/~jylee
http://math.ewha.ac.kr/~jylee/SciComp/sc-crandall/strassen.h
http://stackoverflow.com/questions/9539082/how-to-use-this-c-code-to-multiply-two-matrices-using-strassens-algorithm
일단,
타입 이름이 *matrix가 아니라,
타입 이름은 matrix 인데, 이 타입의 실제 정의는 union _matrix * 인 거고요.
union _matrix 가 유니온인 이유는 (주석문을 보아하니) 그 안에서 d 필드를 써서 접근할 때는 double 의 2차원 배열(더 정확히는, 동적으로 할당된 1차원 배열들을 가리키는 포인터들의 1차원 배열)로 접근할 수도 있고, p 필드를 써서 접근할 때는 또다른 union _matrix 들의 2차원 배열(이 때도 더 정확히는 동적으로 할당된...)로서 접근할 수도 있게 한 거네요. 어떤 때는 실제 원소의 배열이고 어떤 때는 서브행렬의 배열로 쓸 것인지는 나름의 조건(size와 BREAK의 비교)이 있고.
그러면 위 define의 경우,
a12 라고 쓰면, 실제로는 a 변수(이 변수가 가리키는 대상은 적절히 초기화되어 있을 것이고)에 대해 p[1]을 접근하니까, 결과적으로
a 행렬 안에 있는 네 개의 서브행렬 중 두번째 행렬의 포인터를 얻게 되죠.
알고리즘 자체에 대해서야 질문하신 분이 아실 테니까 이 정도면...
좋은 하루 되세요!
댓글 달기