[Bug-gsl] bug in cblas on aix

[Kyle L. Jensen]

- GSL version 1.4 from http://sources.redhat.com/gsl/

- uname -a = "AIX brahmastra 2 5 0026022A4C00 unknown"
  The machine is running aix 5.2L and it's a p670
  8way with power4 processors.

- The error occurs with gcc (version 2.9-aix51-020209)
  and xlc/cc (Version 6).  In both cases, no compiler
  flags except 'CPPFLAGS="-I/usr/local/include"
  LDFLAGS="-L/usr/local/lib"' which were passed to 
  ./configure.

- I find this error only on my aix box, not my linux
  box.  I am completely mystified.  Googling around,
  I found a similar bug that it appears was fixed (at
  least when it builds on a linux box??).  The example
  bug code for that barfs on my machine.  Here's how
  I compiled gsl-1.4 (I used xlc here, but the same
  error occurs with gcc)

  CC="xlc" CPPFLAGS="-I/usr/local/include" LDFLAGS="-L/usr/local/lib"
  CFLAGS="-g" ./configure
  make
  make install

- a piece of C is attached.  here's a transcript of compilation
  and the gdb backtrace (as a side note, compiling the blas.c
  file in doc/examples also works just fine but seg faults on
  execution...both with gcc and xlc)

address@hidden ~/GSL]$ xlc -L /usr/local/lib -I /usr/local/include -lgsl -g 
svd-test.c 
address@hidden ~/GSL]$ !gd
address@hidden ~/GSL]$ gdb ./a.out 
GNU gdb 5.0-aix51-020209
Copyright 2000 Free Software Foundation, Inc.
GDB is free software, covered by the GNU General Public License, and you are
welcome to change it and/or distribute copies of it under certain conditions.
Type "show copying" to see the conditions.
There is absolutely no warranty for GDB.  Type "show warranty" for details.
This GDB was configured as "powerpc-ibm-aix5.1.0.0"...
(gdb) run
Starting program: /home/kljensen/GSL/./a.out 
[tcsetpgrp failed in terminal_inferior: Not owner]
gets to here

Program received signal SIGSEGV, Segmentation fault.
0xd2d09450 in cblas_ddot ()
(gdb) bt
#0  0xd2d09450 in cblas_ddot ()
#1  0xd2d08154 in gsl_blas_ddot (X=0x2ff22708, Y=0x2ff22720, result=0x2ff226d8) 
at blas.c:100
#2  0xd2dda4b0 in gsl_linalg_SV_decomp_jacobi (A=0x200009b8, Q=0x20000a38, 
S=0x20000ab8) at svd.c:465
#3  0x100004c0 in main () at svd-test.c:28
#4  0x100001dc in __start ()
(gdb) quit
A debugging session is active.
Do you still want to close the debugger?(y or n) y

svd-test.c
Description: Text document