| 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.;
|