Mercurial > octave-nkf

comparison liboctave/fCMatrix.cc @ 10630:832732c2eb80

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
use ID for singularity warning in matrix division
author Jaroslav Hajek <highegg@gmail.com>
date Thu, 13 May 2010 13:47:20 +0200
parents a0728e81ed25
children 2786b7ae3f80
comparison
equal deleted inserted replaced
10629:b7c4954e1c00 10630:832732c2eb80
2692 F77_XFCN (cgelsd, CGELSD, (m, n, nrhs, tmp_data, m, pretval, 2692 F77_XFCN (cgelsd, CGELSD, (m, n, nrhs, tmp_data, m, pretval,
2693 maxmn, ps, rcon, rank, 2693 maxmn, ps, rcon, rank,
2694 work.fortran_vec (), lwork, 2694 work.fortran_vec (), lwork,
2695 prwork, piwork, info)); 2695 prwork, piwork, info));
2696 2696
2697 if (rank < minmn)
2698 (*current_liboctave_warning_handler)
2699 ("zgelsd: rank deficient %dx%d matrix, rank = %d, tol = %e",
2700 m, n, rank, rcon);
2701
2702 if (s.elem (0) == 0.0) 2697 if (s.elem (0) == 0.0)
2703 rcon = 0.0; 2698 rcon = 0.0;
2704 else 2699 else
2705 rcon = s.elem (minmn - 1) / s.elem (0); 2700 rcon = s.elem (minmn - 1) / s.elem (0);
2706 2701
2860 work.fortran_vec (), lwork, 2855 work.fortran_vec (), lwork,
2861 prwork, piwork, info)); 2856 prwork, piwork, info));
2862 2857
2863 if (rank < minmn) 2858 if (rank < minmn)
2864 { 2859 {
2865 if (rank < minmn)
2866 (*current_liboctave_warning_handler)
2867 ("zgelsd: rank deficient %dx%d matrix, rank = %d, tol = %e",
2868 m, n, rank, rcon);
2869
2870 if (s.elem (0) == 0.0) 2860 if (s.elem (0) == 0.0)
2871 rcon = 0.0; 2861 rcon = 0.0;
2872 else 2862 else
2873 rcon = s.elem (minmn - 1) / s.elem (0); 2863 rcon = s.elem (minmn - 1) / s.elem (0);
2874 2864