c++ как возвести матрицу в степень n ?

[proef]

так не работает

int I,J,K;
         float e  = 0;
		
for ( I = 0;I<ISizeOfMatrix;++I){    e = 0;
for ( K = 0;K<ISizeOfMatrix;++K){
for ( J= 0;J<ISizeOfMatrix;++J)
{ e += this->InitMatrix[J][K]*InitMatrix[K][J];}
				
				this->InitMatrix[I][J] = e;}}

как нужно сделать?

[proef]

вот если что программа полностью

// ConsoleApplication68.cpp: определяет точку входа для консольного приложения.
//


#include "stdafx.h"
#include <iostream>
#include <locale.h>
#include <math.h>

class Matrix
{protected:
  float **InitMatrix;
   int ISizeOfMatrix;
   int JSizeOfMatrix;
public:
	/*ввод*/   friend std::istream &operator>>(std::istream &Stream,Matrix &a)   
	{
		for (int NumberI = 0;NumberI<a.ISizeOfMatrix;NumberI++){
			for (int NumberJ = 0;NumberJ<a.JSizeOfMatrix;NumberJ++){
				std::cout<<"["<<NumberI<<"]"<<"["<<NumberJ<<"] = ";
              Stream>>a.InitMatrix[NumberI][NumberJ];
			  }
			

		}
			std::cout<<"______________________________"<<std::endl;
		return Stream;
	
	}
	 /*вывод*/ friend std::ostream &operator<<(std::ostream &Stream,Matrix &a)
	{
		for (int NumberI = 0;NumberI<a.ISizeOfMatrix;NumberI++)
		{
			for (int NumberJ = 0;NumberJ<a.JSizeOfMatrix;NumberJ++)
			{
				std::cout<<"["<<NumberI<<"]"<<"["<<NumberJ<<"] = ";
			    Stream<<a.InitMatrix[NumberI][NumberJ];
			}
			Stream<<std::endl;
		
		}
		Stream<<"______________________________"<<std::endl;
		return Stream;
	}
	 /*конструктор*/                             Matrix(int SizeI ,int  SizeJ)
	 {       
		 ISizeOfMatrix = SizeI;
	    JSizeOfMatrix = SizeJ;
           InitMatrix = new float*[ISizeOfMatrix];
		for (int Number = 0;Number<ISizeOfMatrix;Number++)
		{
			InitMatrix[Number] = new float[JSizeOfMatrix];
		}
		  }
	 /*деструктор*/                                                  ~Matrix()
								  { 
									  delete[] this->InitMatrix;
								  }
	 /*+*/                                  Matrix& operator+(Matrix &Element)
	{ 
		std::cout<<"сложение"<<std::endl;
		
		 if((this->ISizeOfMatrix != Element.ISizeOfMatrix) ||(this->JSizeOfMatrix != Element.JSizeOfMatrix)){
             std::cout << "Different sizes" << std::endl;
             return *this;}else {
       
         for (int NumberI = 0;NumberI<this->ISizeOfMatrix; NumberI++){
             for (int NumberJ = 0;NumberJ<this->JSizeOfMatrix;NumberJ++){
                this->InitMatrix[NumberI][NumberJ] += Element.InitMatrix[NumberI][NumberJ];
             }
		 }
		 return *this;
			 }
     
     }
/*возведение в степень*/                       Matrix & operator()(int Degree)
	{
	if (this->ISizeOfMatrix !=this->JSizeOfMatrix){ std::cout<<"матрица не квадратная"<<std::endl;return *this;}
	else
	{    int I,J,K;
         float e  = 0;
		
		for ( I = 0;I<ISizeOfMatrix;++I){e = 0;for ( K = 0;K<ISizeOfMatrix;++K){
		for ( J= 0;J<ISizeOfMatrix;++J){ 
			e += this->InitMatrix[I][K]*InitMatrix[K][J];}
		this->InitMatrix[I][J] = e;}}

									 return *this;
	}
		}
	/*приравнивание*/                       Matrix& operator=(Matrix &Element)
	{
		for (int NumberI = 0;NumberI<this->ISizeOfMatrix;NumberI++){
		
			for (int NumberJ = 0;NumberJ<this->JSizeOfMatrix;NumberJ++)
			{
			this->InitMatrix[NumberI][NumberJ] = Element.InitMatrix[NumberI][NumberJ];
			}}
			
		return *this;}
	/*обращение по индексу*/            float &operator()(int SizeI,int SizeJ)
	{
	 return this->InitMatrix[SizeI][SizeJ];} 
};

int main(int argCount, char** arguments)
{
	setlocale(2,"rus");
	Matrix a(3,3);
	std::cout<<"матрица A"<<std::endl;
	std::cin>>a;
	std::cout<<a;
	Matrix b(3,3);
	std::cout<<"Матрица В"<<std::endl;
	std::cin>>b;
    std::cout<<"a(1,0)+b(1,0) ="; 
    float q = a(1,0)+b(1,0);  std::cout<<q<<std::endl;
    std::cout<<a(2);
    std::cout<<a+b;
 
   system("pause");
	return 0;
}

[proef]

в интернете то толком и нету такого почти (я не нашёл )
может кто знает?

[proef]

теперь только первую строку правильно делает

 int I,J,K;
         float e  = 0.0 ;
		//for (int CurrDegr = 0;CurrDegr<Degree;CurrDegr++){
		for ( I = 0;I<ISizeOfMatrix;++I)
			for ( J = 0;J<ISizeOfMatrix;++J){
			e = 0.0;
		for ( K= 0;K<ISizeOfMatrix;++K)
			 e += InitMatrix[J][K]*InitMatrix[K][J];
		       
		
		
	
		InitMatrix[I][J] = e;
		
		
		
		
		}

[BDA]

Сохраняйте результаты расчета во временную матрицу.

tmp = new float*[ISizeOfMatrix];
for (int Number = 0;Number<ISizeOfMatrix;Number++)
    tmp[Number] = new float[JSizeOfMatrix];
for (I = 0;I<ISizeOfMatrix;++I)
    for (J = 0;J<JSizeOfMatrix;++J) {
        tmp[I][J] = 0;
        for ( K= 0;K<ISizeOfMatrix;++K)
            tmp[I][J] += InitMatrix[I][K]*InitMatrix[K][J];
    }
for (I = 0;I<ISizeOfMatrix;++I)
    delete[] InitMatrix[I];
delete[] InitMatrix;
InitMatrix = tmp;

Кстати, у Вас неправильный деструктор. Количество delete должно совпадать с количество new.

[proef]

Спасибо
сначала тоже делал так но почему - то выдавало
"
Необработанное исключение по адресу 0x00C3889C в ConsoleApplication68.exe: 0xC0000005: нарушение прав доступа при чтении по адресу 0x00008127."
решил переделать так

[proef]

cейчас в любом случае в квадрат возводит

	float** tmp = new float*[ISizeOfMatrix];
	int I,J,K;

	for (int Number = 0;Number<ISizeOfMatrix;++Number)
	{tmp[Number] = new float[ISizeOfMatrix];}

	for (int CurrDeg = 0;CurrDeg<Degree;CurrDeg++){
	for (int I = 0;I<ISizeOfMatrix;++I)
	{
	for (int J = 0;J<ISizeOfMatrix;++J)
	{
	tmp[I][J] = 0;
	for (K = 0;K<ISizeOfMatrix;++K)
	{
     tmp[I][J] += InitMatrix[I][K]*InitMatrix[K][J];
	}
	}
	}
	
	}
	InitMatrix = tmp;

[BDA]

Не тестировал:

float** tmp1 = new float*[ISizeOfMatrix];
float** tmp2 = new float*[ISizeOfMatrix];
float** tmp3 = NULL;
for (int I = 0;I<ISizeOfMatrix;++I) {
    tmp1[I] = new float[ISizeOfMatrix];
    tmp2[I] = new float[ISizeOfMatrix];
}
for (int I = 0;I<ISizeOfMatrix;++I)
    for (int J = 0;J<ISizeOfMatrix;++J)
        tmp2[I][J] = InitMatrix[I][J];
for (int CurrDeg = 1;CurrDeg<Degree;CurrDeg++) {
    tmp3 = tmp2;
    tmp2 = tmp1;
    tmp1 = tmp3;
    for (int I = 0;I<ISizeOfMatrix;++I)
        for (int J = 0;J<ISizeOfMatrix;++J) {
            tmp2[I][J] = 0;
            for (K = 0;K<ISizeOfMatrix;++K)
                tmp2[I][J] += tmp1[I][K]*InitMatrix[K][J];
            }
}
for (int I = 0;I<ISizeOfMatrix;++I) {
    delete[] tmp1[I];
    delete[] InitMatrix[I];
}
delete[] tmp1;
delete[] InitMatrix;
InitMatrix = tmp2;

[proef]

работает

Спасибо
сначала не понял что написано но потом разобрался)

[BDA]

Пожалуйста. Если будет у Вас время, можете реализовать так называемое быстрое возведение в степень.