Re: Grauenhafter Fehler in MMA unter Mac OSX?

[Michael Gamer <mgamer@XXXXXXX.com>]

Hallo Martin, 

    ich habe versucht das nachzuvollziehen. Meiner Meinung nach liegt das vielleicht an der Art und Weise wie "Table" 
Ausdrücke berechnet. Hier mal in aller Kürze mein kurzer Vergleich (Auslagern der Berechnung in eine Funktion geht also 
zumindest):


In[52]:= m1000 = IdentityMatrix[1000];
In[53]:= d10 = Table[m1000[[i]].m1000[[All, i]], {i, 1, 10}]
Out[53]= {1, 1, 1, 1, 1, 1, 1, 1, 1, 1}

In[54]:= func[i_] := m1000[[i]].m1000[[All, i]];

In[55]:= d1000 = Table[func[i], {i, 1, 1000}];
In[56]:= dohnefunc = Table[(m1000[[i]].m1000[[All, i]]), {i, 1, 1000}];

In[57]:= d1000[[;; 10]]
Out[57]= {1, 1, 1, 1, 1, 1, 1, 1, 1, 1}

In[58]:= dohnefunc[[;; 10]]
Out[58]= {123887622, 1, 16384, 1999, 33574528, 1, 71726648, 1, 33809760, 1}




(Notebook erzeugt auf einen MacBookPro mit MMA 7.0.1.0)



Viele Grüße

m.g.


Michael Gamer
Zur Römerbrücke 19
63456 Hanau 

Telefon:        06181 663363
Fax:            06181 663363
Mail:   mg@XXXXXXX.de

home:   www.michaelgamer.de
photos:         www.mg-photo.de

Am 06.02.2010 um 15:10 schrieb Martin Heimann:

> 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
> 
> 
>