On Sat, 05 Dec 2009 11:55:59 +0100,
J.F <j...@p...onet.pl> wrote:

>>> chciałbym prosić o pomoc. Zajmuję się MES-em i jak wiadomo, uzyskuje
>
>>Co masz na myśli pisząc duże? Takie 1k równań matlab
>>rozwiązuje w sensownym czasie jako pełne.
>
> Mowisz o zagadnieniach typu MES, czy ogolnych ?

W zasadzie to jest dla ilustracji. Przy macierzy pełnej
ma tylko znaczenie czy jest symetryczna i czy wymaga
wyboru elementu wiodącego. Elementy zerowe i tak są
przetwarzane. Tak więc dotyczy to zarówno MES jak i
ogólnych, o ile uwzględnimy rozciągliwość pojęcia
,,sensowny czas''.

>
>>Takie około 100k to według mnie tylko LLT.
>>Powyżej 1M równań zaczynają się schody, LLT produkuje
>>wypełnienia i krytyczna jest renumeracja.
>
> 1M rownan to 1T wspolczynnikow. Trzeba duuuuzo pamieci :-)
>

Dokładnie 64TB pamięci. Dlatego nie da się większych
niż choćby owe 1k równań rozwiązywać układu bez
uwzględnienia zer w macierzy. Dla MEB w ogóle rząd
zadania typu 10k jest poza zasięgiem (macierz pełna) i
dlatego MES króluje (rzadka, pasmowa macierz
sztywności), 99.99% elementów macierzy to zera.

> Czy mi sie wydaje czy problemy dokladnosci numerycznej zaczynaja sie
> juz przy dziesiatkach rownan ?
>

To bardzo zależy od konkretnego przypadku. Jeśli
weźmiesz najprostsze zadanie pręta rozciąganego i
włożysz choćby 1000 elementów _równej długości_ to się
bardzo ładnie policzy i to z idealną dokładością.

Jeśli te sam pręt podzielisz w ten sposób, że każdy
następny element będzie dwa razy krótszy (1, 0.5, 0.25,
0.125 itd. ), to gdzieś między 32 a 64 elementem (nie
chce mi się dokładnie liczyć) _każdy_ solver się
wyłoży. Układ będzie po prostu osobliwy. A zanim się
wyłoży (czyli np. dla 30 elementów) rozwiązanie
ucieknie ci w kosmos.

To, kiedy zaczynają się problemy numeryczne nie zależy
tylko od ilości równań (problemy mogą być już przy dwu
równaniach), ale też od wielu innych czynników. Nie ma
prostej recepty.

Czołem

amx

--
adres w rot13
Nyrxfnaqre Znghfmnx r...@b...cy