[Getfem-users] Re: bug in getfem

[Renard Yves]

Selon Eugene Kashdan <address@hidden>:

> Hi,
> 
> There is a bug in the export of  matrices into to the HB format:
> 
> The header file line 3:
> Col. 57-70 -- number of elemental matrix entrances (zero in the case of
> assembled matrices).
> 
> This bug leads to the problems with SuperLU (the program starts to read the
> file, but never ends).

Ok, I added a zero there. I changed the line of the function writeHB_mat_double
in gmm_inoutput.h : 

fprintf(out_file,"%3s%11s%14d%14d%14d\n",Type,"          ", M, N, nz);

to

fprintf(out_file,"%3s%11s%14d%14d%14d%14d\n",Type,"          ", M, N, nz, 0);


I commited it in the SVN tree at gna.org.

Thank you for this bug report.

Yves.