Skip to content

Instantly share code, notes, and snippets.

@certik

certik/gist:830fc3bc015e66c0976249df8c07756b

Created July 25, 2021 09:39
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
Star You must be signed in to star a gist
Embed
What would you like to do?
Learn more about clone URLs
Download ZIP
Raw
gistfile1.txt
$ python
Python 3.9.6 | packaged by conda-forge | (default, Jul 11 2021, 03:35:11)
[Clang 11.1.0 ] on darwin
Type "help", "copyright", "credits" or "license" for more information.
>>> import scipy
>>> scipy.__file__
'/Users/ondrej/repos/scipy/scipy/__init__.py'
>>> scipy.test()
============================= test session starts ==============================
platform darwin -- Python 3.9.6, pytest-6.2.4, py-1.10.0, pluggy-0.13.1
rootdir: /Users/ondrej/repos/scipy, configfile: pytest.ini
plugins: cov-2.12.1, xdist-2.3.0, forked-1.3.0
collected 46069 items / 11135 deselected / 34934 selected
scipy/_build_utils/tests/test_scipy_version.py . [ 0%]
scipy/_lib/tests/test__gcutils.py ...... [ 0%]
scipy/_lib/tests/test__pep440.py ......... [ 0%]
scipy/_lib/tests/test__testutils.py .. [ 0%]
scipy/_lib/tests/test__threadsafety.py .. [ 0%]
scipy/_lib/tests/test__util.py ...................... [ 0%]
scipy/_lib/tests/test_bunch.py ............................ [ 0%]
scipy/_lib/tests/test_ccallback.py .... [ 0%]
scipy/_lib/tests/test_deprecation.py . [ 0%]
scipy/_lib/tests/test_import_cycles.py . [ 0%]
scipy/_lib/tests/test_tmpdirs.py ... [ 0%]
scipy/cluster/tests/test_disjoint_set.py ..................... [ 0%]
scipy/cluster/tests/test_hierarchy.py .................................. [ 0%]
......................................................... [ 0%]
scipy/cluster/tests/test_vq.py ......................... [ 0%]
scipy/constants/tests/test_codata.py ....... [ 0%]
scipy/constants/tests/test_constants.py ... [ 0%]
scipy/fft/_pocketfft/tests/test_basic.py ............................... [ 0%]
........................................................................ [ 0%]
........................................................................ [ 1%]
........................................................................ [ 1%]
........................................................................ [ 1%]
........................................................................ [ 1%]
..................................................................... [ 1%]
scipy/fft/_pocketfft/tests/test_real_transforms.py ..................... [ 2%]
........................................................................ [ 2%]
........................................................................ [ 2%]
........................................................................ [ 2%]
........................................................................ [ 2%]
........................................................................ [ 3%]
........................................................................ [ 3%]
........................................................................ [ 3%]
........................................................................ [ 3%]
........................................................................ [ 3%]
........................................................................ [ 4%]
........................................................................ [ 4%]
........................................................................ [ 4%]
........................................................................ [ 4%]
........................................................................ [ 4%]
........................................................................ [ 5%]
........................................................................ [ 5%]
........................................................................ [ 5%]
........................................................................ [ 5%]
........................................................................ [ 5%]
........................................................................ [ 6%]
........................................................................ [ 6%]
........................................................................ [ 6%]
........................................................................ [ 6%]
........................................................................ [ 6%]
........................................................................ [ 7%]
........................................................................ [ 7%]
....................................................... [ 7%]
scipy/fft/tests/test_backend.py ........................................ [ 7%]
.... [ 7%]
scipy/fft/tests/test_fft_function.py . [ 7%]
scipy/fft/tests/test_fftlog.py ........................................ [ 7%]
scipy/fft/tests/test_helper.py .................. [ 7%]
scipy/fft/tests/test_multithreading.py ................................. [ 7%]
....................... [ 7%]
scipy/fft/tests/test_numpy.py .......................................... [ 8%]
........................................................................ [ 8%]
............ [ 8%]
scipy/fft/tests/test_real_transforms.py ................................ [ 8%]
........................................................................ [ 8%]
........................................................................ [ 8%]
........................................................................ [ 9%]
........................................................................ [ 9%]
........................................................................ [ 9%]
........................................................................ [ 9%]
........................................................................ [ 9%]
........................................................................ [ 10%]
........................................................................ [ 10%]
........................................................................ [ 10%]
........................................................................ [ 10%]
........................................................................ [ 10%]
........................................................................ [ 11%]
........................................................................ [ 11%]
........................................................................ [ 11%]
........................................................................ [ 11%]
........................................................................ [ 11%]
........................................................................ [ 12%]
........................................................................ [ 12%]
........................................................................ [ 12%]
........................................................................ [ 12%]
........................................................................ [ 12%]
........................................................................ [ 13%]
........................................................................ [ 13%]
........................................................................ [ 13%]
........................................................................ [ 13%]
........................................................................ [ 14%]
........................................................................ [ 14%]
........................................................................ [ 14%]
........................................................................ [ 14%]
........................................................................ [ 14%]
........................................................................ [ 15%]
[ 15%]
scipy/fftpack/tests/test_basic.py ......x............................... [ 15%]
........................................................................ [ 15%]
........................................................................ [ 15%]
........................................................................ [ 15%]
.......................................... [ 15%]
scipy/fftpack/tests/test_helper.py .... [ 15%]
scipy/fftpack/tests/test_import.py . [ 15%]
scipy/fftpack/tests/test_pseudo_diffs.py ............................. [ 15%]
scipy/fftpack/tests/test_real_transforms.py ............................ [ 16%]
........................................................................ [ 16%]
........................................................................ [ 16%]
.................................................... [ 16%]
scipy/integrate/_ivp/tests/test_ivp.py ........................ [ 16%]
scipy/integrate/_ivp/tests/test_rk.py .......... [ 16%]
scipy/integrate/tests/test__quad_vec.py .................... [ 16%]
scipy/integrate/tests/test_banded_ode_solvers.py F [ 16%]
scipy/integrate/tests/test_bvp.py ................ [ 16%]
scipy/integrate/tests/test_integrate.py ..F....F........................ [ 16%]
FFFF........... [ 16%]
scipy/integrate/tests/test_odeint_jac.py .. [ 16%]
scipy/integrate/tests/test_quadpack.py ................................ [ 17%]
scipy/integrate/tests/test_quadrature.py ...................... [ 17%]
scipy/interpolate/tests/test_bsplines.py ............................... [ 17%]
........................................x....................... [ 17%]
scipy/interpolate/tests/test_fitpack.py ................s.... [ 17%]
scipy/interpolate/tests/test_fitpack2.py ............................... [ 17%]
....................... [ 17%]
scipy/interpolate/tests/test_interpnd.py ..................... [ 17%]
scipy/interpolate/tests/test_interpolate.py ............................ [ 17%]
........................................................................ [ 17%]
......................................................... [ 18%]
scipy/interpolate/tests/test_ndgriddata.py ............ [ 18%]
scipy/interpolate/tests/test_pade.py .... [ 18%]
scipy/interpolate/tests/test_polyint.py .........................F...... [ 18%]
..................... [ 18%]
scipy/interpolate/tests/test_rbf.py .......... [ 18%]
scipy/interpolate/tests/test_rbfinterp.py .............................. [ 18%]
........................................................................ [ 18%]
........................................................................ [ 18%]
........ [ 18%]
scipy/interpolate/tests/test_regression.py . [ 18%]
scipy/io/arff/tests/test_arffread.py ........................... [ 18%]
scipy/io/harwell_boeing/tests/test_fortran_format.py ........... [ 18%]
scipy/io/harwell_boeing/tests/test_hb.py .. [ 18%]
scipy/io/matlab/tests/test_byteordercodes.py .. [ 18%]
scipy/io/matlab/tests/test_mio.py ...................................... [ 19%]
................... [ 19%]
scipy/io/matlab/tests/test_mio5_utils.py ...... [ 19%]
scipy/io/matlab/tests/test_mio_funcs.py . [ 19%]
scipy/io/matlab/tests/test_mio_utils.py .. [ 19%]
scipy/io/matlab/tests/test_miobase.py . [ 19%]
scipy/io/matlab/tests/test_pathological.py .. [ 19%]
scipy/io/matlab/tests/test_streams.py ........... [ 19%]
scipy/io/tests/test_fortran.py ........... [ 19%]
scipy/io/tests/test_idl.py ............................................. [ 19%]
..................... [ 19%]
scipy/io/tests/test_mmio.py ............................................ [ 19%]
.............................. [ 19%]
scipy/io/tests/test_netcdf.py ......................... [ 19%]
scipy/io/tests/test_paths.py ........... [ 19%]
scipy/io/tests/test_wavfile.py ....................... [ 19%]
scipy/linalg/tests/test_basic.py ......................................F [ 19%]
.....................s...............F..F.F............................. [ 20%]
.s............................ [ 20%]
scipy/linalg/tests/test_blas.py ...............F........................ [ 20%]
........... [ 20%]
scipy/linalg/tests/test_build.py s [ 20%]
scipy/linalg/tests/test_cython_blas.py ...... [ 20%]
scipy/linalg/tests/test_cython_lapack.py .. [ 20%]
scipy/linalg/tests/test_decomp.py ........x..F......FF.F..............FF [ 20%]
..FF......FF.........................................F.F.F.F.F.F.F.F.F.F [ 20%]
.F.F.F.F.F.F.............................................FF............. [ 20%]
........................................................................ [ 21%]
......FFF..FFFFFF..F..s.x............................... [ 21%]
scipy/linalg/tests/test_decomp_cholesky.py ................ [ 21%]
scipy/linalg/tests/test_decomp_cossin.py ..............FF......FF....... [ 21%]
.......FF..............FF......FF..............FFF.......... [ 21%]
scipy/linalg/tests/test_decomp_ldl.py ..F.F [ 21%]
scipy/linalg/tests/test_decomp_polar.py .. [ 21%]
scipy/linalg/tests/test_decomp_update.py ............................... [ 21%]
........................................................................ [ 21%]
........................................................................ [ 22%]
........................................................................ [ 22%]
........................................................................ [ 22%]
........................................................................ [ 22%]
........................................................................ [ 22%]
........................................................................ [ 23%]
........................................................................ [ 23%]
.......................................................... [ 23%]
scipy/linalg/tests/test_fblas.py ....................................... [ 23%]
........................................................................ [ 23%]
.................. [ 23%]
scipy/linalg/tests/test_interpolative.py ...... [ 23%]
scipy/linalg/tests/test_lapack.py ........F............................. [ 24%]
.........................................s...........................F.. [ 24%]
..s..................................................................... [ 24%]
........................................................................ [ 24%]
........................................................................ [ 24%]
........................................................................ [ 25%]
........................................................................ [ 25%]
........................................................................ [ 25%]
........................................................................ [ 25%]
........................................................................ [ 25%]
........................................................................ [ 26%]
........................................................................ [ 26%]
........................................................................ [ 26%]
........................................................................ [ 26%]
........................................................................ [ 26%]
........................................................................ [ 27%]
........................................................................ [ 27%]
........................................................................ [ 27%]
........................................................................ [ 27%]
........................................................................ [ 27%]
........................................................................ [ 28%]
........................................................................ [ 28%]
........................................................................ [ 28%]
........................................................................ [ 28%]
...FF..FF............................................................... [ 28%]
.................FF..................................................... [ 29%]
.................................................FF [ 29%]
scipy/linalg/tests/test_matfuncs.py .................................... [ 29%]
.X............... [ 29%]
scipy/linalg/tests/test_matmul_toeplitz.py ... [ 29%]
scipy/linalg/tests/test_misc.py . [ 29%]
scipy/linalg/tests/test_procrustes.py .......... [ 29%]
scipy/linalg/tests/test_sketches.py ..... [ 29%]
scipy/linalg/tests/test_solve_toeplitz.py ......x [ 29%]
scipy/linalg/tests/test_solvers.py ..FF..... [ 29%]
scipy/linalg/tests/test_special_matrices.py ............................ [ 29%]
........................................................................ [ 29%]
................... [ 29%]
scipy/misc/tests/test_common.py ... [ 29%]
scipy/misc/tests/test_doccer.py ..... [ 29%]
scipy/ndimage/tests/test_c_api.py ... [ 29%]
scipy/ndimage/tests/test_datatypes.py .. [ 29%]
scipy/ndimage/tests/test_filters.py .................................... [ 30%]
........................................................................ [ 30%]
........................................................................ [ 30%]
........................................................................ [ 30%]
........................................................................ [ 30%]
........................................................................ [ 31%]
........................................................................ [ 31%]
........................................................................ [ 31%]
........................................................................ [ 31%]
........................................................................ [ 31%]
........................................................................ [ 32%]
........................................................................ [ 32%]
........................................................................ [ 32%]
........................................................................ [ 32%]
........................................................................ [ 32%]
........................................................................ [ 33%]
........................................................................ [ 33%]
........................................................................ [ 33%]
........................................................................ [ 33%]
........................................................................ [ 33%]
........................................................................ [ 34%]
...................................................................... [ 34%]
scipy/ndimage/tests/test_fourier.py .................................... [ 34%]
........................................ [ 34%]
scipy/ndimage/tests/test_interpolation.py .............................. [ 34%]
........................................................................ [ 34%]
........................................................................ [ 35%]
........................................................................ [ 35%]
........................................................................ [ 35%]
........................................................................ [ 35%]
........................................................................ [ 35%]
........................................................................ [ 36%]
........................................................................ [ 36%]
........................................................................ [ 36%]
........................................................................ [ 36%]
........................................................................ [ 36%]
........................................................................ [ 37%]
........................................................................ [ 37%]
...................... [ 37%]
scipy/ndimage/tests/test_measurements.py ............................... [ 37%]
........................................................................ [ 37%]
................ [ 37%]
scipy/ndimage/tests/test_morphology.py ................................. [ 37%]
........................................................................ [ 38%]
........................................................................ [ 38%]
........................................................................ [ 38%]
........................................................................ [ 38%]
........................................................................ [ 38%]
........................................................................ [ 39%]
........................................................................ [ 39%]
........................................................................ [ 39%]
........................................................................ [ 39%]
........................................................................ [ 39%]
........................................................................ [ 40%]
.......... [ 40%]
scipy/ndimage/tests/test_splines.py .................. [ 40%]
scipy/odr/tests/test_odr.py ................s [ 40%]
scipy/optimize/_trustregion_constr/tests/test_canonical_constraint.py .. [ 40%]
.... [ 40%]
scipy/optimize/_trustregion_constr/tests/test_projections.py .......... [ 40%]
scipy/optimize/_trustregion_constr/tests/test_qp_subproblem.py ......... [ 40%]
.............. [ 40%]
scipy/optimize/_trustregion_constr/tests/test_report.py .. [ 40%]
scipy/optimize/tests/test__basinhopping.py ............................. [ 40%]
[ 40%]
scipy/optimize/tests/test__differential_evolution.py ................... [ 40%]
................................ [ 40%]
scipy/optimize/tests/test__dual_annealing.py ........................... [ 40%]
............... [ 40%]
scipy/optimize/tests/test__linprog_clean_inputs.py ............. [ 40%]
scipy/optimize/tests/test__numdiff.py ............................... [ 40%]
scipy/optimize/tests/test__remove_redundancy.py ........................ [ 40%]
............................................ [ 41%]
scipy/optimize/tests/test__root.py .... [ 41%]
scipy/optimize/tests/test__shgo.py ..................................... [ 41%]
........... [ 41%]
scipy/optimize/tests/test__spectral.py .... [ 41%]
scipy/optimize/tests/test_cobyla.py .... [ 41%]
scipy/optimize/tests/test_constraint_conversion.py ........ [ 41%]
scipy/optimize/tests/test_constraints.py .......... [ 41%]
scipy/optimize/tests/test_cython_optimize.py ..... [ 41%]
scipy/optimize/tests/test_differentiable_functions.py .......... [ 41%]
scipy/optimize/tests/test_hessian_update_strategy.py .... [ 41%]
scipy/optimize/tests/test_lbfgsb_hessinv.py .. [ 41%]
scipy/optimize/tests/test_lbfgsb_setulb.py . [ 41%]
scipy/optimize/tests/test_least_squares.py ............................. [ 41%]
........................................................................ [ 41%]
............. [ 41%]
scipy/optimize/tests/test_linear_assignment.py ...................... [ 41%]
scipy/optimize/tests/test_linesearch.py ........... [ 41%]
scipy/optimize/tests/test_linprog.py ................................... [ 41%]
.........................................sss............................ [ 42%]
...............................................ss....................... [ 42%]
..................................................s..s.ss............... [ 42%]
............................................................s........... [ 42%]
..............................................................s.X....... [ 42%]
...................................................................s.sX. [ 43%]
........................................................................ [ 43%]
...sssss................................................................ [ 43%]
..............ssss...................................................... [ 43%]
....................s................................................... [ 43%]
.....................................s.................................. [ 44%]
..... [ 44%]
scipy/optimize/tests/test_lsq_common.py .......... [ 44%]
scipy/optimize/tests/test_lsq_linear.py .................. [ 44%]
scipy/optimize/tests/test_minimize_constrained.py ................X [ 44%]
scipy/optimize/tests/test_minpack.py ................................... [ 44%]
............................. [ 44%]
scipy/optimize/tests/test_nnls.py .. [ 44%]
scipy/optimize/tests/test_nonlin.py ............................... [ 44%]
scipy/optimize/tests/test_optimize.py ............X.............X....... [ 44%]
......X.............X................................................... [ 44%]
....x................................................................... [ 45%]
................. [ 45%]
scipy/optimize/tests/test_quadratic_assignment.py .............. [ 45%]
scipy/optimize/tests/test_regression.py ... [ 45%]
scipy/optimize/tests/test_slsqp.py ..................................... [ 45%]
[ 45%]
scipy/optimize/tests/test_tnc.py .................. [ 45%]
scipy/optimize/tests/test_trustregion.py ........ [ 45%]
scipy/optimize/tests/test_trustregion_exact.py .......... [ 45%]
scipy/optimize/tests/test_trustregion_krylov.py ..... [ 45%]
scipy/optimize/tests/test_zeros.py .................................. [ 45%]
scipy/signal/tests/test_array_tools.py ...... [ 45%]
scipy/signal/tests/test_bsplines.py .............. [ 45%]
scipy/signal/tests/test_cont2discrete.py ........................ [ 45%]
scipy/signal/tests/test_dltisys.py ................................ [ 45%]
scipy/signal/tests/test_filter_design.py ............................... [ 45%]
........................................................................ [ 45%]
........................................................................ [ 46%]
...... [ 46%]
scipy/signal/tests/test_fir_filter_design.py ........................... [ 46%]
...... [ 46%]
scipy/signal/tests/test_ltisys.py .F.................................... [ 46%]
............................................................. [ 46%]
scipy/signal/tests/test_max_len_seq.py .. [ 46%]
scipy/signal/tests/test_peak_finding.py ................................ [ 46%]
................ [ 46%]
scipy/signal/tests/test_result_type.py ..... [ 46%]
scipy/signal/tests/test_savitzky_golay.py ............ [ 46%]
scipy/signal/tests/test_signaltools.py ................................. [ 46%]
........................................................................ [ 47%]
........................................................................ [ 47%]
........................................................................ [ 47%]
........................................................................ [ 47%]
........................................................................ [ 47%]
........................................................................ [ 48%]
........................................................................ [ 48%]
........................................................................ [ 48%]
........................................................................ [ 48%]
........................................................................ [ 48%]
........................................................................ [ 49%]
........................................................................ [ 49%]
........................................................................ [ 49%]
........................................................................ [ 49%]
........................................................................ [ 49%]
........................................................................ [ 50%]
........................................................................ [ 50%]
........................................................................ [ 50%]
........................................................................ [ 50%]
........................................................................ [ 50%]
........................................................................ [ 51%]
........................................................................ [ 51%]
........................................................................ [ 51%]
...............................................ss....................... [ 51%]
........................................................................ [ 52%]
........... [ 52%]
scipy/signal/tests/test_spectral.py .................................... [ 52%]
........................................................................ [ 52%]
.. [ 52%]
scipy/signal/tests/test_upfirdn.py ..................................... [ 52%]
........................................................................ [ 52%]
........................................................................ [ 52%]
........................................................................ [ 53%]
....................................... [ 53%]
scipy/signal/tests/test_waveforms.py .................................. [ 53%]
scipy/signal/tests/test_wavelets.py ....... [ 53%]
scipy/signal/tests/test_windows.py ..................................... [ 53%]
........ [ 53%]
scipy/sparse/csgraph/tests/test_connected_components.py ...... [ 53%]
scipy/sparse/csgraph/tests/test_conversions.py ... [ 53%]
scipy/sparse/csgraph/tests/test_flow.py ................................ [ 53%]
........... [ 53%]
scipy/sparse/csgraph/tests/test_graph_laplacian.py .... [ 53%]
scipy/sparse/csgraph/tests/test_matching.py ....................s [ 53%]
scipy/sparse/csgraph/tests/test_reordering.py ... [ 53%]
scipy/sparse/csgraph/tests/test_shortest_path.py ....................... [ 53%]
....... [ 53%]
scipy/sparse/csgraph/tests/test_spanning_tree.py . [ 53%]
scipy/sparse/csgraph/tests/test_traversal.py .... [ 53%]
scipy/sparse/linalg/dsolve/tests/test_linsolve.py .s.s.s.s.ss.......s... [ 53%]
................. [ 53%]
scipy/sparse/linalg/eigen/arpack/tests/test_arpack.py .F..FFFFF......... [ 53%]
. [ 53%]
scipy/sparse/linalg/eigen/lobpcg/tests/test_lobpcg.py .......F...... [ 53%]
scipy/sparse/linalg/eigen/tests/test_svds.py .........F.. [ 54%]
scipy/sparse/linalg/isolve/tests/test_gcrotmk.py ....... [ 54%]
scipy/sparse/linalg/isolve/tests/test_iterative.py ..................... [ 54%]
.....xxxx................. [ 54%]
scipy/sparse/linalg/isolve/tests/test_lgmres.py ........ [ 54%]
scipy/sparse/linalg/isolve/tests/test_lsmr.py .............. [ 54%]
scipy/sparse/linalg/isolve/tests/test_lsqr.py ..... [ 54%]
scipy/sparse/linalg/isolve/tests/test_minres.py ....... [ 54%]
scipy/sparse/linalg/isolve/tests/test_utils.py . [ 54%]
scipy/sparse/linalg/tests/test_expm_multiply.py ................. [ 54%]
scipy/sparse/linalg/tests/test_interface.py .............. [ 54%]
scipy/sparse/linalg/tests/test_matfuncs.py ............................. [ 54%]
.... [ 54%]
scipy/sparse/linalg/tests/test_norm.py ....... [ 54%]
scipy/sparse/linalg/tests/test_onenormest.py ss.s.. [ 54%]
scipy/sparse/linalg/tests/test_pydata_sparse.py ssssssssssssssssssssssss [ 54%]
ssssssssssssssssssssssssssss [ 54%]
scipy/sparse/tests/test_base.py ........................................ [ 54%]
......................................s................................. [ 54%]
........................................................................ [ 55%]
..............................................s......................... [ 55%]
...........................................s.ssssss..................... [ 55%]
.................................................s...................... [ 55%]
............................sss.ss........ssssss........................ [ 55%]
..............................................s......................... [ 56%]
.........................sss.s...........ssssss......................... [ 56%]
............................................s......sssss.sssssssssssssss [ 56%]
sssssssssssssssssssss.............x........ssssss....................... [ 56%]
..............................................s......sssss.sssssssssssss [ 57%]
sssssssssssssssssssssss..sss.s.x.s...................................... [ 57%]
......................................s......sssss.sssssssssssssssssssss [ 57%]
sssssssssssssss.......s...........xx.x.................................. [ 57%]
..........................................s............................. [ 57%]
.........................................sss..xx........................ [ 58%]
...................................................s.................... [ 58%]
.............................................sss..xx.................... [ 58%]
....................................................s......sssss.sssssss [ 58%]
sssssssssssssssssssssssssssss..................xx.x.sssxx.ssssss........ [ 58%]
............................................................s......sssss [ 59%]
.ssssssssssssssssssssssssssssssssssss............x......sss..........ss. [ 59%]
ss.................s......s...sssssssssssssssss.....sss.ss.....s..s.x... [ 59%]
.............s..ss.................sssssxxs..sssss.s.s............s..... [ 59%]
..ss.ss.............s.....s..ssssssssssssssssss....ssss.s..s...s..s.xsss [ 59%]
..............ss.ss.................ssss..s..sssss.s.s...............s.. [ 60%]
..................................s....................................s [ 60%]
........................................s..............s................ [ 60%]
s.........................ss.............s.............................. [ 60%]
........s........................................s..............s....... [ 60%]
.........s......s.ss.ss.........s....s..ssssssssssssssssss....ssss.ss.s. [ 61%]
..s..s.xss.........sss..ss.ss................ssss...s..sssss.s.s........ [ 61%]
....s......s........s................ss....................s.....s.s.... [ 61%]
..s..ss.........sss....s..................................s............. [ 61%]
..s.....s.......................ss.....................s.......s.......s [ 61%]
..s...s.........sss...s....................................s............ [ 62%]
..s.......ss.ss.................s......s...sssssssssssssssss.....sss.ss. [ 62%]
....s..s.x................s..ss.................sssssxxs..sssss.s.s..... [ 62%]
.......s.......ss.ss.............s.....s..ssssssssssssssssss....ssss.s.. [ 62%]
s...s..s.xsss..............ss.ss.................ssss..s..sssss.s.s..... [ 62%]
..........s....................................s........................ [ 63%]
............s........................................s..............s... [ 63%]
.............s.........................ss.............s................. [ 63%]
.....................s........................................s......... [ 63%]
.....s................s......s.ss.ss.........s....s..ssssssssssssssssss. [ 64%]
...ssss.ss.s...s..s.xss.........sss..ss.ss................ssss...s..ssss [ 64%]
s.s.s............s......s........s................ss.................... [ 64%]
s.....s.s......s..ss.........sss....s..................................s [ 64%]
...............s.....s.......................ss.....................s... [ 64%]
....s.......s..s...s.........sss...s.................................... [ 65%]
s..............s.......ss.ss.................s......s...ssssssssssssssss [ 65%]
s.....sss.ss.....s..s.x................s..ss.................sssssxxs..s [ 65%]
ssss.s.s............s.......ss.ss.............s.....s..sssssssssssssssss [ 65%]
s....ssss.s..s...s..s.xsss..............ss.ss.................ssss..s..s [ 65%]
ssss.s.s...............s....................................s........... [ 66%]
.........................s........................................s..... [ 66%]
.........s................s.........................ss.............s.... [ 66%]
..................................s..................................... [ 66%]
...s..............s................s......s.ss.ss.........s....s..ssssss [ 66%]
ssssssssssss....ssss.ss.s...s..s.xss.........sss..ss.ss................s [ 67%]
sss...s..sssss.s.s............s......s........s................ss....... [ 67%]
.............s.....s.s......s..ss.........sss....s...................... [ 67%]
............s...............s.....s.......................ss............ [ 67%]
.........s.......s.......s..s...s.........sss...s....................... [ 67%]
.............s..............s.......ss.ss.................s......s...sss [ 68%]
ssssssssssssss.....sss.ss.....s..s.x................s..ss............... [ 68%]
..sssssxxs..sssss.s.s............s.......ss.ss.............s.....s..ssss [ 68%]
ssssssssssssss....ssss.s..s...s..s.xsss..............ss.ss.............. [ 68%]
...ssss..s..sssss.s.s...............s................................... [ 68%]
.s....................................s................................. [ 69%]
.......s..............s................s.........................ss..... [ 69%]
........s......................................s........................ [ 69%]
................s..............s................s......s.ss.ss.........s [ 69%]
....s..ssssssssssssssssss....ssss.ss.s...s..s.xss.........sss..ss.ss.... [ 69%]
............ssss...s..sssss.s.s............s......s........s............ [ 70%]
....ss....................s.....s.s......s..ss.........sss....s......... [ 70%]
.........................s...............s.....s.......................s [ 70%]
s.....................s.......s.......s..s...s.........sss...s.......... [ 70%]
..........................s..............s.......ss.ss.................s [ 71%]
......s...sssssssssssssssss.....sss.ss.....s..s.x................s..ss.. [ 71%]
...............sssssxxs..sssss.s.s............s.......ss.ss............. [ 71%]
s.....s..ssssssssssssssssss....ssss.s..s...s..s.xsss..............ss.ss. [ 71%]
................ssss..s..sssss.s.s...............s...................... [ 71%]
..............s....................................s.................... [ 72%]
....................s..............s................s................... [ 72%]
......ss.............s......................................s........... [ 72%]
.............................s..............s................s......s.ss [ 72%]
.ss.........s....s..ssssssssssssssssss....ssss.ss.s...s..s.xss.........s [ 72%]
ss..ss.ss................ssss...s..sssss.s.s............s......s........ [ 73%]
s................ss....................s.....s.s......s..ss.........sss. [ 73%]
...s..................................s...............s.....s........... [ 73%]
............ss.....................s.......s.......s..s...s.........sss. [ 73%]
..s....................................s..............s. [ 73%]
scipy/sparse/tests/test_construct.py ........................... [ 73%]
scipy/sparse/tests/test_csc.py ........... [ 73%]
scipy/sparse/tests/test_csr.py ........ [ 73%]
scipy/sparse/tests/test_extract.py .. [ 73%]
scipy/sparse/tests/test_matrix_io.py ...... [ 74%]
scipy/sparse/tests/test_sparsetools.py ...ss.. [ 74%]
scipy/sparse/tests/test_spfuncs.py ... [ 74%]
scipy/sparse/tests/test_sputils.py ............. [ 74%]
scipy/spatial/tests/test__plotutils.py ... [ 74%]
scipy/spatial/tests/test__procrustes.py ...... [ 74%]
scipy/spatial/tests/test_distance.py ................................... [ 74%]
........................................................................ [ 74%]
....................................................... [ 74%]
scipy/spatial/tests/test_hausdorff.py ............. [ 74%]
scipy/spatial/tests/test_kdtree.py ..................................... [ 74%]
........................................................................ [ 74%]
........................................................................ [ 75%]
................................................... [ 75%]
scipy/spatial/tests/test_qhull.py ...................................... [ 75%]
........................................................................ [ 75%]
.......................................................... [ 75%]
scipy/spatial/tests/test_slerp.py ...................................... [ 75%]
........................................... [ 75%]
scipy/spatial/tests/test_spherical_voronoi.py .......................... [ 76%]
........................................................................ [ 76%]
...................... [ 76%]
scipy/spatial/transform/tests/test_rotation.py ......................... [ 76%]
................................................................... [ 76%]
scipy/spatial/transform/tests/test_rotation_groups.py .................. [ 76%]
........................................................................ [ 76%]
........................................................................ [ 77%]
...................................................................... [ 77%]
scipy/spatial/transform/tests/test_rotation_spline.py ....... [ 77%]
scipy/special/tests/test_basic.py ...................................... [ 77%]
..............x...X..................................................... [ 77%]
........................................................................ [ 77%]
........................................................................ [ 78%]
........................................................................ [ 78%]
........................................................................ [ 78%]
............................... [ 78%]
scipy/special/tests/test_bdtr.py ........................ [ 78%]
scipy/special/tests/test_boxcox.py ........ [ 78%]
scipy/special/tests/test_cdflib.py .. [ 78%]
scipy/special/tests/test_cosine_distr.py ............................... [ 78%]
[ 78%]
scipy/special/tests/test_cython_special.py ............................. [ 78%]
.......................xxxxx..xx..xx.................................... [ 78%]
........................................................................ [ 79%]
xxxx...........xxxx.................................... [ 79%]
scipy/special/tests/test_data.py ....................................... [ 79%]
........................................................................ [ 79%]
............x.x...............................................x..x...x.. [ 79%]
..................... [ 79%]
scipy/special/tests/test_digamma.py ... [ 79%]
scipy/special/tests/test_ellip_harm.py ..... [ 79%]
scipy/special/tests/test_erfinv.py ............. [ 79%]
scipy/special/tests/test_exponential_integrals.py ........ [ 80%]
scipy/special/tests/test_faddeeva.py ................. [ 80%]
scipy/special/tests/test_gamma.py .. [ 80%]
scipy/special/tests/test_gammainc.py ............................. [ 80%]
scipy/special/tests/test_hyp2f1.py .xx.xxxxxx..x.x...x..FF....F......F.. [ 80%]
.x..xx [ 80%]
scipy/special/tests/test_hypergeometric.py ............................. [ 80%]
......................... [ 80%]
scipy/special/tests/test_kolmogorov.py ................................. [ 80%]
.... [ 80%]
scipy/special/tests/test_lambertw.py ... [ 80%]
scipy/special/tests/test_log_softmax.py .......... [ 80%]
scipy/special/tests/test_loggamma.py ...... [ 80%]
scipy/special/tests/test_logit.py ...... [ 80%]
scipy/special/tests/test_logsumexp.py .......... [ 80%]
scipy/special/tests/test_mpmath.py ................... [ 80%]
scipy/special/tests/test_nan_inputs.py ................................. [ 80%]
........................................................................ [ 80%]
........................................................................ [ 81%]
................................................. [ 81%]
scipy/special/tests/test_ndtr.py .... [ 81%]
scipy/special/tests/test_ndtri_exp.py ............ [ 81%]
scipy/special/tests/test_orthogonal.py ........................... [ 81%]
scipy/special/tests/test_orthogonal_eval.py ............................ [ 81%]
.................................................................. [ 81%]
scipy/special/tests/test_owens_t.py .... [ 81%]
scipy/special/tests/test_pcf.py .. [ 81%]
scipy/special/tests/test_pdtr.py .......... [ 81%]
scipy/special/tests/test_precompute_expn_asy.py s [ 81%]
scipy/special/tests/test_precompute_gammainc.py .s.s [ 81%]
scipy/special/tests/test_round.py .. [ 81%]
scipy/special/tests/test_sf_error.py ....... [ 81%]
scipy/special/tests/test_sici.py .. [ 81%]
scipy/special/tests/test_spence.py .. [ 81%]
scipy/special/tests/test_spfun_stats.py .... [ 81%]
scipy/special/tests/test_sph_harm.py . [ 81%]
scipy/special/tests/test_spherical_bessel.py ........................... [ 81%]
.......... [ 81%]
scipy/special/tests/test_trig.py .... [ 81%]
scipy/special/tests/test_wright_bessel.py .............................. [ 82%]
........................................................................ [ 82%]
........................................................................ [ 82%]
........................................................................ [ 82%]
............................................................xxxxxxxxxxxx [ 82%]
xxx............... [ 82%]
scipy/special/tests/test_wrightomega.py .......... [ 82%]
scipy/special/tests/test_zeta.py ..... [ 82%]
scipy/stats/tests/test_binned_statistic.py ............................. [ 83%]
............ [ 83%]
scipy/stats/tests/test_bootstrap.py ........................x..x..x..... [ 83%]
................................ [ 83%]
scipy/stats/tests/test_contingency.py ........... [ 83%]
scipy/stats/tests/test_continuous_basic.py ....F..................F..... [ 83%]
...................F.........................F..F..s....xx.............. [ 83%]
................F........................F.............................. [ 83%]
.s.....................................F...ss........................... [ 84%]
F....................................s.................................. [ 84%]
.........................FF............................................. [ 84%]
........................................................................ [ 84%]
.............................F.FFFF..................................... [ 84%]
........................................................................ [ 85%]
........................................................................ [ 85%]
........................................................................ [ 85%]
........................................................................ [ 85%]
........................................................................ [ 85%]
........................................................................ [ 86%]
.................................FF..............................FF..... [ 86%]
........................................................................ [ 86%]
............................................................ [ 86%]
scipy/stats/tests/test_crosstab.py ............ [ 86%]
scipy/stats/tests/test_discrete_basic.py .............F................. [ 86%]
...F.................................................................... [ 86%]
........................................................................ [ 87%]
........................................................................ [ 87%]
........................................................................ [ 87%]
......... [ 87%]
scipy/stats/tests/test_discrete_distns.py .............................. [ 87%]
..................F..... [ 87%]
scipy/stats/tests/test_distributions.py ................................ [ 87%]
........................................................................ [ 88%]
........................................................................ [ 88%]
........................................................................ [ 88%]
........................................................................ [ 88%]
s............................................................s.......... [ 88%]
............s.......................F..F................................ [ 89%]
........................................................................ [ 89%]
........................................................................ [ 89%]
.................s.....s................................................ [ 89%]
........................................................................ [ 89%]
........................................................................ [ 90%]
........................................................................ [ 90%]
........s...............s................ [ 90%]
scipy/stats/tests/test_entropy.py .................................... [ 90%]
scipy/stats/tests/test_fit.py .. [ 90%]
scipy/stats/tests/test_hypotests.py .................................... [ 90%]
........................................................................ [ 90%]
............................................ [ 91%]
scipy/stats/tests/test_kdeoth.py ..........ss......ss....sssssssssssssss [ 91%]
s..ss......ss......ss......ss....ssssssssssssssss..ss......ss......ss... [ 91%]
...ss....ssssssssssssssss..ss......ss......ss......ss....sssssssssssssss [ 91%]
s..ss......ss......ss......ss....ssssssssssssssss..ss......ss......ss... [ 91%]
...ss....ssssssssssssssss..ss......ss......ss......ss....sssssssssssssss [ 91%]
s..ss......ss......ss......ss....ssssssssssssssss..ss......ss......ss... [ 92%]
...ss....ssssssssssssssss..ss......ss......ss......ss....sssssssssssssss [ 92%]
s..ss......ss......ss......ss....ssssssssssssssss..ss......ss......ss... [ 92%]
...ss....ssssssssssssssss..ss......ss....sssssssssssssssssssssssssssssss [ 92%]
ssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssss [ 92%]
ssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssss [ 93%]
ssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssss [ 93%]
ssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssss [ 93%]
ssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssss [ 93%]
ssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssss [ 94%]
ssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssss [ 94%]
sssssssssssssssssssssssssssssssssssssssss..ss......ss....sssssssssssssss [ 94%]
s..ss......ss......ss......ss....ssssssssssssssss..ss......ss......ss... [ 94%]
...ss....ssssssssssssssss..ss......ss......ss......ss....sssssssssssssss [ 94%]
s..ss......ss......ss......ss....ssssssssssssssss..ss......ss......ss... [ 95%]
...ss....ssssssssssssssss..ss......ss......ss......ss....sssssssssssssss [ 95%]
s..ss......ss......ss......ss....ssssssssssssssss..ss......ss......ss... [ 95%]
...ss....ssssssssssssssss..ss......ss......ss......ss....sssssssssssssss [ 95%]
s..ss......ss......ss......ss....ssssssssssssssss..ss......ss......ss... [ 95%]
...ss....ssssssssssssssss..ss......ss......s... [ 95%]
scipy/stats/tests/test_morestats.py .................................... [ 96%]
........................................................................ [ 96%]
...................xxx.................................................. [ 96%]
........................................................................ [ 96%]
.............................. [ 96%]
scipy/stats/tests/test_mstats_basic.py ......s..s....................... [ 96%]
........................................................................ [ 97%]
........................ [ 97%]
scipy/stats/tests/test_mstats_extras.py ......... [ 97%]
scipy/stats/tests/test_multivariate.py ................................. [ 97%]
........................................................................ [ 97%]
.................................................................... [ 97%]
scipy/stats/tests/test_qmc.py ..............................s.s.s.s.s.ss [ 97%]
ss............................s.s.s.s.s.s.sss................ [ 97%]
scipy/stats/tests/test_rank.py ....................... [ 98%]
scipy/stats/tests/test_relative_risk.py ............... [ 98%]
scipy/stats/tests/test_stats.py ........................................ [ 98%]
........................................................................ [ 98%]
........................................................................ [ 98%]
........................................................................ [ 98%]
........................................................................ [ 99%]
........................................................................ [ 99%]
........................................................................ [ 99%]
.............x......................x.....xx......................ss.... [ 99%]
.s..........xx..x......F...F.sx......................................... [ 99%]
................................F............... [ 99%]
scipy/stats/tests/test_tukeylambda_stats.py ... [100%]
=================================== FAILURES ===================================
___________________________ test_banded_ode_solvers ____________________________
scipy/integrate/tests/test_banded_ode_solvers.py:218: in test_banded_ode_solvers
check_complex(idx, "zvode", meth, use_jac, with_jac, banded)
a = array([[-0.6+0.3j , 0. +0.j , 0. +0.j , 0. +0.j , 0. +0.j ],
[ 0. +0.j , -0.5+0.25j, 0. +0.j , 0. +... 0. +0.j , 0. +0.j , -0.9+0.45j, 0. +0.j ],
[ 0. +0.j , 0. +0.j , 0. +0.j , 0. +0.j , -0.7+0.35j]])
a_complex = array([[-0.6+0.3j , 0.1-0.05j, 0. +0.j , 0. +0.j , 0. +0.j ],
[ 0.2-0.1j , -0.5+0.25j, 0.9-0.45j, 0. +... 0.3-0.15j, -0.1+0.05j, -0.9+0.45j, -0.3+0.15j],
[ 0. +0.j , 0. +0.j , 0.1-0.05j, 0.1-0.05j, -0.7+0.35j]])
a_complex_diag = array([[-0.6+0.3j , 0. +0.j , 0. +0.j , 0. +0.j , 0. +0.j ],
[ 0. +0.j , -0.5+0.25j, 0. +0.j , 0. +... 0. +0.j , 0. +0.j , -0.9+0.45j, 0. +0.j ],
[ 0. +0.j , 0. +0.j , 0. +0.j , 0. +0.j , -0.7+0.35j]])
a_real = array([[-0.6, 0.1, 0. , 0. , 0. ],
[ 0.2, -0.5, 0.9, 0. , 0. ],
[ 0.1, 0.1, -0.4, 0.1, 0. ],
[ 0. , 0.3, -0.1, -0.9, -0.3],
[ 0. , 0. , 0.1, 0.1, -0.7]])
a_real_diag = array([[-0.6, 0. , 0. , 0. , 0. ],
[ 0. , -0.5, 0. , 0. , 0. ],
[ 0. , 0. , -0.4, 0. , 0. ],
[ 0. , 0. , 0. , -0.9, 0. ],
[ 0. , 0. , 0. , 0. , -0.7]])
a_real_lower = array([[-0.6, 0. , 0. , 0. , 0. ],
[ 0.2, -0.5, 0. , 0. , 0. ],
[ 0.1, 0.1, -0.4, 0. , 0. ],
[ 0. , 0.3, -0.1, -0.9, 0. ],
[ 0. , 0. , 0.1, 0.1, -0.7]])
a_real_upper = array([[-0.6, 0.1, 0. , 0. , 0. ],
[ 0. , -0.5, 0.9, 0. , 0. ],
[ 0. , 0. , -0.4, 0.1, 0. ],
[ 0. , 0. , 0. , -0.9, -0.3],
[ 0. , 0. , 0. , 0. , -0.7]])
banded = False
check_complex = <function test_banded_ode_solvers.<locals>.check_complex at 0x1683b6a60>
check_real = <function test_banded_ode_solvers.<locals>.check_real at 0x168405f70>
complex_matrices = [array([[-0.6+0.3j , 0.1-0.05j, 0. +0.j , 0. +0.j , 0. +0.j ],
[ 0.2-0.1j , -0.5+0.25j, 0.9-0.45j, 0. ... 0. +0.j , 0. +0.j , -0.9+0.45j, 0. +0.j ],
[ 0. +0.j , 0. +0.j , 0. +0.j , 0. +0.j , -0.7+0.35j]])]
complex_solutions = [(array([1.+1.j, 2.+1.j, 3.+1.j, 4.+1.j, 5.+1.j]), array([0. , 0.25, 0.5 , 0.75, 1. ]), array([[1. +1.j ...1j, 1.02529209+0.88779126j,
1.83770279+1.05647438j, 1.2875341 +1.07346963j,
2.16211541+1.31786874j]]))]
idx = 0
meth = 'bdf'
p = [['bdf', 'adams'], [False, True], [False, True], [False, True]]
real_matrices = [array([[-0.6, 0.1, 0. , 0. , 0. ],
[ 0.2, -0.5, 0.9, 0. , 0. ],
[ 0.1, 0.1, -0.4, 0.1, 0. ],
...
[ 0. , 0. , -0.4, 0. , 0. ],
[ 0. , 0. , 0. , -0.9, 0. ],
[ 0. , 0. , 0. , 0. , -0.7]])]
real_solutions = [(array([1, 2, 3, 4, 5]), array([0. , 0.25, 0.5 , 0.75, 1. ]), array([[1. +5.66269967e-17j, 2. +1.0916....37457856, 2.22245466, 2.03662568, 2.95777682],
[0.54881164, 1.21306132, 2.01096014, 1.62627864, 2.48292652]]))]
solver = 'lsoda'
t_exact = array([0. , 0.25, 0.5 , 0.75, 1. ])
use_jac = False
with_jac = False
y0 = array([1.+1.j, 2.+1.j, 3.+1.j, 4.+1.j, 5.+1.j])
y_exact = array([[1. +1.j , 2. +1.j ,
3. +1.j , 4. +1.j ,
...471j, 1.02529209+0.88779126j,
1.83770279+1.05647438j, 1.2875341 +1.07346963j,
2.16211541+1.31786874j]])
scipy/integrate/tests/test_banded_ode_solvers.py:201: in check_complex
t, y = _solve_linear_sys(a, y0,
a = array([[-0.6+0.3j , 0.1-0.05j, 0. +0.j , 0. +0.j , 0. +0.j ],
[ 0.2-0.1j , -0.5+0.25j, 0.9-0.45j, 0. +... 0.3-0.15j, -0.1+0.05j, -0.9+0.45j, -0.3+0.15j],
[ 0. +0.j , 0. +0.j , 0.1-0.05j, 0.1-0.05j, -0.7+0.35j]])
banded = False
complex_matrices = [array([[-0.6+0.3j , 0.1-0.05j, 0. +0.j , 0. +0.j , 0. +0.j ],
[ 0.2-0.1j , -0.5+0.25j, 0.9-0.45j, 0. ... 0. +0.j , 0. +0.j , -0.9+0.45j, 0. +0.j ],
[ 0. +0.j , 0. +0.j , 0. +0.j , 0. +0.j , -0.7+0.35j]])]
complex_solutions = [(array([1.+1.j, 2.+1.j, 3.+1.j, 4.+1.j, 5.+1.j]), array([0. , 0.25, 0.5 , 0.75, 1. ]), array([[1. +1.j ...1j, 1.02529209+0.88779126j,
1.83770279+1.05647438j, 1.2875341 +1.07346963j,
2.16211541+1.31786874j]]))]
idx = 0
meth = 'bdf'
solver = 'zvode'
t_exact = array([0. , 0.25, 0.5 , 0.75, 1. ])
use_jac = False
with_jac = False
y0 = array([1.+1.j, 2.+1.j, 3.+1.j, 4.+1.j, 5.+1.j])
y_exact = array([[1. +1.j , 2. +1.j ,
3. +1.j , 4. +1.j ,
...753j, 3.49020872+1.06056879j,
2.42592262+1.12330981j, 0.93474257+0.9849586j ,
2.57287977+1.38244112j]])
scipy/integrate/tests/test_banded_ode_solvers.py:103: in _solve_linear_sys
r.integrate(r.t + dt)
a = array([[-0.6+0.3j , 0.1-0.05j, 0. +0.j , 0. +0.j , 0. +0.j ],
[ 0.2-0.1j , -0.5+0.25j, 0.9-0.45j, 0. +... 0.3-0.15j, -0.1+0.05j, -0.9+0.45j, -0.3+0.15j],
[ 0. +0.j , 0. +0.j , 0.1-0.05j, 0.1-0.05j, -0.7+0.35j]])
banded = False
dt = 0.25
lband = None
method = 'bdf'
r = <scipy.integrate._ode.ode object at 0x1683c5820>
solver = 'zvode'
t = [0]
t0 = 0
tend = 1.0
uband = None
use_jac = False
with_jacobian = False
y = [array([1.+1.j, 2.+1.j, 3.+1.j, 4.+1.j, 5.+1.j])]
y0 = array([1.+1.j, 2.+1.j, 3.+1.j, 4.+1.j, 5.+1.j])
scipy/integrate/_ode.py:433: in integrate
self._y, self.t = mth(self.f, self.jac or (lambda: None),
mth = <bound method vode.run of <scipy.integrate._ode.zvode object at 0x1683c5b50>>
relax = False
self = <scipy.integrate._ode.ode object at 0x1683c5820>
step = False
t = 0.25
scipy/integrate/_ode.py:1013: in run
warnings.warn('{:s}: {:s}'.format(self.__class__.__name__,
E UserWarning: zvode: Excess work done on this call. (Perhaps wrong MF.)
args = (<function _linear_func at 0x11d4e4e50>, <function ode.integrate.<locals>.<lambda> at 0x16840f040>, array([1.+1.j, 2.+1.j, 3.+1.j, 4.+1.j, 5.+1.j]), 0, 0.25, 1e-09, ...)
f = <function _linear_func at 0x11d4e4e50>
f_params = (array([[-0.6+0.3j , 0.1-0.05j, 0. +0.j , 0. +0.j , 0. +0.j ],
[ 0.2-0.1j , -0.5+0.25j, 0.9-0.45j, 0. ...0.3-0.15j, -0.1+0.05j, -0.9+0.45j, -0.3+0.15j],
[ 0. +0.j , 0. +0.j , 0.1-0.05j, 0.1-0.05j, -0.7+0.35j]]),)
istate = -1
jac = <function ode.integrate.<locals>.<lambda> at 0x16840f040>
jac_params = (array([[-0.6+0.3j , 0.1-0.05j, 0. +0.j , 0. +0.j , 0. +0.j ],
[ 0.2-0.1j , -0.5+0.25j, 0.9-0.45j, 0. ...0.3-0.15j, -0.1+0.05j, -0.9+0.45j, -0.3+0.15j],
[ 0. +0.j , 0. +0.j , 0.1-0.05j, 0.1-0.05j, -0.7+0.35j]]),)
self = <scipy.integrate._ode.zvode object at 0x1683c5b50>
t = 0.0014324127418710484
t0 = 0
t1 = 0.25
unexpected_istate_msg = 'Unexpected istate=-1'
y0 = array([1.+1.j, 2.+1.j, 3.+1.j, 4.+1.j, 5.+1.j])
y1 = array([0.99906963+0.99957098j, 2.00314966+0.99949701j,
2.99921202+1.00021471j, 3.99241394+1.0020111j ,
4.99563249+1.00129052j])
______________________________ TestOde.test_zvode ______________________________
scipy/integrate/tests/test_integrate.py:107: in test_zvode
self._do_problem(problem, 'zvode', 'adams')
problem = <scipy.integrate.tests.test_integrate.SimpleOscillator object at 0x16842cfa0>
problem_cls = <class 'scipy.integrate.tests.test_integrate.SimpleOscillator'>
self = <scipy.integrate.tests.test_integrate.TestOde object at 0x168410580>
scipy/integrate/tests/test_integrate.py:80: in _do_problem
z = ig.integrate(problem.stop_t)
f = <function TestODEClass._do_problem.<locals>.<lambda> at 0x168405e50>
ig = <scipy.integrate._ode.ode object at 0x16842cd00>
integrator = 'zvode'
integrator_params = {}
jac = None
method = 'adams'
problem = <scipy.integrate.tests.test_integrate.SimpleOscillator object at 0x16842cfa0>
self = <scipy.integrate.tests.test_integrate.TestOde object at 0x168410580>
scipy/integrate/_ode.py:433: in integrate
self._y, self.t = mth(self.f, self.jac or (lambda: None),
mth = <bound method vode.run of <scipy.integrate._ode.zvode object at 0x16842c370>>
relax = False
self = <scipy.integrate._ode.ode object at 0x16842cd00>
step = False
t = 1.09
scipy/integrate/_ode.py:1013: in run
warnings.warn('{:s}: {:s}'.format(self.__class__.__name__,
E UserWarning: zvode: Excess work done on this call. (Perhaps wrong MF.)
args = (<function TestODEClass._do_problem.<locals>.<lambda> at 0x168405e50>, <function ode.integrate.<locals>.<lambda> at 0x1684598b0>, array([1. +0.j, 0.1+0.j]), 0.0, 1.09, 1.0000000000000002e-06, ...)
f = <function TestODEClass._do_problem.<locals>.<lambda> at 0x168405e50>
f_params = ()
istate = -1
jac = <function ode.integrate.<locals>.<lambda> at 0x1684598b0>
jac_params = ()
self = <scipy.integrate._ode.zvode object at 0x16842c370>
t = 0.03249684834867783
t0 = 0.0
t1 = 1.09
unexpected_istate_msg = 'Unexpected istate=-1'
y0 = array([1. +0.j, 0.1+0.j])
y1 = array([ 1.0011345 -0.00211519j, -0.03010698-0.00011959j])
__________________________ TestOde.test_concurrent_ok __________________________
scipy/integrate/tests/test_integrate.py:172: in test_concurrent_ok
assert_allclose(r.y, 0.1)
E AssertionError:
E Not equal to tolerance rtol=1e-07, atol=0
E
E Mismatched elements: 1 / 1 (100%)
E Max absolute difference: 0.02059729
E Max relative difference: 0.20597291
E x: array([0.1-0.020597j])
E y: array(0.1)
f = <function TestOde.test_concurrent_ok.<locals>.<lambda> at 0x168459ee0>
k = 0
r = <scipy.integrate._ode.ode object at 0x168471dc0>
r2 = <scipy.integrate._ode.ode object at 0x168471fd0>
self = <scipy.integrate.tests.test_integrate.TestOde object at 0x168471730>
sol = 'zvode'
__________________ TestZVODECheckParameterUse.test_no_params ___________________
scipy/integrate/tests/test_integrate.py:602: in test_no_params
self._check_solver(solver)
self = <scipy.integrate.tests.test_integrate.TestZVODECheckParameterUse object at 0x16846f310>
solver = <scipy.integrate._ode.ode object at 0x16846fdc0>
scipy/integrate/tests/test_integrate.py:597: in _check_solver
solver.integrate(pi)
ic = [1.0, 0.0]
self = <scipy.integrate.tests.test_integrate.TestZVODECheckParameterUse object at 0x16846f310>
solver = <scipy.integrate._ode.ode object at 0x16846fdc0>
scipy/integrate/_ode.py:433: in integrate
self._y, self.t = mth(self.f, self.jac or (lambda: None),
mth = <bound method vode.run of <scipy.integrate._ode.zvode object at 0x16846fd90>>
relax = False
self = <scipy.integrate._ode.ode object at 0x16846fdc0>
step = False
t = 3.141592653589793
scipy/integrate/_ode.py:1013: in run
warnings.warn('{:s}: {:s}'.format(self.__class__.__name__,
E UserWarning: zvode: Excess work done on this call. (Perhaps wrong MF.)
args = (<function f at 0x11d640ee0>, <function jac at 0x11d625280>, array([1.+0.j, 0.+0.j]), 0.0, 3.141592653589793, 1e-07, ...)
f = <function f at 0x11d640ee0>
f_params = ()
istate = -1
jac = <function jac at 0x11d625280>
jac_params = ()
self = <scipy.integrate._ode.zvode object at 0x16846fd90>
t = 0.012148892003676552
t0 = 0.0
t1 = 3.141592653589793
unexpected_istate_msg = 'Unexpected istate=-1'
y0 = array([1.+0.j, 0.+0.j])
y1 = array([ 0.9999262 -7.37968809e-05j, -0.01214859+2.98852130e-07j])
_______________ TestZVODECheckParameterUse.test_one_scalar_param _______________
scipy/integrate/tests/test_integrate.py:610: in test_one_scalar_param
self._check_solver(solver)
omega = 1.0
self = <scipy.integrate.tests.test_integrate.TestZVODECheckParameterUse object at 0x1684c2f10>
solver = <scipy.integrate._ode.ode object at 0x1684c2c10>
scipy/integrate/tests/test_integrate.py:597: in _check_solver
solver.integrate(pi)
ic = [1.0, 0.0]
self = <scipy.integrate.tests.test_integrate.TestZVODECheckParameterUse object at 0x1684c2f10>
solver = <scipy.integrate._ode.ode object at 0x1684c2c10>
scipy/integrate/_ode.py:433: in integrate
self._y, self.t = mth(self.f, self.jac or (lambda: None),
mth = <bound method vode.run of <scipy.integrate._ode.zvode object at 0x1684c2b50>>
relax = False
self = <scipy.integrate._ode.ode object at 0x1684c2c10>
step = False
t = 3.141592653589793
scipy/integrate/_ode.py:1013: in run
warnings.warn('{:s}: {:s}'.format(self.__class__.__name__,
E UserWarning: zvode: Excess work done on this call. (Perhaps wrong MF.)
args = (<function f1 at 0x11d625310>, <function jac1 at 0x11d6253a0>, array([1.+0.j, 0.+0.j]), 0.0, 3.141592653589793, 1e-07, ...)
f = <function f1 at 0x11d625310>
f_params = (1.0,)
istate = -1
jac = <function jac1 at 0x11d6253a0>
jac_params = (1.0,)
self = <scipy.integrate._ode.zvode object at 0x1684c2b50>
t = 0.012148892003676552
t0 = 0.0
t1 = 3.141592653589793
unexpected_istate_msg = 'Unexpected istate=-1'
y0 = array([1.+0.j, 0.+0.j])
y1 = array([ 0.9999262 -7.37968809e-05j, -0.01214859+2.98852130e-07j])
______________ TestZVODECheckParameterUse.test_two_scalar_params _______________
scipy/integrate/tests/test_integrate.py:619: in test_two_scalar_params
self._check_solver(solver)
omega1 = 1.0
omega2 = 1.0
self = <scipy.integrate.tests.test_integrate.TestZVODECheckParameterUse object at 0x1684fa0a0>
solver = <scipy.integrate._ode.ode object at 0x1684faf10>
scipy/integrate/tests/test_integrate.py:597: in _check_solver
solver.integrate(pi)
ic = [1.0, 0.0]
self = <scipy.integrate.tests.test_integrate.TestZVODECheckParameterUse object at 0x1684fa0a0>
solver = <scipy.integrate._ode.ode object at 0x1684faf10>
scipy/integrate/_ode.py:433: in integrate
self._y, self.t = mth(self.f, self.jac or (lambda: None),
mth = <bound method vode.run of <scipy.integrate._ode.zvode object at 0x1684fae50>>
relax = False
self = <scipy.integrate._ode.ode object at 0x1684faf10>
step = False
t = 3.141592653589793
scipy/integrate/_ode.py:1013: in run
warnings.warn('{:s}: {:s}'.format(self.__class__.__name__,
E UserWarning: zvode: Excess work done on this call. (Perhaps wrong MF.)
args = (<function f2 at 0x11d625430>, <function jac2 at 0x11d6254c0>, array([1.+0.j, 0.+0.j]), 0.0, 3.141592653589793, 1e-07, ...)
f = <function f2 at 0x11d625430>
f_params = (1.0, 1.0)
istate = -1
jac = <function jac2 at 0x11d6254c0>
jac_params = (1.0, 1.0)
self = <scipy.integrate._ode.zvode object at 0x1684fae50>
t = 0.012148892003676552
t0 = 0.0
t1 = 3.141592653589793
unexpected_istate_msg = 'Unexpected istate=-1'
y0 = array([1.+0.j, 0.+0.j])
y1 = array([ 0.9999262 -7.37968809e-05j, -0.01214859+2.98852130e-07j])
_________________ TestZVODECheckParameterUse.test_vector_param _________________
scipy/integrate/tests/test_integrate.py:627: in test_vector_param
self._check_solver(solver)
omega = [1.0, 1.0]
self = <scipy.integrate.tests.test_integrate.TestZVODECheckParameterUse object at 0x1684c1df0>
solver = <scipy.integrate._ode.ode object at 0x1684c1b80>
scipy/integrate/tests/test_integrate.py:597: in _check_solver
solver.integrate(pi)
ic = [1.0, 0.0]
self = <scipy.integrate.tests.test_integrate.TestZVODECheckParameterUse object at 0x1684c1df0>
solver = <scipy.integrate._ode.ode object at 0x1684c1b80>
scipy/integrate/_ode.py:433: in integrate
self._y, self.t = mth(self.f, self.jac or (lambda: None),
mth = <bound method vode.run of <scipy.integrate._ode.zvode object at 0x1684c1ee0>>
relax = False
self = <scipy.integrate._ode.ode object at 0x1684c1b80>
step = False
t = 3.141592653589793
scipy/integrate/_ode.py:1013: in run
warnings.warn('{:s}: {:s}'.format(self.__class__.__name__,
E UserWarning: zvode: Excess work done on this call. (Perhaps wrong MF.)
args = (<function fv at 0x11d625550>, <function jacv at 0x11d6255e0>, array([1.+0.j, 0.+0.j]), 0.0, 3.141592653589793, 1e-07, ...)
f = <function fv at 0x11d625550>
f_params = ([1.0, 1.0],)
istate = -1
jac = <function jacv at 0x11d6255e0>
jac_params = ([1.0, 1.0],)
self = <scipy.integrate._ode.zvode object at 0x1684c1ee0>
t = 0.012148892003676552
t0 = 0.0
t1 = 3.141592653589793
unexpected_istate_msg = 'Unexpected istate=-1'
y0 = array([1.+0.j, 0.+0.j])
y1 = array([ 0.9999262 -7.37968809e-05j, -0.01214859+2.98852130e-07j])
_________________________ TestBarycentric.test_append __________________________
scipy/interpolate/tests/test_polyint.py:318: in test_append
P.add_xi(self.xs[3:], self.ys[3:])
P = <scipy.interpolate.polyint.BarycentricInterpolator object at 0x14acc24f0>
self = <scipy.interpolate.tests.test_polyint.TestBarycentric object at 0x14acbca60>
scipy/interpolate/polyint.py:617: in add_xi
self.wi **= -1
E RuntimeWarning: divide by zero encountered in reciprocal
old_n = 3
self = <scipy.interpolate.polyint.BarycentricInterpolator object at 0x14acc24f0>
xi = array([0.5, 1. ])
yi = array([[-1.],
[ 3.]])
____________________ TestSolve.test_simple_her_actuallysym _____________________
scipy/linalg/tests/test_basic.py:579: in test_simple_her_actuallysym
assert_array_almost_equal(dot(a, x), b)
E AssertionError:
E Arrays are not almost equal to 6 decimals
E
E Mismatched elements: 1 / 2 (50%)
E Max absolute difference: 1.03161697
E Max relative difference: 0.
E x: array([0. +1.j , 0.748905+0.709489j])
E y: array([0.+1.j, 0.+0.j])
a = [[2, 3], [3, -5]]
b = [1j, 0]
lower = 1
self = <scipy.linalg.tests.test_basic.TestSolve object at 0x14ae94100>
x = array([ 0.11824818+0.37518248j, -0.07883212+0.08321168j])
____________________ TestLstsq.test_simple_overdet_complex _____________________
scipy/linalg/tests/test_basic.py:1059: in test_simple_overdet_complex
assert_allclose(abs((dot(a, x) - b)**2).sum(axis=0),
E AssertionError:
E Not equal to tolerance rtol=2.98023e-06, atol=2.98023e-06
E driver: gelsd
E Mismatched elements: 1 / 1 (100%)
E Max absolute difference: 6.891783
E Max relative difference: 0.96745825
E x: array(14.015381, dtype=float32)
E y: array(7.123598, dtype=float32)
a = array([[1.+2.j, 2.+0.j],
[4.+0.j, 5.+0.j],
[3.+0.j, 4.+0.j]], dtype=complex64)
a1 = array([[1.+2.j, 2.+0.j],
[4.+0.j, 5.+0.j],
[3.+0.j, 4.+0.j]], dtype=complex64)
b = array([1.+0.j, 2.+4.j, 3.+0.j], dtype=complex64)
b1 = array([1.+0.j, 2.+4.j, 3.+0.j], dtype=complex64)
dtype = <class 'numpy.complex64'>
lapack_driver = 'gelsd'
out = (array([-0.48314652+0.02559727j, 0.94566095+0.08548044j], dtype=complex64), 7.1235976, 2, array([8.517418 , 1.5663979], dtype=float32))
overwrite = True
r = 2
residuals = 7.1235976
self = <scipy.linalg.tests.test_basic.TestLstsq object at 0x14ae821f0>
x = array([-0.48314652+0.02559727j, 0.94566095+0.08548044j], dtype=complex64)
_____________________ TestLstsq.test_random_complex_exact ______________________
scipy/linalg/tests/test_basic.py:1149: in test_random_complex_exact
assert_allclose(
E AssertionError:
E Not equal to tolerance rtol=4.76837e-05, atol=4.76837e-05
E driver: gelsd
E Mismatched elements: 60 / 60 (100%)
E Max absolute difference: 8.837243
E Max relative difference: 338.7324
E x: array([[ 1.023868-0.516378j, 0.93382 -0.630808j, 1.028059-0.395364j],
E [ 1.487115+0.568498j, 1.259903+0.727326j, 1.122598+0.464395j],
E [-2.461137+8.200006j, -2.912361+8.148932j, -2.881072+8.141109j],...
E y: array([[0.616139+0.j, 0.219291+0.j, 0.872524+0.j],
E [0.647865+0.j, 0.473003+0.j, 0.23419 +0.j],
E [0.833822+0.j, 0.098516+0.j, 0.205027+0.j],...
a = array([[5.8303890e+00+1.42998781e+01j, 6.2210876e-01+7.24091470e-01j,
4.3772775e-01+1.86764393e-02j, 7.8535861...505e-01+3.70635867e-01j,
2.2599553e-01+7.59584546e-01j, 8.0830746e+00+6.12367821e+00j]],
dtype=complex64)
a1 = array([[5.8303890e+00+1.42998781e+01j, 6.2210876e-01+7.24091470e-01j,
4.3772775e-01+1.86764393e-02j, 7.8535861...505e-01+3.70635867e-01j,
2.2599553e-01+7.59584546e-01j, 8.0830746e+00+6.12367821e+00j]],
dtype=complex64)
b = array([[0.61613923+0.j, 0.21929118+0.j, 0.87252367+0.j],
[0.64786494+0.j, 0.47300294+0.j, 0.23419032+0.j],
...52694 +0.j, 0.1732783 +0.j, 0.57995284+0.j],
[0.7273775 +0.j, 0.73040235+0.j, 0.4039958 +0.j]], dtype=complex64)
b1 = array([[0.61613923+0.j, 0.21929118+0.j, 0.87252367+0.j],
[0.64786494+0.j, 0.47300294+0.j, 0.23419032+0.j],
...52694 +0.j, 0.1732783 +0.j, 0.57995284+0.j],
[0.7273775 +0.j, 0.73040235+0.j, 0.4039958 +0.j]], dtype=complex64)
dtype = <class 'numpy.complex64'>
i = 0
lapack_driver = 'gelsd'
n = 20
out = (array([[ 0.01294096+0.01694867j, -0.00192093+0.01694867j,
0.01272957+0.01694867j],
[ 0.08152949-0.093...519 , 10.6915245, 8.947641 ,
8.566785 , 8.184578 , 7.360665 , 4.6518607, 3.4158041],
dtype=float32))
overwrite = True
r = 20
self = <scipy.linalg.tests.test_basic.TestLstsq object at 0x14ae82a30>
x = array([[ 0.01294096+0.01694867j, -0.00192093+0.01694867j,
0.01272957+0.01694867j],
[ 0.08152949-0.0939...181478j],
[ 0.02685548-0.3564865j , 0.01576587-0.37385136j,
0.00933838-0.35475665j]], dtype=complex64)
____________________ TestLstsq.test_random_complex_overdet _____________________
scipy/linalg/tests/test_basic.py:1211: in test_random_complex_overdet
assert_allclose(
E AssertionError:
E Not equal to tolerance rtol=2.98023e-06, atol=2.98023e-06
E driver: gelsd
E Mismatched elements: 45 / 45 (100%)
E Max absolute difference: 0.7172899
E Max relative difference: 149.43454
E x: array([[ 0.062104+6.638021e-01j, 0.057934+6.638021e-01j,
E 0.002216+6.638021e-01j],
E [-0.058405-3.299081e-02j, -0.058758-1.167682e-02j,...
E y: array([[ 0.062104-0.053488j, 0.057934-0.045925j, 0.002216+0.00382j ],
E [ 0.033628-0.025135j, -0.009612+0.005401j, 0.016862-0.020744j],
E [-0.004304+0.002065j, 0.0038 -0.002958j, 0.004777-0.00959j ],...
a = array([[5.8303890e+00+5.71755362e+00j, 6.2210876e-01+9.86747205e-01j,
4.3772775e-01+4.31800842e-01j, 7.8535861...8.2215977e-01+5.79493463e-01j, 6.2796509e-01+6.17560923e-01j,
1.1792306e-01+1.77117482e-01j]], dtype=complex64)
a1 = array([[5.8303890e+00+5.71755362e+00j, 6.2210876e-01+9.86747205e-01j,
4.3772775e-01+4.31800842e-01j, 7.8535861...8.2215977e-01+5.79493463e-01j, 6.2796509e-01+6.17560923e-01j,
1.1792306e-01+1.77117482e-01j]], dtype=complex64)
b = array([[0.8649818 +0.j, 0.7102707 +0.j, 0.17063169+0.j],
[0.85508263+0.j, 0.01985662+0.j, 0.64168376+0.j],
...500082+0.j, 0.02463633+0.j, 0.6935125 +0.j],
[0.42121193+0.j, 0.85834944+0.j, 0.57507354+0.j]], dtype=complex64)
b1 = array([[0.8649818 +0.j, 0.7102707 +0.j, 0.17063169+0.j],
[0.85508263+0.j, 0.01985662+0.j, 0.64168376+0.j],
...500082+0.j, 0.02463633+0.j, 0.6935125 +0.j],
[0.42121193+0.j, 0.85834944+0.j, 0.57507354+0.j]], dtype=complex64)
dtype = <class 'numpy.complex64'>
i = 0
lapack_driver = 'gelsd'
m = 15
n = 20
out = (array([[ 0.06210371+6.6380215e-01j, 0.05793425+6.6380215e-01j,
0.00221577+6.6380215e-01j],
[-0.05840...38003 , 15.058733, 13.859146, 11.219447, 9.0597 , 8.560572,
8.029268, 7.581321, 6.87882 ], dtype=float32))
overwrite = True
r = 15
self = <scipy.linalg.tests.test_basic.TestLstsq object at 0x149481f10>
x = array([[ 0.06210371+6.6380215e-01j, 0.05793425+6.6380215e-01j,
0.00221577+6.6380215e-01j],
[-0.058405...
[ 0.02510021-1.1387596e-01j, 0.0231653 -1.1449404e-01j,
0.02585518-9.7305194e-02j]], dtype=complex64)
________________________ TestFBLAS2Simple.test_syr_her _________________________
scipy/linalg/tests/test_blas.py:317: in test_syr_her
assert_allclose(f(1.0, z, lower=True), resz.T, rtol=rtol)
E AssertionError:
E Not equal to tolerance rtol=1e-07, atol=0
E
E Mismatched elements: 10 / 16 (62.5%)
E Max absolute difference: 53.15072906
E Max relative difference: 2.23606798
E x: array([[ -1. +3.j, 0. +0.j, 0. +0.j, 0. +0.j],
E [ -9. +13.j, -3. +21.j, 0. +0.j, 0. +0.j],
E [-25. +31.j, -15. +43.j, -5. +55.j, 0. +0.j],
E [-49. +57.j, -35. +73.j, -21. +89.j, -7.+105.j]], dtype=complex64)
E y: array([[ -3. +4.j, 0. +0.j, 0. +0.j, 0. +0.j],
E [ -5. +10.j, -7. +24.j, 0. +0.j, 0. +0.j],
E [ -7. +16.j, -9. +38.j, -11. +60.j, 0. +0.j],
E [ -9. +22.j, -11. +52.j, -13. +82.j, -15.+112.j]])
a = array([[ 1., 2., 3., 4.],
[ 0., 4., 6., 8.],
[ 0., 0., 9., 12.],
[ 0., 0., 0., 16.]])
b = array([[ 3., 6., 9., 12.],
[ 0., 12., 18., 24.],
[ 0., 0., 27., 36.],
[ 0., 0., 0., 48.]])
f = <fortran object>
p = 'c'
rehz = array([[ 5.+0.j, 11.+2.j, 17.+4.j, 23.+6.j],
[ 0.+0.j, 25.+0.j, 39.+2.j, 53.+4.j],
[ 0.+0.j, 0.+0.j, 61.+0.j, 83.+2.j],
[ 0.+0.j, 0.+0.j, 0.+0.j, 113.+0.j]])
rehz_reverse = array([[113.+0.j, 83.-2.j, 53.-4.j, 23.-6.j],
[ 0.+0.j, 61.+0.j, 39.-2.j, 17.-4.j],
[ 0.+0.j, 0.+0.j, 25.+0.j, 11.-2.j],
[ 0.+0.j, 0.+0.j, 0.+0.j, 5.+0.j]])
resx = array([[ 1., 2., 3., 4.],
[ 0., 4., 6., 8.],
[ 0., 0., 9., 12.],
[ 0., 0., 0., 16.]])
resx_reverse = array([[16., 12., 8., 4.],
[ 0., 9., 6., 3.],
[ 0., 0., 4., 2.],
[ 0., 0., 0., 1.]])
resz = array([[ -3. +4.j, -5. +10.j, -7. +16.j, -9. +22.j],
[ 0. +0.j, -7. +24.j, -9. +38.j, -11. +52.j],
[ 0. +0.j, 0. +0.j, -11. +60.j, -13. +82.j],
[ 0. +0.j, 0. +0.j, 0. +0.j, -15.+112.j]])
resz_reverse = array([[-15.+112.j, -13. +82.j, -11. +52.j, -9. +22.j],
[ 0. +0.j, -11. +60.j, -9. +38.j, -7. +16.j],
[ 0. +0.j, 0. +0.j, -7. +24.j, -5. +10.j],
[ 0. +0.j, 0. +0.j, 0. +0.j, -3. +4.j]])
rtol = 1e-07
self = <scipy.linalg.tests.test_blas.TestFBLAS2Simple object at 0x14949abb0>
w = array([0.+0.j, 1.+2.j, 0.+0.j, 0.+0.j, 3.+4.j, 0.+0.j, 0.+0.j, 5.+6.j,
0.+0.j, 0.+0.j, 7.+8.j, 0.+0.j])
x = array([1., 2., 3., 4.])
y = array([0. , 0.5, 1. , 1.5, 2. , 2.5, 3. , 3.5, 4. , 4.5, 5. , 5.5, 6. ,
6.5, 7. , 7.5, 8. ])
z = array([1.+2.j, 3.+4.j, 5.+6.j, 7.+8.j])
__________________________ TestEig.test_make_eigvals ___________________________
scipy/linalg/tests/test_decomp.py:351: in test_make_eigvals
self._check_gen_eig(A, B)
A = array([[0.36488598+0.56809865j, 0.61539618+0.86912739j,
0.07538124+0.43617342j],
[0.36882401+0.80214764... 0.65137814+0.70426097j],
[0.39720258+0.70458131j, 0.78873014+0.21879211j,
0.31683612+0.92486763j]])
B = array([[0.44214076+0.56143308j, 0.90931596+0.32966845j,
0.05980922+0.50296683j],
[0.18428708+0.11189432... 0.67488094+0.56594464j],
[0.59462478+0.00676406j, 0.53331016+0.61744171j,
0.04332406+0.91212289j]])
self = <scipy.linalg.tests.test_decomp.TestEig object at 0x149468bb0>
scipy/linalg/tests/test_decomp.py:235: in _check_gen_eig
assert_allclose(val1[:, i], val2[:, i],
E AssertionError:
E Not equal to tolerance rtol=1e-13, atol=1e-13
E
E array([[0.36488598+0.56809865j, 0.61539618+0.86912739j,
E 0.07538124+0.43617342j],
E [0.36882401+0.80214764j, 0.9331401 +0.14376682j,
E 0.65137814+0.70426097j],
E [0.39720258+0.70458131j, 0.78873014+0.21879211j,
E 0.31683612+0.92486763j]])
E array([[0.44214076+0.56143308j, 0.90931596+0.32966845j,
E 0.05980922+0.50296683j],
E [0.18428708+0.11189432j, 0.04735528+0.60719371j,
E 0.67488094+0.56594464j],
E [0.59462478+0.00676406j, 0.53331016+0.61744171j,
E 0.04332406+0.91212289j]])
E Mismatched elements: 3 / 3 (100%)
E Max absolute difference: 0.90239291
E Max relative difference: 0.30065
E x: array([ 0.666376-2.307768j, -0.150503-2.942555j, 0.3752 -2.735541j])
E y: array([0.087702-3.000192j, 0.331766-2.40094j , 0.403462-2.903364j])
A = array([[0.36488598+0.56809865j, 0.61539618+0.86912739j,
0.07538124+0.43617342j],
[0.36882401+0.80214764... 0.65137814+0.70426097j],
[0.39720258+0.70458131j, 0.78873014+0.21879211j,
0.31683612+0.92486763j]])
B = array([[0.44214076+0.56143308j, 0.90931596+0.32966845j,
0.05980922+0.50296683j],
[0.18428708+0.11189432... 0.67488094+0.56594464j],
[0.59462478+0.00676406j, 0.53331016+0.61744171j,
0.04332406+0.91212289j]])
B0 = array([[0.44214076+0.56143308j, 0.90931596+0.32966845j,
0.05980922+0.50296683j],
[0.18428708+0.11189432... 0.67488094+0.56594464j],
[0.59462478+0.00676406j, 0.53331016+0.61744171j,
0.04332406+0.91212289j]])
i = 0
msg = '\narray([[0.36488598+0.56809865j, 0.61539618+0.86912739j,\n 0.07538124+0.43617342j],\n [0.36882401+0.802... 0.67488094+0.56594464j],\n [0.59462478+0.00676406j, 0.53331016+0.61744171j,\n 0.04332406+0.91212289j]])'
self = <scipy.linalg.tests.test_decomp.TestEig object at 0x149468bb0>
val1 = array([[ 0.66637578-2.30776823j, -0.01908767+0.00947095j,
0.09765692-0.15772331j],
[-0.15050348-2.9425... -0.55883967-0.15461035j],
[ 0.37520026-2.735541j , 0.14129369-0.05101752j,
-0.1945099 -0.13107334j]])
val2 = array([[ 8.77023558e-02-3.00019153e+00j, 4.65650672e-03-2.64720113e-03j,
3.87012365e-04-2.61067454e-04j],
...],
[ 4.03462148e-01-2.90336431e+00j, -1.57577078e-02+5.65717628e-03j,
-5.46545800e-04+2.63856650e-03j]])
vr = array([[-0.43040566-0.43040566j, -0.59839128-0.59839128j,
0.54544764+0.54544764j],
[-0.36063072-0.3606... -0.32737811-0.32737811j],
[-0.42976326-0.42976326j, 0.20434679+0.20434679j,
-0.3087239 -0.3087239j ]])
w = array([[ 2.67726000e+00+0.00000000e+00j, 2.83445479e-02-1.01134379e-15j,
-2.49703427e-03+2.88834243e-03j],
[ 1.99500344e+00+0.00000000e+00j, 5.88559241e-01+0.00000000e+00j,
1.16801448e+00+0.00000000e+00j]])
wt = array([[ 2.67726000e+00+0.00000000e+00j, 2.83445479e-02-1.01134379e-15j,
-2.49703427e-03+2.88834243e-03j],
[ 1.99500344e+00+0.00000000e+00j, 5.88559241e-01+0.00000000e+00j,
1.16801448e+00+0.00000000e+00j]])
__________________________ TestEigBanded.test_zhbevd ___________________________
scipy/linalg/tests/test_decomp.py:502: in test_zhbevd
assert_array_almost_equal(abs(evec_), abs(self.evec_herm_lin))
E AssertionError:
E Arrays are not almost equal to 6 decimals
E
E Mismatched elements: 100 / 100 (100%)
E Max absolute difference: 1.58968705
E Max relative difference: 6.82730566
E x: array([[0.359877, 0.031817, 0.638579, 0.030717, 0.772035, 0.227172,
E 0.696111, 0.244334, 0.462393, 1.014976],
E [0.159978, 0.440427, 0.222026, 0.754656, 0.076127, 0.741219,...
E y: array([[0.254471, 0.022498, 0.451543, 0.02172 , 0.545911, 0.160635,
E 0.492225, 0.17277 , 0.326962, 0.165873],
E [0.113121, 0.311429, 0.156995, 0.533622, 0.05383 , 0.524121,...
evec = array([[-0.25447137-0.2544714j , -0.0224981 -0.0224981j ,
-0.45154314-0.45154325j, 0.0217204 +0.0217204j ,
...,
-0.4922247 +0.4922247j , 0.17277077-0.17276965j,
-0.32696153+0.32696153j, -0.14312147+0.16587337j]])
evec_ = array([[-0.25447137-0.2544714j , -0.0224981 -0.0224981j ,
-0.45154314-0.45154325j, 0.0217204 +0.0217204j ,
...,
-0.4922247 +0.4922247j , 0.17277077-0.17276965j,
-0.32696153+0.32696153j, -0.14312147+0.16587337j]])
info = 0
self = <scipy.linalg.tests.test_decomp.TestEigBanded object at 0x14ae13af0>
w = array([-4.6749823 , -4.49981592, -3.43214201, -2.84350424, -1.66090423,
-0.86217766, 0.43023792, 0.52994106, 2.64944971, 4.36389767])
__________________________ TestEigBanded.test_zhbevx ___________________________
scipy/linalg/tests/test_decomp.py:513: in test_zhbevx
assert_array_almost_equal(abs(evec_), abs(self.evec_herm_lin))
E AssertionError:
E Arrays are not almost equal to 6 decimals
E
E Mismatched elements: 100 / 100 (100%)
E Max absolute difference: 0.84910218
E Max relative difference: 5.11897831
E x: array([[3.598769e-01, 3.181712e-02, 6.385785e-01, 3.071729e-02,
E 7.720354e-01, 2.271717e-01, 6.961108e-01, 2.443340e-01,
E 4.623934e-01, 1.014976e+00],...
E y: array([[0.254471, 0.022498, 0.451543, 0.02172 , 0.545911, 0.160635,
E 0.492225, 0.17277 , 0.326962, 0.165873],
E [0.113121, 0.311429, 0.156995, 0.533622, 0.05383 , 0.524121,...
N = 10
evec = array([[-0.25447137-2.54471398e-01j, -0.0224981 -2.24981016e-02j,
-0.45154314-4.51543251e-01j, 0.0217204 +2.1... +2.56791415e-30j, 0. +1.18038862e-11j,
0. +0.00000000e+00j, 0. +9.99983105e-01j]])
evec_ = array([[-0.25447137-2.54471398e-01j, -0.0224981 -2.24981016e-02j,
-0.45154314-4.51543251e-01j, 0.0217204 +2.1... +2.56791415e-30j, 0. +1.18038862e-11j,
0. +0.00000000e+00j, 0. +9.99983105e-01j]])
ifail = array([0, 0, 0, 0, 0, 0, 0, 0, 0, 0], dtype=int32)
info = 0
num = 10
self = <scipy.linalg.tests.test_decomp.TestEigBanded object at 0x14ae13d90>
w = array([-4.6749823 , -4.49981592, -3.43214201, -2.84350424, -1.66090423,
-0.86217766, 0.43023792, 0.52994106, 2.64944971, 4.36389767])
________________________ TestEigBanded.test_eig_banded _________________________
scipy/linalg/tests/test_decomp.py:568: in test_eig_banded
assert_array_almost_equal(abs(evec_herm_), abs(self.evec_herm_lin))
E AssertionError:
E Arrays are not almost equal to 6 decimals
E
E Mismatched elements: 100 / 100 (100%)
E Max absolute difference: 1.58968705
E Max relative difference: 6.82730566
E x: array([[0.359877, 0.031817, 0.638579, 0.030717, 0.772035, 0.227172,
E 0.696111, 0.244334, 0.462393, 1.014976],
E [0.159978, 0.440427, 0.222026, 0.754656, 0.076127, 0.741219,...
E y: array([[0.254471, 0.022498, 0.451543, 0.02172 , 0.545911, 0.160635,
E 0.492225, 0.17277 , 0.326962, 0.165873],
E [0.113121, 0.311429, 0.156995, 0.533622, 0.05383 , 0.524121,...
evec_herm = array([[-0.25447137-0.2544714j , -0.0224981 -0.0224981j ,
-0.45154314-0.45154325j, 0.0217204 +0.0217204j ,
...,
-0.4922247 +0.4922247j , 0.17277077-0.17276965j,
-0.32696153+0.32696153j, -0.14312147+0.16587337j]])
evec_herm_ = array([[-0.25447137-0.2544714j , -0.0224981 -0.0224981j ,
-0.45154314-0.45154325j, 0.0217204 +0.0217204j ,
...,
-0.4922247 +0.4922247j , 0.17277077-0.17276965j,
-0.32696153+0.32696153j, -0.14312147+0.16587337j]])
evec_sym = array([[ 0.16587337, 0.32696153, 0.17276965, -0.4922247 , 0.16063466,
0.54591149, 0.0217204 , 0.45154314...6153, -0.17276965, -0.4922247 , 0.16063466,
-0.54591149, -0.0217204 , 0.45154314, -0.0224981 , -0.25447137]])
evec_sym_ = array([[ 0.16587337, 0.32696153, 0.17276965, -0.4922247 , 0.16063466,
0.54591149, 0.0217204 , 0.45154314...6153, -0.17276965, -0.4922247 , 0.16063466,
-0.54591149, -0.0217204 , 0.45154314, -0.0224981 , -0.25447137]])
self = <scipy.linalg.tests.test_decomp.TestEigBanded object at 0x149453dc0>
w_herm = array([-4.6749823 , -4.49981592, -3.43214201, -2.84350424, -1.66090423,
-0.86217766, 0.43023792, 0.52994106, 2.64944971, 4.36389767])
w_sym = array([-4.36389767, -2.64944971, -0.52994106, -0.43023792, 0.86217766,
1.66090423, 2.84350424, 3.43214201, 4.49981592, 4.6749823 ])
_____________ TestEigh.test_various_drivers_standard[ev-complex64] _____________
scipy/linalg/tests/test_decomp.py:851: in test_various_drivers_standard
assert_allclose(a @ v - (v * w), 0.,
E AssertionError:
E Not equal to tolerance rtol=0, atol=0.000119209
E
E Mismatched elements: 400 / 400 (100%)
E Max absolute difference: 9.878377
E Max relative difference: inf
E x: array([[-1.024111e-01-1.129053e-01j, -9.103125e-02-3.533325e-01j,
E 1.826477e-01-3.730512e-02j, 2.729376e-01+2.975390e-01j,
E -1.050508e-01+2.498122e-01j, 3.369594e-01-3.583230e-02j,...
E y: array(0.)
a = array([[0.6445826 +0.00000000e+00j, 0.32437274+9.31576341e-02j,
0.61983734+6.31575584e-02j, 0.32045603+1.76963...j, 0.41159117-4.80479747e-02j,
0.68176913+2.78819323e-01j, 0.0289948 +0.00000000e+00j]],
dtype=complex64)
driver = 'ev'
dtype_ = <class 'numpy.complex64'>
self = <scipy.linalg.tests.test_decomp.TestEigh object at 0x14942eee0>
v = array([[-0.0739059 -7.39059001e-02j, -0.30122298-3.01222980e-01j,
0.14191413+1.41914129e-01j, 0.17931043+1.7...-0.05785928-5.78592792e-02j,
0.01692364+9.84401166e-01j, 0.21668212+2.16682121e-01j]],
dtype=complex64)
w = array([-2.3175747 , -1.6421384 , -1.5190917 , -1.3529652 , -1.0369561 ,
-0.91279966, -0.8082594 , -0.59898156, ... 0.9107243 , 1.1070666 ,
1.2276475 , 1.6616762 , 2.0099874 , 2.3414361 , 10.686417 ],
dtype=float32)
____________ TestEigh.test_various_drivers_standard[ev-complex128] _____________
scipy/linalg/tests/test_decomp.py:851: in test_various_drivers_standard
assert_allclose(a @ v - (v * w), 0.,
E AssertionError:
E Not equal to tolerance rtol=0, atol=2.22045e-13
E
E Mismatched elements: 400 / 400 (100%)
E Max absolute difference: 4.96917367
E Max relative difference: inf
E x: array([[-2.983289e-02-2.611873e-02j, -2.292372e-01-1.617507e-01j,
E 2.889755e-01+2.519307e-01j, 6.053284e-02+5.959726e-01j,
E 4.928402e-02-3.507084e-02j, -8.426184e-02-3.250729e-01j,...
E y: array(0.)
a = array([[0.37492896+0.j , 0.45935991-0.18855491j,
0.82462976-0.20893167j, 0.20393551+0.34284233j,
...255j,
0.56518081-0.23093414j, 0.29804915+0.15078988j,
0.40743044-0.35303982j, 0.15557592+0.j ]])
driver = 'ev'
dtype_ = <class 'numpy.complex128'>
self = <scipy.linalg.tests.test_decomp.TestEigh object at 0x14ae38df0>
v = array([[ 9.12043145e-02+9.12043145e-02j, -2.28230214e-01-2.28230214e-01j,
3.26174267e-01+3.26174267e-01j, 2....0736e-01j, -2.91865337e-02-2.91865732e-02j,
-6.23002395e-02+6.33701440e-01j, 1.82498824e-01+1.82498824e-01j]])
w = array([-2.32053699, -1.92537705, -1.69232193, -1.33074409, -1.0528957 ,
-0.76211562, -0.54819113, -0.24444331, ...028108, 0.68326401, 0.79165356, 0.93821424,
1.32300894, 1.51791058, 1.76330908, 2.3285108 , 10.21982614])
____________ TestEigh.test_various_drivers_standard[evd-complex64] _____________
scipy/linalg/tests/test_decomp.py:851: in test_various_drivers_standard
assert_allclose(a @ v - (v * w), 0.,
E AssertionError:
E Not equal to tolerance rtol=0, atol=0.000119209
E
E Mismatched elements: 20 / 400 (5%)
E Max absolute difference: 0.17769675
E Max relative difference: inf
E x: array([[-1.341105e-07+1.483355e-02j, 3.576279e-07+2.980232e-08j,
E 2.980232e-07+2.980232e-07j, 2.980232e-08+8.940697e-08j,
E -7.450581e-08-3.576279e-07j, -8.195639e-08+2.980232e-08j,...
E y: array(0.)
a = array([[0.93991697+0.00000000e+00j, 0.39091605+1.65591612e-01j,
0.26158276-4.73397635e-02j, 0.37927064-1.44070...j, 0.16194409+1.27978563e-01j,
0.4030868 -2.19227582e-01j, 0.38564578+0.00000000e+00j]],
dtype=complex64)
driver = 'evd'
dtype_ = <class 'numpy.complex64'>
self = <scipy.linalg.tests.test_decomp.TestEigh object at 0x149436fd0>
v = array([[-1.07141308e-01-0.10309255j, -1.86949298e-01-0.1869493j ,
-3.60869735e-01-0.36086974j, -1.99490339e-01... 6.41181841e-02-0.40363362j,
-3.97245467e-01+0.00978629j, 1.36577025e-01+0.17734522j]],
dtype=complex64)
w = array([-2.2366235e+00, -2.1211147e+00, -1.6599302e+00, -1.3560791e+00,
-1.0950742e+00, -8.0291504e-01, -5.99710...21013e+00, 1.2699474e+00,
1.6432656e+00, 1.9971005e+00, 2.2034056e+00, 9.8387175e+00],
dtype=float32)
____________ TestEigh.test_various_drivers_standard[evd-complex128] ____________
scipy/linalg/tests/test_decomp.py:851: in test_various_drivers_standard
assert_allclose(a @ v - (v * w), 0.,
E AssertionError:
E Not equal to tolerance rtol=0, atol=2.22045e-13
E
E Mismatched elements: 20 / 400 (5%)
E Max absolute difference: 3.91507493
E Max relative difference: inf
E x: array([[ 5.551115e-16+1.072247e-01j, -2.775558e-16-5.551115e-16j,
E 3.053113e-16+4.996004e-16j, -1.110223e-16-2.775558e-16j,
E -2.775558e-16+5.551115e-17j, -1.804112e-16-2.914335e-16j,...
E y: array(0.)
a = array([[0.80326946+0.j , 0.59936528+0.05300101j,
0.95832685+0.0462222j , 0.55680251+0.17499864j,
...294j,
0.60840317-0.22914422j, 0.59827221-0.00135214j,
0.26349336-0.23383477j, 0.60505255+0.j ]])
driver = 'evd'
dtype_ = <class 'numpy.complex128'>
self = <scipy.linalg.tests.test_decomp.TestEigh object at 0x149436160>
v = array([[-7.20684005e-02+9.15398123e-01j, 1.51473516e-01+1.51473516e-01j,
-8.83182172e-02-8.83182172e-02j, 2....5449e-01j, -3.57786128e-01-2.07193101e-01j,
1.11792605e-01+1.76430821e-01j, 2.37966225e-01+2.04786508e-01j]])
w = array([-2.1933784 , -1.94651785, -1.5997267 , -1.57196017, -0.97569209,
-0.88542367, -0.71470547, -0.3696185 , ...027221, 0.58746605, 0.83714941, 1.03364908,
1.26355445, 1.46203797, 1.7543098 , 2.06181122, 10.41503829])
____________ TestEigh.test_various_drivers_standard[evx-complex64] _____________
scipy/linalg/tests/test_decomp.py:851: in test_various_drivers_standard
assert_allclose(a @ v - (v * w), 0.,
E AssertionError:
E Not equal to tolerance rtol=0, atol=0.000119209
E
E Mismatched elements: 400 / 400 (100%)
E Max absolute difference: 13.744765
E Max relative difference: inf
E x: array([[-3.651766e-01+1.175622e+01j, -6.825393e-02-5.268404e-01j,
E 3.930492e-01+1.110123e-01j, -4.181564e-04-3.892285e-02j,
E 1.544308e-01+3.418234e-02j, -2.727565e-01-1.666747e-01j,...
E y: array(0.)
a = array([[0.9203447 +0.j , 0.2483745 -0.03459179j,
0.6773846 +0.12101765j, 0.5054712 -0.07335676j,
...704623 +0.23423882j, 0.63310504+0.28264895j,
0.38397732+0.3441984j , 0.07380515+0.j ]], dtype=complex64)
driver = 'evx'
dtype_ = <class 'numpy.complex64'>
self = <scipy.linalg.tests.test_decomp.TestEigh object at 0x14945c7f0>
v = array([[-0.0429727 +1.00964117e+00j, -0.14481622-1.44816220e-01j,
0.26343495+2.63434947e-01j, -0.00761492-7.6... 0.04240745+4.24074531e-02j,
-0.17201637-1.72016367e-01j, 0.2040696 +2.04069585e-01j]],
dtype=complex64)
w = array([-2.2410452 , -1.9704098 , -1.5324193 , -1.3974583 , -1.078885 ,
-0.72428995, -0.64166534, -0.3583088 , ... 0.9276952 , 1.150788 ,
1.3305326 , 1.5351841 , 1.7345213 , 2.1826546 , 10.071262 ],
dtype=float32)
____________ TestEigh.test_various_drivers_standard[evx-complex128] ____________
scipy/linalg/tests/test_decomp.py:851: in test_various_drivers_standard
assert_allclose(a @ v - (v * w), 0.,
E AssertionError:
E Not equal to tolerance rtol=0, atol=2.22045e-13
E
E Mismatched elements: 400 / 400 (100%)
E Max absolute difference: 7.80396382
E Max relative difference: inf
E x: array([[ 3.311647e-01+4.732797e-01j, -4.936677e-01-4.542667e-01j,
E -3.088349e-01-9.947181e-02j, -2.556849e-02+1.133738e-02j,
E 2.515436e-03-5.583135e-02j, -5.922909e-02-4.024031e-01j,...
E y: array(0.)
a = array([[0.26914684+0.00000000e+00j, 0.29146127-4.10195053e-01j,
0.95705696-3.48075063e-01j, 0.87025858+1.15636...33202859-3.21654094e-01j, 0.59800323+1.94424008e-01j,
0.55331016+4.14019946e-01j, 0.25478618+0.00000000e+00j]])
driver = 'evx'
dtype_ = <class 'numpy.complex128'>
self = <scipy.linalg.tests.test_decomp.TestEigh object at 0x14942fdf0>
v = array([[ 3.30891400e-01+3.30891400e-01j, -4.17037287e-01-4.17037287e-01j,
-1.07798794e-01-1.07798794e-01j, -1....3780e-02j, -1.89615754e-01-1.89615754e-01j,
-4.41563754e-02+8.01356524e-01j, 1.85068789e-01+1.85068789e-01j]])
w = array([-2.27173441e+00, -2.03989104e+00, -1.62632289e+00, -1.41799595e+00,
-1.07636252e+00, -8.08232864e-01, -6...9124e-01, 1.21517278e+00, 1.39670038e+00,
1.80985410e+00, 1.96426791e+00, 2.36322476e+00, 1.00990438e+01])
_________________ TestEigh.test_eigh[6-F-True-True-True-None] __________________
scipy/linalg/tests/test_decomp.py:892: in test_eigh
assert_allclose(diag1_, w, rtol=0., atol=atol)
E AssertionError:
E Not equal to tolerance rtol=0, atol=0.0001
E
E Mismatched elements: 6 / 6 (100%)
E Max absolute difference: 0.46809623
E Max relative difference: 1.0000427
E x: array([-0.485601, -0.312531, -0.074911, 0.219079, 0.5228 , 0.936172],
E dtype=float32)
E y: array([-0.2428 +0.j, -0.156266+0.j, -0.037455+0.j, 0.109539+0.j,
E 0.2614 +0.j, 0.468076+0.j], dtype=complex64)
a = array([[0.94189006+0.j, 0.3177541 +0.j, 0.48921168+0.j, 0.29285413+0.j,
0.62844056+0.j, 0.11751232+0.j],
...1751232+0.j, 0.47093302+0.j, 0.41464418+0.j, 0.5263543 +0.j,
0.6622156 +0.j, 0.56324327+0.j]], dtype=complex64)
atol = 0.0001
b = array([[3.8814993 +0.j, 0.7532733 +0.j, 0.5667193 +0.j, 0.91612804+0.j,
0.18487613+0.j, 0.43611383+0.j],
...3611383+0.j, 0.39468226+0.j, 0.4902285 +0.j, 0.3416944 +0.j,
0.246093 +0.j, 4.1263437 +0.j]], dtype=complex64)
diag1_ = array([-0.48560095, -0.31253147, -0.0749108 , 0.21907865, 0.5228004 ,
0.9361725 ], dtype=float32)
diag_ = array([-0.60781467, -0.1966563 , 0.19238439, 0.44767025, 0.61062247,
2.7847323 ], dtype=float32)
dim = 6
dtype_ = 'F'
eigvals = None
lower = True
overwrite = True
self = <scipy.linalg.tests.test_decomp.TestEigh object at 0x14942d3a0>
turbo = True
w = array([-0.24280049+0.j, -0.15626569+0.j, -0.03745542+0.j, 0.10953932+0.j,
0.2614001 +0.j, 0.46807626+0.j], dtype=complex64)
z = array([[ 0.02692915+0.02692915j, -0.021659 -0.021659j ,
0.19803704+0.19803704j, 0.1290862 +0.1290862j ,
....33830327j, -0.27845508-0.27845508j,
-0.15060899-0.15060899j, -0.1726871 -0.1711923j ]],
dtype=complex64)
_________________ TestEigh.test_eigh[6-F-True-True-False-None] _________________
scipy/linalg/tests/test_decomp.py:892: in test_eigh
assert_allclose(diag1_, w, rtol=0., atol=atol)
E AssertionError:
E Not equal to tolerance rtol=0, atol=0.0001
E
E Mismatched elements: 6 / 6 (100%)
E Max absolute difference: 0.32953268
E Max relative difference: 1.0000012
E x: array([-0.541171, -0.24827 , -0.16693 , 0.01691 , 0.242608, 0.149311],
E dtype=float32)
E y: array([-0.270585+0.j, -0.124142+0.j, -0.083465+0.j, 0.008455+0.j,
E 0.121304+0.j, 0.478843+0.j], dtype=complex64)
a = array([[0.24046203+0.j, 0.7203457 +0.j, 0.51699066+0.j, 0.72209704+0.j,
0.46539277+0.j, 0.29535893+0.j],
...9535893+0.j, 0.8227661 +0.j, 0.7007535 +0.j, 0.4971041 +0.j,
0.5259398 +0.j, 0.06233544+0.j]], dtype=complex64)
atol = 0.0001
b = array([[4.3692007 +0.j, 0.5208768 +0.j, 0.4924582 +0.j, 0.6876096 +0.j,
0.7258201 +0.j, 0.33202532+0.j],
...3202532+0.j, 0.48081577+0.j, 0.39746702+0.j, 0.4460055 +0.j,
0.4292883 +0.j, 4.4085026 +0.j]], dtype=complex64)
diag1_ = array([-0.54117084, -0.2482701 , -0.16692963, 0.01691022, 0.24260822,
0.14931066], dtype=float32)
diag_ = array([-0.44902846, -0.16516139, 0.14273405, 0.51106656, 0.8109742 ,
3.2629323 ], dtype=float32)
dim = 6
dtype_ = 'F'
eigvals = None
lower = True
overwrite = True
self = <scipy.linalg.tests.test_decomp.TestEigh object at 0x14942df40>
turbo = False
w = array([-0.27058533+0.j, -0.12414176+0.j, -0.08346476+0.j, 0.00845512+0.j,
0.12130409+0.j, 0.47884333+0.j], dtype=complex64)
z = array([[ 0.23658736+0.23658854j, -0.05110337-0.04964942j,
0.2828052 +0.28268144j, -0.27289012-0.2729052j ,
....13521922j, 0.20471877+0.2047004j ,
0.02705213+0.02705213j, 0.15438505-0.02807979j]],
dtype=complex64)
_________________ TestEigh.test_eigh[6-F-True-False-True-None] _________________
scipy/linalg/tests/test_decomp.py:892: in test_eigh
assert_allclose(diag1_, w, rtol=0., atol=atol)
E AssertionError:
E Not equal to tolerance rtol=0, atol=0.0001
E
E Mismatched elements: 6 / 6 (100%)
E Max absolute difference: 0.54828775
E Max relative difference: 1.1236883
E x: array([-0.478913, -0.163269, -0.036441, 0.180954, 0.371176, 1.036224],
E dtype=float32)
E y: array([-0.239456+0.j, -0.081634+0.j, -0.018221+0.j, 0.090477+0.j,
E 0.185588+0.j, 0.487936+0.j], dtype=complex64)
a = array([[0.14704211+0.j, 0.89636266+0.j, 0.6964647 +0.j, 0.45620486+0.j,
0.08488958+0.j, 0.5465375 +0.j],
...465375 +0.j, 0.4186488 +0.j, 0.64523673+0.j, 0.3573673 +0.j,
0.44854876+0.j, 0.13988845+0.j]], dtype=complex64)
atol = 0.0001
b = array([[3.7057762 +0.j, 0.5976716 +0.j, 0.7405352 +0.j, 0.38507336+0.j,
0.18701029+0.j, 0.75921696+0.j],
...5921696+0.j, 0.07421535+0.j, 0.699367 +0.j, 0.57520294+0.j,
0.1353806 +0.j, 4.269169 +0.j]], dtype=complex64)
diag1_ = array([-0.4789126 , -0.16326874, -0.03644104, 0.1809543 , 0.3711757 ,
1.0362235 ], dtype=float32)
diag_ = array([-0.69455516, -0.35041702, -0.06705053, 0.13802083, 0.39961326,
2.6240501 ], dtype=float32)
dim = 6
dtype_ = 'F'
eigvals = None
lower = False
overwrite = True
self = <scipy.linalg.tests.test_decomp.TestEigh object at 0x14940e460>
turbo = True
w = array([-0.23945613+0.j, -0.08163433+0.j, -0.01822052+0.j, 0.09047712+0.j,
0.18558776+0.j, 0.48793578+0.j], dtype=complex64)
z = array([[ 0.452138 +0.452138j , 0.13338469+0.13338473j,
0.15458602+0.15457486j, -0.01659884-0.01659893j,
....09108649j, -0.15876791-0.15876791j,
0.08394218+0.08394222j, 0.11395784+0.49494687j]],
dtype=complex64)
________________ TestEigh.test_eigh[6-F-True-False-False-None] _________________
scipy/linalg/tests/test_decomp.py:892: in test_eigh
assert_allclose(diag1_, w, rtol=0., atol=atol)
E AssertionError:
E Not equal to tolerance rtol=0, atol=0.0001
E
E Mismatched elements: 6 / 6 (100%)
E Max absolute difference: 0.27752507
E Max relative difference: 1.0000007
E x: array([-0.55505 , -0.213289, 0.045369, 0.183871, 0.273931, 0.396452],
E dtype=float32)
E y: array([-0.277525+0.j, -0.106645+0.j, 0.022685+0.j, 0.091935+0.j,
E 0.136965+0.j, 0.453675+0.j], dtype=complex64)
a = array([[0.4506064 +0.j, 0.38382778+0.j, 0.81977063+0.j, 0.26001203+0.j,
0.53073037+0.j, 0.19195576+0.j],
...9195576+0.j, 0.59642404+0.j, 0.4405579 +0.j, 0.5363573 +0.j,
0.48572755+0.j, 0.21342213+0.j]], dtype=complex64)
atol = 0.0001
b = array([[3.5921583 +0.j, 0.14073387+0.j, 0.620896 +0.j, 0.46062562+0.j,
0.2517855 +0.j, 0.49150977+0.j],
...9150977+0.j, 0.5499142 +0.j, 0.54769075+0.j, 0.5205428 +0.j,
0.3506582 +0.j, 4.38437 +0.j]], dtype=complex64)
diag1_ = array([-0.55505013, -0.21328928, 0.04536902, 0.1838706 , 0.2739307 ,
0.396452 ], dtype=float32)
diag_ = array([-0.856958 , -0.4250862 , -0.15582177, -0.01355576, 0.23185635,
2.9385388 ], dtype=float32)
dim = 6
dtype_ = 'F'
eigvals = None
lower = False
overwrite = True
self = <scipy.linalg.tests.test_decomp.TestEigh object at 0x14940ea90>
turbo = False
w = array([-0.27752507+0.j, -0.1066446 +0.j, 0.02268453+0.j, 0.09193532+0.j,
0.13696538+0.j, 0.45367455+0.j], dtype=complex64)
z = array([[-0.24175027-0.24175027j, 0.13737267+0.13737267j,
-0.19659264-0.19659264j, 0.35927033+0.35927033j,
....01595859j, -0.26116204-0.26116204j,
-0.08484881-0.08484881j, -0.10266025+0.35609558j]],
dtype=complex64)
_________________ TestEigh.test_eigh[6-F-False-True-True-None] _________________
scipy/linalg/tests/test_decomp.py:892: in test_eigh
assert_allclose(diag1_, w, rtol=0., atol=atol)
E AssertionError:
E Not equal to tolerance rtol=0, atol=0.0001
E
E Mismatched elements: 6 / 6 (100%)
E Max absolute difference: 0.49825108
E Max relative difference: 1.0000006
E x: array([-0.44047 , -0.188464, 0.038233, 0.278568, 0.352137, 0.996502],
E dtype=float32)
E y: array([-0.220235+0.j, -0.094232+0.j, 0.019117+0.j, 0.139284+0.j,
E 0.176069+0.j, 0.498251+0.j], dtype=complex64)
a = array([[0.882543 +0.j, 0.60269994+0.j, 0.3709318 +0.j, 0.5176232 +0.j,
0.5265752 +0.j, 0.3746794 +0.j],
...746794 +0.j, 0.68633467+0.j, 0.63627464+0.j, 0.4105417 +0.j,
0.69689256+0.j, 0.76453674+0.j]], dtype=complex64)
atol = 0.0001
b = array([[4.4315524 +0.j, 0.41940913+0.j, 0.14003824+0.j, 0.5052338 +0.j,
0.65201926+0.j, 0.96257436+0.j],
...6257436+0.j, 0.9009159 +0.j, 0.21293898+0.j, 0.35946187+0.j,
0.5895454 +0.j, 3.6991673 +0.j]], dtype=complex64)
diag1_ = array([-0.44046974, -0.18846393, 0.03823304, 0.2785678 , 0.35213748,
0.9965023 ], dtype=float32)
diag_ = array([-0.49769592, -0.42596987, -0.07926966, 0.35028592, 0.5796132 ,
2.511098 ], dtype=float32)
dim = 6
dtype_ = 'F'
eigvals = None
lower = True
overwrite = False
self = <scipy.linalg.tests.test_decomp.TestEigh object at 0x14ac61dc0>
turbo = True
w = array([-0.22023499+0.j, -0.09423202+0.j, 0.01911653+0.j, 0.13928385+0.j,
0.17606875+0.j, 0.4982512 +0.j], dtype=complex64)
z = array([[ 0.02304962+0.02304962j, 0.19031185+0.19031185j,
0.00609896+0.00609896j, 0.26555067+0.26555067j,
....38324526j, -0.13731356-0.13731356j,
-0.23705967-0.23705967j, -0.20846404-0.2084637j ]],
dtype=complex64)
________________ TestEigh.test_eigh[6-F-False-True-False-None] _________________
scipy/linalg/tests/test_decomp.py:892: in test_eigh
assert_allclose(diag1_, w, rtol=0., atol=atol)
E AssertionError:
E Not equal to tolerance rtol=0, atol=0.0001
E
E Mismatched elements: 6 / 6 (100%)
E Max absolute difference: 0.22122192
E Max relative difference: 1.0000068
E x: array([-0.442444, -0.201545, 0.017607, 0.16953 , 0.357098, 0.721418],
E dtype=float32)
E y: array([-0.221222+0.j, -0.100773+0.j, 0.008804+0.j, 0.084765+0.j,
E 0.178549+0.j, 0.525691+0.j], dtype=complex64)
a = array([[0.22101411+0.j, 0.33586907+0.j, 0.63221896+0.j, 0.53261286+0.j,
0.24382536+0.j, 0.35554168+0.j],
...5554168+0.j, 0.50479794+0.j, 0.7114129 +0.j, 0.64866424+0.j,
0.73286283+0.j, 0.47943765+0.j]], dtype=complex64)
atol = 0.0001
b = array([[3.853972 +0.j, 0.6447499 +0.j, 0.48021883+0.j, 0.4430055 +0.j,
0.3101057 +0.j, 0.365514 +0.j],
...65514 +0.j, 0.46312022+0.j, 0.6176184 +0.j, 0.4594969 +0.j,
0.60431015+0.j, 4.433647 +0.j]], dtype=complex64)
diag1_ = array([-0.44244403, -0.20154493, 0.01760732, 0.16953002, 0.3570978 ,
0.7214176 ], dtype=float32)
diag_ = array([-0.7402807 , -0.37563217, 0.06704548, 0.12726681, 0.40285635,
3.2228992 ], dtype=float32)
dim = 6
dtype_ = 'F'
eigvals = None
lower = True
overwrite = False
self = <scipy.linalg.tests.test_decomp.TestEigh object at 0x1494072b0>
turbo = False
w = array([-0.2212221 +0.j, -0.10077266+0.j, 0.00880363+0.j, 0.08476513+0.j,
0.17854875+0.j, 0.5256907 +0.j], dtype=complex64)
z = array([[-0.17789046-0.17789167j, -0.21407811-0.21442269j,
-0.3454858 -0.3451161j , -0.25057662-0.25060064j,
....32072693j, -0.1055024 -0.10548996j,
-0.04871838-0.04871838j, 0.14641598-0.0470806j ]],
dtype=complex64)
________________ TestEigh.test_eigh[6-F-False-False-True-None] _________________
scipy/linalg/tests/test_decomp.py:892: in test_eigh
assert_allclose(diag1_, w, rtol=0., atol=atol)
E AssertionError:
E Not equal to tolerance rtol=0, atol=0.0001
E
E Mismatched elements: 6 / 6 (100%)
E Max absolute difference: 0.47770208
E Max relative difference: 1.0000014
E x: array([-0.111534, -0.045158, 0.016268, 0.269463, 0.379338, 0.955574],
E dtype=float32)
E y: array([-0.055767+0.j, -0.022579+0.j, 0.008134+0.j, 0.134731+0.j,
E 0.189669+0.j, 0.477872+0.j], dtype=complex64)
a = array([[0.6452997 +0.j, 0.5964542 +0.j, 0.5424254 +0.j, 0.5148393 +0.j,
0.57869583+0.j, 0.58155155+0.j],
...8155155+0.j, 0.38993406+0.j, 0.45339584+0.j, 0.44374526+0.j,
0.26262194+0.j, 0.35480455+0.j]], dtype=complex64)
atol = 0.0001
b = array([[3.587556 +0.j, 0.35207134+0.j, 0.3048609 +0.j, 0.28332126+0.j,
0.22228493+0.j, 0.38722858+0.j],
...8722858+0.j, 0.7595627 +0.j, 0.65089786+0.j, 0.46203086+0.j,
0.33025602+0.j, 3.6631072 +0.j]], dtype=complex64)
diag1_ = array([-0.11153392, -0.04515787, 0.01626768, 0.26946256, 0.37933812,
0.9555737 ], dtype=float32)
diag_ = array([-0.54924893, -0.46740282, -0.08804972, 0.21839237, 0.5904868 ,
3.1221344 ], dtype=float32)
dim = 6
dtype_ = 'F'
eigvals = None
lower = False
overwrite = False
self = <scipy.linalg.tests.test_decomp.TestEigh object at 0x149407610>
turbo = True
w = array([-0.05576696+0.j, -0.02257896+0.j, 0.00813383+0.j, 0.13473132+0.j,
0.18966901+0.j, 0.4778716 +0.j], dtype=complex64)
z = array([[-0.28307393-0.28307393j, -0.24174953-0.24174953j,
-0.26248735-0.26248735j, 0.01137766+0.01137766j,
....33739296j, 0.11381758+0.11381758j,
-0.05777605-0.05777605j, -0.12800555-0.12734008j]],
dtype=complex64)
________________ TestEigh.test_eigh[6-F-False-False-False-None] ________________
scipy/linalg/tests/test_decomp.py:892: in test_eigh
assert_allclose(diag1_, w, rtol=0., atol=atol)
E AssertionError:
E Not equal to tolerance rtol=0, atol=0.0001
E
E Mismatched elements: 6 / 6 (100%)
E Max absolute difference: 0.5541795
E Max relative difference: 1.0005932
E x: array([-0.419803, -0.325263, -0.18648 , -0.031227, 0.369174, 1.10803 ],
E dtype=float32)
E y: array([-0.209902+0.j, -0.162631+0.j, -0.09324 +0.j, -0.015613+0.j,
E 0.184587+0.j, 0.553851+0.j], dtype=complex64)
a = array([[0.8223416 +0.j, 0.7209296 +0.j, 0.17268606+0.j, 0.4822605 +0.j,
0.5141309 +0.j, 0.91035676+0.j],
...1035676+0.j, 0.42631674+0.j, 0.44917846+0.j, 0.57719296+0.j,
0.73331213+0.j, 0.13215497+0.j]], dtype=complex64)
atol = 0.0001
b = array([[3.9971466 +0.j, 0.4122611 +0.j, 0.6195147 +0.j, 0.27250645+0.j,
0.37227485+0.j, 0.27615562+0.j],
...7615562+0.j, 0.77285135+0.j, 0.6765339 +0.j, 0.56973404+0.j,
0.21664503+0.j, 3.7920017 +0.j]], dtype=complex64)
diag1_ = array([-0.41980323, -0.325263 , -0.18648006, -0.03122663, 0.36917424,
1.1080304 ], dtype=float32)
diag_ = array([-0.72434413, -0.25737715, 0.12102426, 0.35222927, 0.89393467,
2.8990715 ], dtype=float32)
dim = 6
dtype_ = 'F'
eigvals = None
lower = False
overwrite = False
self = <scipy.linalg.tests.test_decomp.TestEigh object at 0x14ac88610>
turbo = False
w = array([-0.20990156+0.j, -0.16263144+0.j, -0.09324001+0.j, -0.01561319+0.j,
0.18458709+0.j, 0.55385095+0.j], dtype=complex64)
z = array([[ 0.20411697+0.20411697j, -0.15473579-0.15473579j,
-0.00154966-0.00154966j, 0.05437683+0.05437683j,
....15470162j, 0.157751 +0.157751j ,
-0.11361784-0.11361784j, -0.17066842-0.16511473j]],
dtype=complex64)
_________________ TestEigh.test_eigh[6-D-True-True-True-None] __________________
scipy/linalg/tests/test_decomp.py:892: in test_eigh
assert_allclose(diag1_, w, rtol=0., atol=atol)
E AssertionError:
E Not equal to tolerance rtol=0, atol=1e-11
E
E Mismatched elements: 6 / 6 (100%)
E Max absolute difference: 0.46721899
E Max relative difference: 1.
E x: array([-0.208976, 0.012785, 0.053826, 0.312004, 0.48156 , 0.934438])
E y: array([-0.104488+0.j, 0.006393+0.j, 0.026913+0.j, 0.156002+0.j,
E 0.24078 +0.j, 0.467219+0.j])
a = array([[0.71151329+0.j, 0.3093607 +0.j, 0.52717286+0.j, 0.46442694+0.j,
0.40183626+0.j, 0.58061673+0.j],
...0.j],
[0.58061673+0.j, 0.12754241+0.j, 0.24559284+0.j, 0.3463392 +0.j,
0.85971392+0.j, 0.78403115+0.j]])
atol = 1e-11
b = array([[3.53282904+0.j, 0.56861146+0.j, 0.47723147+0.j, 0.58576209+0.j,
0.66667433+0.j, 0.56435813+0.j],
...0.j],
[0.56435813+0.j, 0.5475766 +0.j, 0.3108049 +0.j, 0.73431137+0.j,
0.61806327+0.j, 4.14883347+0.j]])
diag1_ = array([-0.20897559, 0.01278527, 0.05382603, 0.31200375, 0.48156019,
0.93443799])
diag_ = array([-0.40828611, -0.26025893, 0.14991251, 0.31064419, 0.67324649,
2.68692907])
dim = 6
dtype_ = 'D'
eigvals = None
lower = True
overwrite = True
self = <scipy.linalg.tests.test_decomp.TestEigh object at 0x149439ee0>
turbo = True
w = array([-0.10448779+0.j, 0.00639263+0.j, 0.02691302+0.j, 0.15600187+0.j,
0.2407801 +0.j, 0.467219 +0.j])
z = array([[ 0.10179151+0.10179151j, -0.19182457-0.19182457j,
-0.28036969-0.28036969j, 0.38349856+0.38349856j,
...,
-0.00571041-0.00571041j, -0.04042499-0.04042499j,
0.3445809 +0.3445809j , -0.14721246-0.14723448j]])
_________________ TestEigh.test_eigh[6-D-True-True-False-None] _________________
scipy/linalg/tests/test_decomp.py:892: in test_eigh
assert_allclose(diag1_, w, rtol=0., atol=atol)
E AssertionError:
E Not equal to tolerance rtol=0, atol=1e-11
E
E Mismatched elements: 6 / 6 (100%)
E Max absolute difference: 0.48914139
E Max relative difference: 1.
E x: array([-0.440268, -0.108758, -0.045163, 0.209923, 0.304468, 0.978286])
E y: array([-0.220134+0.j, -0.054379+0.j, -0.022581+0.j, 0.104961+0.j,
E 0.152234+0.j, 0.489144+0.j])
a = array([[0.83184019+0.j, 0.44850066+0.j, 0.48224936+0.j, 0.93117523+0.j,
0.7562725 +0.j, 0.61425982+0.j],
...0.j],
[0.61425982+0.j, 0.10912077+0.j, 0.19157452+0.j, 0.4631397 +0.j,
0.41174307+0.j, 0.36281503+0.j]])
atol = 1e-11
b = array([[4.44894656+0.j, 0.38731322+0.j, 0.53842472+0.j, 0.40216099+0.j,
0.2844065 +0.j, 0.48004812+0.j],
...0.j],
[0.48004812+0.j, 0.9264048 +0.j, 0.40715136+0.j, 0.50295856+0.j,
0.697232 +0.j, 4.43435686+0.j]])
diag1_ = array([-0.44026798, -0.10875835, -0.04516259, 0.2099229 , 0.30446826,
0.97828587])
diag_ = array([-0.82087746, -0.36496181, -0.16815352, 0.37702657, 0.79248229,
2.88464801])
dim = 6
dtype_ = 'D'
eigvals = None
lower = True
overwrite = True
self = <scipy.linalg.tests.test_decomp.TestEigh object at 0x149439040>
turbo = False
w = array([-0.22013399+0.j, -0.05437918+0.j, -0.0225813 +0.j, 0.10496145+0.j,
0.15223413+0.j, 0.48914448+0.j])
z = array([[ 0.18728137+0.18728137j, -0.01264387-0.01264387j,
-0.30275018-0.30275018j, 0.02387516+0.02387516j,
...,
0.36218632+0.36218632j, -0.04173572-0.04173572j,
-0.30788209-0.30788209j, -0.04906664-0.0492157j ]])
_________________ TestEigh.test_eigh[6-D-True-False-True-None] _________________
scipy/linalg/tests/test_decomp.py:892: in test_eigh
assert_allclose(diag1_, w, rtol=0., atol=atol)
E AssertionError:
E Not equal to tolerance rtol=0, atol=1e-11
E
E Mismatched elements: 6 / 6 (100%)
E Max absolute difference: 0.45220174
E Max relative difference: 1.
E x: array([-0.430458, -0.343752, -0.126716, 0.142294, 0.225441, 0.904403])
E y: array([-0.215229+0.j, -0.171876+0.j, -0.063358+0.j, 0.071147+0.j,
E 0.11272 +0.j, 0.452202+0.j])
a = array([[0.20165145+0.j, 0.57751939+0.j, 0.46180506+0.j, 0.54345016+0.j,
0.50644849+0.j, 0.3592125 +0.j],
...0.j],
[0.3592125 +0.j, 0.43640971+0.j, 0.45532387+0.j, 0.95986398+0.j,
0.27387325+0.j, 0.45295321+0.j]])
atol = 1e-11
b = array([[4.29417849+0.j, 0.86871958+0.j, 0.47298458+0.j, 0.76512999+0.j,
0.75171527+0.j, 0.3393883 +0.j],
...0.j],
[0.3393883 +0.j, 0.76192083+0.j, 0.1011976 +0.j, 0.80656604+0.j,
0.92261296+0.j, 3.48071444+0.j]])
diag1_ = array([-0.4304576 , -0.3437524 , -0.12671572, 0.14229393, 0.22544075,
0.9044035 ])
diag_ = array([-0.71297136, -0.50447194, -0.25125495, -0.08989238, 0.30587581,
2.84549194])
dim = 6
dtype_ = 'D'
eigvals = None
lower = False
overwrite = True
self = <scipy.linalg.tests.test_decomp.TestEigh object at 0x14941b550>
turbo = True
w = array([-0.2152288 +0.j, -0.1718762 +0.j, -0.06335786+0.j, 0.07114697+0.j,
0.11272037+0.j, 0.45220175+0.j])
z = array([[-0.01940712-0.01940712j, -0.28205045-0.28205045j,
-0.34314316-0.34314316j, -0.20374353-0.20374353j,
...,
0.12946989+0.12946989j, 0.3215643 +0.3215643j ,
-0.17511595-0.17511595j, -0.20654279-0.20652544j]])
________________ TestEigh.test_eigh[6-D-True-False-False-None] _________________
scipy/linalg/tests/test_decomp.py:892: in test_eigh
assert_allclose(diag1_, w, rtol=0., atol=atol)
E AssertionError:
E Not equal to tolerance rtol=0, atol=1e-11
E
E Mismatched elements: 6 / 6 (100%)
E Max absolute difference: 0.21458406
E Max relative difference: 1.
E x: array([-0.429168, -0.226661, -0.057259, 0.03991 , 0.368466, 0.453574])
E y: array([-0.214584+0.j, -0.11333 +0.j, -0.028629+0.j, 0.019955+0.j,
E 0.184233+0.j, 0.534686+0.j])
a = array([[0.84857284+0.j, 0.54979222+0.j, 0.61389136+0.j, 0.71715856+0.j,
0.71611654+0.j, 0.14062217+0.j],
...0.j],
[0.14062217+0.j, 0.28443538+0.j, 0.38469069+0.j, 0.09659502+0.j,
0.75176294+0.j, 0.47611267+0.j]])
atol = 1e-11
b = array([[3.56819358+0.j, 0.09731193+0.j, 0.76233493+0.j, 0.61009798+0.j,
0.50667643+0.j, 0.22415343+0.j],
...0.j],
[0.22415343+0.j, 0.5987129 +0.j, 0.79354542+0.j, 0.35780076+0.j,
0.52505517+0.j, 4.03412453+0.j]])
diag1_ = array([-0.42916812, -0.2266605 , -0.05725854, 0.03990991, 0.36846629,
0.45357373])
diag_ = array([-0.45961248, -0.27479489, 0.03087143, 0.18251432, 0.37059721,
3.67821052])
dim = 6
dtype_ = 'D'
eigvals = None
lower = False
overwrite = True
self = <scipy.linalg.tests.test_decomp.TestEigh object at 0x14937ef10>
turbo = False
w = array([-0.21458406+0.j, -0.11333025+0.j, -0.02862927+0.j, 0.01995495+0.j,
0.18423314+0.j, 0.534686 +0.j])
z = array([[-0.07062765-0.07062764j, -0.15635187-0.15635196j,
0.15145832+0.15145832j, 0.37445072+0.37445072j,
...,
0.31546351+0.31546351j, -0.01766316-0.01766316j,
0.35818957+0.35818957j, 0.04611843+0.51517196j]])
_________________ TestEigh.test_eigh[6-D-False-True-True-None] _________________
scipy/linalg/tests/test_decomp.py:892: in test_eigh
assert_allclose(diag1_, w, rtol=0., atol=atol)
E AssertionError:
E Not equal to tolerance rtol=0, atol=1e-11
E
E Mismatched elements: 6 / 6 (100%)
E Max absolute difference: 0.26617536
E Max relative difference: 1.
E x: array([-0.532351, -0.32691 , -0.032601, 0.087917, 0.275316, 0.862828])
E y: array([-0.266175+0.j, -0.163455+0.j, -0.0163 +0.j, 0.043958+0.j,
E 0.137658+0.j, 0.608115+0.j])
a = array([[0.19726992+0.j, 0.74086873+0.j, 0.247372 +0.j, 0.36777871+0.j,
0.74380484+0.j, 0.6093195 +0.j],
...0.j],
[0.6093195 +0.j, 0.4378835 +0.j, 0.7100433 +0.j, 0.34691623+0.j,
0.83602676+0.j, 0.58308823+0.j]])
atol = 1e-11
b = array([[3.9406136 +0.j, 0.54713909+0.j, 0.13773987+0.j, 0.23951298+0.j,
0.48738011+0.j, 0.69066462+0.j],
...0.j],
[0.69066462+0.j, 0.5388175 +0.j, 0.64299084+0.j, 0.29212669+0.j,
0.45186728+0.j, 4.04272396+0.j]])
diag1_ = array([-0.53235071, -0.32691028, -0.03260072, 0.08791695, 0.27531566,
0.86282839])
diag_ = array([-0.67218586, -0.07246483, -0.0063822 , 0.43931345, 0.87416187,
3.23329803])
dim = 6
dtype_ = 'D'
eigvals = None
lower = True
overwrite = False
self = <scipy.linalg.tests.test_decomp.TestEigh object at 0x14937e790>
turbo = True
w = array([-0.26617536+0.j, -0.16345514+0.j, -0.01630036+0.j, 0.04395847+0.j,
0.13765783+0.j, 0.60811538+0.j])
z = array([[-1.77013896e-01-1.77013896e-01j, 3.42370878e-01+3.42370878e-01j,
-1.95549287e-01-1.95549287e-01j, 2....5464e-01j, -1.64962022e-01-1.64962022e-01j,
2.10487466e-01+2.10487466e-01j, 1.36618340e-01+3.10923431e-01j]])
________________ TestEigh.test_eigh[6-D-False-True-False-None] _________________
scipy/linalg/tests/test_decomp.py:892: in test_eigh
assert_allclose(diag1_, w, rtol=0., atol=atol)
E AssertionError:
E Not equal to tolerance rtol=0, atol=1e-11
E
E Mismatched elements: 6 / 6 (100%)
E Max absolute difference: 0.53240441
E Max relative difference: 1.57960175
E x: array([-0.869454, -0.277809, 0.139136, 0.220895, 0.452656, 0.861342])
E y: array([-0.33705 +0.j, -0.138904+0.j, 0.069568+0.j, 0.110448+0.j,
E 0.226328+0.j, 0.430671+0.j])
a = array([[0.27261062+0.j, 0.04205495+0.j, 0.51513428+0.j, 0.13874065+0.j,
0.44880645+0.j, 0.81168894+0.j],
...0.j],
[0.81168894+0.j, 0.87849289+0.j, 0.06339601+0.j, 0.24770736+0.j,
0.58057842+0.j, 0.22153806+0.j]])
atol = 1e-11
b = array([[3.85261617+0.j, 0.35970308+0.j, 0.58669934+0.j, 0.34498049+0.j,
0.32272017+0.j, 0.6520558 +0.j],
...0.j],
[0.6520558 +0.j, 0.49537484+0.j, 0.6639669 +0.j, 0.61580777+0.j,
0.7446846 +0.j, 3.50357128+0.j]])
diag1_ = array([-0.86945419, -0.27780894, 0.13913612, 0.22089507, 0.45265554,
0.8613416 ])
diag_ = array([-0.80106543, -0.59442604, -0.00617281, 0.15573334, 0.4056252 ,
3.01405565])
dim = 6
dtype_ = 'D'
eigvals = None
lower = True
overwrite = False
self = <scipy.linalg.tests.test_decomp.TestEigh object at 0x14ac6e3d0>
turbo = False
w = array([-0.33704977+0.j, -0.13890447+0.j, 0.06956806+0.j, 0.11044753+0.j,
0.22632777+0.j, 0.4306708 +0.j])
z = array([[ 0.27253022+0.6251033j , 0.14597259+0.14597259j,
-0.2026133 -0.2026133j , -0.32725165-0.32725165j,
...,
0.07314678+0.07314679j, -0.0589045 -0.05890449j,
0.34644687+0.34644687j, -0.14780682-0.14780682j]])
________________ TestEigh.test_eigh[6-D-False-False-True-None] _________________
scipy/linalg/tests/test_decomp.py:892: in test_eigh
assert_allclose(diag1_, w, rtol=0., atol=atol)
E AssertionError:
E Not equal to tolerance rtol=0, atol=1e-11
E
E Mismatched elements: 6 / 6 (100%)
E Max absolute difference: 0.46400199
E Max relative difference: 1.
E x: array([-0.47534 , -0.277137, -0.075012, 0.075452, 0.183818, 0.020822])
E y: array([-0.23767 +0.j, -0.138569+0.j, -0.037506+0.j, 0.037726+0.j,
E 0.091909+0.j, 0.484824+0.j])
a = array([[0.16360209+0.j, 0.74914928+0.j, 0.30434003+0.j, 0.57516177+0.j,
0.52795367+0.j, 0.72242201+0.j],
...0.j],
[0.72242201+0.j, 0.46992984+0.j, 0.60366644+0.j, 0.66528242+0.j,
0.30804763+0.j, 0.37233133+0.j]])
atol = 1e-11
b = array([[3.4803405 +0.j, 0.65967944+0.j, 0.65800516+0.j, 0.34342059+0.j,
0.52096372+0.j, 0.40808915+0.j],
...0.j],
[0.40808915+0.j, 0.12017123+0.j, 0.55640349+0.j, 0.50448368+0.j,
0.26358134+0.j, 3.94508881+0.j]])
diag1_ = array([-0.47534025, -0.27713705, -0.07501199, 0.07545225, 0.18381825,
0.0208216 ])
diag_ = array([-0.47864621, 0.09167204, 0.24625594, 0.42640824, 0.5104482 ,
3.23696316])
dim = 6
dtype_ = 'D'
eigvals = None
lower = False
overwrite = False
self = <scipy.linalg.tests.test_decomp.TestEigh object at 0x14ac6e8b0>
turbo = True
w = array([-0.23767012+0.j, -0.13856853+0.j, -0.03750599+0.j, 0.03772612+0.j,
0.09190912+0.j, 0.48482359+0.j])
z = array([[ 0.43073616+0.43073616j, -0.14501967-0.14501967j,
-0.05811003-0.05811003j, 0.16165742+0.16165743j,
...,
-0.39225138-0.39225139j, -0.10643115-0.10643115j,
0.1075131 +0.1075131j , 0.19312526+0.51353835j]])
________________ TestEigh.test_eigh[6-D-False-False-False-None] ________________
scipy/linalg/tests/test_decomp.py:892: in test_eigh
assert_allclose(diag1_, w, rtol=0., atol=atol)
E AssertionError:
E Not equal to tolerance rtol=0, atol=1e-11
E
E Mismatched elements: 6 / 6 (100%)
E Max absolute difference: 0.5927337
E Max relative difference: 1.
E x: array([-0.3406 , -0.109432, -0.048392, 0.114588, 0.43081 , 1.185471])
E y: array([-0.1703 +0.j, -0.054716+0.j, -0.024196+0.j, 0.057294+0.j,
E 0.215405+0.j, 0.592737+0.j])
a = array([[0.66935268+0.j, 0.77780905+0.j, 0.55409121+0.j, 0.61948759+0.j,
0.61420033+0.j, 0.55203713+0.j],
...0.j],
[0.55203713+0.j, 0.64195696+0.j, 0.64273712+0.j, 0.83521969+0.j,
0.34542958+0.j, 0.41701111+0.j]])
atol = 1e-11
b = array([[4.1987233 +0.j, 0.30511017+0.j, 0.44010466+0.j, 0.7037236 +0.j,
0.31646723+0.j, 0.57026433+0.j],
...0.j],
[0.57026433+0.j, 0.40125308+0.j, 0.73838623+0.j, 0.30892016+0.j,
0.06773588+0.j, 4.12817453+0.j]])
diag1_ = array([-0.34059977, -0.10943169, -0.04839208, 0.11458819, 0.4308101 ,
1.18547053])
diag_ = array([-1.04998803, -0.49886791, -0.01795759, 0.26587977, 0.49271342,
3.33323287])
dim = 6
dtype_ = 'D'
eigvals = None
lower = False
overwrite = False
self = <scipy.linalg.tests.test_decomp.TestEigh object at 0x14935b580>
turbo = False
w = array([-0.17029989+0.j, -0.05471584+0.j, -0.02419604+0.j, 0.05729409+0.j,
0.21540505+0.j, 0.59273683+0.j])
z = array([[-0.045203 -0.045203j , -0.26777574-0.26777574j,
0.2606524 +0.2606524j , -0.24979522-0.24979522j,
...,
0.16915973+0.16915973j, 0.28749646+0.28749646j,
-0.01672724-0.01672724j, 0.15803141+0.15789633j]])
______________________ TestSVD_GESVD.test_simple_complex _______________________
scipy/linalg/tests/test_decomp.py:1111: in test_simple_complex
assert_array_almost_equal(u.conj().T @ u, eye(u.shape[1]))
E AssertionError:
E Arrays are not almost equal to 6 decimals
E
E Mismatched elements: 9 / 9 (100%)
E Max absolute difference: 0.97054768
E Max relative difference: 0.97054768
E x: array([[ 1.863167-1.286402e-17j, 0.252408-5.630535e-02j,
E 0.024053+4.787066e-04j],
E [ 0.252408+5.630535e-02j, 0.915003-1.847309e-17j,...
E y: array([[1., 0., 0.],
E [0., 1., 0.],
E [0., 0., 1.]])
a = [[1, 2, 3], [1, 2j, 3], [2, 5, 6]]
full_matrices = True
s = array([9.34187884e+00, 2.39359557e+00, 6.05073700e-16])
self = <scipy.linalg.tests.test_decomp.TestSVD_GESVD object at 0x14935bee0>
u = array([[-0.39828431-0.37357977j, -0.09163318-0.09133861j,
0.90369977+0.90369977j],
[-0.29651332-0.3668... -0.11583439-0.11583439j],
[-0.8576312 -0.77905014j, 0.15293019+0.15386718j,
-0.39393269-0.39393269j]])
vh = array([[-0.25798451+0.00000000e+00j, -0.57043932-9.50732193e-02j,
-0.77395353-4.44968492e-18j],
[-0.182...9496108e-16j],
[-0.9486833 -0.00000000e+00j, -0. +6.24500451e-17j,
0.31622777-4.49574451e-17j]])
______________________ TestSVD_GESVD.test_random_complex _______________________
scipy/linalg/tests/test_decomp.py:1127: in test_random_complex
assert_array_almost_equal(u.conj().T @ u,
E AssertionError:
E Arrays are not almost equal to 6 decimals
E
E Mismatched elements: 375 / 400 (93.8%)
E Max absolute difference: 19.52939442
E Max relative difference: 19.52939442
E x: array([[ 8.586701e-01+1.046603e-18j, 6.314249e-03+2.182309e-08j,
E 2.876587e-03+2.722746e-05j, -4.915836e-02-2.165741e-15j,
E 5.267591e-03+2.317472e-16j, -2.283455e-02-1.240581e-15j,...
E y: array([[1., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
E 0., 0., 0., 0.],
E [0., 1., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,...
a = array([[0.19151945+0.86498179j, 0.62210877+0.71027069j,
0.43772774+0.17063169j, 0.78535858+0.85508266j,
...856j, 0.60667522+0.8112992j ,
0.82215979+0.8727098j , 0.62796507+0.66598823j,
0.11792306+0.58878655j]])
full_matrices = True
i = 0
m = 15
n = 20
s = array([12.86143638, 2.94843302, 2.55481777, 2.3477689 , 2.11704004,
2.09284945, 1.8806212 , 1.70856509, 1.42069197, 1.32153852,
1.04277929, 0.93685169, 0.78054485, 0.71154826, 0.38778732])
self = <scipy.linalg.tests.test_decomp.TestSVD_GESVD object at 0x14ac68820>
u = array([[-1.65378616e-01-1.65378616e-01j, 1.16847751e-01+1.16848109e-01j,
2.29235212e-02+2.29707931e-02j, -3....0115e-01j, -8.46257457e-02+1.85156223e-01j,
1.85181013e-02-1.78106928e-02j, 4.70380837e-01+4.10219065e-02j]])
vh = array([[-0.25763384+0.j , -0.24905824+0.02672578j,
-0.26160805-0.00707567j, -0.2333506 +0.02518936j,
..., -0.18131298+0.0098072j ,
-0.06945942-0.10975804j, 0.12818036+0.03224759j,
0.21116326+0.40193129j]])
____________________________ TestQZ.test_qz_complex ____________________________
scipy/linalg/tests/test_decomp.py:2102: in test_qz_complex
assert_array_almost_equal(Q @ AA @ Z.conj().T, A)
E AssertionError:
E Arrays are not almost equal to 6 decimals
E
E Mismatched elements: 25 / 25 (100%)
E Max absolute difference: 1.24268613
E Max relative difference: 1.04389533
E x: array([[-0.004922+0.003162j, 0.03643 -0.042272j, -0.004731+0.007782j,
E 0.007437-0.010201j, -0.017222+0.015339j],
E [-0.003225-0.000939j, 0.017724-0.004525j, -0.002029+0.001507j,...
E y: array([[0.929616+0.729689j, 0.316376+0.994015j, 0.183919+0.676874j,
E 0.20456 +0.790823j, 0.567725+0.170914j],
E [0.595545+0.026849j, 0.964515+0.80037j , 0.653177+0.903723j,...
A = array([[0.92961609+0.72968908j, 0.31637555+0.99401459j,
0.18391881+0.67687371j, 0.20456028+0.79082252j,
...906j, 0.71745362+0.46810575j,
0.46759901+0.4593732j , 0.32558468+0.70950978j,
0.43964461+0.17805301j]])
AA = array([[-1.23486358e-01-7.15399961e-15j, -1.81158310e-01+7.82707232e-15j,
3.58634296e-04-5.29225448e-04j, -8....j,
0.00000000e+00+0.00000000e+00j, 0.00000000e+00+0.00000000e+00j,
4.87840651e-02-6.69110860e-02j]])
B = array([[0.53144988+0.22563761j, 0.16774223+0.24497444j,
0.76881392+0.7928007j , 0.92817055+0.49517241j,
...536j, 0.49541476+0.33605516j,
0.35161653+0.82638043j, 0.71423037+0.89810064j,
0.50392912+0.0427153j ]])
BB = array([[ 0.00161829+0.00000000e+00j, -0.06737054+9.71445147e-17j,
0.00277447-4.09419490e-03j, 0.03008124+1.8...00000000e+00j,
0. +0.00000000e+00j, 0. +0.00000000e+00j,
0.78879555+0.00000000e+00j]])
Q = array([[-0.0124049 -0.0124049j , -0.08098177-0.08098177j,
-0.00857385-0.00857385j, 0.00053118+0.00053118j,
..., -0.13537673-0.13537673j,
-0.01433514-0.01433514j, 0.00358038+0.00358038j,
-0.09321653-0.09321653j]])
Z = array([[ 7.69175751e-05+9.66024136e-17j, 1.28001480e-02+2.56045185e-15j,
-2.95299146e-03+4.35763742e-03j, 1....j,
-7.15035354e-03+1.05515538e-02j, 2.66907399e-01+4.51503378e-08j,
3.83963303e-01-5.26635112e-01j]])
n = 5
self = <scipy.linalg.tests.test_decomp.TestQZ object at 0x14ac639a0>
___________________________ TestQZ.test_qz_complex64 ___________________________
scipy/linalg/tests/test_decomp.py:2114: in test_qz_complex64
assert_array_almost_equal(Q @ AA @ Z.conj().T, A, decimal=5)
E AssertionError:
E Arrays are not almost equal to 5 decimals
E
E Mismatched elements: 25 / 25 (100%)
E Max absolute difference: 1.3936292
E Max relative difference: 2.9828136
E x: array([[ 0.07049+0.40556j, 0.06169+0.15344j, 0.00871+0.00982j,
E 0.00099+0.09088j, -0.03962+0.14536j],
E [ 0.36104-0.14102j, -0.03418-0.17164j, -0.05379-0.05545j,...
E y: array([[0.92962+0.72969j, 0.31638+0.99401j, 0.18392+0.67687j,
E 0.20456+0.79082j, 0.56773+0.17091j],
E [0.59554+0.02685j, 0.96451+0.80037j, 0.65318+0.90372j,...
A = array([[0.9296161 +0.72968906j, 0.31637555+0.99401456j,
0.18391882+0.6768737j , 0.20456028+0.7908225j ,
...0.46810576j,
0.467599 +0.4593732j , 0.32558468+0.7095098j ,
0.4396446 +0.178053j ]], dtype=complex64)
AA = array([[-3.17302674e-01-3.1730002e-01j, -1.10790467e+00+5.9604645e-08j,
-9.15605366e-01-5.3644180e-07j, 2.068...00000000e+00+0.0000000e+00j, 0.00000000e+00+0.0000000e+00j,
-1.09895885e-01-6.9619828e-01j]], dtype=complex64)
B = array([[0.5314499 +0.2256376j , 0.16774222+0.24497443j,
0.7688139 +0.7928007j , 0.92817056+0.4951724j ,
...0.33605516j,
0.35161653+0.82638043j, 0.71423036+0.8981006j ,
0.50392914+0.0427153j ]], dtype=complex64)
BB = array([[ 0. +0.0000000e+00j, -0.37612104+0.0000000e+00j,
-0.15505296-8.1956387e-08j, 0.53143716-3.2782...
0. +0.0000000e+00j, 0. +0.0000000e+00j,
0.3922476 +0.0000000e+00j]], dtype=complex64)
Q = array([[-0.02583549-0.13805246j, -0.0003476 -0.00034761j,
0.00028997+0.00028999j, 0.01062595+0.01062595j,
...0810727j,
0.00789832+0.00789865j, 0.15082175+0.15082175j,
0.23858337+0.23858342j]], dtype=complex64)
Z = array([[ 0.17593992+1.7593990e-01j, 0.73558056-2.9802322e-08j,
0.37611586+2.0861626e-07j, -0.8283269 +5.3644...
0.06691903+1.8626451e-08j, -0.6479421 +4.1723251e-07j,
-0.13334364-8.4474075e-01j]], dtype=complex64)
n = 5
self = <scipy.linalg.tests.test_decomp.TestQZ object at 0x1497dcdc0>
________________________ TestQZ.test_qz_double_complex _________________________
scipy/linalg/tests/test_decomp.py:2125: in test_qz_double_complex
AA, BB, Q, Z = qz(A, B, output='complex')
A = array([[0.92961609, 0.31637555, 0.18391881, 0.20456028, 0.56772503],
[0.5955447 , 0.96451452, 0.6531771 , 0.748... 0.80981255, 0.87217591, 0.9646476 , 0.72368535],
[0.64247533, 0.71745362, 0.46759901, 0.32558468, 0.43964461]])
B = array([[0.72968908, 0.99401459, 0.67687371, 0.79082252, 0.17091426],
[0.02684928, 0.80037024, 0.90372254, 0.024... 0.50022275, 0.81018941, 0.09596853, 0.21895004],
[0.25871906, 0.46810575, 0.4593732 , 0.70950978, 0.17805301]])
n = 5
self = <scipy.linalg.tests.test_decomp.TestQZ object at 0x1497dc790>
scipy/linalg/_decomp_qz.py:258: in qz
result, _ = _qz(A, B, output=output, lwork=lwork, sort=sort,
A = array([[0.92961609, 0.31637555, 0.18391881, 0.20456028, 0.56772503],
[0.5955447 , 0.96451452, 0.6531771 , 0.748... 0.80981255, 0.87217591, 0.9646476 , 0.72368535],
[0.64247533, 0.71745362, 0.46759901, 0.32558468, 0.43964461]])
B = array([[0.72968908, 0.99401459, 0.67687371, 0.79082252, 0.17091426],
[0.02684928, 0.80037024, 0.90372254, 0.024... 0.50022275, 0.81018941, 0.09596853, 0.21895004],
[0.25871906, 0.46810575, 0.4593732 , 0.70950978, 0.17805301]])
check_finite = True
lwork = None
output = 'complex'
overwrite_a = False
overwrite_b = False
sort = None
scipy/linalg/_decomp_qz.py:127: in _qz
warnings.warn("The QZ iteration failed. (a,b) are not in Schur "
E scipy.linalg.misc.LinAlgWarning: The QZ iteration failed. (a,b) are not in Schur form, but ALPHAR(j), ALPHAI(j), and BETA(j) should be correct for J=4,...,N
A = array([[0.92961609, 0.31637555, 0.18391881, 0.20456028, 0.56772503],
[0.5955447 , 0.96451452, 0.6531771 , 0.748... 0.80981255, 0.87217591, 0.9646476 , 0.72368535],
[0.64247533, 0.71745362, 0.46759901, 0.32558468, 0.43964461]])
B = array([[0.72968908, 0.99401459, 0.67687371, 0.79082252, 0.17091426],
[0.02684928, 0.80037024, 0.90372254, 0.024... 0.50022275, 0.81018941, 0.09596853, 0.21895004],
[0.25871906, 0.46810575, 0.4593732 , 0.70950978, 0.17805301]])
a1 = array([[0.92961609+0.j, 0.31637555+0.j, 0.18391881+0.j, 0.20456028+0.j,
0.56772503+0.j],
[0.5955447 +0.... 0.72368535+0.j],
[0.64247533+0.j, 0.71745362+0.j, 0.46759901+0.j, 0.32558468+0.j,
0.43964461+0.j]])
a_m = 5
a_n = 5
b1 = array([[0.72968908+0.j, 0.99401459+0.j, 0.67687371+0.j, 0.79082252+0.j,
0.17091426+0.j],
[0.02684928+0.... 0.21895004+0.j],
[0.25871906+0.j, 0.46810575+0.j, 0.4593732 +0.j, 0.70950978+0.j,
0.17805301+0.j]])
b_m = 5
b_n = 5
check_finite = True
gges = <fortran object>
info = 5
lwork = 165
output = 'complex'
overwrite_a = True
overwrite_b = True
result = (array([[-2.34524812e-01-2.34524833e-01j, 5.96844304e-01+5.97497505e-01j,
3.16520812e-01+3.16520812e-01j, -1...6647458-0.0664754j ,
-0.00333645-0.00333645j, 0.01643819+0.0165189j ,
-0.21114442-0.21115542j]]), ...)
sfunction = <function _qz.<locals>.<lambda> at 0x149360af0>
sort = None
typa = 'D'
typb = 'D'
______________________________ TestOrdQZ.test_lhp ______________________________
scipy/linalg/tests/test_decomp.py:2366: in test_lhp
self.check_all('lhp')
self = <scipy.linalg.tests.test_decomp.TestOrdQZ object at 0x1497edd00>
scipy/linalg/tests/test_decomp.py:2363: in check_all
self.check(Ai, Bi, sort, *reti)
Ai = array([[-21.1 -22.5j , 53.5 -50.5j , -34.5 +127.5j , 7.5 +0.5j ],
[ -0.46 -7.78j, -3.5 -37.5j , -15.5....7 -17.1j , -68.5 +12.5j , -7.5 -3.5j ],
[ 5.5 +4.4j , 14.4 +43.3j , -32.5 -46.j , -19. -32.5j ]])
Bi = array([[ 1. -5.j , 1.6+1.2j, -3. +0.j , 0. -1.j ],
[ 0.8-0.6j, 0. -5.j , -4. +3.j , -2.4-3.2j],
[ 1. +0.j , 2.4+1.8j, -4. -5.j , 0. -3.j ],
[ 0. +1.j , -1.8+2.4j, 0. -4.j , 4. -5.j ]])
ret = ((array([[ -31.1141381 -6.46012214e-13j, -10.8963297 -4.14218588e-09j,
-157.68589005+1.84081236e-08j, 10.33...]]), array([1.+0.j, 0.+0.j]), array([1., 0.]), array([[1., 0.],
[0., 1.]]), array([[1., 0.],
[0., 1.]])))
reti = (array([[ -31.1141381 -6.46012214e-13j, -10.8963297 -4.14218588e-09j,
-157.68589005+1.84081236e-08j, 10.335...8402-1.38261738e-15j, -0.07905359+3.10199171e-12j,
0.26549587-3.32887046e-11j, 0.75392979-7.45269858e-01j]]))
self = <scipy.linalg.tests.test_decomp.TestOrdQZ object at 0x1497edd00>
sort = 'lhp'
scipy/linalg/tests/test_decomp.py:2316: in check
assert_array_almost_equal(Q @ Q.T.conj(), Id)
E AssertionError:
E Arrays are not almost equal to 6 decimals
E
E Mismatched elements: 16 / 16 (100%)
E Max absolute difference: 0.98869617
E Max relative difference: 0.98869617
E x: array([[ 0.011304-2.726199e-19j, 0.009174-4.493415e-18j,
E 0.010538-1.450496e-18j, 0.005477-1.979775e-17j],
E [ 0.009174+4.636163e-18j, 0.012032+1.067203e-19j,...
E y: array([[1., 0., 0., 0.],
E [0., 1., 0., 0.],
E [0., 0., 1., 0.],
E [0., 0., 0., 1.]])
A = array([[-21.1 -22.5j , 53.5 -50.5j , -34.5 +127.5j , 7.5 +0.5j ],
[ -0.46 -7.78j, -3.5 -37.5j , -15.5....7 -17.1j , -68.5 +12.5j , -7.5 -3.5j ],
[ 5.5 +4.4j , 14.4 +43.3j , -32.5 -46.j , -19. -32.5j ]])
AA = array([[ -31.1141381 -6.46012214e-13j, -10.8963297 -4.14218588e-09j,
-157.68589005+1.84081236e-08j, 10.3350...0.00000000e+00j, 0. +0.00000000e+00j,
0. +0.00000000e+00j, 22.19705256-2.19420886e+01j]])
B = array([[ 1. -5.j , 1.6+1.2j, -3. +0.j , 0. -1.j ],
[ 0.8-0.6j, 0. -5.j , -4. +3.j , -2.4-3.2j],
[ 1. +0.j , 2.4+1.8j, -4. -5.j , 0. -3.j ],
[ 0. +1.j , -1.8+2.4j, 0. -4.j , 4. -5.j ]])
BB = array([[ 0.83578309+0.00000000e+00j, -4.65649438-2.87034546e-10j,
-4.33106291+4.03125777e-10j, -1.77606997+1.7... +0.00000000e+00j, 0. +0.00000000e+00j,
0. +0.00000000e+00j, 5.05710627+0.00000000e+00j]])
Id = array([[1., 0., 0., 0.],
[0., 1., 0., 0.],
[0., 0., 1., 0.],
[0., 0., 0., 1.]])
Q = array([[ 7.46840448e-02+7.46840448e-02j, -5.77557803e-03-5.77557803e-03j,
4.80710267e-03+4.80710267e-03j, -4....2355e-05j, -1.50548800e-03-1.50548800e-03j,
3.29492283e-04+3.29492283e-04j, -6.49405396e-01-6.49405396e-01j]])
Z = array([[ 0.05099366-3.58927444e-16j, -0.55491431-2.82923687e-11j,
-0.26510151+1.83312152e-11j, 0.00117845-1.1...38402-1.38261738e-15j, -0.07905359+3.10199171e-12j,
0.26549587-3.32887046e-11j, 0.75392979-7.45269858e-01j]])
alpha = array([-31.1141381 -6.46012214e-13j, -7.22439206+5.68310534e-10j,
11.80413294+3.74449419e-15j, 22.19705256-2.19420886e+01j])
beta = array([0.83578309+0.j, 3.14087286+0.j, 0.09375939+0.j, 5.05710627+0.j])
self = <scipy.linalg.tests.test_decomp.TestOrdQZ object at 0x1497edd00>
sort = 'lhp'
______________________________ TestOrdQZ.test_rhp ______________________________
scipy/linalg/tests/test_decomp.py:2369: in test_rhp
self.check_all('rhp')
self = <scipy.linalg.tests.test_decomp.TestOrdQZ object at 0x1497ed5b0>
scipy/linalg/tests/test_decomp.py:2363: in check_all
self.check(Ai, Bi, sort, *reti)
Ai = array([[-21.1 -22.5j , 53.5 -50.5j , -34.5 +127.5j , 7.5 +0.5j ],
[ -0.46 -7.78j, -3.5 -37.5j , -15.5....7 -17.1j , -68.5 +12.5j , -7.5 -3.5j ],
[ 5.5 +4.4j , 14.4 +43.3j , -32.5 -46.j , -19. -32.5j ]])
Bi = array([[ 1. -5.j , 1.6+1.2j, -3. +0.j , 0. -1.j ],
[ 0.8-0.6j, 0. -5.j , -4. +3.j , -2.4-3.2j],
[ 1. +0.j , 2.4+1.8j, -4. -5.j , 0. -3.j ],
[ 0. +1.j , -1.8+2.4j, 0. -4.j , 4. -5.j ]])
ret = ((array([[134.39477527+4.26325641e-14j, -8.94635726+1.86775391e+01j,
69.27716232-6.76779276e+00j, 53.836549...]]), array([1.+0.j, 0.+0.j]), array([1., 0.]), array([[1., 0.],
[0., 1.]]), array([[1., 0.],
[0., 1.]])))
reti = (array([[134.39477527+4.26325641e-14j, -8.94635726+1.86775391e+01j,
69.27716232-6.76779276e+00j, 53.8365497...5001+1.38777878e-16j, 0.06082859-1.73544744e-01j,
-0.59212176+3.79578056e-01j, 0.57441682-5.35857893e-01j]]))
self = <scipy.linalg.tests.test_decomp.TestOrdQZ object at 0x1497ed5b0>
sort = 'rhp'
scipy/linalg/tests/test_decomp.py:2316: in check
assert_array_almost_equal(Q @ Q.T.conj(), Id)
E AssertionError:
E Arrays are not almost equal to 6 decimals
E
E Mismatched elements: 16 / 16 (100%)
E Max absolute difference: 0.98869617
E Max relative difference: 0.98869617
E x: array([[ 0.011304-3.547364e-19j, 0.009174-4.922254e-18j,
E 0.010538-2.343044e-18j, 0.005477-1.931867e-17j],
E [ 0.009174+5.530040e-18j, 0.012032+1.565478e-19j,...
E y: array([[1., 0., 0., 0.],
E [0., 1., 0., 0.],
E [0., 0., 1., 0.],
E [0., 0., 0., 1.]])
A = array([[-21.1 -22.5j , 53.5 -50.5j , -34.5 +127.5j , 7.5 +0.5j ],
[ -0.46 -7.78j, -3.5 -37.5j , -15.5....7 -17.1j , -68.5 +12.5j , -7.5 -3.5j ],
[ 5.5 +4.4j , 14.4 +43.3j , -32.5 -46.j , -19. -32.5j ]])
AA = array([[134.39477527+4.26325641e-14j, -8.94635726+1.86775391e+01j,
69.27716232-6.76779276e+00j, 53.83654971... +0.00000000e+00j, 0. +0.00000000e+00j,
0. +0.00000000e+00j, -7.2906963 +5.73525227e-10j]])
B = array([[ 1. -5.j , 1.6+1.2j, -3. +0.j , 0. -1.j ],
[ 0.8-0.6j, 0. -5.j , -4. +3.j , -2.4-3.2j],
[ 1. +0.j , 2.4+1.8j, -4. -5.j , 0. -3.j ],
[ 0. +1.j , -1.8+2.4j, 0. -4.j , 4. -5.j ]])
BB = array([[ 1.06748818+0.j , -1.51349933+0.6342589j ,
2.50124443+0.68327563j, 5.93919865-2.79324658j],
...,
[ 0. +0.j , 0. +0.j ,
0. +0.j , 3.16969926+0.j ]])
Id = array([[1., 0., 0., 0.],
[0., 1., 0., 0.],
[0., 0., 1., 0.],
[0., 0., 0., 1.]])
Q = array([[-7.43612701e-02-7.43612701e-02j, -3.70941403e-04-4.54338691e-03j,
-5.87340373e-03+6.80560633e-03j, -1....6865e-18j, -7.14045418e-01+2.51211141e-01j,
-3.82062145e-01-7.77407639e-03j, -2.69248399e-01-2.27958306e-01j]])
Z = array([[ 0.06688086+2.77555756e-17j, 0.05007278-1.10312020e-01j,
-0.41418762+9.10066489e-03j, -0.41835419+1.2...05001+1.38777878e-16j, 0.06082859-1.73544744e-01j,
-0.59212176+3.79578056e-01j, 0.57441682-5.35857893e-01j]])
alpha = array([134.39477527+4.26325641e-14j, 4.98648107-4.92920441e+00j,
-12.05432554-2.48689958e-13j, -7.2906963 +5.73525227e-10j])
beta = array([1.06748818+0.j, 1.13605915+0.j, 0.3238014 +0.j, 3.16969926+0.j])
self = <scipy.linalg.tests.test_decomp.TestOrdQZ object at 0x1497ed5b0>
sort = 'rhp'
______________________________ TestOrdQZ.test_iuc ______________________________
scipy/linalg/tests/test_decomp.py:2372: in test_iuc
self.check_all('iuc')
self = <scipy.linalg.tests.test_decomp.TestOrdQZ object at 0x1497dee80>
scipy/linalg/tests/test_decomp.py:2363: in check_all
self.check(Ai, Bi, sort, *reti)
Ai = array([[-21.1 -22.5j , 53.5 -50.5j , -34.5 +127.5j , 7.5 +0.5j ],
[ -0.46 -7.78j, -3.5 -37.5j , -15.5....7 -17.1j , -68.5 +12.5j , -7.5 -3.5j ],
[ 5.5 +4.4j , 14.4 +43.3j , -32.5 -46.j , -19. -32.5j ]])
Bi = array([[ 1. -5.j , 1.6+1.2j, -3. +0.j , 0. -1.j ],
[ 0.8-0.6j, 0. -5.j , -4. +3.j , -2.4-3.2j],
[ 1. +0.j , 2.4+1.8j, -4. -5.j , 0. -3.j ],
[ 0. +1.j , -1.8+2.4j, 0. -4.j , 4. -5.j ]])
ret = ((array([[134.39477527+4.26325641e-14j, -7.48252835-1.22568622e-13j,
89.4978431 -7.04037702e-09j, -10.570688...]]), array([1.+0.j, 0.+0.j]), array([1., 0.]), array([[1., 0.],
[0., 1.]]), array([[1., 0.],
[0., 1.]])))
reti = (array([[134.39477527+4.26325641e-14j, -7.48252835-1.22568622e-13j,
89.4978431 -7.04037702e-09j, -10.5706889...5001+1.38777878e-16j, 0.13130749+2.80331314e-15j,
0.06726816-5.29180866e-12j, 0.75392979-7.45269858e-01j]]))
self = <scipy.linalg.tests.test_decomp.TestOrdQZ object at 0x1497dee80>
sort = 'iuc'
scipy/linalg/tests/test_decomp.py:2316: in check
assert_array_almost_equal(Q @ Q.T.conj(), Id)
E AssertionError:
E Arrays are not almost equal to 6 decimals
E
E Mismatched elements: 16 / 16 (100%)
E Max absolute difference: 0.98869617
E Max relative difference: 0.98869617
E x: array([[ 0.011304-3.505470e-19j, 0.009174-5.035516e-18j,
E 0.010538-2.263648e-18j, 0.005477-1.979775e-17j],
E [ 0.009174+5.557735e-18j, 0.012032+1.067203e-19j,...
E y: array([[1., 0., 0., 0.],
E [0., 1., 0., 0.],
E [0., 0., 1., 0.],
E [0., 0., 0., 1.]])
A = array([[-21.1 -22.5j , 53.5 -50.5j , -34.5 +127.5j , 7.5 +0.5j ],
[ -0.46 -7.78j, -3.5 -37.5j , -15.5....7 -17.1j , -68.5 +12.5j , -7.5 -3.5j ],
[ 5.5 +4.4j , 14.4 +43.3j , -32.5 -46.j , -19. -32.5j ]])
AA = array([[134.39477527+4.26325641e-14j, -7.48252835-1.22568622e-13j,
89.4978431 -7.04037702e-09j, -10.57068897... +0.00000000e+00j, 0. +0.00000000e+00j,
0. +0.00000000e+00j, 22.19705256-2.19420886e+01j]])
B = array([[ 1. -5.j , 1.6+1.2j, -3. +0.j , 0. -1.j ],
[ 0.8-0.6j, 0. -5.j , -4. +3.j , -2.4-3.2j],
[ 1. +0.j , 2.4+1.8j, -4. -5.j , 0. -3.j ],
[ 0. +1.j , -1.8+2.4j, 0. -4.j , 4. -5.j ]])
BB = array([[ 1.06748818+0.00000000e+00j, -1.18208167-2.42028619e-14j,
6.47769615-5.09570209e-10j, 2.14942102-2.1... +0.00000000e+00j, 0. +0.00000000e+00j,
0. +0.00000000e+00j, 5.05710627+0.00000000e+00j]])
Id = array([[1., 0., 0., 0.],
[0., 1., 0., 0.],
[0., 0., 1., 0.],
[0., 0., 0., 1.]])
Q = array([[-7.43612701e-02-7.43612701e-02j, 7.01704469e-03+7.01704469e-03j,
-7.43872879e-03-7.43872879e-03j, -4....6865e-18j, -9.01568861e-04-9.01568861e-04j,
-1.25368789e-03-1.25368789e-03j, -6.49405396e-01-6.49405396e-01j]])
Z = array([[ 0.06688086+2.77555756e-17j, 0.03340062+7.07767178e-16j,
-0.61255237+4.81867324e-11j, 0.00117845-1.1...05001+1.38777878e-16j, 0.13130749+2.80331314e-15j,
0.06726816-5.29180866e-12j, 0.75392979-7.45269858e-01j]])
alpha = array([134.39477527+4.26325641e-14j, -3.50777004-7.28306304e-14j,
-5.62833156+4.42755610e-10j, 22.19705256-2.19420886e+01j])
beta = array([1.06748818+0.j, 0.09422517+0.j, 2.44697045+0.j, 5.05710627+0.j])
self = <scipy.linalg.tests.test_decomp.TestOrdQZ object at 0x1497dee80>
sort = 'iuc'
______________________________ TestOrdQZ.test_ouc ______________________________
scipy/linalg/tests/test_decomp.py:2375: in test_ouc
self.check_all('ouc')
self = <scipy.linalg.tests.test_decomp.TestOrdQZ object at 0x1497efc70>
scipy/linalg/tests/test_decomp.py:2363: in check_all
self.check(Ai, Bi, sort, *reti)
Ai = array([[-21.1 -22.5j , 53.5 -50.5j , -34.5 +127.5j , 7.5 +0.5j ],
[ -0.46 -7.78j, -3.5 -37.5j , -15.5....7 -17.1j , -68.5 +12.5j , -7.5 -3.5j ],
[ 5.5 +4.4j , 14.4 +43.3j , -32.5 -46.j , -19. -32.5j ]])
Bi = array([[ 1. -5.j , 1.6+1.2j, -3. +0.j , 0. -1.j ],
[ 0.8-0.6j, 0. -5.j , -4. +3.j , -2.4-3.2j],
[ 1. +0.j , 2.4+1.8j, -4. -5.j , 0. -3.j ],
[ 0. +1.j , -1.8+2.4j, 0. -4.j , 4. -5.j ]])
ret = ((array([[134.39477527+4.26325641e-14j, -7.48252835-1.22568622e-13j,
89.4978431 -7.04037702e-09j, -10.570688...]]), array([1.+0.j, 0.+0.j]), array([1., 0.]), array([[1., 0.],
[0., 1.]]), array([[1., 0.],
[0., 1.]])))
reti = (array([[134.39477527+4.26325641e-14j, -7.48252835-1.22568622e-13j,
89.4978431 -7.04037702e-09j, -10.5706889...5001+1.38777878e-16j, 0.13130749+2.80331314e-15j,
0.06726816-5.29180866e-12j, 0.75392979-7.45269858e-01j]]))
self = <scipy.linalg.tests.test_decomp.TestOrdQZ object at 0x1497efc70>
sort = 'ouc'
scipy/linalg/tests/test_decomp.py:2316: in check
assert_array_almost_equal(Q @ Q.T.conj(), Id)
E AssertionError:
E Arrays are not almost equal to 6 decimals
E
E Mismatched elements: 16 / 16 (100%)
E Max absolute difference: 0.98869617
E Max relative difference: 0.98869617
E x: array([[ 0.011304-3.505470e-19j, 0.009174-5.035516e-18j,
E 0.010538-2.263648e-18j, 0.005477-1.979775e-17j],
E [ 0.009174+5.557735e-18j, 0.012032+1.067203e-19j,...
E y: array([[1., 0., 0., 0.],
E [0., 1., 0., 0.],
E [0., 0., 1., 0.],
E [0., 0., 0., 1.]])
A = array([[-21.1 -22.5j , 53.5 -50.5j , -34.5 +127.5j , 7.5 +0.5j ],
[ -0.46 -7.78j, -3.5 -37.5j , -15.5....7 -17.1j , -68.5 +12.5j , -7.5 -3.5j ],
[ 5.5 +4.4j , 14.4 +43.3j , -32.5 -46.j , -19. -32.5j ]])
AA = array([[134.39477527+4.26325641e-14j, -7.48252835-1.22568622e-13j,
89.4978431 -7.04037702e-09j, -10.57068897... +0.00000000e+00j, 0. +0.00000000e+00j,
0. +0.00000000e+00j, 22.19705256-2.19420886e+01j]])
B = array([[ 1. -5.j , 1.6+1.2j, -3. +0.j , 0. -1.j ],
[ 0.8-0.6j, 0. -5.j , -4. +3.j , -2.4-3.2j],
[ 1. +0.j , 2.4+1.8j, -4. -5.j , 0. -3.j ],
[ 0. +1.j , -1.8+2.4j, 0. -4.j , 4. -5.j ]])
BB = array([[ 1.06748818+0.00000000e+00j, -1.18208167-2.42028619e-14j,
6.47769615-5.09570209e-10j, 2.14942102-2.1... +0.00000000e+00j, 0. +0.00000000e+00j,
0. +0.00000000e+00j, 5.05710627+0.00000000e+00j]])
Id = array([[1., 0., 0., 0.],
[0., 1., 0., 0.],
[0., 0., 1., 0.],
[0., 0., 0., 1.]])
Q = array([[-7.43612701e-02-7.43612701e-02j, 7.01704469e-03+7.01704469e-03j,
-7.43872879e-03-7.43872879e-03j, -4....6865e-18j, -9.01568861e-04-9.01568861e-04j,
-1.25368789e-03-1.25368789e-03j, -6.49405396e-01-6.49405396e-01j]])
Z = array([[ 0.06688086+2.77555756e-17j, 0.03340062+7.07767178e-16j,
-0.61255237+4.81867324e-11j, 0.00117845-1.1...05001+1.38777878e-16j, 0.13130749+2.80331314e-15j,
0.06726816-5.29180866e-12j, 0.75392979-7.45269858e-01j]])
alpha = array([134.39477527+4.26325641e-14j, -3.50777004-7.28306304e-14j,
-5.62833156+4.42755610e-10j, 22.19705256-2.19420886e+01j])
beta = array([1.06748818+0.j, 0.09422517+0.j, 2.44697045+0.j, 5.05710627+0.j])
self = <scipy.linalg.tests.test_decomp.TestOrdQZ object at 0x1497efc70>
sort = 'ouc'
______________________________ TestOrdQZ.test_ref ______________________________
scipy/linalg/tests/test_decomp.py:2386: in test_ref
self.check_all(sort)
self = <scipy.linalg.tests.test_decomp.TestOrdQZ object at 0x1495d0b20>
sort = <function TestOrdQZ.test_ref.<locals>.sort at 0x149353940>
scipy/linalg/tests/test_decomp.py:2363: in check_all
self.check(Ai, Bi, sort, *reti)
Ai = array([[-21.1 -22.5j , 53.5 -50.5j , -34.5 +127.5j , 7.5 +0.5j ],
[ -0.46 -7.78j, -3.5 -37.5j , -15.5....7 -17.1j , -68.5 +12.5j , -7.5 -3.5j ],
[ 5.5 +4.4j , 14.4 +43.3j , -32.5 -46.j , -19. -32.5j ]])
Bi = array([[ 1. -5.j , 1.6+1.2j, -3. +0.j , 0. -1.j ],
[ 0.8-0.6j, 0. -5.j , -4. +3.j , -2.4-3.2j],
[ 1. +0.j , 2.4+1.8j, -4. -5.j , 0. -3.j ],
[ 0. +1.j , -1.8+2.4j, 0. -4.j , 4. -5.j ]])
ret = ((array([[134.39477527+4.26325641e-14j, -7.48252835-1.22568622e-13j,
89.4978431 -7.04037702e-09j, -10.570688...]]), array([1.+0.j, 0.+0.j]), array([1., 0.]), array([[1., 0.],
[0., 1.]]), array([[1., 0.],
[0., 1.]])))
reti = (array([[134.39477527+4.26325641e-14j, -7.48252835-1.22568622e-13j,
89.4978431 -7.04037702e-09j, -10.5706889...5001+1.38777878e-16j, 0.13130749+2.80331314e-15j,
0.06726816-5.29180866e-12j, 0.75392979-7.45269858e-01j]]))
self = <scipy.linalg.tests.test_decomp.TestOrdQZ object at 0x1495d0b20>
sort = <function TestOrdQZ.test_ref.<locals>.sort at 0x149353940>
scipy/linalg/tests/test_decomp.py:2316: in check
assert_array_almost_equal(Q @ Q.T.conj(), Id)
E AssertionError:
E Arrays are not almost equal to 6 decimals
E
E Mismatched elements: 16 / 16 (100%)
E Max absolute difference: 0.98869617
E Max relative difference: 0.98869617
E x: array([[ 0.011304-3.505470e-19j, 0.009174-5.035516e-18j,
E 0.010538-2.263648e-18j, 0.005477-1.979775e-17j],
E [ 0.009174+5.557735e-18j, 0.012032+1.067203e-19j,...
E y: array([[1., 0., 0., 0.],
E [0., 1., 0., 0.],
E [0., 0., 1., 0.],
E [0., 0., 0., 1.]])
A = array([[-21.1 -22.5j , 53.5 -50.5j , -34.5 +127.5j , 7.5 +0.5j ],
[ -0.46 -7.78j, -3.5 -37.5j , -15.5....7 -17.1j , -68.5 +12.5j , -7.5 -3.5j ],
[ 5.5 +4.4j , 14.4 +43.3j , -32.5 -46.j , -19. -32.5j ]])
AA = array([[134.39477527+4.26325641e-14j, -7.48252835-1.22568622e-13j,
89.4978431 -7.04037702e-09j, -10.57068897... +0.00000000e+00j, 0. +0.00000000e+00j,
0. +0.00000000e+00j, 22.19705256-2.19420886e+01j]])
B = array([[ 1. -5.j , 1.6+1.2j, -3. +0.j , 0. -1.j ],
[ 0.8-0.6j, 0. -5.j , -4. +3.j , -2.4-3.2j],
[ 1. +0.j , 2.4+1.8j, -4. -5.j , 0. -3.j ],
[ 0. +1.j , -1.8+2.4j, 0. -4.j , 4. -5.j ]])
BB = array([[ 1.06748818+0.00000000e+00j, -1.18208167-2.42028619e-14j,
6.47769615-5.09570209e-10j, 2.14942102-2.1... +0.00000000e+00j, 0. +0.00000000e+00j,
0. +0.00000000e+00j, 5.05710627+0.00000000e+00j]])
Id = array([[1., 0., 0., 0.],
[0., 1., 0., 0.],
[0., 0., 1., 0.],
[0., 0., 0., 1.]])
Q = array([[-7.43612701e-02-7.43612701e-02j, 7.01704469e-03+7.01704469e-03j,
-7.43872879e-03-7.43872879e-03j, -4....6865e-18j, -9.01568861e-04-9.01568861e-04j,
-1.25368789e-03-1.25368789e-03j, -6.49405396e-01-6.49405396e-01j]])
Z = array([[ 0.06688086+2.77555756e-17j, 0.03340062+7.07767178e-16j,
-0.61255237+4.81867324e-11j, 0.00117845-1.1...05001+1.38777878e-16j, 0.13130749+2.80331314e-15j,
0.06726816-5.29180866e-12j, 0.75392979-7.45269858e-01j]])
alpha = array([134.39477527+4.26325641e-14j, -3.50777004-7.28306304e-14j,
-5.62833156+4.42755610e-10j, 22.19705256-2.19420886e+01j])
beta = array([1.06748818+0.j, 0.09422517+0.j, 2.44697045+0.j, 5.05710627+0.j])
self = <scipy.linalg.tests.test_decomp.TestOrdQZ object at 0x1495d0b20>
sort = <function TestOrdQZ.test_ref.<locals>.sort at 0x149353940>
______________________________ TestOrdQZ.test_cef ______________________________
scipy/linalg/tests/test_decomp.py:2397: in test_cef
self.check_all(sort)
self = <scipy.linalg.tests.test_decomp.TestOrdQZ object at 0x1495d0cd0>
sort = <function TestOrdQZ.test_cef.<locals>.sort at 0x149353f70>
scipy/linalg/tests/test_decomp.py:2363: in check_all
self.check(Ai, Bi, sort, *reti)
Ai = array([[-21.1 -22.5j , 53.5 -50.5j , -34.5 +127.5j , 7.5 +0.5j ],
[ -0.46 -7.78j, -3.5 -37.5j , -15.5....7 -17.1j , -68.5 +12.5j , -7.5 -3.5j ],
[ 5.5 +4.4j , 14.4 +43.3j , -32.5 -46.j , -19. -32.5j ]])
Bi = array([[ 1. -5.j , 1.6+1.2j, -3. +0.j , 0. -1.j ],
[ 0.8-0.6j, 0. -5.j , -4. +3.j , -2.4-3.2j],
[ 1. +0.j , 2.4+1.8j, -4. -5.j , 0. -3.j ],
[ 0. +1.j , -1.8+2.4j, 0. -4.j , 4. -5.j ]])
ret = ((array([[134.39477527+4.26325641e-14j, -7.48252835-1.22568622e-13j,
89.4978431 -7.04037702e-09j, -10.570688...]]), array([1.+0.j, 0.+0.j]), array([1., 0.]), array([[1., 0.],
[0., 1.]]), array([[1., 0.],
[0., 1.]])))
reti = (array([[134.39477527+4.26325641e-14j, -7.48252835-1.22568622e-13j,
89.4978431 -7.04037702e-09j, -10.5706889...5001+1.38777878e-16j, 0.13130749+2.80331314e-15j,
0.06726816-5.29180866e-12j, 0.75392979-7.45269858e-01j]]))
self = <scipy.linalg.tests.test_decomp.TestOrdQZ object at 0x1495d0cd0>
sort = <function TestOrdQZ.test_cef.<locals>.sort at 0x149353f70>
scipy/linalg/tests/test_decomp.py:2316: in check
assert_array_almost_equal(Q @ Q.T.conj(), Id)
E AssertionError:
E Arrays are not almost equal to 6 decimals
E
E Mismatched elements: 16 / 16 (100%)
E Max absolute difference: 0.98869617
E Max relative difference: 0.98869617
E x: array([[ 0.011304-3.505470e-19j, 0.009174-5.035516e-18j,
E 0.010538-2.263648e-18j, 0.005477-1.979775e-17j],
E [ 0.009174+5.557735e-18j, 0.012032+1.067203e-19j,...
E y: array([[1., 0., 0., 0.],
E [0., 1., 0., 0.],
E [0., 0., 1., 0.],
E [0., 0., 0., 1.]])
A = array([[-21.1 -22.5j , 53.5 -50.5j , -34.5 +127.5j , 7.5 +0.5j ],
[ -0.46 -7.78j, -3.5 -37.5j , -15.5....7 -17.1j , -68.5 +12.5j , -7.5 -3.5j ],
[ 5.5 +4.4j , 14.4 +43.3j , -32.5 -46.j , -19. -32.5j ]])
AA = array([[134.39477527+4.26325641e-14j, -7.48252835-1.22568622e-13j,
89.4978431 -7.04037702e-09j, -10.57068897... +0.00000000e+00j, 0. +0.00000000e+00j,
0. +0.00000000e+00j, 22.19705256-2.19420886e+01j]])
B = array([[ 1. -5.j , 1.6+1.2j, -3. +0.j , 0. -1.j ],
[ 0.8-0.6j, 0. -5.j , -4. +3.j , -2.4-3.2j],
[ 1. +0.j , 2.4+1.8j, -4. -5.j , 0. -3.j ],
[ 0. +1.j , -1.8+2.4j, 0. -4.j , 4. -5.j ]])
BB = array([[ 1.06748818+0.00000000e+00j, -1.18208167-2.42028619e-14j,
6.47769615-5.09570209e-10j, 2.14942102-2.1... +0.00000000e+00j, 0. +0.00000000e+00j,
0. +0.00000000e+00j, 5.05710627+0.00000000e+00j]])
Id = array([[1., 0., 0., 0.],
[0., 1., 0., 0.],
[0., 0., 1., 0.],
[0., 0., 0., 1.]])
Q = array([[-7.43612701e-02-7.43612701e-02j, 7.01704469e-03+7.01704469e-03j,
-7.43872879e-03-7.43872879e-03j, -4....6865e-18j, -9.01568861e-04-9.01568861e-04j,
-1.25368789e-03-1.25368789e-03j, -6.49405396e-01-6.49405396e-01j]])
Z = array([[ 0.06688086+2.77555756e-17j, 0.03340062+7.07767178e-16j,
-0.61255237+4.81867324e-11j, 0.00117845-1.1...05001+1.38777878e-16j, 0.13130749+2.80331314e-15j,
0.06726816-5.29180866e-12j, 0.75392979-7.45269858e-01j]])
alpha = array([134.39477527+4.26325641e-14j, -3.50777004-7.28306304e-14j,
-5.62833156+4.42755610e-10j, 22.19705256-2.19420886e+01j])
beta = array([1.06748818+0.j, 0.09422517+0.j, 2.44697045+0.j, 5.05710627+0.j])
self = <scipy.linalg.tests.test_decomp.TestOrdQZ object at 0x1495d0cd0>
sort = <function TestOrdQZ.test_cef.<locals>.sort at 0x149353f70>
____________________ TestOrdQZWorkspaceSize.test_decompose _____________________
scipy/linalg/tests/test_decomp.py:2474: in test_decompose
_ = ordqz(A, B, sort=lambda alpha, beta: alpha < beta,
A = array([[0.3400652 +0.j, 0.40679654+0.j, 0.03785227+0.j, ...,
0.46569912+0.j, 0.34000689+0.j, 0.18197288+0.j],
...
[0.35160158+0.j, 0.41201777+0.j, 0.68358587+0.j, ...,
0.50537399+0.j, 0.08104583+0.j, 0.39046819+0.j]])
B = array([[0.01923255+0.j, 0.58262855+0.j, 0.18256915+0.j, ...,
0.98463685+0.j, 0.60828429+0.j, 0.56107299+0.j],
...
[0.3774488 +0.j, 0.22516928+0.j, 0.08706877+0.j, ...,
0.9015293 +0.j, 0.18932018+0.j, 0.78947954+0.j]])
N = 202
_ = (array([[-4.07674921e+00, -3.65758548e-01, 6.84469640e-03, ...,
1.01177726e+00, -5.76916826e+00, 2.12013153...28, -0.10160696],
[-0.02071906, -0.03215941, -0.04681432, ..., 0.01070236,
0.07517671, 0.03073703]]))
ddtype = <class 'numpy.complex128'>
self = <scipy.linalg.tests.test_decomp.TestOrdQZWorkspaceSize object at 0x1497ddd30>
scipy/linalg/_decomp_qz.py:353: in ordqz
(AA, BB, _, *ab, Q, Z, _, _), typ = _qz(A, B, output=output, sort=None,
A = array([[0.3400652 +0.j, 0.40679654+0.j, 0.03785227+0.j, ...,
0.46569912+0.j, 0.34000689+0.j, 0.18197288+0.j],
...
[0.35160158+0.j, 0.41201777+0.j, 0.68358587+0.j, ...,
0.50537399+0.j, 0.08104583+0.j, 0.39046819+0.j]])
B = array([[0.01923255+0.j, 0.58262855+0.j, 0.18256915+0.j, ...,
0.98463685+0.j, 0.60828429+0.j, 0.56107299+0.j],
...
[0.3774488 +0.j, 0.22516928+0.j, 0.08706877+0.j, ...,
0.9015293 +0.j, 0.18932018+0.j, 0.78947954+0.j]])
check_finite = True
output = 'complex'
overwrite_a = False
overwrite_b = False
sort = <function TestOrdQZWorkspaceSize.test_decompose.<locals>.<lambda> at 0x149353670>
scipy/linalg/_decomp_qz.py:127: in _qz
warnings.warn("The QZ iteration failed. (a,b) are not in Schur "
E scipy.linalg.misc.LinAlgWarning: The QZ iteration failed. (a,b) are not in Schur form, but ALPHAR(j), ALPHAI(j), and BETA(j) should be correct for J=98,...,N
A = array([[0.3400652 +0.j, 0.40679654+0.j, 0.03785227+0.j, ...,
0.46569912+0.j, 0.34000689+0.j, 0.18197288+0.j],
...
[0.35160158+0.j, 0.41201777+0.j, 0.68358587+0.j, ...,
0.50537399+0.j, 0.08104583+0.j, 0.39046819+0.j]])
B = array([[0.01923255+0.j, 0.58262855+0.j, 0.18256915+0.j, ...,
0.98463685+0.j, 0.60828429+0.j, 0.56107299+0.j],
...
[0.3774488 +0.j, 0.22516928+0.j, 0.08706877+0.j, ...,
0.9015293 +0.j, 0.18932018+0.j, 0.78947954+0.j]])
a1 = array([[0.3400652 +0.j, 0.40679654+0.j, 0.03785227+0.j, ...,
0.46569912+0.j, 0.34000689+0.j, 0.18197288+0.j],
...
[0.35160158+0.j, 0.41201777+0.j, 0.68358587+0.j, ...,
0.50537399+0.j, 0.08104583+0.j, 0.39046819+0.j]])
a_m = 202
a_n = 202
b1 = array([[0.01923255+0.j, 0.58262855+0.j, 0.18256915+0.j, ...,
0.98463685+0.j, 0.60828429+0.j, 0.56107299+0.j],
...
[0.3774488 +0.j, 0.22516928+0.j, 0.08706877+0.j, ...,
0.9015293 +0.j, 0.18932018+0.j, 0.78947954+0.j]])
b_m = 202
b_n = 202
check_finite = True
gges = <fortran object>
info = 99
lwork = 6666
output = 'complex'
overwrite_a = False
overwrite_b = False
result = (array([[ 2.52544948e+00+2.52544986e+00j, -3.15241315e-01-3.15241315e-01j,
3.33983793e-01+3.33983793e-01j, ..... -1.82661396e-10-1.82661396e-10j, -3.33609739e-08-3.33609739e-08j,
-4.46818968e-02-4.46818968e-02j]]), ...)
sfunction = <function _qz.<locals>.<lambda> at 0x1493535e0>
sort = None
typa = 'D'
typb = 'D'
______________________ test_cossin[True-4-2-2-complex64] _______________________
scipy/linalg/tests/test_decomp_cossin.py:39: in test_cossin
assert_allclose(x, u @ cs @ vh, rtol=0., atol=m*1e3*np.finfo(dtype_).eps)
E AssertionError:
E Not equal to tolerance rtol=0, atol=0.000476837
E
E Mismatched elements: 16 / 16 (100%)
E Max absolute difference: 0.984457
E Max relative difference: 43.62083
E x: array([[ 0.247539+0.212893j, -0.097059+0.043832j, 0.495105+0.6285j ,
E 0.305295+0.385632j],
E [-0.378363-0.106404j, -0.07387 +0.060195j, 0.551738-0.08671j ,...
E y: array([[ 0.247539+0.250358j, -0.078446-0.101413j, -0.40529 +1.026551j,
E 0.196049+0.111241j],
E [-0.378363-0.376737j, -0.056753-0.073382j, -0.062781+0.191695j,...
cs = array([[ 0.5111532 , 0. , 0.8594896 , 0. ],
[ 0. , 0.14221044, 0. , 0.9898364... , 0.5111532 , 0. ],
[-0. , -0.98983645, 0. , 0.14221044]],
dtype=float32)
dtype_ = <class 'numpy.complex64'>
m = 4
p = 2
q = 2
swap_sign = True
u = array([[ 0.4789903 +0.48450702j, 0.68741924+0.68741924j,
0. +0.j , 0. +0.j ],
....j , 0. +0.j ,
-0.55586 -0.5809152j , -0.81093484-0.81093484j]],
dtype=complex64)
vh = array([[ 0.99958915+0.j , 0.02842978+0.00364315j,
0. +0.j , 0. +0.j ],
....j , 0. +0.j ,
0.3809979 +0.8601899j , -0.31575269+0.12334777j]],
dtype=complex64)
x = array([[ 0.24753878+0.21289341j, -0.09705902+0.04383202j,
0.49510536+0.6285002j , 0.30529493+0.3856325j ],
....0961457j , -0.76420313-0.24125303j,
-0.08124338-0.188753j , -0.19499378+0.13314982j]],
dtype=complex64)
______________________ test_cossin[True-4-2-2-complex128] ______________________
scipy/linalg/tests/test_decomp_cossin.py:39: in test_cossin
assert_allclose(x, u @ cs @ vh, rtol=0., atol=m*1e3*np.finfo(dtype_).eps)
E AssertionError:
E Not equal to tolerance rtol=0, atol=8.88178e-13
E
E Mismatched elements: 16 / 16 (100%)
E Max absolute difference: 0.98445686
E Max relative difference: 43.6209201
E x: array([[ 0.247539+0.212893j, -0.097059+0.043832j, 0.495105+0.6285j ,
E 0.305295+0.385632j],
E [-0.378363-0.106405j, -0.07387 +0.060195j, 0.551738-0.08671j ,...
E y: array([[ 0.247539+0.250358j, -0.078446-0.101413j, -0.405289+1.026551j,
E 0.196049+0.111241j],
E [-0.378363-0.376737j, -0.056753-0.073382j, -0.062781+0.191695j,...
cs = array([[ 0.51115322, 0. , 0.85948961, 0. ],
[ 0. , 0.14221044, 0. , 0.9898364... [-0.85948961, -0. , 0.51115322, 0. ],
[-0. , -0.98983645, 0. , 0.14221044]])
dtype_ = <class 'numpy.complex128'>
m = 4
p = 2