cmatrixmultiply2222.cpp

来自「dfb laaser c++ novel」· C++ 代码 · 共 28 行

/*  本程序用以计算两个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];
}