[Scipy-tickets] [SciPy] #709: ValueError: array dimensions are not compatible for copy

SciPy scipy-tickets@scipy....
Sat Dec 6 09:08:49 CST 2008


#709: ValueError: array dimensions are not compatible for copy
--------------------------+-------------------------------------------------
 Reporter:  nils          |        Owner:  somebody
     Type:  defect        |       Status:  new     
 Priority:  high          |    Milestone:  0.7.0   
Component:  scipy.linalg  |      Version:  devel   
 Severity:  blocker       |   Resolution:          
 Keywords:                |  
--------------------------+-------------------------------------------------
Comment (by pv):

 In the end, I forgot to send the mail :/
 Sent now, and CC'd to scipy-dev.

 Anyway, this is starting to look like a subtle compiler/optimization bug:
 {{{
 # lapack 3.2, with only -g flag in make.inc OPTS
 $ gfortran --version
 4.3.2-1ubuntu11
 $ make clean
 ...
 $ make lib
 ...
 $ gfortran -o 2dof 2dof.f lapack_LINUX.a -lblas-3
 $ ./2dof
 Checking for missing eigenvalue. If you see no
 further output, you don't have this problem.
 }}}
 No bug there. But,
 {{{
 # lapack 3.2, with -O2 flag in make.inc OPTS
 $ make clean
 ...
 $ make lib
 $ gfortran -o 2dof 2dof.f lapack_LINUX.a -lblas-3
 $ ./2dof
  Checking for missing eigenvalue. If you see no
  further output, you don't have this problem.
  Missing (?) complex pair for omega =  0.85858585858585856
            1   0.0000000000000000        0.0000000000000000
 0.0000000000000000
            2   0.0000000000000000      -1.56357833001198695E-016
 4.05220561966457282E-017
            3   0.0000000000000000       5.05904065833670526E-016
 2.36247653384591416E-016
            4   0.0000000000000000      -5.12056639374098328E-016
 2.39120789141678003E-016
  Missing (?) complex pair for omega =   2.6767676767676769
            1   0.0000000000000000        1.0350059973275267
 0.18232312052566751
            2   0.0000000000000000       -1.0350059973275267
 0.18232312052566751
            3   0.0000000000000000        0.0000000000000000
 0.0000000000000000
            4   0.0000000000000000      -1.21185328887430225E-016
 3.74917111245487433E-016
 }}}
 Lapack 3.1.1 has -O3 as the default, and also this produces an error
 (whereas -g doesn't). Ubuntu compiles by default with -O3.

 On another machine, 64-bit instead of 32-bit, with gfortran 4.3.1, I
 always get an error independent of optimization flags.

-- 
Ticket URL: <http://scipy.org/scipy/scipy/ticket/709#comment:9>
SciPy <http://www.scipy.org/>
SciPy is open-source software for mathematics, science, and engineering.


More information about the Scipy-tickets mailing list