[C++] Behöver hjälp med matriser!

Tja det effektiva sättet att göra det är att endera använda ovan nämnda bibliotek, direkt använda SSE/AVX instrincts (inte att rekommendera) eller utnyttja någon teknik som "Cilk+ array notation".

Vill du göra det "för hand" med standard C++ är artikeln för Cilk+ ändå en bra referens då man visar logisk ekvivalent C++. Fördelen med Cilk+ metoden är dels att den är kompaktare och dels att den gör det möjligt för kompilatorn att lägga ut SSE/AVX instruktioner, något som är extremt svårt med "vanlig" C++.

Andra "roll your own" varianter där du står för logiken är

• Intel® SPMD Program Compiler

• OpenCL (går då att köra både på CPU med SSE/AVX och på GPU)

• samt CUDA (går att köra på Nvidia GPU och x86 CPUer med SSE/AVX via CUDA-x86)

Edit: för att svara på om C-arrayer eller C++ vector<> är bäst, finns ingen anledning att köra med C-arrayer i C++ då vector<> är lika snabbt för själva matrisoperationerna.

Tänk bara på att den inre loopen måste iterera över den inre vector<> delen. Detta så att du gör minnesaccesserna i följd i minnet, något som ökar den spatiala lokaliteten på data-accesserna (högre CPU-cache hit rate -> snabbare).

Edit 2: Intel må ha köpt CilkArts som gör Cilk+, men tekniken finns numera implementerad i gcc (och är endera på gång för clang eller redan på plats). Så går att använda gratis.