/*  本程序用以計算兩個2×2的復矩陣的乘積。c_f3[2][2]=c_f2[2][2]×c_f1[2][2]
使用方法：
主程序中聲明：  extern void CMatrMuti(complex<long double> [2][2],complex<long double> [2][2],complex<long double> [2][2]);

調用方式： CMatrMuti(c_f3,c_f2,c_f1); 三個矩陣可以相同可以為相同矩陣；
*/

#include <iostream>
#include <complex>
using namespace std;

void CMatrMuti2222(complex<long double> [2][2],complex<long double> [2][2],complex<long double> [2][2]);

void CMatrMuti2222(complex<long double>c_f3[2][2],complex<long double>c_f2[2][2],complex<long double> c_f1[2][2])
{
	int i,j;
	complex<long double> c_t3[2][2],c_t2[2][2],c_t1[2][2];
	for(i=0;i<2;i++)
		for(j=0;j<2;j++)
		{
			c_t2[i][j]=c_f2[i][j];
			c_t1[i][j]=c_f1[i][j];
		}
	c_f3[0][0]=c_t2[0][0]*c_t1[0][0]+c_t2[0][1]*c_t1[1][0];
	c_f3[0][1]=c_t2[0][0]*c_t1[0][1]+c_t2[0][1]*c_t1[1][1];
	c_f3[1][0]=c_t2[1][0]*c_t1[0][0]+c_t2[1][1]*c_t1[1][0];
	c_f3[1][1]=c_t2[1][0]*c_t1[0][1]+c_t2[1][1]*c_t1[1][1];
}