The R Project SVN R

Rev

Rev 68947 | Rev 85063 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 68947 Rev 76006
Line 1... Line 1...
1
/*
1
/*
2
 *  R : A Computer Language for Statistical Data Analysis
2
 *  R : A Computer Language for Statistical Data Analysis
3
 *  Copyright (C) 2000-2015 The R Core Team
3
 *  Copyright (C) 2000-2019 The R Core Team
4
 *
4
 *
5
 *  This program is free software; you can redistribute it and/or modify
5
 *  This program is free software; you can redistribute it and/or modify
6
 *  it under the terms of the GNU General Public License as published by
6
 *  it under the terms of the GNU General Public License as published by
7
 *  the Free Software Foundation; either version 2 of the License, or
7
 *  the Free Software Foundation; either version 2 of the License, or
8
 *  (at your option) any later version.
8
 *  (at your option) any later version.
Line 869... Line 869...
869
	    theta, &rr, &dr, &stp, &dtd);
869
	    theta, &rr, &dr, &stp, &dtd);
870
/*     Form the upper half of the pds T = theta*SS + L*D^(-1)*L'; */
870
/*     Form the upper half of the pds T = theta*SS + L*D^(-1)*L'; */
871
/*	  Store T in the upper triangular of the array wt; */
871
/*	  Store T in the upper triangular of the array wt; */
872
/*	  Cholesky factorize T to J*J' with */
872
/*	  Cholesky factorize T to J*J' with */
873
/*	     J' stored in the upper triangular of wt. */
873
/*	     J' stored in the upper triangular of wt. */
874
    formt(m, &wt[wt_offset], &sy[sy_offset], &ss[ss_offset], &col, &theta, &
874
    formt(m, &wt[wt_offset], &sy[sy_offset], &ss[ss_offset], &col, &theta,
875
	    info);
875
	  &info);
876
    if (info != 0) {
876
    if (info != 0) {
877
/*	    nonpositive definiteness in Cholesky factorization; */
-
 
878
/*	    refresh the lbfgs memory and restart the iteration. */
-
 
879
	if (iprint >= 0)
877
	if (iprint >= 0)
880
	    Rprintf("%s\n%s\n",
878
	    Rprintf("%s\n%s\n",
881
		    "Nonpositive definiteness in Cholesky factorization in formk;",
879
		    "Nonpositive definiteness in Cholesky factorization in formt();",
882
		    "   refresh the lbfgs memory and restart the iteration.");
880
		    "   refresh the lbfgs memory and restart the iteration.");
883
	info = 0;
881
	info = 0;
884
	col = 0;
882
	col = 0;
885
	head = 1;
883
	head = 1;
886
	theta = 1.;
884
	theta = 1.;