X-Git-Url: https://pintos-os.org/cgi-bin/gitweb.cgi?a=blobdiff_plain;f=lib%2Flinreg%2Fsweep.c;h=000a5e3e0f740ed6b26957ccb2ee5703649db1e6;hb=6ed74d8699a65ef63f04ccef33ab1c8600d6eb9f;hp=8805840e2c4f0261af1cc3df92267919a20d626d;hpb=009350385d97d9fba2b33ea68cfe8b1aeebf1529;p=pspp-builds.git
diff --git a/lib/linreg/sweep.c b/lib/linreg/sweep.c
index 8805840e..000a5e3e 100644
--- a/lib/linreg/sweep.c
+++ b/lib/linreg/sweep.c
@@ -1,30 +1,25 @@
-/* lib/linreg/sweep.c
+/* PSPP - a program for statistical analysis.
+ Copyright (C) 2005, 2009 Free Software Foundation, Inc.
- Copyright (C) 2005 Free Software Foundation, Inc.
- Written by Jason H Stover.
+ This program is free software: you can redistribute it and/or modify
+ it under the terms of the GNU General Public License as published by
+ the Free Software Foundation, either version 3 of the License, or
+ (at your option) any later version.
- This program is free software; you can redistribute it and/or modify
- it under the terms of the GNU General Public License as published by
- the Free Software Foundation; either version 2 of the License, or (at
- your option) any later version.
+ This program is distributed in the hope that it will be useful,
+ but WITHOUT ANY WARRANTY; without even the implied warranty of
+ MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
+ GNU General Public License for more details.
- This program is distributed in the hope that it will be useful, but
- WITHOUT ANY WARRANTY; without even the implied warranty of
- MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
- General Public License for more details.
-
- You should have received a copy of the GNU General Public License
- along with this program; if not, write to the Free Software
- Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA
- 02111-1307, USA.
- */
+ You should have received a copy of the GNU General Public License
+ along with this program. If not, see . */
/*
Find the least-squares estimate of b for the linear model:
Y = Xb + Z
- where Y is an n-by-1 column vector, X is an n-by-p matrix of
+ where Y is an n-by-1 column vector, X is an n-by-p matrix of
independent variables, b is a p-by-1 vector of regression coefficients,
and Z is an n-by-1 normally-distributed random vector with independent
identically distributed components with mean 0.
@@ -45,6 +40,8 @@
Springer. 1998. ISBN 0-387-98542-5.
*/
+#include
+
#include "sweep.h"
/*
@@ -58,7 +55,7 @@
-- --
X refers to the design matrix and Y to the vector of dependent
- observations. reg_sweep sweeps on the diagonal elements of
+ observations. reg_sweep sweeps on the diagonal elements of
X'X.
The matrix A is assumed to be symmetric, so the sweep operation is
@@ -95,7 +92,7 @@ reg_sweep (gsl_matrix * A)
for (j = i; j < A->size2; j++)
{
/*
- Use only the upper triangle of A.
+ Use only the upper triangle of A.
*/
if (j < k)
{