#pragma once
#include "afxtempl.h"
#include "FeVector.h"
class CFeMarix
{
public:
	CFeMarix(void);
	CFeMarix(long m_lNIn, CArray<long,long>& m_lMarixBandIn);
	CFeMarix(long m_lNIn, CArray<long,long>& m_lMarixBandIn,double** data);
	CFeMarix(CFeMarix&	a);
	void Init(long m_lNIn,CArray<long,long>& m_lMarixBandIn);
	~CFeMarix(void);
public:
	CFeMarix& operator= (CFeMarix&);
	CFeMarix& operator+= (CFeMarix&);
	CFeMarix& operator-= (CFeMarix&);
	CFeMarix& operator*= (double);
	CFeVector operator* (CFeVector&);
	long Judge(long i);
public:
	long m_lN;
	CArray<long,long> m_lMarixBand;
	double** data;
	double** ic_data;
public:
	void Set(CFeMarix& a);
	void Set_E(long i, long j, double a);
	double Get_E(long i, long j);
	void Add_E(long i, long j, double a);
public:
	void Save(void);
public:
	void CreateData(long m_lNin, CArray<long,long>& m_lMarixBandIn, double** datades);
	void ICC();
	CFeVector ICCINV(CFeVector& R);
public:
	void CholeskyEquation(CFeVector& X, CFeVector& B);
	void Zero(void);
};