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 14:22 | Komunita

    Dle plánu byl vývoj Firefoxu přesunut z Mercurialu na Git. Oficiální repozitář se zdrojovými kódy je na GitHubu.

    Ladislav Hagara | Komentářů: 2
    dnes 04:33 | Bezpečnostní upozornění

    V terminálovém multiplexoru GNU Screen byly nalezeny a v upstreamu ve verzi 5.0.1 už opraveny bezpečnostních chyby CVE-2025-23395, CVE-2025-46802, CVE-2025-46803, CVE-2025-46804 a CVE-2025-46805. Podrobnosti na blogu SUSE Security Teamu.

    Ladislav Hagara | Komentářů: 31
    včera 19:33 | Bezpečnostní upozornění

    Training Solo (Paper, GitHub) je nejnovější bezpečnostní problém procesorů Intel s eIBRS a některých procesorů ARM. Intel vydal opravnou verzi 20250512 mikrokódů pro své procesory.

    Ladislav Hagara | Komentářů: 0
    včera 11:44 | Nová verze

    Byla vydána nová verze 25.05.11 svobodného multiplatformního video editoru Shotcut (Wikipedie) postaveného nad multimediálním frameworkem MLT. Nejnovější Shotcut je již vedle zdrojových kódů k dispozici také ve formátech AppImage, Flatpak a Snap.

    Ladislav Hagara | Komentářů: 0
    včera 11:11 | Nová verze

    Svobodný elektronický platební systém GNU Taler (Wikipedie, cgit) byl vydán ve verzi 1.0. GNU Taler chrání soukromí plátců a zároveň zajišťuje, aby byl příjem viditelný pro úřady. S vydáním verze 1.0 byl systém spuštěn ve Švýcarsku.

    Ladislav Hagara | Komentářů: 10
    včera 00:55 | Pozvánky

    Spolek OpenAlt zve příznivce otevřených řešení a přístupu na 209. brněnský sraz, který proběhne tento pátek 16. května od 18:00 ve studentském klubu U Kachničky na Fakultě informačních technologií Vysokého učení technického na adrese Božetěchova 2/1. Jelikož se Brno stalo jedním z hlavních míst, kde se vyvíjí open source knihovna OpenSSL, tentokrát se OpenAlt komunita potká s komunitou OpenSSL. V rámci srazu Anton Arapov z OpenSSL

    … více »
    Ladislav Hagara | Komentářů: 0
    včera 00:22 | Komunita

    GNOME Foundation má nového výkonného ředitele. Po deseti měsících skončil dočasný výkonný ředitel Richard Littauer. Vedení nadace převzal Steven Deobald.

    Ladislav Hagara | Komentářů: 8
    10.5. 15:00 | Zajímavý článek

    Byl publikován přehled vývoje renderovacího jádra webového prohlížeče Servo (Wikipedie) za uplynulé dva měsíce. Servo zvládne už i Gmail. Zakázány jsou příspěvky generované pomocí AI.

    Ladislav Hagara | Komentářů: 29
    9.5. 17:22 | Nová verze

    Raspberry Pi Connect, tj. oficiální služba Raspberry Pi pro vzdálený přístup k jednodeskovým počítačům Raspberry Pi z webového prohlížeče, byla vydána v nové verzi 2.5. Nejedná se už o beta verzi.

    Ladislav Hagara | Komentářů: 6
    9.5. 15:22 | Komunita

    Google zveřejnil seznam 1272 projektů (vývojářů) od 185 organizací přijatých do letošního, již jednadvacátého, Google Summer of Code. Plánovaným vylepšením v grafických a multimediálních aplikacích se věnuje článek na Libre Arts.

    Ladislav Hagara | Komentářů: 0
    Jaký filesystém primárně používáte?
     (57%)
     (1%)
     (8%)
     (22%)
     (4%)
     (2%)
     (2%)
     (1%)
     (0%)
     (3%)
    Celkem 601 hlasů
     Komentářů: 26, poslední 8.5. 09:58
    Rozcestník

    Dotaz: Ukazatele v C, pole doublu

    25.12.2010 21:56 Preston Pigglesworth
    Ukazatele v C, pole doublu
    Přečteno: 383×
    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: 72 | 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: 72 | 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: 72 | 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.