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:11 | IT novinky

    Řada vestavěných počítačových desek a vývojových platforem NVIDIA Jetson se rozrostla o NVIDIA Jetson Thor. Ve srovnání se svým předchůdcem NVIDIA Jetson Orin nabízí 7,5krát vyšší výpočetní výkon umělé inteligence a 3,5krát vyšší energetickou účinnost. Softwarový stack NVIDIA JetPack 7 je založen na Ubuntu 24.04 LTS.

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

    Národní úřad pro kybernetickou a informační bezpečnost (NÚKIB) spolu s NSA a dalšími americkými úřady upozorňuje (en) na čínského aktéra Salt Typhoon, který kompromituje sítě po celém světě.

    Ladislav Hagara | Komentářů: 17
    včera 16:33 | IT novinky

    Společnost Framework Computer představila (YouTube) nový výkonnější Framework Laptop 16. Rozhodnou se lze například pro procesor Ryzen AI 9 HX 370 a grafickou kartu NVIDIA GeForce RTX 5070.

    Ladislav Hagara | Komentářů: 1
    včera 14:22 | IT novinky

    Google oznamuje, že na „certifikovaných“ zařízeních s Androidem omezí instalaci aplikací (včetně „sideloadingu“) tak, že bude vyžadovat, aby aplikace byly podepsány centrálně registrovanými vývojáři s ověřenou identitou. Tato politika bude implementována během roku 2026 ve vybraných zemích (jihovýchodní Asie, Brazílie) a od roku 2027 celosvětově.

    Fluttershy, yay! | Komentářů: 7
    včera 13:11 | Nová verze

    Byla vydána nová verze 21.1.0, tj. první stabilní verze z nové řady 21.1.x, překladačové infrastruktury LLVM (Wikipedie). Přehled novinek v poznámkách k vydání: LLVM, Clang, LLD, Extra Clang Tools a Libc++.

    Ladislav Hagara | Komentářů: 0
    včera 05:11 | Komunita

    Alyssa Anne Rosenzweig v příspěvku na svém blogu oznámila, že opustila Asahi Linux a nastoupila do Intelu. Místo Apple M1 a M2 se bude věnovat architektuře Intel Xe-HPG.

    Ladislav Hagara | Komentářů: 17
    26.8. 12:55 | IT novinky

    EU chce (pořád) skenovat soukromé zprávy a fotografie. Návrh "Chat Control" by nařídil skenování všech soukromých digitálních komunikací, včetně šifrovaných zpráv a fotografií.

    Ladislav Hagara | Komentářů: 59
    26.8. 12:11 | Nová verze

    Byly publikovány fotografie a všechny videozáznamy z Python konference PyCon US 2025 proběhlé v květnu.

    Ladislav Hagara | Komentářů: 0
    26.8. 11:55 | IT novinky

    Společnost xAI a sociální síť X amerického miliardáře Elona Muska zažalovaly firmy Apple a OpenAI. Viní je z nezákonné konspirace s cílem potlačit konkurenci v oblasti umělé inteligence (AI).

    Ladislav Hagara | Komentářů: 3
    26.8. 05:44 | Nová verze

    Byla vydána nová verze 9.16 z Debianu vycházející linuxové distribuce DietPi pro (nejenom) jednodeskové počítače. Přehled novinek v poznámkách k vydání.

    Ladislav Hagara | Komentářů: 0
    Pro otevření více webových stránek ve webovém prohlížečí používám
     (78%)
     (9%)
     (3%)
     (4%)
     (5%)
     (1%)
    Celkem 101 hlasů
     Komentářů: 9, poslední dnes 11:53
    Rozcestník

    Interface mezi C++ a F77

    6.3.2008 10:35 | Přečteno: 1175× | Výběrový blog | poslední úprava: 10.6.2008 17:33

    Jednoduchy příklad interface mezi C++ a f77

    Minulý víkend jsem za po dlouhé době potřeboval udělat interface mezi kódem ve fortranu a kódem v C++. Jelikož jsem to již dlouho nedělal, chvíli mi trvalo, než jsem si vzpomněl, jak to udělat (neměl jsem přístup na net; při této příležitosti jsem si zas uvědomil, že bez netu jsem jako bez ruky :-( - tedy co se programování týče). Poté, co se mi to zadařilo, jsem se rozhodnul udělat si jednoduchý prográmek pro přístě, až zas budu bez netu a budu vzpomínat, jak se to dělá. Vše níže uvedené platí pro kompilátor gcc.

    Volání fortranovské rutiny z C++ je vcelku jednoduché, stačí funkci deklarovat jako externí, název je stejný jako ve fortranu, ale přídáme podtržítko. Případný parameter předáme pomocí jeho adresy:

    extern "C" {void vypisparam_(int* cislo);}
    ...
    vypisparam_(&mojeCislo);
    

    Obvykle ale potřebujeme víc. Velmi často je potřeba nějak přistupovat k proměnným ve fortranovských common blocích. Mějme následně definovaný common blok a jednoduchý fotranovský prográmek:
    test.inc:

           DOUBLE PRECISION dp1,dp2
           INTEGER int1
           COMMON/MYTEST/dp1,dp2,int1
    
    test.f:
          PROGRAM TEST
          INCLUDE 'test.inc'
          CALL NASTAV
          CALL VYPIS
          CALL VYPISPARAM(20)
          CALL VYPIS
          END
    
          SUBROUTINE NASTAV
          INCLUDE 'test.inc'
          PRINT*,'Rutina Nastav - nastavuji parametry'
          dp1 = 1.D0
          dp2 = 3.5D0
          int1 = 10
          END
    
          SUBROUTINE VYPIS
          INCLUDE 'test.inc'
          PRINT*,"Rutina Vypis:"
          PRINT*,'dp1 = ',dp1
          PRINT*,'dp2 = ',dp2
          PRINT*,'int1 = ',int1
          END
    
          SUBROUTINE VYPISPARAM(a)
          INCLUDE 'test.inc'
          INTEGER a
          dp1 = a
          dp2 = a
          int1 = a
          PRINT*,'Rutina VypisParam:'
          PRINT*,'a=',a
          END
    
    Common blok bude v C++ reprezentovat struktura. Proměnné v common bloku jsou uloženy v paměti za sebou. Inicializace tedy probíhá tak, že se ukazateli na strukturu v C++ předá adresa první proměnné common bloku z fortranu. Je proto nutné, aby proměnné ve struktuře byly stejného typu a ve stejném pořadí jako ve fortranu (pro vícerozměrná pole pak ješte je nutno dát pozor indexy). O vrácení adresy common bloku se nám postará jednoduchá fortranovská funkce:
          FUNCTION common_block_address(common_block_name)
          INCLUDE 'test.inc'
          CHARACTER*(*) common_block_name
          INTEGER common_block_address
          INTEGER aaadress
          IF(common_block_name.EQ.'MYTEST')THEN
             common_block_address = adress(dp1)
          ELSE
             PRINT*,'Neexistujici common block'
          ENDIF
          RETURN
          END
    
    Jednoduchý prográmek v C++, který volá rutiny fortranovského kódu a pracuje s proměnnými v common bloku, by mohl vypadat takto (v C++ by tedy bylo vhodnejší, abychom byli objektoví:-), vytvořit objekt, který by reprezentoval fortanovský kód):
    #include <iostream>
    
    double loc_dp1,loc_dp2;
    int loc_int1;
    
    
    struct MyTest_t{
      double dp1,dp2;
      int int1;
    };
    
    extern "C" {
      void nastav_();
      void vypis_();
      void vypisparam_(int* cislo);
    }
    
    extern "C" void* adress_(void* var){
      return var;
    }
    
    extern "C" void* common_block_address_(char*,int len);
    
    MyTest_t* fMyTest;
    void InitCommonBlock(){
      fMyTest = (MyTest_t*)common_block_address_("MYTEST",6);
    }
    
    void LocalParToFor(){
        fMyTest->dp1 = loc_dp1;
        fMyTest->dp2 = loc_dp2;
        fMyTest->int1 = loc_int1;
    }
    
    void ForToLocalPar(){
        loc_dp1 = fMyTest->dp1;
        loc_dp2 = fMyTest->dp2;
        loc_int1 = fMyTest->int1;
    }
    
    int main(){
      int aaa = 100;
      InitCommonBlock();
      nastav_();
      vypis_();
      ForToLocalPar();
      std::cout << "loc_dp1 je " << loc_dp1 << std::endl;
      std::cout << "loc_dp2 je " << loc_dp2 << std::endl;
      std::cout << "loc_int1 je " << loc_int1 << std::endl;
      vypisparam_(&aaa);
      loc_dp1 = 4.5;
      loc_dp2 = 8.3;
      loc_int1 = 15;
      LocalParToFor();
      vypis_();
      return 0;
    }
    
    Vše kompilujeme pomocí g++ a poté slikujeme dohromady společně s knihovnou gfortan, u mě např. takto:
    g++ common_block_address.o interface.o test.o -L/usr/lib/gcc/i386-redhat-linux/4.3.0/ -lgfortran        

    Hodnocení: 100 %

            špatnédobré        

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

    Komentáře

    Vložit další komentář

    6.3.2008 11:14 anicka | blog: ze_zivota
    Rozbalit Rozbalit vše Re: Interface mezi C++ a F77
    Hezky clanek a pekne osklive dejavu. Staci nad clankem jen trochu privrit oci a predstavim si zdrojaky octave... :-)
    ^D
    6.3.2008 13:55 Let_Me_Be | skóre: 20 | blog: cat /proc/idea/current | Brno
    Rozbalit Rozbalit vše Re: Interface mezi C++ a F77
    Sice nechapu co to dela ve vyberu, ale jinak dobry zapisek.
    Linked in profil - Můj web - Nemůžete vyhrát hádku s blbcem. Nejdřív vás stáhne na svoji úroveň a pak ubije zkušenostmi.
    ISSN 1214-1267   www.czech-server.cz
    © 1999-2015 Nitemedia s. r. o. Všechna práva vyhrazena.