The R Project SVN R

Rev

Rev 89016 | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 89016 Rev 89125
Line 1116... Line 1116...
1116 
		warning(_("the matrix is either rank-deficient or not positive definite"));
 1116 
		warning(_("the matrix is either rank-deficient or not positive definite"));
1117 
	    else
 1117 
	    else
1118 
		error(_("argument %d of Lapack routine %s had invalid value"),
 1118 
		error(_("argument %d of Lapack routine %s had invalid value"),
1119 
		      -info, "dpstrf");
 1119 
		      -info, "dpstrf");
1120 
	}
 1120 
	}
- 
 
 1121 
 	/* zero remaining upper triangle if rank < m */
- 
 
 1122 
	for (int j = rank ; j < m ; j++)
- 
 
 1123 
            for (int i = rank ; i <= j ; i++)
- 
 
 1124 
                REAL(ans) [i + N * j] = 0. ;
1121 
	setAttrib(ans, install("pivot"), piv);
 1125 
	setAttrib(ans, install("pivot"), piv);
1122 
	SEXP s_rank = install("rank");
 1126 
	SEXP s_rank = install("rank");
1123 
	setAttrib(ans, s_rank, ScalarInteger(rank));
 1127 
	setAttrib(ans, s_rank, ScalarInteger(rank));
1124 
	SEXP cn, dn = getAttrib(ans, R_DimNamesSymbol);
 1128 
	SEXP cn, dn = getAttrib(ans, R_DimNamesSymbol);
1125 
	if (!isNull(dn) && !isNull(cn = VECTOR_ELT(dn, 1))) {
 1129 
	if (!isNull(dn) && !isNull(cn = VECTOR_ELT(dn, 1))) {