bugGNU Scientific Library - Bugs: bug #24162, solve_cyc_tridiag_nonsym does not...

 
 

You are not allowed to post comments on this tracker with your current authentification level.

bug #24162: solve_cyc_tridiag_nonsym does not respect solution vector's stride

Submitted by:  Paul Accisano <shiningmasamune>
Submitted on:  Fri 29 Aug 2008 05:30:10 PM UTC  
 
Category: Runtime errorSeverity: 4 - Important
Operating System: Status: Fixed
Assigned to: -Deleted Account- <bjg>Open/Closed: Closed
Release: 1.11

Sat 30 Aug 2008 01:26:10 PM UTC, comment #1:

Thanks for the bug report, I've committed your fix to the repository and extended the tests to check non-unit strides.

http://git.savannah.gnu.org/gitweb/?p=gsl.git;a=commit;h=456f944ac4fcf8b730a28b7cac13f7a84f29a9f1

-Deleted Account- <bjg>In charge of this item.
Fri 29 Aug 2008 05:30:10 PM UTC, original submission:

The function solve_cyc_tridiag_nonsym in linalg_tridiag.c does not respect the stride for the x (solution) vector. Line 403 reads:

x[N-1] = zb[N-1]/alpha[N-1];

when it should be

x[x_stride * (N-1)] = zb[N-1]/alpha[N-1];

In addition, Line 423 reads:

x[i] -= vx/(1 + vw)*w[i];

when it should be

x[x_stride * i] -= vx/(1 + vw)*w[i];

Those are the only two instances I could find. After fixing them, the function appears to operate correctly.

Paul Accisano <shiningmasamune>

 

No files currently attached

 

Depends on the following items: None found

Items that depend on this one: None found

 

Carbon-Copy List
  • -unavailable- added by shiningmasamune (Submitted the item)
  •  

    Please enter the title of George Orwell's famous dystopian book (it's a date):

     

     

    Follow 4 latest changes.

    Date Changed By Updated Field Previous Value => Replaced By
    Sat 30 Aug 2008 01:26:10 PM UTCbjgStatusNone=>Fixed
      Assigned toNone=>bjg
      Open/ClosedOpen=>Closed
    Sat 30 Aug 2008 01:26:09 PM UTCbjgSeverity3 - Normal=>4 - Important

    Back to the top


    Powered by Savane 3.1-cleanup