abclinuxu.cz AbcLinuxu.cz itbiz.cz ITBiz.cz HDmag.cz HDmag.cz abcprace.cz AbcPráce.cz
AbcLinuxu hledá autory!
Inzerujte na AbcPráce.cz od 950 Kč
Rozšířené hledání
×
    dnes 16:11 | Humor

    Win8DE je desktopové prostředí pro Wayland, inspirované nechvalně proslulým uživatelským rozhraním Metro z Windows 8. Nabízí dlaždicové rozhraní s velkými tlačítky a jednoduchou navigací, optimalizované pro dotyková zařízení. Cílem projektu je přetvořit design operačního systému Windows 8 do funkčního a minimalistického rozhraní vhodného pro každodenní použití na Linuxu.

    NUKE GAZA! 🎆 | Komentářů: 7
    dnes 14:33 | Nová verze

    Laboratoře CZ.NIC vydaly Datovku 4.28.0 a Mobilní Datovku 2.6.0. Hlavní novinkou je ukládání rozpracovaných datových zpráv do konceptů. Datovka je svobodné multiplatformní aplikace pro přístup k datovým schránkám a k trvalému uchovávání datových zpráv v lokální databázi.

    Ladislav Hagara | Komentářů: 0
    dnes 14:22 | Zajímavý projekt

    Unix Pipe Game je vzdělávací karetní hra zaměřená na děti a rodiče, která děti učí používat unixové příkazy prostřednictvím interaktivních úkolů. Klíčovým prvkem hry je využití symbolu | pro pipeline neboli 'rouru', který umožňuje propojit výstupy a vstupy jednotlivých unixových příkazů, v tomto případě vytištěných na kartičkách. Předpokládá se, že rodič má alespoň nějaké povědomí o unixových příkazech a jejich provazování pomocí |.

    … více »
    NUKE GAZA! 🎆 | Komentářů: 3
    dnes 13:55 | Zajímavý software

    PCIem je linuxový framework, který vytváří virtuální zařízení PCIe pomocí technik, které umožňují hostitelskému operačnímu systému rozpoznat tyto syntetické 'neexistující' karty jako fyzické zařízení přítomné na sběrnici. Framework PCIem je primárně zamýšlen jako pomůcka pro vývoj a testování ovladačů bez nutnosti použít skutečný hardware. Dle tvrzení projektu si fungování PCIem můžeme představit jako MITM (Man-in-the-Middle), který se nachází mezi ovladači a kernelem.

    NUKE GAZA! 🎆 | Komentářů: 0
    dnes 05:33 | Bezpečnostní upozornění

    Byla nalezena vážná bezpečnostní chyba v telnetd z balíčku GNU InetUtils. Týká se verzí GNU InetUtils od 1.9.3 z 12. května 2015 až po aktuální 2.7 z 14. prosince 2025. Útočník může obejít autentizaci a získat root přístup, jelikož telnetd nekontroluje předaný obsah proměnné prostředí USER a pokud obsahuje "-f root"…

    Ladislav Hagara | Komentářů: 5
    včera 17:22 | Komunita

    Stanislav Aleksandrov předložil patch rozšiřující KWin (KDE Plasma) na 3D virtuální desktopové prostředí (videoukázka v mp4).

    Ladislav Hagara | Komentářů: 5
    včera 05:55 | IT novinky

    Digg (Wikipedie), "místo, kde můžete sdílet a objevovat to nejlepší z internetu – a nejen to", je zpět. Ve veřejné betě.

    Ladislav Hagara | Komentářů: 4
    včera 01:44 | Komunita

    Po .deb balíčcích Mozilla nově poskytuje také .rpm balíčky Firefoxu Nightly.

    Ladislav Hagara | Komentářů: 9
    včera 00:33 | IT novinky

    Vývojové prostředí IntelliJ IDEA slaví 25. narozeniny (YouTube).

    Ladislav Hagara | Komentářů: 0
    19.1. 23:55 | IT novinky

    Vedení společnosti NVIDIA údajně povolilo použití milionů knih ze známého 'warez' archivu Anna's Archive k výcviku umělé inteligence, ačkoliv vědělo, že archiv tyto knihy nezískal legální cestou. Žaloba, ve které se objevují i citace interních dokumentů společnosti NVIDIA, tvrdí, že NVIDIA přímo kontaktovala Anna's Archive a požadovala vysokorychlostní přístup k datům knihovny.

    NUKE GAZA! 🎆 | Komentářů: 24
    Které desktopové prostředí na Linuxu používáte?
     (18%)
     (5%)
     (0%)
     (10%)
     (20%)
     (3%)
     (5%)
     (2%)
     (11%)
     (37%)
    Celkem 557 hlasů
     Komentářů: 16, poslední včera 17:04
    Rozcestník

    Dotaz: Ukazatele v C, pole doublu

    25.12.2010 21:56 Preston Pigglesworth
    Ukazatele v C, pole doublu
    Přečteno: 408×
    Přílohy:
    Dobry den, nevim si rady z nasledujicim problemem v C. Z knihovny lapack mam funkci pro reseni soustavy linearnich rovnic _dgesv, volam ji presne podle prikladu z http://software.intel.com/sites/products/documentation/hpc/mkl/lapack/mkl_lapack_examples/dgesv_ex.c.htm

    Napsal jsem mensi program na jeji odzkouseni, viz priloha lapack2.c -- funguje OK

    Kdyz ji ale pouzivam tam kde ji potrebuju, tak jsou vysledky pokazde spatne, krome prvniho prvku. Napada me chyba v ukazatelich, ale podle clanku, co jsem procital tady i jinde, by to melo byt v poradku.Zkousel jsem i zapisovat do souboru jednotlive sekvence vypoctu.. vsechno je v poradku dokud nedojde na _dgesv. Ta sice spravne zapise vysledk do vstupniho pole B, ale s naprosto spatnymi cisly! Zkopiroval jsem obsah souboru B.txt a S.txt do zkouseciho programu a obdrzel jsem spravny vysledek.[overeny GNU/octave]

    Muzete mi prosim vysvetlit, jaky je z pohledu te funkce rozdil v tom jestli ji dam odkaz pole, pevne dane velikosti s na pevno naplnenymi cisly, nebo odkaz dynamicky alokovane pole?


    Řešení dotazu:


    Odpovědi

    25.12.2010 22:33 Michal Kubeček | skóre: 71 | Luštěnice
    Rozbalit Rozbalit vše Re: Ukazatele v C, pole doublu
    Nezkoumal jsem to nijak podrobně, ale problém by mohl být spíš v těch ostatních parametrech, které nejsou pointery na double. V jedné ukázce totiž u všech předáváte pointery na int, zatímco v druhé pointery na long.
    25.12.2010 22:56 Dr. Eddy | skóre: 9 | blog: glog | České Budějovice
    Rozbalit Rozbalit vše Re: Ukazatele v C, pole doublu
    no, kdyz uz nejsou srozumitelne nazvy promennych, tak bys aspon ten kod mohl okomentovat. aspon nejak zhruba, abychom vedeli, co se kde deje. A velmi prakticke by bylo napsat, co je ktera promenna (vstupni / vystupni data, ridici promenne atd.).
    25.12.2010 23:46 Preston Pigglesworth
    Rozbalit Rozbalit vše Re: Ukazatele v C, pole doublu
    dgesv je subroutina na reseni soustavy rovnic A*X=b , pritom pri volani dgesv( &n, &nrhs, a, &lda, ipiv, b, &ldb, &info ) je A predstavovano pomoci pole a. b je na VSTUPU predstavovano pomoci pole b. Na VYSTUPU se reseni X zapisuje do pole b! [Myslim, ze takovy komentar tam nekdeje]. n je pocet radku a sloupcu ctvercove matice A. nrhs=1 - pocet pravych stran ,lda=ldb=n ipiv jsou koeficienty z LU rozkladu, info je navratova hodnota 0= spravne ukonceni
    25.12.2010 22:58 luky
    Rozbalit Rozbalit vše Re: Ukazatele v C, pole doublu
    Po letmom pohlade na oba programy som si vsimol jeden rozdiel vo funkciach LVP_solveGE() (zo suboru LVP_solver1.c) a funkciou LVP_solveG() (zo suboru lapack2.c). V jednom su lokalne premenne typu long int a v druhom int. Mozno to bude v tom.
    25.12.2010 23:35 Preston Pigglesworth
    Rozbalit Rozbalit vše Re: Ukazatele v C, pole doublu
    Příloha:
    Dekuji za reakce, omlouvam se, ze sem vas zmatl tema ukazkama, zkousel jsem, jestli nemuze byt problem v tom, ze fortranovsky integer odpovida Cckovskemu long int.. V tomto problem neni, je to osetreno pres lapacke.h

    Funkce dgesv( &n, &nrhs, a, &lda, ipiv, b, &ldb, &info ) da spravny vysledek pokazde, kdyz jsou pole a a b zavedeny jako double a[121]={...}; double b[11]={...};

    Kdyz se ale snazim ty pole zavest dynamicky, tak to pokazde vypocita uplne nesmysly

    Zjednodusil jsem ukazku,jak to jen slo.. kdyz nahradim volani LVP_solveGE(n,S,B); za LVP_solveGE(n,as,bs); tak se dostavi spravny vysledek.

    P.S.Omlouvam se, ze to vypada tak neprehledne, ale to je proste vypocet..Resp. to je neco, co "proste musi byt" v C, i kdyz je to v C zbytecne, spis az nevhodne

    26.12.2010 00:52 Dr. Eddy | skóre: 9 | blog: glog | České Budějovice
    Rozbalit Rozbalit vše Re: Ukazatele v C, pole doublu
    Funkce dgesv( &n, &nrhs, a, &lda, ipiv, b, &ldb, &info ) da spravny vysledek pokazde, kdyz jsou pole a a b zavedeny jako double a[121]={...}; double b[11]={...}; Kdyz se ale snazim ty pole zavest dynamicky, tak to pokazde vypocita uplne nesmysly
    To bude ten problem, protoze je rozdil mezi dynamicky a staticky alokovanym polem. Funkce totiz ocekava staticke pole, ktere by melo byt alokovane jako souvisly blok v pameti, zatimco dynamicky alokovane pole velikosti [n][m] je n poli velikosti m ruzne po pameti. Takze by asi bylo vhnodne ty matice definovat jako jednorozmerne pole a program prepsat.
    26.12.2010 01:05 Dr. Eddy | skóre: 9 | blog: glog | České Budějovice
    Rozbalit Rozbalit vše Re: Ukazatele v C, pole doublu
    jo, pardon, blbe koukam a pak blbe kecam... jsem uz unaveny...
    26.12.2010 10:54 Michal Kubeček | skóre: 71 | Luštěnice
    Rozbalit Rozbalit vše Re: Ukazatele v C, pole doublu

    Dost by pomohlo, kdyby ta ukázka šla přeložit bez spousty souborů, které vy máte, ale my se o jejich obsahu můžeme jen dohadovat. Pro začátek např. zkuste alokovat dynamicky pole se stejným obsahem jako mají as a bs a zavolat funkci na ně. Pokud bude výsledek stejný jako u staticky alokovaného (což je více než pravděpodobné), problém je jinde, než kde se ho snažíte hledat.

    Ještě mne napadá: nikde nekontrolujete návratové hodnoty malloc(), nemůže být problém v tom, že selže alokace? Jak velká je ta vaše matice?

    26.12.2010 11:46 Preston Pigglesworth
    Rozbalit Rozbalit vše Re: Ukazatele v C, pole doublu
    Příloha:
    Upravil jsem solver1.c tak aby ho bylo mozne spustit bez tech veci okolo, kompiluju to pomoci gcc LVP_solver1.c lapack_LINUX.a lapacke.a blas_LINUX.a tmglib_LINUX.a libgfortran.so.3.0.0 Dal bych jsem i ty knihovny, ale maji dohromady 30MB.

    Pro začátek např. zkuste alokovat dynamicky pole se stejným obsahem jako mají as a bs a zavolat funkci na ně.

    To ted jsem prave vyzkousel pomoci tech statickych poli as, bs. prekopiroval jsem jejich obsah do B a S (pred tim, nez sem volal _dgesv) pomoci for (i=0; i<n; i++) { *(B + i) =bs[i]; } for (i=0; i<n*n; i++) { *(S + i) =as[i]; } Vysledek je kupodivu v poradku!

    O to vic tomu ted ale nerozumim, protoze ty cisla v tech statickych polich as, bs jsem ziskal primo z tohoto programu tim, ze jsem to pole ihned po naplneni for cykly zapsal to textoveho souboru

    26.12.2010 13:57 Preston Pigglesworth
    Rozbalit Rozbalit vše Re: Ukazatele v C, pole doublu
    Vypada to, ze problem je nekde v tom plneni matic ( ////FILL MATRIX // ) : kdyz jsem presunul

    for (i=0; i<n; i++) { *(B + i) =bs[i]; } for (i=0; i<n*n; i++) { *(S + i) =as[i]; } Pred vsechno to naplnovani, tak to zase vyhodilo spatny vysledek..

    jeste jsem za kazde to plneni matice napsal jeji vypis pomoci for (i=0; i<n; i++) {printf("%20.20f \n",*(B+i) ); } ... vsechno je v poradku,vstupy do _dgesv jsou stejne, jako v pripade pouziti statickeho pole
    26.12.2010 23:00 Michal Kubeček | skóre: 71 | Luštěnice
    Rozbalit Rozbalit vše Re: Ukazatele v C, pole doublu

    Hodnoty ve statických polích nejsou úplně stejné jako ty počítané, vypadá to, že jsou zaokrouhlené na šest desetinných míst. A pokud místo zkopírování statického pole dám

      for (i=0; i<n; i++) {
        B[i] = 1E-6 * lround(1E6 * B[i]);
      }
      for (i=0; i<n; i++) {
        for (j=0; j<n; j++) {
          S[i+j*n] = 1E-6 * lround(1E6 * S[i+j*n]);
        }
      }
    

    dostanu stejné výsledky jako se statickým polem. Takže můj tip je, že celý problém je v (ne)stabilitě té soustavy, tj. že i malá změna v koeficientech nebo pravé straně může způsobit (relativně) velkou změnu řešení. Ale nechce se mi počítat vlastní čísla, abych si to ověřil.

    27.12.2010 12:29 Preston Pigglesworth
    Rozbalit Rozbalit vše Re: Ukazatele v C, pole doublu
    Dekuji vsem za pomoc! S tim zaokrouhlovanim mate pravdu, _dgesv pocita spravne. Nakonec jsem nasel chybu ve spatnem umisteni jedne zavorky pri plneni matice S, vzniklo to behem kopirovani z matlabovskeho prototypu. Protoze me ale v minulosti jeden vypocet pocital nesmysly diky ukazatelum, podezrival jsem spatno spravu pameti, za to se omlouvam.. Kazdopadne jeste jednou: diky za pomoc!
    26.12.2010 09:46 Preston Pigglesworth
    Rozbalit Rozbalit vše Re: Ukazatele v C, pole doublu
    Jeste jsem zkousel pouzit jinou knihovnu (intel MKL), vysledek stale stejne spatny, takze knihovnou to neni.

    Založit nové vláknoNahoru

    Tiskni Sdílej: Linkuj Jaggni to Vybrali.sme.sk Google Del.icio.us Facebook

    ISSN 1214-1267   www.czech-server.cz
    © 1999-2015 Nitemedia s. r. o. Všechna práva vyhrazena.