[Drongo]

Нахождение обратной матрицы методом расширенной матрицы

Код:

//Нахождение обратной матрицы методом расширенной матрицы
#include <vcl.h>
#include <conio.h>
#include <iostream.h>
#include <math.h>
#pragma hdrstop
#pragma argsused

const int dim1=20;
const int dim2=40;
float**a,**a1;
int n;
//vvod matrici
void MatrIn()
{
   int i,
       j;

   for(i = 0; i < n; i++){
     for(j = 0; j < n; j++){
        cout<<"vvedite cherez probel elementi metrici[ "<<i<<" ][ "<<j<<" ] = ";
        cin>>a[i][j];
      }
    }
}
//obraschenie matrici
void invMatr(float e)
{
   int i,
       j,
       k,
       n1;
    float f,
          f1;
    float matr[20][40];
    n1 = 2 * n;

    for(int i = 0; i < n; i++)
      for(j = 0; j < n1; j++)
         if(j < n)
           matr[i][j] = a[i][j];
         else if(j == (n + 1))
           matr[i][j] = 1;
         else
           matr[i][j] = 0;

    for(int i=0;i<n;i++){
       k = i;
       f1 = matr[i][j];
       while(fabs(f1) < e){
          if(k < n){
            k++;
            f1 = matr[k][i];
           }
          else{
            cout<<"matrica virojdena";
            getch();
           }
       for(j = i; j < n1; j++){
         f = matr[k][j];
         matr[k][j] = matr[i][j];
         matr[i][j] = f / f1;
       }

       for(k = 0; k < n; k++)
         if(k != i)
           for(j = n1 - 1; j >= i; j--)
              matr[k][j] = matr[k][j] - matr[i][j] * matr[k][j];
     }

       for(i = 0; i < n; i++)
          for(j = 0; j < n; j++)
            a[i][j] = matr[i][j + n];
  }
}
//vivod obratnoy matrici
void MatrOut()
{
   int i,
       j;

   for(i = 0; i < n; i++)
      for(j = 0; j < n; j++)
         cout<<a1[i][j]<<" ";
}
//glavnaya proga
void main()
{
   int e;

   cout<<"vvedite razm matrici: ";
   cin>>n;

   a = new float* [n];

   MatrIn();

   a1 = new float* [n];

   for(int i = 0; i < n; i++)
      a1[i] = new float[n];

   cout<<"vvedite tochnost' vichisleniya: ";
   cin>>e;
   invMatr(e);
   cout<<"obratnaya matrica: ";
   MatrOut();
   getch();
}
//--------------------------------------------------------------------------

Код:

Вопрос, что делает эта строка кода:
a=new float*[n];  // эта строка
MatrIn();
a1=new float*[n];  // И эта строка

А так же, что это? Указатель на указатель?

Код:

float **a,**a1;

Ошибка появляется после ввода данных, на этапе

Код:

a=new float*[n];  // эта строка

А так всё комплируется и до этой строки работает.

[D.Y.]

Вопрос, что делает эта строка кода:
a=new float*[n]; // эта строка
MatrIn();
a1=new float*[n]; // И эта строка
это создается новая матрица типу float

[Drongo]

D.Y.,

это создается новая матрица типу float »

Динамическое выделение

Код:

float *array = new float[size];

[Drongo]

D.Y., Тогда, что делает эта строка в главной функции main()?

Код:

 for(int i = 0; i < n; i++)
      a1[i] = new int[n];

У меня просто у самого знаний нехватает.

Эти строки не используются в программе

Код:

const int dim1=20;
const int dim2=40;

В общем, ошибка где-то здесь:

Код:

...
   a1 = new float* [n];

   for(int i = 0; i < n; i++)
      a1[i] = new float[n];
...

[Admiral]

D.Y., в первом посте неплохо бы было указать сразу тему с другого форума есть рабочая прога(содрана с книги) не могу понять под какую среду писалась так сказать доложить о проделанных попытках по преодолении сложности, здесь бы дублей по догадкам не было.
А книга именно так называется, может всё же полное название Языки Си С++ для решения инженерных и экономических задач?

if (j!=i) a[k][i]-=b[j]*c[k];} for (i=0;i<n;i++){ »

Выделенная скобка закрывает функцию int InvMatrG(int n, float(*a)[dim], float e,float d) всё что дальше в кодовом космосе находится.

[Drongo]

Admiral,

оффтоп

бы было указать сразу тему с другого форума »

А как вы его нашли там?

D.Y.,

оффтоп

D.Y., так поставь себе задачу: выбрось этот жуткий, кривой код и напиши свой, красивый и работающий

и

Придется так и сделать но я хотел не напрягаясь содрать

Не обессудьте, но таки придётся написать самому. Хоть бы сказали, что вам его подправили здесь на OSZone, а то: "я подправил." Всё-таки первое впечатление

В общем, я точно уверен, что вы не писали эту программу »

правильное.

[Admiral]

Drongo, уточнял книгу и главу в которой сабж и вышел на ту тему.
http://www.google.com/search?hl=ru&c...1%83%D0%BA&lr=

[Drongo]

Код:

...
void MatrIn()
{
   int i,
       j;

   for(i = 0; i < n; i++){
     for(j = 0; j < n; j++){
        cout<<"vvedite cherez probel elementi metrici[ "<<i<<" ][ "<<j<<" ] = ";
        cin>>a[i][j];
      }
    }
}
...

Как я понял в первом куске кода выполняется присвоение массиву размером n * n если n равна 2, то получается 4, в то время, как в куске кода ниже, память выделяется только под n = 2

Код:

...
cin>>n;

   a = new float* [n];

   MatrIn();

   a1 = new float* [n];

   for(int i = 0; i < n; i++)
      a1[i] = new float[n];
...

Хотя если честно я не знаю, проверил этот способ, но всё равно ошибка.

[D.Y.]

Код:

а не моглибы вы прокомментировать(что в какой строке делается) программу я взял у паренька но он кудато пропал и я не пойму что в ней творится особенно в конце
было задание решение СЛАУ методом Крамера
#include <iostream.h>  
#include <stdio.h> 
#include <math.h> 
double DET (unsigned int columns);
bool PROVCOL(unsigned int num, char colNum);
unsigned int Column1(unsigned int num, char colNum);
unsigned int Column0(unsigned int num, char colNum);
double **Matrix;
int n;
main()
{
     int i, j; 
     double mainDET; 
     cout<<"Vvedite kolichestvo uravneniy "<<endl;
     cin>>n;
     while (n > 10 || n < 1)
     {
          cout<<"Vu ne mogete vvesti bol'she 10 ili men'she 1 uravneniya, vvedite esche raz ";
          cin>>n;
     }
     
     Matrix = (double**) new double[n];
     for (i=0; i<= n-1; i++)
     {Matrix[i] = new double[n+1];}
     for (i = 0; i <= n-1; i++)
               for (j = 0; j <= n; j++)
               {
                    cout<<"Vvedite element matricu s koordinatami ("<<i<<","<<j<<") ";
                    cin>>Matrix[i][j];
               }
     cout<<"Vasha matrica:\n";
     for (i = 0; i <= n-1; i++)
     {for (j = 0; j <= n; j++)
          cout<<Matrix[i][j]<<"\t";
          cout<<endl;
     }
     mainDET = DET(pow(2,n));
     if (mainDET == 0)
     {
		 cout<<"Net resheniy!\n";
     }
     else
     {
      printf ("Korni sistemu:\n");
               for (i=0; i <= n-1; i++)
               printf ("%.3f ", -pow(-1,n+i)*DET(pow(2,i))/mainDET);
          printf ("\n");
     }
	 cout<<"Detrmenant sistemu raven: "<<mainDET<<endl;
     for (i = 0; i <= (n-1); i++) 
     delete Matrix[i];
     delete Matrix;
     cout<<"Vuchisleniya zakonchenu. Spasibo chto vubrali nas!\n";
     return 0;
}


double DET (unsigned int columns)
{double det = 0;
    int i, numCols = 0, lastFalse, cntr = 0;
     for (i=0; i<=n; i++) 
     {
	  if (PROVCOL(columns, i)) 
      numCols++;
      else lastFalse = i;
     }
    if ((n - numCols) == 0) det = Matrix[n-1][lastFalse];
     else
         for (i = 0; i <= n; i++)
               if (!PROVCOL(columns, i))
               {
					columns = Column1(columns, i);
                    det += pow(-1,(cntr))*Matrix[numCols-1][i]*DET(columns);
                    columns = Column0(columns, i);
                    cntr++;
               }
     return det;
}
bool PROVCOL(unsigned int num, char colNum)
{
     _asm
     {
          xor eax, eax; 
          xor ebx, ebx; 
          mov bl, colNum; 
          bt num, ebx; 
          jnc none; 
          inc eax; 
          none:
     }
}
unsigned int Column1(unsigned int num, char colNum)
{
     _asm
     {
          xor eax, eax; 
          mov al, colNum; 
          bts num, eax; 
     }
     return num; 
}
unsigned int Column0(unsigned int num, char colNum)
{
     _asm
     {
          xor eax, eax; 
          mov al, colNum; 
          btr num, eax; 
     }
     return num; 
}

[Drongo]

D.Y.,

в этих строках

Код:

bool PROVCOL(unsigned int num, char colNum)
{
     _asm
     {
          xor eax, eax; 
          xor ebx, ebx; 
          mov bl, colNum; 
          bt num, ebx; 
          jnc none; 
          inc eax; 
          none:
     }
}
unsigned int Column1(unsigned int num, char colNum)
{
     _asm
     {
          xor eax, eax; 
          mov al, colNum; 
          bts num, eax; 
     }
     return num; 
}
unsigned int Column0(unsigned int num, char colNum)
{
     _asm
     {
          xor eax, eax; 
          mov al, colNum; 
          btr num, eax; 
     }
     return num;

Ассемблерная вставка _asm

В остальном, чужой код комментировать трудно. Простой ещё можно, а тут, поможет только автор.
Вот тут добавь строку cin>>i; иначе будет программа закрываться после выполнения

Код:

...
 delete Matrix;
     cout<<"Vuchisleniya zakonchenu. Spasibo chto vubrali nas!\n";
     cin>>i;
     return 0;
...