Vielleicht hat ja schon einmal jemand von euch mit cublas gearbeitet, und weiß, welchen Anfängerfehler ich mache.

Ich sitze momentan an dem Problem, dass ich mit ziemlich großen Matrizen (1000*1000 Elementen) komplexe Berechnungen durchführen muss. Momentan ist meine Ausführungsgeschwindigkeit bei ca 70h pro Experiment, weswegen ich mich entschlossen habe, es einmal mit cublas zu probieren und meinen alten Code darauf umzustellen.

Leider scheint das alles nicht so einfach zu sein, wie ich es mir vorstelle. Folgender vereinfachter Code:

void updateStatistics(const Array<double>& s)
{
//allokation der Vektoren
//m_dimensions=1000=dimension der Eingabedaten  
cublasAlloc(m_dimensions,sizeof(double),(void**)&m_theta);
cublasAlloc(m_dimensions,sizeof(double),(void**)&m_state);
//auf 0 setzen
cublasDscal(m_dimensions,0,m_theta,1);

//...

//Kopieren der Eingabe in die Grafikkarte
//elemvec liefert das interne Array...
cublasSetVector(m_dimensions,sizeof(double),s.elemvec(),1,m_state,1);

//skalarprodukt für double Precision
double delta=cublasDdot(m_dimensions,m_theta,1,m_nextState,1);
std::cout<<cublasGetError()<<" "<<delta<<" ";//kein Fehler, delta=NaN
}

Wie schon in den Kommentaren steht, wird delta=nan. Die Allokation scheint zu klappen, zumindest haben die Pointer am Ende Werte und es wird auch kein Fehler ausgegeben. Ich bin grad ziemlich ratlos, weil ich glaube alles richtig gemacht zu haben. Und an der Eingabe habe ich seit meinem letzten (korrekten) Experiment nichts verändert.

Weiß jemand, wa sich da falsch gemacht habe? In der Doku steht leider auch nicht, wann NaN auftritt. dass 0 raus kommt ist zu dem zeitpunkt okay. Und die Eingabe besteht nur aus 0 und 1 Werten. Von daher kann ich auch dort NaN ausschließen.

cublasDscal(m_dimensions,0,m_theta,1);

Ist zwar toll in der Theorie, aber wenn die Werte da drin nach allokation des Speichers NaN sind, dann kann ich die mit 0 multiplizieren soviel ich will, die bleiben NaN.

Lass ich mal stehen für alle, die mit cublas mal arbeiten wollen(man ist das Schnell!!!!)