eGospodarka.pl
eGospodarka.pl poleca

  • 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



Podziel się

Poleć ten post znajomemu poleć

Wydrukuj ten post drukuj


Następne wpisy z tego wątku

Najnowsze wątki z tej grupy


Najnowsze wątki

Szukaj w grupach

Eksperci egospodarka.pl

1 1 1

Wpisz nazwę miasta, dla którego chcesz znaleźć jednostkę ZUS.

Wzory dokumentów

Bezpłatne wzory dokumentów i formularzy.
Wyszukaj i pobierz za darmo: