Mercurial > hg > octave-lyh

comparison liboctave/CSparse.cc @ 5877:50d43cdbec80

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
[project @ 2006-07-01 21:52:05 by dbateman]
author dbateman
date Sat, 01 Jul 2006 21:52:06 +0000
parents 565d0cd4d9d0
children d73ffe42f2c8
comparison
equal deleted inserted replaced
5876:565d0cd4d9d0 5877:50d43cdbec80
780 780
781 // get A(m,m) 781 // get A(m,m)
782 if (typ == MatrixType::Upper) 782 if (typ == MatrixType::Upper)
783 colUp = cidx(j+1) - 1; 783 colUp = cidx(j+1) - 1;
784 else 784 else
785 colUp = cidx(j) - 1; 785 colUp = cidx(j);
786 Complex pivot = data(colUp); 786 Complex pivot = data(colUp);
787 if (pivot == 0. || colUp != j) 787 if (pivot == 0. || ridx(colUp) != j)
788 { 788 {
789 (*current_liboctave_error_handler) 789 (*current_liboctave_error_handler)
790 ("division by zero"); 790 ("division by zero");
791 goto inverse_singular; 791 goto inverse_singular;
792 } 792 }
808 // get A(m,m) 808 // get A(m,m)
809 octave_idx_type colUp; 809 octave_idx_type colUp;
810 if (typ == MatrixType::Upper) 810 if (typ == MatrixType::Upper)
811 colUp = cidx(i+1) - 1; 811 colUp = cidx(i+1) - 1;
812 else 812 else
813 colUp = cidx(i) - 1; 813 colUp = cidx(i);
814 Complex pivot = data(colUp); 814 Complex pivot = data(colUp);
815 if (pivot == 0. || colUp != i) 815 if (pivot == 0. || ridx(colUp) != i)
816 { 816 {
817 (*current_liboctave_error_handler) ("division by zero"); 817 (*current_liboctave_error_handler) ("division by zero");
818 goto inverse_singular; 818 goto inverse_singular;
819 } 819 }
820 820
876 // get A(m,m) 876 // get A(m,m)
877 Complex pivot; 877 Complex pivot;
878 if (typ == MatrixType::Permuted_Upper) 878 if (typ == MatrixType::Permuted_Upper)
879 pivot = data(cidx(jidx+1) - 1); 879 pivot = data(cidx(jidx+1) - 1);
880 else 880 else
881 pivot = data(cidx(jidx) - 1); 881 pivot = data(cidx(jidx));
882 if (pivot == 0.) 882 if (pivot == 0.)
883 { 883 {
884 (*current_liboctave_error_handler) 884 (*current_liboctave_error_handler)
885 ("division by zero"); 885 ("division by zero");
886 goto inverse_singular; 886 goto inverse_singular;
892 // get A(m,m) 892 // get A(m,m)
893 octave_idx_type colUp; 893 octave_idx_type colUp;
894 if (typ == MatrixType::Permuted_Upper) 894 if (typ == MatrixType::Permuted_Upper)
895 colUp = cidx(perm[iidx]+1) - 1; 895 colUp = cidx(perm[iidx]+1) - 1;
896 else 896 else
897 colUp = cidx(perm[iidx]) - 1; 897 colUp = cidx(perm[iidx]);
898 898
899 Complex pivot = data(colUp); 899 Complex pivot = data(colUp);
900 if (pivot == 0.) 900 if (pivot == 0.)
901 { 901 {
902 (*current_liboctave_error_handler) 902 (*current_liboctave_error_handler)