Re: octave with enable-64

[Jaroslav Hajek]

On Fri, Jun 5, 2009 at 12:04 PM, Alexander
Barth<address@hidden> wrote:
> On Wed, Jun 3, 2009 at 4:08 PM, Alexander
> Barth<address@hidden> wrote:
>> Hi all,
>> I would like to compile octave with --enable-64 to create matrices
>> larger than 2GB.
>>
>> Does somebody has experience in compiling BLAS, LAPACK and SuiteSparse
>> so that this option can be used?
>>
>> So far I tried to compile the reference BLAS and LAPACK with the
>> gfortran option -fdefault-integer-8 and SuiteSparse with the option
>> -DLP64 (I got this from UMFPACK's User Guide). However the example
>> program in UMFPACK (SuiteSparse/UMFPACK/Demo/umfpack_simple) failed at
>> the call of dgemv (a BLAS function).
>>
>> I have the same error with gotoBLAS 1.26 (using BINARY64=1 INTERFACE64=1).
>>
>> Any insight would be greatly appreciated!
>>
>> Thanks
>> Alex
>>
>> My system and software versions:
>> gcc/gfortran 4.2.4
>> octave 3.0.5
>> metis 4.0.1
>> SuiteSparse 3.4.0
>> lapack 3.1.1
>> Ubuntu 8.04 / Xeon 64-bit / 16 GB RAM
>>
>
> It seems that compiling suitesparse with the flags " -DLP64
> -D'LONGBLAS=long int' -D'LONG=long int' " works. It compiles correctly
> the demo program in UMFPACK and the building of octave works too.
>
> In octave, splu does also works:
>
>>> A = sprandn(2000,2000,.1);
>>> [L,U] = splu(A);
>>> d = L*U - A; max(abs(d(:)))
> ans = Compressed Column Sparse (rows = 1, cols = 1, nnz = 1)
>
>  (1, 1) ->  6.8834e-14
>
>
> Cheers,
> Alex

This is nice. Since you're likely a pioneer in this regard, maybe
you'd want to write a short info on the Octave wiki?



-- 
RNDr. Jaroslav Hajek
computing expert & GNU Octave developer
Aeronautical Research and Test Institute (VZLU)
Prague, Czech Republic
url: www.highegg.matfyz.cz