// dot product
template <class T>
T Vector<T>::operator*(const Vector& v) const
{
T ret_val=0;
for(std::size_t i=0; i < m_Dim; i++) ret_val += m_data[i] * v.m_data[i]; // 86.61%
return ret_val;
}

ich denke mal, da irgendwas zu optimieren ist hoffnungslos, oder?
ich könnte per hand die schleife aufrollen,aber da pfusche ich imho dem
compiler dazwischen

m_data ist ein pures c-array vom type T, inline dürfte auch nichts bringen,
denn 86.61% laufen in dieser Zeile ab, und die kosten für die gesamten Aufrufe sind nur unbedeutend höher : 86.77

vorschläge sind willkommen ...

profiling mit kcachegrind ist leider nur im debug-modus möglich, also weiss ich auch nicht direkt, wie der gcc (oder icc, ich benutze beide) sowas optimiert, und wie man da noch was rausholen kann.

aber immerhin schön, das ich nix optimieren muss, wenn nicht diese eine zeile

je nach größe der vektoren könnte man hier mit ner art caching nachhelfen. die vector klasse muss sich dazu möglicherweise merken, ob sie sich geändert hat.

wichtig ist hier allerdings dann die frage, wie gross so ein vektor ist. wenn m_dim klein ist, wird man auch mit caching nicht viel erreichen. dann hilft nur die aufrufe zu reduzieren.

Das sind dann in der Klasse darüber 500-1000 aufrufe dieser Funktion aber über
openMP parallelisiert.

also eigentlich ist das ganze eine multiplikation einer
(500-1000)x(500-1000) Matrix mit einem Vektor, die Matrix ist voll besetzt

Gruß
Don06