Обратная неквадратная матрица, ayudame por favor! )

[lecume]

Всем добрый вечер!
Ребят, может у кого-то был опыт нахождения обратной матрицы для имеющейся неквадратной? Для квадратной код есть, все работает. Могу выложить, но, думаю, для неквадратной подход совсем другой.
Подскажите, пожалуйста!

int gaus_ostatok; // дефект линейного преобразования
float gaus_deter; //  в случае успешного обращения детерминант
float gaus_minved; //  минимальный ведущий элемент, можно использовать для оценки точности

#define _NaN() (float&)*"Implementation dependent"; // заполнитель для неразрешенных измерений //  при невозможности обращения

inline static  void  invers(float m[100][100], int N)
   {
     int rn[100], cn[100];
     int j, k;
     for(j = N; j--; )  rn[j] = cn[j] = j;
     gaus_minved = 1e99;
     gaus_deter = 1;
     for (gaus_ostatok = N; gaus_ostatok; gaus_ostatok--)
     {  
		 int jved, kved;
		 float vved = -1,t; // поиск ведущего
		 for (j = N; j--; )
		 {  
			if(~rn[j])
			for (k = N; k--;)
			  if(~cn[k])
			  if(vved < fabs(m[j][k]))  vved = fabs(m[j][k]), jved = j, kved = k;
		 }        
		 if(gaus_minved > vved)  gaus_minved = vved;
		 gaus_deter*=m[jved][kved];
		 if(vved < 1e-99)
		 {  
			for (j = N; j--; )
			{  
				if(~rn[j]) for (k = N; k--; ) m[j][k] = _NaN();
				if(~cn[j]) for(k = N; k--; ) m[k][j] = _NaN();
			}
		 return;
		 }
        int jt = rn[jved], kt = cn[kved];  // перестановки
        for (j = N; j--; )  t = m[kt][j], m[kt][j] = m[jved][j], m[jved][j] = t;
        for (j = N; j--; )  t = m[j][jt], m[j][jt] = m[j][kved], m[j][kved] = t;
        rn[jved] = rn[kt];
        cn[kved] = cn[jt];
        rn[kt] = cn[jt] = -1;
        vved = m[kt][jt];  
	    m[kt][jt] = 1;
        for (j=N;j--;)
        {  
		   if(j == kt)  continue;
// if (fabs (vved) < FLT_EPSILON) {}
           float mul = m[j][jt] / vved;
           m[j][jt] = 0;
           for (k = N; k--; ) m[j][k] -= m[kt][k] * mul; // самый внутренний цикл ровно N^3 операций
        }
        for (k = N; k--; ) m[kt][k] /= vved;
     }
   }

[ex.cluz]

Посмотрите функцию ReverseMatrix(MPtr : MatrixPtr) : MatrixPtr;
в теме: http://www.programmersforum.ru/showt...t=70310&page=3
посты #25 - #28.
Там, правда на Паскале, но если с ним знакомы - переведете на Си.

[lecume]

Спасибо большое! Я поиском, честно, пользовалась
С Паскалем незнакома совершенно. Для знающего человека много времени займет перевести это на Си/С++?
P. S. С не квадратными матрицами модуль хорошо работает?

[ex.cluz]

Я на нем прогу по стат.анализу делал (реализовывал метод множественной линейной регрессии для обработки данных, а там, если мне не изменяет память, обращение идет произведения матриц, а при произведении получается квадратная).

А вообще, кто теорию читать будет?

Для неквадратных матриц и вырожденных матриц обратных матриц не существует.

[Vago]

Сообщение от ex.cluz

А вообще, кто теорию читать будет?

Очень верное замечание!

Сообщение от ex.cluz

Для неквадратных матриц и вырожденных матриц обратных матриц не существует.

Для неквадратных матриц, при выполнении определённых условий, существуют понятия левой [псевдо-]обратной (left [pseudo-]inverse ) и правой [псевдо-]обратной (right [pseudo-]inverse) матрицы. Метод Мура-Пенроуза (Moore-Penrose method) это...

[lecume]

Точно-точно! Простите, 2-ой курс был давно
А программа эта все равно пригодится.

[Vago]

Armadillo.
A pesar de su título Español, todo es en idioma Inglés