Liebe Kollegen,
Beim Berechnen von einzelnen Diagonal-Elementen eines Matrizenprodukts mit Hilfe der Skalarproduktfunktion "." bin ich
auf einen grauenhaften Fehler gestossen. Reduziert auf ein ganz triviales Beispiel zeigt er sich wie folgt:
In[1]:= m=IdentityMatrix[10000];
In[2]:= d10=Table[m[[i]] . m[[All,i]],{i,10}];
In[3]:= d1000=Table[m[[i]] . m[[All,i]],{i,1000}];
In[4]:= d10
Out[4]= {1, 1, 1, 1, 1, 1, 1, 1, 1, 1}
In[5]:= d1000[[;;10]]
Out[5]= {0, 0, 0, 0, 0, 0, 0, 0, 0, 0}
Die Diagonalelemente müssten doch alle 1 geben - wenn ich nur die ersten 10 berechne (d10), dann ergibt sich dies; wenn
ich die ersten 1000 berechne, erhalte ich Schrott. Wie ich leider erfahren musste, tritt das Problem nicht nur bei
dieser trivialen Identitätsmatrix auf.
Der Fehler tritt konsistent in zwei verschiedenen Mma Installationen auf zwei verschiedenen Mac's auf (Version 7.0.1),
auf einer Linux Installation hingegen nicht.
Frage: kann dies jemand von Euch auf einer Mac OSX - Installation auch nachvollziehen?
Gruss,
Martin
----------------------------------------------------------------------------
Max-Planck-Institute for Biogeochemistry, PF 100164, D-07701 Jena, Germany
Street Address: Beutenberg Campus, Hans-Knoell-Straße 10, D-07745 Jena
Office: +49-3641-57-6350/6301
Mobile No: +49-151-12035946
Home: +49-3641-618247
Fax.: +49-3641-57-7300
Skype: mheimann
Email: martin.heimann@XXXXXXX.de,
office.bgc-systems@XXXXXXX.de
Web: http://www.bgc-jena.mpg.de/~martin.heimann