eGospodarka.pl

eGospodarka.plGrupypl.comp.programming › Dokładny podział integera na k jednakowych integerów
Ilość wypowiedzi w tym wątku: 4

  • 1. Data: 2018-01-30 18:39:39
    Temat: Dokładny podział integera na k jednakowych integerów
    Od: Borneq <b...@a...hidden.pl>

    To coś co ma wiele wspólnego przy powiększaniu rastra metodą Nearest
    przy niecałkowitym powiększeniu.

    Metoda naiwna:
    void divide_float(int n, int k)
    {
    double interval = double(n)/k;
    //double
    for (int i=0; i<k; i++)
    {
    printf("%d %f %f\n",i,i*interval, (i+1)*interval);
    }
    }

    chcę działać na integerach ale i tak mam kłopot z zaokrągleniami:
    void divide_int(int n, int k)
    {
    double interval = double(n)/k;
    //int
    uint32_t intPartPos, fracPartPos;
    uint32_t intPartDelta,fracPartDelta,oldFrac;
    intPartPos = 0; fracPartPos = 0;
    intPartDelta = floor(interval);
    fracPartDelta = floor((interval - intPartDelta) * 65536 * 65536);
    for (int i = 0; i<k; i++)
    {
    intPartPos += intPartDelta;
    oldFrac = fracPartPos;
    fracPartPos += fracPartDelta;
    if (oldFrac > fracPartPos) // when LongWord overflow
    intPartPos++;
    printf("%d %d\n", i, intPartPos);
    }
    }

    divide(10, 3) zwraca 3,6,9 zamiast 3,6,10 dlatego że w ostatnim jest
    błąd rzędu 2^-32 i to wystarczy by było źle.


  • 2. Data: 2018-01-30 18:54:25
    Temat: Re: Dokładny podział integera na k jednakowych integerów
    Od: Borneq <b...@a...hidden.pl>

    W dniu 30.01.2018 o 18:39, Borneq pisze:
    >     fracPartDelta = floor((interval - intPartDelta) * 65536 * 65536);

    fracPartDelta = ceil((interval - intPartDelta) * 65536 * 65536);

    ceil zamiast floor tutaj pomaga ale to niezbyt eleganckie rozwiązanie,
    trzeba nie używać interval ale wymiernego n/k, jak?


  • 3. Data: 2018-01-30 18:56:26
    Temat: Re: Dokładny podział integera na k jednakowych integerów
    Od: Borneq <b...@a...hidden.pl>

    W dniu 30.01.2018 o 18:54, Borneq pisze:
    > W dniu 30.01.2018 o 18:39, Borneq pisze:
    >>      fracPartDelta = floor((interval - intPartDelta) * 65536 * 65536);
    >
    > fracPartDelta = ceil((interval - intPartDelta) * 65536 * 65536);
    >
    > ceil zamiast floor tutaj pomaga ale to niezbyt eleganckie rozwiązanie,
    > trzeba nie używać interval ale wymiernego n/k, jak?
    ceil - nie!
    a może to ma coś wspólnego z algorytmem Bresenhama?


  • 4. Data: 2018-01-30 19:35:49
    Temat: Re: Dokładny podział integera na k jednakowych integerów
    Od: Borneq <b...@a...hidden.pl>

    W dniu 30.01.2018 o 18:56, Borneq pisze:
    >> ceil zamiast floor tutaj pomaga ale to niezbyt eleganckie rozwiązanie,
    >> trzeba nie używać interval ale wymiernego n/k, jak?
    > ceil - nie!

    Rozwiązanie:
    Chodzi o przydział wątków, więc mogę "po chamsku" wstawić ostatnią wartość:
    vector<uint32_t> getRanges(int n, int k, int mintask)
    {
    double interval = double(n) / k;
    if (interval<mintask)
    {
    k = n / mintask;
    interval = double(n) / k;
    }
    vector<uint32_t> result;
    result.push_back(0);
    uint32_t intPartPos, fracPartPos;
    uint32_t intPartDelta, fracPartDelta, oldFrac;
    intPartPos = 0; fracPartPos = 0;
    intPartDelta = (uint32_t)floor(interval);
    fracPartDelta = (uint32_t)floor((interval - intPartDelta) * 65536 * 65536);
    for (int i = 0; i<k - 1; i++)
    {
    intPartPos += intPartDelta;
    oldFrac = fracPartPos;
    fracPartPos += fracPartDelta;
    if (oldFrac > fracPartPos) // when LongWord overflow
    intPartPos++;
    result.push_back(intPartPos);
    }
    result.push_back(n);
    return result;
    }


    mintask to liczba , ile minimalnie obiektami może zająć się wątek, aby
    nie było nagle 1 czy 0.

strony : [ 1 ]



Szukaj w grupach

Szukaj w grupach

REKLAMA

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:

Strona korzysta z plików cookies w celu realizacji usług i zgodnie z Polityką Prywatności. Możesz określić warunki przechowywania lub dostępu do plików cookies w Twojej przeglądarce.