Gianluca:<div>Checking Superlu source code SuperLU_4.3/SRC/dgssvx.c,</div><div>this input format is required by superlu.</div><div><br></div><div>Hong</div><div><br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hong,<br>
<br>
thanks for your answer. I do not know if my problem request 64 bit support.<br>
<br>
My understanding is that I can preallocate the working array for<br>
superlu with the PETSc database options -mat_superlu_lwork. Following<br>
PETSc manual, this has to be specified in bytes.<br>
<br>
>From my point of view, the problem is that I cannot pass to PETSc an<br>
integer greater than <a href="tel:2.147.483.647" value="+12147483647">2.147.483.647</a> (that is, LONG_MAX), thus I am<br>
limited to allocating roughly 2 GB of memory. If my factored matrix is<br>
bigger, I have to rely on the mallocs/copy.<br>
<br>
Let me know if it is not clear yet.<br>
<br>
Gianlu<br>
<div class="HOEnZb"><div class="h5"><br>
<br>
<br>
<br>
<br>
On 3 July 2012 15:53, Hong Zhang <<a href="mailto:hzhang@mcs.anl.gov">hzhang@mcs.anl.gov</a>> wrote:<br>
> Gianluca:<br>
><br>
> Do you mean 64-bit support for superlu?<br>
> I do not understand your request, thus forward your request to superlu<br>
> developer.<br>
><br>
> Hong<br>
><br>
>> Dear all,<br>
>><br>
>> I am trying to use superlu as solver for a large, sparse matrix, and I<br>
>> would like to use -mat_superlu_lwork to speed up the computation. The<br>
>> problem I encounter is that my lwork size I would need is greated than<br>
>> the maximum value for an object of type long int and is not correctly<br>
>> read by the code.<br>
>><br>
>> Is there a workaround to that? Does using another solver (mumps<br>
>> maybe?) solve this problem and, if so, which option should I use?<br>
>><br>
>> Thanks in advance<br>
>><br>
>> Gianluca<br>
>><br>
>> --<br>
>> "[Je pense que] l'homme est un monde qui vaut des fois les mondes et<br>
>> que les plus ardentes ambitions sont celles qui ont eu l'orgueil de<br>
>> l'Anonymat" -- Non omnibus, sed mihi et tibi<br>
>> Amedeo Modigliani<br>
><br>
><br>
<br>
<br>
<br>
--<br>
"[Je pense que] l'homme est un monde qui vaut des fois les mondes et<br>
que les plus ardentes ambitions sont celles qui ont eu l'orgueil de<br>
l'Anonymat" -- Non omnibus, sed mihi et tibi<br>
Amedeo Modigliani<br>
</div></div></blockquote></div><br></div>