Переполнение

[impersonalis]

Есть небольшая либа для решения матричной системы.
Всё работает нормально, пока при попытке инвертировать матрицу порядка 7 на 7, ф-ция аварийно не завершает работу, т.к. получает детерминант=NaN (хотя он чуть превышает 10^6).
Само собой разумеется, в остальных случаях всй считается.

Const TMatrixMaxIndex2D=99
Type TMatrix
	Field Buffer$[TMatrixMaxIndex2D]
	Field StringSize
	Field CSize
End Type

Function TMatrix_Create.TMatrix(size,size2=-1)
	M.TMatrix=New TMatrix
	M\StringSize=size
	M\CSize=size2
	If M\CSize=-1 M\CSize=size
	Return M
End Function

Function TMatrix_IN(M.TMatrix,SNumber,CNumber,Z$)
	Local Index2d=(SNumber-1)*M\StringSize+(CNumber-1)
	M\Buffer[Index2d]=Z$
End Function

Function TMatrix_OUT$(M.TMatrix,SNumber,CNumber)
	Local Index2d=(SNumber-1)*M\StringSize+(CNumber-1)
	Return M\Buffer[Index2d]
End Function

Function TMatrix_Delete(M.TMatrix)
	Delete M
End Function

Function TMatrix_Determinant#(M.TMatrix)
	
	SS=M\StringSize
	CS=M\CSize
	If SS=1 And CS=1
 *Return Float(M\buffer[0])
	EndIf
	Local DET#=0
	For i=1 To SS
 *MINOR.TMATRIX=TMatrix_minor(M,1,i)
 *DET=DET+(-1)^(i+1)*TMatrix_Determinant(MINOR)*Float(TMatrix_OUT(M,1,i))
 *TMatrix_Delete(MINOR)
	Next
	Return Float(DET)
End Function

Function TMatrix_minor.TMatrix(M.TMatrix,NS,NC)
	minSS=M\StringSize-1
	minCS=M\CSize-1
	Minor.TMatrix=TMATRIX_Create(minSS,minCS)
	Local Smin=1
	For S=1 To M\CSize
 *Local Cmin=1
 *For C=1 To M\StringSize
 *	If C=NC
 * *Cmin=Cmin-1
 *	Else
 * *TMatrix_IN(Minor,Smin,Cmin,TMatrix_OUT(M,S,C))
 *	EndIf
 *	Cmin=Cmin+1
 *Next
 *If S=NS
 *	Smin=Smin-1
 *EndIf
 *Smin=Smin+1
	Next
	Return Minor
End Function

вот - выдрал необходимые куски кода.
Собсна вопрос - можно как-нить уменьшить кол-во рекурсий, не сильно загромождая код?

лутче скажи что собствено код делает... ну я в математике матриц не очень... всештаки только 9 класс

[impersonalis]

проблема отпала

[impersonalis]

Проблему с нехваткой места в 4-байтовых переменных сначала хотел решить эмуляцией длинных чисел на основе массива коээфициентов. Но код получился слишком громоздкий (мне была необходима поддержка положительных/отрицательных дробных чисел.
Поэтому на данный момент пишу систему округлений вычислений, т.к. точность некритична а вот скокрость работы на эмуляции упадёт больше чем на округениях с вынесением степени десятки.

[pax]

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

[pax]

И вообще лучше и быстрее всего решать системы лин. уравнений методом Жордано-Гаусса... никаких переполнений переменных не будет...

'Процедура преобразования матрицы методом Жордана-Гаусса
Public Sub MAT(A() As Double, I1 As Integer)
 * If A(I1, I1) = 0 Then *'Проверка на "0" по главной диагонали
 * *If I1 = N Then R = True: Exit Sub
 * * For i = I1 + 1 To N
 * * *If A(i, I1) <> 0 Then
 * * * For j = 0 To N + 1
 * * * *b = A(I1, j): A(I1, j) = A(i, j): A(i, j) = b
 * * * Next j
 * * * If A(I1, I1) <> 0 Then Exit For
 * * *End If
 * * Next i
 * * *If A(I1, I1) = 0 Then R = 1: Exit Sub
 * * *If N < NNN Then
 * * * Form2.Text1 = Form2.Text1 + "МЕНЯЕМ СТРОКИ МЕСТАМИ" + vbCrLf: VIVOD A()
 * * *End If
 * End If
 * * *For j = N + 1 To I1 Step -1 'Получение единицы по главной диагонали
 * * * *A(I1, j) = A(I1, j) / A(I1, I1)
 * * *Next j
 * * For i = 0 To N
 * * *If i <> I1 Then
 * * * For j = N + 1 To I1 Step -1
 * * * *A(i, j) = A(i, j) - A(I1, j) * A(i, I1)
 * * * Next j
 * * *End If
 * * Next i
End Sub

Преобразований таких N, т.е. число переменных (цикл for I1=0 to N-1, где I1 - номер итерации матрица A() - имеет размер [N,N+1])

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

VB рулез, С++ отстой

[impersonalis]

1)"И вообще лучше и быстрее всего решать системы лин. уравнений методом Жордано-Гаусса... никаких переполнений переменных не будет..." - это типа, того , вводная, млин!
2)"Инвертированную матрицу можно найти намного быстрее используя дополнительную единичную путем преобразования обеих матриц методом Жордано-Гаусса (точнее методом пямоугольника, а еще точнее смесью)" - неизящно + потеря точности
3)"VB рулез, С++ отстой" - а на это я и обидеться могу. Тем блее не понял : причём здесь С++? На С++ это всё реализовано очень даже.

[pax]

2)"Инвертированную матрицу можно найти намного быстрее используя дополнительную единичную путем преобразования обеих матриц методом Жордано-Гаусса (точнее методом пямоугольника, а еще точнее смесью)" - неизящно + потеря точности

какие тут могут быть разговоры о точности с 4-х байтовыми переменными? зато никакого переполнения, да и при расчете определителей высших порядков у тебя всегда будет необходимость понизить степень, что ведет к понижению точности, или сделать что-то извращенное а не изящное, а этим способом можно алгоритм на одном экране уписать...

3)"VB рулез, С++ отстой" - а на это я и обидеться могу. Тем блее не понял : причём здесь С++? На С++ это всё реализовано очень даже.

извини, но я просто для себя решил недавно, что VB намного круче C++... просто дело в том, что я в будующем не собираюсь в одиночку изобретать супер программы, а для обычных прог VB просто идеальный язык (это мое мнение... можешь меня забросать булыжниками ), к томуже он круче C# по производительности, и не уступает по возможностям... если я и буду использовать C++ то только для математических рассчетов...

[impersonalis]

Я в будующем не собираюсь писать обычные программы (их и так навалом),а вот для супер-прог - С++ самое оно (это мое мнение... можешь меня забросать булыжниками =/ ), к тому же он круче VB по производительности и возможностям... если я буду использовать VB - то даже не знаю на кой фиг (вариант : готовые компонеты/библиотеки ?) ....