lineinnercell.cpp

qrdrls算法

#include "lineInnerCell.h"

lineInnerCell()
{
	Win=0;
	for(i=0;i<M;i++)
	{
		Zin[i][M-1]=0;
	}
	
}

int lineInnerCell::caculation()
{
	int i,k;
	for(i=M-1;i>=0;i++)
		for(k=M-1;k>=i,k++)
		{
			Zout[k][k-1]=Zin[k][k]+r[i][k].conjugate()*Win[k].conjugate();
			Zin[k][k-1]=Zout[k][k-1];
		}
}

complex2 lineInnerCell::get_r()
{
	return r;
}
complex2 lineInnerCell::get_p()
{
	return p;
}
complex2 lineInnerCell::get_Wout()
{
	return Wout;
}

complex2 lineInnerCell::get_Zout()
{
	return Zout;
}