On Thursday 15 May 2008, William Stein wrote: > > Btw. I don't have access to Magma 2.14 which I believe to be the fastest > > in linear algebra over GF(2). In version 2.14 they added SSE2 support > > too. So if anybody could compare the new M4RI with Magma 2.14 I'd be > > happy to hear about the results. I don't know what else to compare with > > except Magma. > > Here are magma V2.14-9 timings on sage.math: [EMAIL PROTECTED]:~/magma64/bin$ ./magma Magma V2.14-9 Wed May 14 2008 17:33:23 on sage [Seed = 793568142] Type ? for help. Type <Ctrl>-D to quit. > A := RandomMatrix(GF(2),10^4,10^4); > B := RandomMatrix(GF(2),10^4,10^4); > time C := A*B; Time: 4.020
> time C := A*B; Time: 3.980 > time C := A*B; Time: 3.990 > A := RandomMatrix(GF(2),2^14,2^14); > B := RandomMatrix(GF(2),2^14,2^14); > time C := A*B; Time: 15.450 > time C := A*B; Time: 15.420 This doesn't seem much different from 2.13, which is good for my purposes. Once I fixed the OSX issue mentioned below could you run the benchmark on the C2D? I want to see if they can make better use of the 4MB L2 cache. > ---- > > I tried your new code up at #3204 under OS X and get this: > > sage: A = random_matrix(GF(2),10^4,10^4) > sage: B = random_matrix(GF(2),10^4,10^4) > sage: time C = A._multiply_strassen(B,cutoff=3200) > sage.bin(39971) malloc: *** error for object 0xb95c010: Non-aligned > pointer being freed (2) > *** set a breakpoint in malloc_error_break to debug > sage.bin(39971) malloc: *** error for object 0x79c9c10: Non-aligned > pointer being freed (2) > *** set a breakpoint in malloc_error_break to debug > sage.bin(39971) malloc: *** error for object 0x7465a00: > non-page-aligned, non-allocated pointer being freed > *** set a breakpoint in malloc_error_break to debug > sage.bin(39971) malloc: *** error for object 0x79ca610: Non-aligned > pointer being freed (2) > *** set a breakpoint in malloc_error_break to debug > ... > CPU times: user 10.29 s, sys: 0.26 s, total: 10.55 s > Wall time: 16.31 > > Maybe you're doing something wrong? Likely, though my Valgrind log is clean under Linux. I'll try under OSX. Martin -- name: Martin Albrecht _pgp: http://pgp.mit.edu:11371/pks/lookup?op=get&search=0x8EF0DC99 _www: http://www.informatik.uni-bremen.de/~malb _jab: [EMAIL PROTECTED] --~--~---------~--~----~------------~-------~--~----~ To post to this group, send email to sage-devel@googlegroups.com To unsubscribe from this group, send email to [EMAIL PROTECTED] For more options, visit this group at http://groups.google.com/group/sage-devel URLs: http://www.sagemath.org -~----------~----~----~----~------~----~------~--~---
