-
Path: news-archive.icm.edu.pl!newsfeed.gazeta.pl!opal.futuro.pl!not-for-mail
From: AMX <r...@b...cy>
Newsgroups: pl.sci.inzynieria
Subject: Re: Duuuuuuuuuuuuże układy równań
Date: Sat, 5 Dec 2009 16:04:33 +0000 (UTC)
Organization: Home Access
Lines: 61
Message-ID: <s...@a...within.nat>
References: <f...@j...googlegroups.com>
<s...@a...within.nat>
<q...@4...com>
NNTP-Posting-Host: 81.219.244.142
Mime-Version: 1.0
Content-Type: text/plain; charset=iso-8859-2
Content-Transfer-Encoding: 8bit
X-Trace: opal.futuro.pl 1260029073 28584 81.219.244.142 (5 Dec 2009 16:04:33 GMT)
X-Complaints-To: u...@f...pl
NNTP-Posting-Date: Sat, 5 Dec 2009 16:04:33 +0000 (UTC)
User-Agent: slrn/pre0.9.9-111 (Linux)
Xref: news-archive.icm.edu.pl pl.sci.inzynieria:25627
[ ukryj nagłówki ]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
Następne wpisy z tego wątku
- 14.12.09 10:22 Jacek_P
Najnowsze wątki z tej grupy
- Filtr do pompy ruskiej
- Wyważanie kół rowerowych
- Belka
- Precyzyjne cięcie opony samochodowej
- Nieparzyste dmuchanie
- Klej "samopoziomujący"
- Kocioł CO po raz kolejny
- zapora Stronie Slaskie
- powodz
- Nie atom tylko fotowoltanika i elektroliza
- Test samoodkręcania nakrętek
- Budowlańcy pomóżcie
- wyciskanie/odlewanie hdpe. Co ma sens?
- Pomysł na czujnik przeciążenia siłownika.
- dokładne poziomowanie
Najnowsze wątki
- 2024-11-13 Filtr do pompy ruskiej
- 2024-11-12 Gdzie kosz?
- 2024-11-13 elektrycznie
- 2024-11-12 Jebane kurwa, kurwy.
- 2024-11-13 karta parkingowa
- 2024-11-13 Wl/Wyl (On/Off) bialy/niebieski
- 2024-11-12 I3C
- 2024-11-13 Kraków => DevOps Engineer (Junior or Regular level) <=
- 2024-11-13 Łódź => Senior SAP HANA Developer <=
- 2024-11-13 Zabrze => Senior PHP Symfony Developer <=
- 2024-11-13 Karlino => Konsultant wewnętrzny SAP (FI/CO) <=
- 2024-11-13 Kraków => QA Inżynier <=
- 2024-11-13 Żerniki => Dyspozytor Międzynarodowy <=
- 2024-11-13 Warszawa => Analityk Biznesowo-Systemowy <=
- 2024-11-13 Lublin => Delphi Programmer <=