Re: [Mplapack-devel] Bug in Rgerfs

[Maho N. <ch...@ma...>]

Hi Fletcher,

committed, thanks!

Best regards,
 Nakata Maho

From: "Fletcher, John P" <j.p...@as...>
Subject: Bug in Rgerfs
Date: Wed, 06 Apr 2011 09:18:41 +0100

> Nakata Maho
> 
> I have several things to send you and will start with this one.
> 
> I have been testing Rgerfs for double and dd and did not get the results I get with LAPACK.
> 
> I think the reason is that some of the accesses in Rgerfs to ipiv, B and X start at [1] and not [0], and also access to work is wrongly indexed so that one more variable is accessed at the end.
> 
> I attach two test routines.  solve_double_refine.cpp shows it not working and solve_dd_refine.cpp shows reindexing the ipiv, B and X arrays - I copied them to new arrays - and extra work.  This works except that the error results are wrong - file attached.
> 
> The other file is the correct results for the same problem using the routine from LAPACK.
> 
> I hope this helps.
> 
> John Fletcher
> 
> P.S.  I will send other things later.
> 
> Dr John P. Fletcher Tel: (44) 121 204 3389 (direct line), FAX: (44) 121 204 3678
> Chemical Engineering and Applied Chemistry (CEAC),
> Associate Dean - External Relations,
> School of Engineering and Applied Science (EAS),
> Aston University, Aston Triangle, BIRMINGHAM B4 7ET  U.K.
> 
>