-
Path: news-archive.icm.edu.pl!agh.edu.pl!news.agh.edu.pl!newsfeed2.atman.pl!newsfeed.
atman.pl!.POSTED!not-for-mail
From: bartekltg <b...@g...com>
Newsgroups: pl.comp.programming
Subject: Re: RSM i spline
Date: Thu, 15 Nov 2012 18:26:49 +0100
Organization: ATMAN - ATM S.A.
Lines: 58
Message-ID: <k838kr$lr8$1@node1.news.atman.pl>
References: <509ee300$0$26682$65785112@news.neostrada.pl>
<k7olf5$rpm$1@news.task.gda.pl> <k7oo6p$3ut$1@news.task.gda.pl>
<50a082a2$0$1301$65785112@news.neostrada.pl>
<k7qgii$cqo$1@news.task.gda.pl>
<f...@g...com>
<k7ujqc$2gh$1@node1.news.atman.pl> <k7ukdi$1nb$1@news.task.gda.pl>
<50a37a59$0$1313$65785112@news.neostrada.pl>
<k7vvq2$62t$1@mx1.internetia.pl> <k8153f$duq$1@node2.news.atman.pl>
<k822k1$fpc$1@news.task.gda.pl> <k82c89$lg6$1@news.task.gda.pl>
<k82en7$ckq$1@zeus.man.szczecin.pl> <k82frk$3m3$1@news.task.gda.pl>
NNTP-Posting-Host: 144-mi3-6.acn.waw.pl
Mime-Version: 1.0
Content-Type: text/plain; charset=UTF-8; format=flowed
Content-Transfer-Encoding: 8bit
X-Trace: node1.news.atman.pl 1353000411 22376 85.222.69.144 (15 Nov 2012 17:26:51
GMT)
X-Complaints-To: u...@a...pl
NNTP-Posting-Date: Thu, 15 Nov 2012 17:26:51 +0000 (UTC)
User-Agent: Mozilla/5.0 (Windows NT 6.1; WOW64; rv:16.0) Gecko/20121026
Thunderbird/16.0.2
In-Reply-To: <k82frk$3m3$1@news.task.gda.pl>
Xref: news-archive.icm.edu.pl pl.comp.programming:200950
[ ukryj nagłówki ]W dniu 2012-11-15 11:23, AK pisze:
>
>> Oczywiście wszystko się zmienia, jeżeli mnożenie przez wagi jest pod
>> pętlą (a tak jest np. w implementacji Simpsona z Wiki -
>
> Chyba tylko jakis _kompletny nieuk_ mnozylby nie wiadomo po co pod petla
> :)).
Ojtam, ja nie tylko mnożyłem pod pętlą, ale i budowałem
wielkie wektory tych wag - do zastosowań bojowych zupełnie
bez sensu;)
Roman W w swoim rozwiązaniu wstawił mnożenie pod pętlą.
Jak całkujesz skomplikowaną funkcję zawierajacą
takie kwiatki: http://mathworld.wolfram.com/Polylogarithm.html
to dodatkowe mnożenie niewiele zmienia.
Dziwi mnie ta gałąź. Sławek się burzy, a wjego kodzie nie
było mnożenia pod pętlą:
double integrate_simpson(double x[], double y[], int n)
{
int i;
double odd = 0., even = 0.;
for(i = 1 ; i <= n ; i+=2) odd += y[i];
for(i = 2 ; i <= n ; i+=2) even += y[i];
return (2.0*odd-y[1]-y[n] + 4.0*even)*(x[3]-x[1])/6.0;
}
Jak ktoś nie wierzy kompilatorowi, można wprost mu kazać:
//n parzyste!
double even=-0.5*y[0];
double odd = 0.0;
for(i = 1 ; i <= n ; i+=2)
{
odd += y[i];
even += y[i-1]
}
return (2.0*even + 4.0*odd + y[n])*stale
Albo od razu
if (n%2==0) return (2.0*even + 4.0*odd + y[n])*stale;
else return (2.0*even + 4.0*odd + a*y[n] + b*y[n-1] + c*[n-2] )*stale;
gdzie a,b,c to współczynniki, o których pisałem wczoraj.
pzdr
bartekltg
Następne wpisy z tego wątku
- 15.11.12 18:31 slawek
- 15.11.12 19:22 AK
- 15.11.12 19:27 AK
- 15.11.12 19:45 AK
- 15.11.12 19:48 slawek
- 15.11.12 19:58 slawek
- 15.11.12 20:06 AK
- 15.11.12 20:08 slawek
- 15.11.12 20:24 slawek
- 15.11.12 20:30 AK
- 15.11.12 20:33 AK
- 15.11.12 20:46 AK
- 15.11.12 20:53 R.e.m.e.K
- 15.11.12 21:58 AK
- 16.11.12 01:00 bartekltg
Najnowsze wątki z tej grupy
- TCL - problem z escape ostatniego \ w nawiasach {}
- Nauka i Praca Programisty C++ w III Rzeczy (pospolitej)
- testy-wyd-sort - Podsumowanie
- Tworzenie Programów Nieuprzywilejowanych Opartych Na Wtyczkach
- Do czego nadaje się QDockWidget z bibl. Qt?
- Bibl. Qt jest sztucznie ograniczona - jest nieprzydatna do celów komercyjnych
- Co sciaga kretynow
- AEiC 2024 - Ada-Europe conference - Deadlines Approaching
- Jakie są dobre zasady programowania programów opartych na wtyczkach?
- sprawdzanie słów kluczowych dot. zła
- Re: W czym sie teraz pisze programy??
- Re: (PDF) Surgical Pathology of Non-neoplastic Gastrointestinal Diseases by Lizhi Zhang
- CfC 28th Ada-Europe Int. Conf. Reliable Software Technologies
- Młodzi programiści i tajna policja
- Ada 2022 Language Reference Manual to be Published by Springer
Najnowsze wątki
- 2024-11-11 Wyważanie kół rowerowych
- 2024-11-11 Kosz, gdzie??
- 2024-11-11 Coraz mniej ludzi robi prawo jazdy
- 2024-11-11 Opole => SAP HANA Cloud Dev / Data Engineer <=
- 2024-11-11 Warszawa => Spedytor Międzynarodowy <=
- 2024-11-11 Lublin => Senior PHP Developer <=
- 2024-11-11 Marki => Senior PHP Symfony Developer <=
- 2024-11-11 Chrzanów => Team Lead / Tribe Lead FrontEnd <=
- 2024-11-11 Gliwice => Specjalista ds. public relations <=
- 2024-11-11 Gdańsk => Kierownik Działu Spedycji Międzynarodowej <=
- 2024-11-11 Gdańsk => Head of International Freight Forwarding Department <=
- 2024-11-11 Warszawa => Sales Development Representative (in German) <=
- 2024-11-11 Marsz niepodległości
- 2024-11-08 Belka
- 2024-11-09 pierdolec na punkcie psa